結果
| 問題 | 
                            No.1693 Invasion
                             | 
                    
| コンテスト | |
| ユーザー | 
                             brthyyjp
                         | 
                    
| 提出日時 | 2021-10-03 13:28:31 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 163 ms / 2,000 ms | 
| コード長 | 786 bytes | 
| コンパイル時間 | 210 ms | 
| コンパイル使用メモリ | 82,448 KB | 
| 実行使用メモリ | 92,288 KB | 
| 最終ジャッジ日時 | 2024-07-21 07:28:25 | 
| 合計ジャッジ時間 | 3,651 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 21 | 
ソースコード
N = 10**6
mod = 998244353
fac = [1]*(N+1)
finv = [1]*(N+1)
for i in range(N):
    fac[i+1] = fac[i] * (i+1) % mod
finv[-1] = pow(fac[-1], mod-2, mod)
for i in reversed(range(N)):
    finv[i] = finv[i+1] * (i+1) % mod
def cmb1(n, r, mod):
    if r <0 or r > n:
        return 0
    r = min(r, n-r)
    return fac[n] * finv[r] * finv[n-r] % mod
n, m = map(int, input().split())
A = list(map(int, input().split()))
INF = 10**18
dp = [INF]*(m+1)
# dp[i]: iをAの要素の重複を許す和で表すときの使う要素の最小値
dp[0] = 0
for i in range(m+1):
    for a in A:
        if i+a <= m:
            dp[i+a] = min(dp[i+a], dp[i]+1)
#print(dp)
ans = 0
for k in range(m+1):
    if dp[k] == INF:
        continue
    ans += cmb1(m-dp[k], k-dp[k], mod)
    ans %= mod
print(ans)
            
            
            
        
            
brthyyjp