結果
問題 |
No.1474 かさまJ
|
ユーザー |
|
提出日時 | 2021-04-09 23:15:53 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,258 ms / 2,500 ms |
コード長 | 1,148 bytes |
コンパイル時間 | 634 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 357,236 KB |
最終ジャッジ日時 | 2024-06-25 07:06:29 |
合計ジャッジ時間 | 8,534 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 17 |
ソースコード
MAX = 10**5 MOD = 10**9+7 fact = [0]*MAX #fact[i]: i! inv = [0]*MAX #inv[i]: iの逆元 finv = [0]*MAX #finv[i]: i!の逆元 fact[0] = 1 fact[1] = 1 finv[0] = 1 finv[1] = 1 inv[1] = 1 for i in range(2, MAX): fact[i] = fact[i-1]*i%MOD inv[i] = MOD-inv[MOD%i]*(MOD//i)%MOD finv[i] = finv[i-1]*inv[i]%MOD def C(n, r): if n<r: return 0 if n<0 or r<0: return 0 return fact[n]*(finv[r]*finv[n-r]%MOD)%MOD N, Mp, Mq, L = map(int, input().split()) S = list(map(int, input().split())) dp = [[[0]*(Mq+1) for _ in range(N+1)] for _ in range(N+1)] dp[0][0][0] = 1 for i in range(N): for j in range(i+1): acc = [0] for k in range(Mq+1): acc.append((acc[-1]+dp[i][j][k])%MOD) dp[i+1][j][k] += dp[i][j][k] dp[i+1][j][k] %= MOD for k in range(Mq+1): dp[i+1][j+1][k] += (acc[k]-acc[max(k-S[i], 0)])%MOD dp[i+1][j+1][k] %= MOD ans = 0 for i in range(N+1): for j in range(Mq+1): p = i*L-j if Mp>=p: ans += dp[N][i][j]*C(N-1+Mp-p, Mp-p) ans %= MOD print(ans)