結果
| 問題 | 
                            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 = 998244353
N,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]] += 1
COMB = [0 for j in range(K+2)]
comb = 1
for i in range(N-1):
    comb *= (M+N-1-i) * pow(i+1,mod-2,mod) % mod
    comb %= mod
COMB[0] = comb
for j in range(1,K+2):
    COMB[j] = COMB[j-1] * (M+N-1+j) * pow(N-1+j,mod-2,mod) % mod
    COMB[j] %= mod
tmp = [0] * (K+3)
tmp[-1] = 1
res = 0
for i in range(1,K+1):
    for j in range(-K-1,0):
        if not tmp[j]:
            continue
        tmp[j-1] = -j*tmp[j] % mod
        tmp[j] = 0
    for j in range(-K-1,0)[::-1]:
        if not tmp[j]:
            continue
        tmp[j+1] -= tmp[j]
        tmp[j+1] %= mod
    
    for j in range(-K-1,0):
        res += tmp[j] * (COMB[-j] * cnt[i] % mod) % mod
        res %= mod
print(res)