C++の練習を兼ねて, AtCoder Regular Contest 141 の 問題B (Increasing Prefix XOR) を解いてみた.
■感想.
1. 問題Bは, 方針が見えなかったので, 解説を参考に, AC版に到達できたと思う.
2. 苦手な動的計画法の訓練を積めたので, 非常に良かったと思う.
3. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト AtCoder Regular Contest 141 解説 の 各リンク を ご覧下さい.
■C++版プログラム(問題B/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 60 61 62 63 64 65 66 |
// 解き直し. // https://atcoder.jp/contests/arc141/editorial/4022 // C++(GCC 9.2.1) #include <bits/stdc++.h> using namespace std; using LL = long long; using vl = vector<LL>; #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 = 998244353; LL dp[66][66]; int main(){ // 1. 入力情報. LL N, M; scanf("%lld %lld", &N, &M); // 2. 例外. if(N > 60){ puts("0"); return 0; } // 3. 初期化. LL cPow2 = 1; int d = 0; repx(k, 1, 61){ // 設定. dp[1][k] = (cPow2 <= M) ? cPow2 : M; dp[1][k] %= MOD; // 減算. M -= cPow2; // 終了条件. if(M <= 0){ d = k; break; } // 2倍. cPow2 <<= 1LL; } // 4. dp更新. repx(i, 2, (int)N + 1){ repx(k1, 1, d + 1){ repx(k2, k1 + 1, d + 1){ dp[i][k2] += (dp[i - 1][k1] * dp[1][k2] % MOD); dp[i][k2] %= MOD; } } } // 5. 集計. LL ans = 0; repx(k, (int)N, d + 1) ans = (ans + dp[(int)N][k]) % MOD; // 6. 出力. printf("%lld\n", ans); 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 |
[入力例] 2 4 [出力例] 5 ※AtCoderのテストケースより [入力例] 4 4 [出力例] 0 ※AtCoderのテストケースより [入力例] 10 123456789 [出力例] 205695670 ※AtCoderのテストケースより [入力例] 3 14 [出力例] 106 [入力例] 5 123 [出力例] 2564352 [入力例] 3 234 [出力例] 665226 [入力例] 12 34567 [出力例] 340820379 [入力例] 50 12345678987654321 [出力例] 967073402 |
■参照サイト
AtCoder Regular Contest 141