#include <bits/stdc++.h>
using namespace std;
using INT = long long;
int power(int a, int b, int m, int ans = 1) {
  for (; b; b >>= 1, a = 1ll * a * a % m)
    if (b & 1) {
      ans = 1ll * ans * a % m;
    }
  return ans;
}
int dp[301000], a[301000], b[301000];
int s[301000];
int rev[301000];
int main() {
  int n;
  cin >> n;
  int tot = 0;
  rev[0] = 1;
  for (int i = 1; i < 301000; i++) rev[i] = power(i, 998244353 - 2, 998244353);
  for (int i = 1; i <= n; i++) {
    scanf("%d", s + i);
    tot += s[i];
  }
  int flag = 1;
  for (int i = tot - 1; i; i--) {
    int A = (INT)i * rev[tot] % 998244353;
    int B = (INT)(tot - i) * rev[tot] % 998244353 * rev[n - 1] % 998244353;
    int C = (INT)(tot - i) * (n - 2) % 998244353 * rev[tot] % 998244353 *
            rev[n - 1] % 998244353;
    C = (1 - C + 998244353) % 998244353;
    C = (C - (INT)B * a[i + 1] % 998244353 + 998244353) % 998244353;
    B = ((INT)B * b[i + 1] % 998244353 + 1) % 998244353;
    if (C == 0) {
      a[i - 1] = 0;
      b[i - 1] = (INT)power(A, 998244353 - 2, 998244353) * (-B + 998244353) %
                 998244353;
      if (i == 1) dp[0] = b[i - 1], flag = 0;
      i--;
      continue;
    }
    a[i] = (INT)A * power(C, 998244353 - 2, 998244353) % 998244353;
    b[i] = (INT)B * power(C, 998244353 - 2, 998244353) % 998244353;
  }
  if (flag)
    dp[0] =
        (INT)(b[1] + n - 1) *
        power((1 - a[1] + 998244353) % 998244353, 998244353 - 2, 998244353) %
        998244353;
  for (int i = 0; i < tot; i++) {
    int A = (INT)i * rev[tot] % 998244353;
    int B = (INT)(tot - i) * rev[tot] % 998244353 * rev[n - 1] % 998244353;
    int C = (INT)(tot - i) * (n - 2) % 998244353 * rev[tot] % 998244353 *
            rev[n - 1] % 998244353;
    int p = 0;
    if (i) p = dp[i - 1];
    dp[i + 1] = (((INT)dp[i] - (INT)A * p % 998244353 -
                  (INT)C * dp[i] % 998244353 - 1) %
                     998244353 +
                 998244353) %
                998244353;
    dp[i + 1] = (INT)dp[i + 1] * power(B, 998244353 - 2, 998244353) % 998244353;
  }
  int ans = 0;
  for (int i = 1; i <= n; i++) {
    ans += dp[s[i]];
    if (ans >= 998244353) ans -= 998244353;
  }
  ans = ans - (INT)dp[0] * (n - 1) % 998244353;
  ans = (ans % 998244353 + 998244353) % 998244353;
  ans = (INT)ans * power(n, 998244353 - 2, 998244353) % 998244353;
  cout << ans << endl;
}
