結果
問題 | No.2005 Sum of Power Sums |
ユーザー |
|
提出日時 | 2022-02-12 22:30:11 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 918 ms / 2,000 ms |
コード長 | 856 bytes |
コンパイル時間 | 219 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 105,600 KB |
最終ジャッジ日時 | 2024-06-12 05:56:31 |
合計ジャッジ時間 | 6,755 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
mod = 998244353N,M = map(int,input().split())k = list(map(int,input().split()))K = max(k)cnt = [0 for i in range(K+1)]for i in range(N):cnt[k[i]] += 1COMB = [0 for j in range(K+2)]comb = 1for i in range(N-1):comb *= (M+N-1-i) * pow(i+1,mod-2,mod) % modcomb %= modCOMB[0] = combfor j in range(1,K+2):COMB[j] = COMB[j-1] * (M+N-1+j) * pow(N-1+j,mod-2,mod) % modCOMB[j] %= modtmp = [0] * (K+3)tmp[-1] = 1res = 0for i in range(1,K+1):for j in range(-K-1,0):if not tmp[j]:continuetmp[j-1] = -j*tmp[j] % modtmp[j] = 0for j in range(-K-1,0)[::-1]:if not tmp[j]:continuetmp[j+1] -= tmp[j]tmp[j+1] %= modfor j in range(-K-1,0):res += tmp[j] * (COMB[-j] * cnt[i] % mod) % modres %= modprint(res)