C++の練習を兼ねて, AtCoder Regular Contest 110 の 問題D (Binomial Coefficient is Fun) を解いてみた.
■感想.
1. 問題D は, 解答の方針が見えなかったので, 解説を参照して, 実装したところ, ようやくAC版に到達できた.
2. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイト AtCoder Regular Contest 110 解説 の 各リンク を ご覧下さい.
■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 |
// 解き直し. // https://atcoder.jp/contests/arc110/editorial/385 // 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; const int LIMIT = 4040404; LL a[2020], FAC[LIMIT], INV[LIMIT], MNFAC[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 % MOD; } int main(){ // 1. 入力情報. LL N, M, S = 0; scanf("%lld %lld", &N, &M); rep(i, (int)N) scanf("%lld", &a[i]); rep(i, (int)N) S += a[i]; // 2. 解説にある公式を計算. FAC[0] = 1; repx(i, 1, LIMIT) FAC[i] = (LL)i * FAC[i - 1] % MOD; rep(i, LIMIT) INV[i] = mPow(FAC[i], MOD - 2); MNFAC[0] = (M + N) % MOD; repx(i, 1, LIMIT) MNFAC[i] = (M + N + MOD - i) % MOD * MNFAC[i - 1] % MOD; // 3. 出力. int idx = (int)(S + N - 1); LL ans = MNFAC[idx] * INV[idx + 1] % MOD; 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 |
[入力例] 3 5 1 2 1 [出力例] 8 ※AtCoderのテストケースより [入力例] 10 998244353 31 41 59 26 53 58 97 93 23 84 [出力例] 642612171 ※AtCoderのテストケースより [入力例] 12 20201207 785 700 1233 872 580 739 842 520 963 202 390 840 [出力例] 76009999 |