結果
問題 | No.1474 かさまJ |
ユーザー | roaris |
提出日時 | 2021-04-09 23:15:53 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,215 ms / 2,500 ms |
コード長 | 1,148 bytes |
コンパイル時間 | 295 ms |
コンパイル使用メモリ | 86,936 KB |
実行使用メモリ | 357,608 KB |
最終ジャッジ日時 | 2023-09-07 13:02:11 |
合計ジャッジ時間 | 8,848 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 78 ms
78,208 KB |
testcase_01 | AC | 76 ms
78,232 KB |
testcase_02 | AC | 107 ms
79,348 KB |
testcase_03 | AC | 81 ms
78,236 KB |
testcase_04 | AC | 79 ms
78,248 KB |
testcase_05 | AC | 494 ms
194,700 KB |
testcase_06 | AC | 197 ms
100,092 KB |
testcase_07 | AC | 85 ms
79,012 KB |
testcase_08 | AC | 91 ms
78,984 KB |
testcase_09 | AC | 735 ms
236,932 KB |
testcase_10 | AC | 129 ms
86,256 KB |
testcase_11 | AC | 90 ms
78,812 KB |
testcase_12 | AC | 102 ms
79,308 KB |
testcase_13 | AC | 867 ms
249,304 KB |
testcase_14 | AC | 300 ms
128,148 KB |
testcase_15 | AC | 201 ms
101,676 KB |
testcase_16 | AC | 128 ms
88,604 KB |
testcase_17 | AC | 361 ms
139,204 KB |
testcase_18 | AC | 1,049 ms
306,764 KB |
testcase_19 | AC | 1,215 ms
343,792 KB |
testcase_20 | AC | 1,120 ms
357,608 KB |
ソースコード
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)