C++の練習を兼ねて, AtCoder Beginner Contest 156 の 問題F (F – Modularness) を解いてみた.
■感想.
1. 方針が全く見えなかったので, 解説を参照して, 実装した.
2. 何とか実装したものの, 結構苦労した気がする(※本質が理解出来てない気もする).
3. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイトABC 156解説をご覧下さい.
■C++版プログラム(問題F/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 |
// 解き直し. // https://atcoder.jp/contests/abc156/tasks/abc156_f #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--) LL d[5050], dMod[5050]; int main(){ // 1. 入力情報. int k, q; scanf("%d %d", &k, &q); rep(i, k) scanf("%lld", &d[i]); // 2. クエリに回答(解説通り). rep(i, q){ LL n, x, m, s = 0, t = 0; scanf("%lld %lld %lld", &n, &x, &m); // 2-1. d の 要素 を m で 割った余りで置き換える. rep(j, k) dMod[j] = d[j] % m; // 2-2. (a[i] mod m) == (a[i + 1] mod m) となる場合. rep(j, k){ // if(dMod[j] % k == 0){ // 誤答となるっぽい. if(dMod[j] == 0){ LL q = (n - 1) / k; LL r = (n - 1 - j) - q * k; s += (q + (r > 0)); } } // 2-3. (a[i] mod m) > (a[i + 1] mod m) となる場合. // a[n - 1] を m で 割った商 から a[0] を m で 割った商 を 引く. // ex. // k = 5, n = 12 // i = 0 1 2 3 4 5 6 7 8 9 10 11 (= n - 1) // a[0] = x // a[1] = x + 1 * d[0] // a[2] = x + 1 * d[0] + 1 * d[1] // a[3] = x + 1 * d[0] + 1 * d[1] + 1 * d[2] // a[4] = x + 1 * d[0] + 1 * d[1] + 1 * d[2] + 1 * d[3] // a[5] = x + 1 * d[0] + 1 * d[1] + 1 * d[2] + 1 * d[3] + 1 * d[4] // a[6] = x + 2 * d[0] + 1 * d[1] + 1 * d[2] + 1 * d[3] + 1 * d[4] // a[7] = x + 2 * d[0] + 2 * d[1] + 1 * d[2] + 1 * d[3] + 1 * d[4] // a[8] = x + 2 * d[0] + 2 * d[1] + 2 * d[2] + 1 * d[3] + 1 * d[4] // a[9] = x + 2 * d[0] + 2 * d[1] + 2 * d[2] + 2 * d[3] + 1 * d[4] // a[10] = x + 2 * d[0] + 2 * d[1] + 2 * d[2] + 2 * d[3] + 2 * d[4] // a[11] = x + 3 * d[0] + 2 * d[1] + 2 * d[2] + 2 * d[3] + 2 * d[4] t = x; rep(j, k){ LL q = (n - 1) / k; LL r = (n - 1 - j) - q * k; LL u = (q + (r > 0)) * dMod[j]; t += u; } // printf("s=%lld t=%lld n=%lld\n", s, t, n); t /= m; // (a[n - 1] / m) を 計算. t -= (x / m); // (a[0] / m) を マイナス. // 2-4. (n - 1) から, 前項までの結果をマイナス. LL ans = (n - 1) - (s + t); 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 |
[入力例] 3 1 3 1 4 5 3 2 [出力例] 1 ※AtCoderテストケースより [入力例] 7 3 27 18 28 18 28 46 1000000000 1000000000 1 7 1000000000 2 10 1000000000 3 12 [出力例] 224489796 214285714 559523809 ※AtCoderテストケースより [入力例] 10 5 120 20 56 903 15 95 94 285 3997 3816151 46008408 19 1182 936523469 6 958 555662761 22 703 776784387 12 1198 353200718 4 152 828319977 25 1207 181860716 17 1118 291738659 26 192 762764016 29 1018 939777463 28 501 [出力例] 35522228 722433031 426508714 613906058 176600359 |
■参照サイト
AtCoder Beginner Contest 156