#include <bits/stdc++.h>
using namespace std;
int n, A, B, w[105][105], cnt, f[8000], s[8000][205], val[8000], ans[8000],
    p[5], sum;
int tx[] = {-1, 0, 1, 0}, ty[] = {0, -1, 0, 1};
int qpow(int x, int y) {
  int ret = 1;
  while (y) {
    if (y & 1) ret = (long long)ret * x % 1000000007;
    x = (long long)x * x % 1000000007;
    y >>= 1;
  }
  return ret;
}
int gauss() {
  int i, j, k;
  for (i = 1; i <= cnt; i++) {
    for (j = i + 1; j <= cnt && (j - i) <= A; j++)
      if (s[j][A - (j - i)]) {
        int num = (long long)s[j][A - (j - i)] * qpow(s[i][A], 1000000007 - 2) %
                  1000000007;
        for (k = A; k <= B; k++) {
          s[j][k - (j - i)] =
              (s[j][k - (j - i)] - (long long)num * s[i][k] % 1000000007 +
               1000000007) %
              1000000007;
        }
        val[j] = (val[j] - (long long)num * val[i] % 1000000007 + 1000000007) %
                 1000000007;
      }
  }
  for (i = cnt; i >= w[n][n]; i--) {
    ans[i] = (long long)val[i] * qpow(s[i][A], 1000000007 - 2) % 1000000007;
    for (int j = i - 1; j && i - j <= A; j--)
      if (s[j][A + (i - j)])
        val[j] = (val[j] - (long long)ans[i] * s[j][A + (i - j)] % 1000000007 +
                  1000000007) %
                 1000000007;
  }
  return ans[w[n][n]];
}
int main() {
  int i, j, k;
  scanf("%d", &n);
  A = n * 2;
  B = n * 4;
  for (i = 0; i < 4; i++) {
    scanf("%d", p + i);
    sum += p[i];
  }
  for (i = 0; i <= A; i++)
    for (j = 0; j <= A; j++)
      if ((i - n) * (i - n) + (j - n) * (j - n) <= n * n) w[i][j] = ++cnt;
  for (i = 0; i <= A; i++)
    for (j = 0; j <= A; j++)
      if (w[i][j]) {
        for (k = 0; k < 4; k++) {
          int x = i + tx[k], y = j + ty[k];
          if (w[x][y]) s[w[i][j]][w[x][y] - w[i][j] + A] = p[k];
        }
        s[w[i][j]][A] = 1000000007 - sum;
        val[w[i][j]] = 1000000007 - sum;
      }
  printf("%d\n", gauss());
  return 0;
}
