#include #include using namespace std; using mint = atcoder::modint998244353; const int mx = 1<<17; int cnt[mx]; mint dp[mx], fact[mx], ifact[mx], pow2[mx], ipow2[mx]; void table_build(){ fact[0] = 1; for (int i = 1; i < mx; i++){ fact[i] = fact[i-1] * i; } ifact[mx-1] = fact[mx-1].inv(); for (int i = mx-2; i >= 0; i--){ ifact[i] = ifact[i+1] * (i+1); } pow2[0] = ipow2[0] = 1; mint i2 = mint(2).inv(); for (int i = 1; i < mx; i++){ pow2[i] = pow2[i-1] * 2; ipow2[i] = ipow2[i-1] * i2; } } array,19> fs, ex; int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); int n; cin >> n; while (n--){ int x; cin >> x; cnt[x]++; } table_build(); for (int d = 1; d < 19; d++){ fs[d].resize(1< void { if (r - l == 1){ dp[l] += ifact[l]; dp[l] *= pow2[cnt[l]]-1; ans += dp[l] * fact[l]; return ; } int m = (l + r) / 2; sfs(sfs,l,m,k-1); if (k <= 6){ for (int i = l; i < m; i++) for (int j = m; j < r; j++){ dp[j] += dp[i] * ifact[j-i]; } } else { fill(fs[k+1].begin(),fs[k+1].end(),0); for (int i = l; i < m; i++){ fs[k+1][i-l] = dp[i]; } atcoder::internal::butterfly(fs[k+1]); for (int i = 0; i < (1<<(k+1)); i++){ fs[k+1][i] *= ex[k+1][i]; } atcoder::internal::butterfly_inv(fs[k+1]); for (int i = (1<<(k-1)); i < (1<