#include <bits/stdc++.h>
using namespace std;
long long ex[500001], coeff[500001], c[500001];
long long rev[500001];
int n, a[500001], m;
long long bigmod(long long a, long long p) {
  if (!p) return 1;
  long long t = bigmod(a, p >> 1);
  t = t * t % 998244353;
  if (p & 1) t = t * a % 998244353;
  return t;
}
long long dp[500001];
int main() {
  scanf("%d", &n);
  for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
  for (int i = 1; i <= n; i++) m += a[i];
  long long rm = bigmod(m, 998244353 - 2);
  long long bb = bigmod(n - 1, 998244353 - 2);
  ex[0] = (n - 1);
  for (int i = 1; i <= m - 1; i++) {
    long long cf = (i * rm % 998244353 + (m - i) * rm % 998244353 * (n - 2) %
                                             998244353 * bb % 998244353) %
                   998244353;
    cf = (1 - cf + 998244353) % 998244353;
    long long base =
        (1 + i * rm % 998244353 * ex[i - 1] % 998244353) % 998244353;
    ex[i] = base * bigmod(cf, 998244353 - 2) % 998244353;
  }
  for (int i = m - 1; i >= 0; i--) {
    dp[i] = (dp[i + 1] + ex[i]) % 998244353;
  }
  long long ret = 0;
  for (int i = 1; i <= n; i++) ret = (ret + dp[a[i]]) % 998244353;
  ret = (ret + 998244353 - (n - 1) * dp[0] % 998244353) % 998244353;
  cout << ret * bigmod(n, 998244353 - 2) % 998244353 << endl;
}
