C++の練習を兼ねて, AtCoder Beginner Contest 122 の 問題D (We Like AGC) を解いてみた.
■感想.
1. 問題Dは, 方針が見えなかったので, 解説を参考に, AC版に到達できたと思う.
2. 実装に苦労したものの, 苦手な動的計画法の訓練を積めたので, 非常に良かったと思う.
3. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト AtCoder Beginner Contest 122 解説 を ご覧下さい.
■C++版プログラム(問題D／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 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 |
// 解き直し. // https://img.atcoder.jp/abc122/editorial.pdf // C++(GCC 9.2.1) #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; // (i - 3) 文字目, (i - 2) 文字目, (i - 1) 文字目, i 文字目 を 管理. // 文字 A, C, G, T は, 0, 1, 2, 3 に 対応させる. LL dp[101][4][4][4][4]; string sAgc[] = {"AGCA", "GACA", "ACGA", "AAGC", "CAGC", "GAGC", "TAGC", "AGAC", "CGAC", "GGAC", "TGAC", "ACGC", "AGGC", "ATGC", "AGCC", "AGTC", "AGCG", "CACG", "ACGG", "AACG", "CACG", "TACG", "AGCT", "GACT", "ACGT"}; int agc[303]; int main(){ // 1. 入力情報. int N; scanf("%d", &N); // 2. 除外したい部分文字列. // -> 解説にあるように, 特定の文字列が含まれない場合をチェック. // 文字 A, C, G, T は, 0, 1, 2, 3 に 対応させる. map<char, int> m; m['A'] = 0; m['C'] = 1; m['G'] = 2; m['T'] = 3; for(auto &s : sAgc){ int x = 0; x += m[s[0]] * 64; x += m[s[1]] * 16; x += m[s[2]] * 4; x += m[s[3]] * 1; ++agc[x]; } // 3. dp初期化. rep(i, 4) dp[0][0][0][0][i] = 1; rep(i, 4) rep(j, 4) dp[1][0][0][j][i] += dp[0][0][0][0][j]; rep(a, 4){ rep(b, 4){ rep(c, 4){ // 特定の部分文字列を含むか？ // AGC (021). if(c == 0 && b == 2 && a == 1) continue; // GAC (201). if(c == 2 && b == 0 && a == 1) continue; // ACG (012). if(c == 0 && b == 1 && a == 2) continue; // 特定の部分文字列を含まない場合. dp[2][0][c][b][a] += dp[1][0][0][c][b]; } } } // 4. dp更新. repx(i, 3, N){ // i 文字目. rep(a, 4){ // (i - 1) 文字目. rep(b, 4){ // (i - 2) 文字目. rep(c, 4){ // (i - 3) 文字目. rep(d, 4){ // 除外したい部分文字列であるか？ int x = 0; x += d * 64; x += c * 16; x += b * 4; x += a * 1; if(agc[x]) continue; // 前回分反映. rep(e, 4){ dp[i][d][c][b][a] += dp[i - 1][e][d][c][b]; dp[i][d][c][b][a] %= MOD; } } } } } } // 5. 集計. LL ans = 0; rep(a, 4) rep(b, 4) rep(c, 4) rep(d, 4) ans = (ans + dp[N - 1][d][c][b][a]) % 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 |
[入力例] 3 [出力例] 61 ※AtCoderテストケースより [入力例] 4 [出力例] 230 ※AtCoderテストケースより [入力例] 100 [出力例] 388130742 ※AtCoderテストケースより [入力例] 5 [出力例] 865 [入力例] 6 [出力例] 3247 [入力例] 10 [出力例] 643550 [入力例] 12 [出力例] 9058467 [入力例] 77 [出力例] 781042769 |
■参照サイト
AtCoder Beginner Contest 122