#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
const int lim = 1e5;
int a[200007];
long long f[200007];
long long pw(long long x, long long p) {
  long long ans = 1;
  while (p) {
    if (p & 1) ans = ans * x % mod;
    x = x * x % mod, p >>= 1;
  }
  return ans;
}
int main() {
  int n;
  long long s = 0;
  scanf("%d", &n);
  for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
  for (int i = 1; i <= n; i++) s += a[i];
  f[1] = (s - 1) * (s - 1) % mod * pw(s, mod - 2) % mod;
  long long sum = f[1];
  for (int i = 1; i < lim; i++) {
    sum = (sum - (s - 1) * pw(s - i, mod - 2) % mod + mod) % mod;
    f[i + 1] = (f[i] + sum) % mod;
  }
  long long ans = 0;
  for (int i = 1; i <= n; i++) ans = (ans + f[a[i]]) % mod;
  printf("%lld\n", ans);
  return 0;
}
