結果
問題 | No.1474 かさまJ |
ユーザー | Shibuyap |
提出日時 | 2021-04-10 04:29:02 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 164 ms / 2,500 ms |
コード長 | 1,826 bytes |
コンパイル時間 | 2,090 ms |
コンパイル使用メモリ | 200,080 KB |
実行使用メモリ | 204,560 KB |
最終ジャッジ日時 | 2023-09-08 00:44:34 |
合計ジャッジ時間 | 4,168 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 12 ms
19,944 KB |
testcase_01 | AC | 11 ms
19,968 KB |
testcase_02 | AC | 17 ms
44,520 KB |
testcase_03 | AC | 12 ms
19,992 KB |
testcase_04 | AC | 13 ms
24,084 KB |
testcase_05 | AC | 81 ms
111,364 KB |
testcase_06 | AC | 38 ms
100,276 KB |
testcase_07 | AC | 26 ms
96,268 KB |
testcase_08 | AC | 27 ms
97,420 KB |
testcase_09 | AC | 118 ms
166,732 KB |
testcase_10 | AC | 19 ms
38,344 KB |
testcase_11 | AC | 11 ms
18,208 KB |
testcase_12 | AC | 12 ms
20,060 KB |
testcase_13 | AC | 109 ms
161,380 KB |
testcase_14 | AC | 62 ms
151,940 KB |
testcase_15 | AC | 46 ms
145,540 KB |
testcase_16 | AC | 17 ms
28,252 KB |
testcase_17 | AC | 48 ms
83,548 KB |
testcase_18 | AC | 133 ms
180,420 KB |
testcase_19 | AC | 151 ms
198,472 KB |
testcase_20 | AC | 164 ms
204,560 KB |
ソースコード
#include <bits/stdc++.h> #define rep(i,n) for(int i = 0; i < (n); ++i) #define srep(i,s,t) for (int i = s; i < t; ++i) #define drep(i,n) for(int i = (n)-1; i >= 0; --i) using namespace std; typedef long long int ll; typedef pair<int,int> P; #define yn {puts("YES");}else{puts("NO");} #define MAX_N 200005 const int MAX = 510000; const int MOD = 1000000007; long long fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++){ fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } // 二項係数計算 long long COM(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } long long FINV(int n){ if (n < 0) return 0; return finv[n]; } ll dp[41][41][20001]; int main() { COMinit(); ll n, p, q, l; cin >> n >> p >> q >> l; ll s[n] = {}; rep(i,n) cin >> s[i]; dp[0][0][0] = 1; rep(i,n){ rep(j,i+2){ rep(k,q+1){ dp[i+1][j][k] = dp[i][j][k]; } } rep(j,i+1){ ll sum = 0; rep(k,q+1){ dp[i+1][j+1][k] += sum; dp[i+1][j+1][k] %= MOD; sum += dp[i][j][k]; if(k>=s[i]) sum += MOD - dp[i][j][k-s[i]]; sum %= MOD; } } } ll ans = 0; rep(j,n+1){ rep(k,q+1){ int distributed = j*l - k; if(distributed > p) continue; ans += dp[n][j][k] * COM(p-distributed+(n-1), n-1) % MOD; ans %= MOD; } } cout << ans << endl; return 0; }