C++の練習を兼ねて, AtCoder Regular Contest 144 の 問題D (AND OR Equation) を解いてみた.
■感想.
1. 問題D は, 方針が見えなかったので, 解説を参考に, AC版に到達できたと思う.
2. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト AtCoder Regular Contest 144 解説 の 各リンク を ご覧下さい.
■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 |
// 解き直し. // https://atcoder.jp/contests/arc144/editorial/4307 // 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; const int LIMIT = 303030; LL kComb[LIMIT], nComb[LIMIT], mPow2[LIMIT]; // 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. 入力情報. LL N, K; scanf("%lld %lld", &N, &K); // 2. 2 の 冪乗 を 保存. mPow2[0] = 1; rep(i, (int)N + 1) mPow2[i + 1] = (mPow2[i] << 1) % MOD; // 3. 組み合わせ. kComb[0] = K + 1; nComb[0] = 1; repx(i, 1, (int)N + 1){ kComb[i] = kComb[i - 1] % MOD * max(0LL, (K + 1 - i + MOD) % MOD) % MOD * mPow(i + 1, MOD - 2) % MOD; nComb[i] = nComb[i - 1] % MOD * max(0LL, (N + 1 - i + MOD) % MOD) % MOD * mPow(i, MOD - 2) % MOD; } // 4. 集計. LL ans = 0; rep(i, (int)N + 1) ans = (ans + mPow2[i] * nComb[i] % MOD * kComb[i] % MOD) % MOD; // 5. 出力. 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 |
[入力例] 2 1 [出力例] 6 ※AtCoderテストケースより [入力例] 2 2 [出力例] 19 ※AtCoderテストケースより [入力例] 100 123456789123456789 [出力例] 34663745 ※AtCoderテストケースより [入力例] 3 2 [出力例] 33 [入力例] 11 5 [出力例] 71000 [入力例] 2022 1123 [出力例] 560596535 [入力例] 77777 202207162100 [出力例] 836824508 |
■参照サイト
AtCoder Regular Contest 144