結果
問題 | No.1474 かさまJ |
ユーザー | googol_S0 |
提出日時 | 2021-04-09 21:55:14 |
言語 | PyPy3 (7.3.13) |
結果 |
AC
|
実行時間 | 1,413 ms / 2,500 ms |
コード長 | 935 bytes |
コンパイル時間 | 502 ms |
コンパイル使用メモリ | 86,924 KB |
実行使用メモリ | 361,332 KB |
最終ジャッジ日時 | 2023-09-07 11:03:15 |
合計ジャッジ時間 | 10,629 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge14 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 89 ms
83,948 KB |
testcase_01 | AC | 86 ms
83,740 KB |
testcase_02 | AC | 108 ms
84,552 KB |
testcase_03 | AC | 88 ms
84,024 KB |
testcase_04 | AC | 85 ms
83,780 KB |
testcase_05 | AC | 668 ms
202,332 KB |
testcase_06 | AC | 222 ms
101,256 KB |
testcase_07 | AC | 101 ms
84,112 KB |
testcase_08 | AC | 103 ms
83,996 KB |
testcase_09 | AC | 930 ms
242,708 KB |
testcase_10 | AC | 157 ms
93,280 KB |
testcase_11 | AC | 105 ms
84,520 KB |
testcase_12 | AC | 114 ms
85,108 KB |
testcase_13 | AC | 912 ms
260,272 KB |
testcase_14 | AC | 373 ms
130,716 KB |
testcase_15 | AC | 226 ms
102,316 KB |
testcase_16 | AC | 153 ms
92,748 KB |
testcase_17 | AC | 410 ms
148,176 KB |
testcase_18 | AC | 1,175 ms
310,104 KB |
testcase_19 | AC | 1,337 ms
347,304 KB |
testcase_20 | AC | 1,413 ms
361,332 KB |
ソースコード
N,P,Q,L=map(int,input().split()) S=list(map(int,input().split())) mod=10**9+7 def cmb(n,r): if r<0 or r>n: return 0 return (g1[n]*g2[r]*g2[n-r])%mod g1=[1,1] g2=[1,1] inv=[0,1] for i in range(2,50003): g1.append((g1[-1]*i)%mod) inv.append((-inv[mod%i]*(mod//i))%mod) g2.append((g2[-1]*inv[-1])%mod) DP=[[[0]*(Q+2) for j in range(i+2)] for i in range(N+1)] DP[0][0][0]=1 C=[[[0]*(Q+2) for j in range(i+2)] for i in range(N+1)] for i in range(N+1): for j in range(i+1): for k in range(Q+1): C[i][j][k]%=mod DP[i][j][k]+=C[i][j][k] DP[i][j][k]%=mod C[i][j][k+1]+=C[i][j][k] if i==N: break for j in range(i+1): for k in range(Q+1): DP[i+1][j+1][k]+=DP[i][j][k] C[i+1][j][k+1]+=DP[i][j][k] C[i+1][j][min(Q+1,k+S[i]+1)]-=DP[i][j][k] ANS=0 for i in range(N+1): for j in range(Q+1): x=L*(N-i)-j if x<=P: ANS+=cmb(N+P-x-1,P-x)*DP[N][i][j] print(ANS%mod)