C++の練習を兼ねて, AtCoder Regular Contest 166 の 問題C (LU / RD Marking) を解いてみた.
■感想.
1. 問題C は, 方針が見えなかったので, 解説を参考に実装したところ, AC版に到達出来た.
2. 個人的には, 問題C で, Fibonacci数列の復習が出来たので, 良かったとおもう.
3. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト AtCoder Regular Contest 166 解説 の 各リンク を ご覧下さい.
■C++版プログラム(問題C／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 |
// 解き直し. // https://atcoder.jp/contests/arc166/editorial/7184 // C++20(GCC 12.2) #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 = 998244353; LL f[2020202], g[2020202]; // Fermat's little theorem から, 大きな冪乗の計算を行う. // @param a: べき乗したい正整数. // @param b: 指数. // @return: べき乗した計算結果(mod版). LL mPow(LL a, LL b){ LL t = 1; while(b){ if(b & 1) t = (t * a) % MOD; a = a * a % MOD; b >>= 1; } return t; } int main(){ // 1. 入力情報. int T; scanf("%d", &T); // 2. 事前準備. f[0] = f[1] = 1; repx(i, 2, 2020202) f[i] = (f[i - 1] + f[i - 2]) % MOD; g[0] = 1; repx(i, 1, 1010101 - 1) g[i] = (g[i - 1] * f[i * 2]) % MOD; // 3. テストケース. // ex. H = 3, W = 8 // 上辺(左から右): 2, 4, 6, 7, 7, 7, 7, 7 // 右辺(上から下): 6, 4, 2 // -> 解説の数式に照らすと, n = 3, m = 5 が分かる. rep(i, T){ // 3-1. 各テストケース. int H, W; scanf("%d %d", &H, &W); // 3-2. m, n. int m = abs(H - W), n = min(H, W); // 3-3. 計算. LL ans = g[n] * g[n] % MOD * mPow(f[2 * n + 1], (LL)m) % MOD; // 3-4. 出力. 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 |
[入力例] 2 1 1 2 3 [出力例] 4 800 ※AtCoderのテストケースより [入力例] 3 123 456 654 321 1000000 1000000 [出力例] 60549740 298307903 656009181 ※AtCoderのテストケースより [入力例] 10 1 2 2 1 3 4 5 3 3 6 3 7 3 8 5 10 15 25 2023 1111 [出力例] 12 12 354900 7452900 156510900 291995841 142446543 546223419 631651090 451184128 |
■参照サイト
AtCoder Regular Contest 166