結果
問題 | No.1474 かさまJ |
ユーザー |
|
提出日時 | 2021-02-19 13:30:02 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,220 ms / 2,500 ms |
コード長 | 888 bytes |
コンパイル時間 | 454 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 83,328 KB |
最終ジャッジ日時 | 2024-09-16 20:41:04 |
合計ジャッジ時間 | 9,008 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 17 |
ソースコード
N, MP, Mq, L = map(int, input().split()) S = list(map(int, input().split())) MOD = 10 ** 9 + 7 dp = [[0] * (Mq + 1) for k in range(N + 1)] dp[0][Mq] = 1 for i in range(1, min(S[0], Mq) + 1): dp[1][Mq - i] = 1 for i in range(1, N): s = S[i] for k in range(N, 0, -1): summ = 0 for r in range(Mq, -1, -1): dp[k][r] = (dp[k][r] + summ) % MOD summ = (summ + dp[k - 1][r]) % MOD if r + s <= Mq: summ = (summ - dp[k - 1][r + s]) % MOD fac = [1] * (50000) for i in range(1, 50000): fac[i] = (fac[i - 1] * i) % MOD def nCr(n, r): if n < r: return 0 return (fac[n] * pow(fac[n - r], MOD - 2, MOD) * pow(fac[r], MOD - 2, MOD)) % MOD ans = 0 for m in range(Mq + 1): for x in range(N + 1): f = dp[x][Mq - m] ans = (ans + f * nCr(MP + m - L * x + N - 1, N - 1)) % MOD print(ans)