#include <bits/stdc++.h>
const int N = 10005;
const int mod = 1e9 + 7;
using namespace std;
int po(int n, int k) {
  int ans = 1;
  while (k) {
    if (k & 1) ans = 1LL * ans * n % mod;
    n = 1LL * n * n % mod;
    k >>= 1;
  }
  return ans;
}
vector<int> mv[N];
int dx[] = {-1, 0, 1, 0}, dy[] = {0, -1, 0, 1};
int R, a[4], id[101][101], sum, cnt;
int ans[N];
int ck(int i, int j) {
  return (i - 50) * (i - 50) + (j - 50) * (j - 50) <= R * R;
}
void add(int& a, int b) {
  a += b;
  if (a >= mod) a -= mod;
  if (a < 0) a += mod;
}
int nani[] = {666666673, 666666673, 666666674, 666666673, 666666673};
bool ck() {
  for (int i = 0; i < cnt; i++) {
    int sum = 0;
    for (int j = 0; j < cnt; j++) {
      sum = (sum + 1LL * nani[j] * mv[i][j]) % mod;
    }
    if (sum != mv[i][cnt]) return false;
  }
  return true;
}
int main() {
  cin >> R;
  for (int i = 0; i < 4; i++) cin >> a[i], sum += a[i];
  sum = po(sum, mod - 2);
  for (int i = 0; i < 4; i++) a[i] = 1LL * a[i] * sum % mod;
  for (int i = 0; i <= 100; i++)
    for (int j = 0; j <= 100; j++)
      if (ck(i, j)) id[i][j] = cnt++;
  for (int i = 0; i <= 100; i++)
    for (int j = 0; j <= 100; j++) {
      if (ck(i, j)) {
        int num = id[i][j];
        mv[num].resize(cnt + 1);
        mv[num][num] = mv[num][cnt] = mod - 1;
        for (int nid = 0; nid < 4; nid++) {
          int x = i + dx[nid], y = j + dy[nid];
          if (ck(x, y)) mv[num][id[x][y]] = a[nid];
        }
      }
    }
  for (int i = 0; i < cnt; i++) {
    for (int j = 0; j < 120; j++)
      if (i + j < cnt && mv[i + j][i]) {
        swap(mv[i], mv[i + j]);
        break;
      }
    for (int j = 0; j < 120; j++)
      if (j && 0 <= i + j && i + j < cnt && mv[i + j][i]) {
        int mul = 1LL * mv[i + j][i] * po(mv[i][i], mod - 2) % mod;
        for (int k = 0; k < 120; k++) {
          if (i + k < cnt)
            mv[i + j][i + k] = (mv[i + j][i + k] - 1LL * mv[i][i + k] * mul +
                                1LL * mod * mod) %
                               mod;
        }
        mv[i + j][cnt] =
            (mv[i + j][cnt] - 1LL * mv[i][cnt] * mul + 1LL * mod * mod) % mod;
      }
  }
  for (int i = cnt - 1; i >= 0; i--) {
    int sum = mv[i][cnt];
    for (int j = i + 1; j < cnt; j++)
      sum = (sum - 1LL * ans[j] * mv[i][j] + 1LL * mod * mod) % mod;
    ans[i] = 1LL * sum * po(mv[i][i], mod - 2) % mod;
  }
  cout << ans[id[50][50]];
}
