C++の練習を兼ねて, AtCoder Regular Contest 113 の 問題D (Sky Reflector) を解いてみた.
■感想.
1. 問題Dは, 方針が見えなかったので, 解説を参考に実装して, AC版に到達できたので, 良かったと思う.
2. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイト AtCoder Regular Contest 113 解説 の 各リンクを ご覧下さい.
■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 |
// 解き直し. // https://atcoder.jp/contests/arc113/editorial/709 // 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 = 998244353; // 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 % MOD; } int main(){ // 1. 入力情報. LL N, M, K; scanf("%lld %lld %lld", &N, &M, &K); // 2. 例外ケース. // 2-1. N が 1. if(N == 1){ printf("%lld\n", mPow(K, M)); return 0; } // 2-2. M が 1. if(M == 1){ printf("%lld\n", mPow(K, N)); return 0; } // 3. M, N が 2以上の場合. LL ans = 0; repx(x, 1, (int)(K + 1)){ // x の N乗. LL a = mPow((LL)x, N); // (x - 1) の N乗. LL b = mPow((LL)((x - 1 + MOD) % MOD), N); // (K - x + 1) の M乗. LL c = mPow((LL)((K - x + 1 + MOD) % MOD), M); // 集計. ans = (ans + (a - b + MOD) * c) % MOD; } // 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 |
[入力例] 2 2 2 [出力例] 7 ※AtCoderのテストケースより [入力例] 1 1 100 [出力例] 100 ※AtCoderのテストケースより [入力例] 31415 92653 58979 [出力例] 469486242 ※AtCoderのテストケースより [入力例] 3 3 2 [出力例] 15 [入力例] 4 3 2 [出力例] 23 [入力例] 12345 23456 34567 [出力例] 134323523 |
■参照サイト
AtCoder Regular Contest 113