結果
問題 | No.1474 かさまJ |
ユーザー | ぷら |
提出日時 | 2021-04-09 23:52:06 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 407 ms / 2,500 ms |
コード長 | 1,810 bytes |
コンパイル時間 | 3,163 ms |
コンパイル使用メモリ | 167,528 KB |
実行使用メモリ | 272,196 KB |
最終ジャッジ日時 | 2023-09-07 18:50:17 |
合計ジャッジ時間 | 5,547 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 15 ms
9,848 KB |
testcase_01 | AC | 14 ms
9,560 KB |
testcase_02 | AC | 16 ms
10,436 KB |
testcase_03 | AC | 15 ms
9,652 KB |
testcase_04 | AC | 15 ms
9,628 KB |
testcase_05 | AC | 192 ms
115,056 KB |
testcase_06 | AC | 47 ms
32,740 KB |
testcase_07 | AC | 18 ms
12,020 KB |
testcase_08 | AC | 21 ms
15,072 KB |
testcase_09 | AC | 269 ms
166,564 KB |
testcase_10 | AC | 28 ms
16,188 KB |
testcase_11 | AC | 16 ms
9,900 KB |
testcase_12 | AC | 15 ms
9,812 KB |
testcase_13 | AC | 265 ms
175,152 KB |
testcase_14 | AC | 99 ms
63,552 KB |
testcase_15 | AC | 54 ms
36,320 KB |
testcase_16 | AC | 26 ms
14,596 KB |
testcase_17 | AC | 106 ms
65,724 KB |
testcase_18 | AC | 336 ms
223,620 KB |
testcase_19 | AC | 387 ms
259,560 KB |
testcase_20 | AC | 407 ms
272,196 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; int mod = 1000000007; long long dp[45][45][20005]; long long fac[200005], finv[200005], inv[200005]; void COMinit(int m) { fac[0] = fac[1] = finv[0] = finv[1] = inv[1] = 1; for (int i = 2; i < 200005; i++){ fac[i] = fac[i - 1] * i % m; inv[i] = m - inv[m % i] * (m / i) % m; finv[i] = finv[i - 1] * inv[i] % m; } } long long COM(int n, int k,int m){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % m) % m; } int main() { int N,Mp,Mq,L; cin >> N >> Mp >> Mq >> L; vector<int>S(N); for(int i = 0; i < N; i++) { cin >> S[i]; } dp[0][0][0] = 1; for(int i = 0; i < N; i++) { for(int j = 0; j <= N; j++) { long long sum = 0; for(int k = 0; k <= Mq; k++) { dp[i+1][j][k] += dp[i][j][k]; if(dp[i+1][j][k] > mod) { dp[i+1][j][k] -= mod; } if(k > S[i]) { sum += mod-dp[i][j][k-S[i]-1]; if(sum > mod) { sum -= mod; } } dp[i+1][j+1][k] += sum; if(dp[i+1][j+1][k] > mod) { dp[i+1][j+1][k] -= mod; } sum += dp[i][j][k]; if(sum > mod) { sum -= mod; } } } } long long ans = 0; COMinit(mod); for(int i = 0; i <= N; i++) { for(int j = 0; j <= Mq; j++) { int p = Mp-(L*i-j); ans += dp[N][i][j]*COM(p+N-1,N-1,mod)%mod; if(ans > mod) { ans -= mod; } } } cout << ans << endl; }