結果
問題 | No.1474 かさまJ |
ユーザー | ぷら |
提出日時 | 2021-04-09 23:58:42 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 380 ms / 2,500 ms |
コード長 | 1,820 bytes |
コンパイル時間 | 1,527 ms |
コンパイル使用メモリ | 169,000 KB |
実行使用メモリ | 272,212 KB |
最終ジャッジ日時 | 2023-09-07 18:57:26 |
合計ジャッジ時間 | 4,573 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 15 ms
9,528 KB |
testcase_01 | AC | 14 ms
9,564 KB |
testcase_02 | AC | 16 ms
10,424 KB |
testcase_03 | AC | 15 ms
9,532 KB |
testcase_04 | AC | 15 ms
9,644 KB |
testcase_05 | AC | 167 ms
115,112 KB |
testcase_06 | AC | 45 ms
32,976 KB |
testcase_07 | AC | 18 ms
12,088 KB |
testcase_08 | AC | 21 ms
14,792 KB |
testcase_09 | AC | 237 ms
166,560 KB |
testcase_10 | AC | 25 ms
16,272 KB |
testcase_11 | AC | 15 ms
9,848 KB |
testcase_12 | AC | 15 ms
9,816 KB |
testcase_13 | AC | 247 ms
175,212 KB |
testcase_14 | AC | 90 ms
63,416 KB |
testcase_15 | AC | 51 ms
36,388 KB |
testcase_16 | AC | 22 ms
14,616 KB |
testcase_17 | AC | 95 ms
65,716 KB |
testcase_18 | AC | 310 ms
223,548 KB |
testcase_19 | AC | 363 ms
259,564 KB |
testcase_20 | AC | 380 ms
272,212 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; constexpr 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; }