C++の練習を兼ねて, AtCoder Regular Contest 059 の 問題E (キャンディーとN人の子供) を解いてみた.
■感想.
1. 問題E は, 解答方針が, 全く見えなかったので, 解答を参照して, 実装して何とかAC版となった.
2. 2. 苦手なdpの訓練を積めたので, 非常に良かったと思う.
3. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイトARC 059 解説をご覧下さい.
■C++版プログラム(問題E／AC版).
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 |
// 解き直し. // https://img.atcoder.jp/data/arc/059/editorial.pdf #include <bits/stdc++.h> using namespace std; using LL = long long; #define repex(i, a, b, c) for(int i = a; i < b; i += c) #define repx(i, a, b) repex(i, a, b, 1) #define rep(i, n) repx(i, 0, n) #define repr(i, a, b) for(int i = a; i >= b; i--) const LL MOD = 1e9 + 7; LL dp[404][404]; LL a[404], b[404], abPow[404][404], cPow[404][404]; int main(){ // 1. 入力情報. int N, C; scanf("%d %d", &N, &C); rep(i, N) scanf("%lld", &a[i]); rep(i, N) scanf("%lld", &b[i]); // 2. 冪乗 を 保存. rep(i, 404) abPow[i][0] = 1; repx(i, 1, 403){ rep(k, 403){ abPow[i][k + 1] = abPow[i][k] * (LL)i; abPow[i][k + 1] %= MOD; } } // 3. 累積和 を 保存. rep(k, 403) rep(i, 403) cPow[i][k] = abPow[i][k]; rep(k, 403){ rep(i, 403){ cPow[i + 1][k] += cPow[i][k]; cPow[i + 1][k] %= MOD; } } // rep(k, C + 1){ // rep(i, N + 1) printf("k=%d i=%d %lld ", k, i, cPow[i][k]); // puts(""); // } // 4. dp更新. dp[0][0] = 1; rep(i, N + 1){ rep(j, C + 1){ rep(k, j + 1){ dp[i + 1][j] += dp[i][j - k] * (cPow[b[i]][k] + MOD - cPow[a[i] - 1][k]); dp[i + 1][j] %= MOD; } } } // 5. 出力. printf("%lld\n", dp[N][C]); return 0; } |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 |
[入力例] 2 3 1 1 1 1 [出力例] 4 ※AtCoderのテストケースより [入力例] 1 2 1 3 [出力例] 14 ※AtCoderのテストケースより [入力例] 2 3 1 1 2 2 [出力例] 66 ※AtCoderのテストケースより [入力例] 4 8 3 1 4 1 3 1 4 1 [出力例] 421749 ※AtCoderのテストケースより [入力例] 3 100 7 6 5 9 9 9 [出力例] 139123417 ※AtCoderのテストケースより [入力例] 5 6 1 2 3 4 5 2 3 4 5 6 [出力例] 19291985 |
■参照サイト
AtCoder Regular Contest 059