C++の練習を兼ねて, AtCoder Beginner Contest 177 の 問題E (Coprime) を解いてみた.
■感想.
1. 途中ロジック誤りがあったものの, AC版に到達できたので良かったと思う.
2. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイトABC 177 解説をご覧下さい.
■C++版プログラム(問題E/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 |
// C++(GCC 9.2.1) #include <bits/stdc++.h> using namespace std; #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--) string ans[] = {"not coprime", "setwise coprime", "pairwise coprime"}; int a[1010101], b[1010101], c[1010101]; set<int> divisors[1010101]; int main() { // 1. 入力情報取得. int N; scanf("%d", &N); rep(i, N) scanf("%d", &a[i]); // 2. 素因数を取得. repx(i, 2, 1010101){ if(b[i]) continue; // max_12.txt, max_13.txt で, WA版となったことで, // j が i で 割り切れる場合の考慮が漏れていたことが判明し, ロジック修正. // repex(j, i * 2, 1010101, i){ repex(j, i, 1010101, i){ b[j]++; if(j % i == 0) divisors[j].insert(i); } } // 3. "pairwise coprime" であるかの判定. int idx = 2; rep(i, N){ bool ok = false; for(auto &p : divisors[a[i]]){ // 前回以前に出てきた素因数が見つかったら, // "pairwise coprime" でないことが分かるので, 終了. if(!c[p]){ c[p]++; }else{ idx--; ok = !ok; break; } } if(ok) break; } // 4. GCD(A[1], ..., A[N]) の 計算. // -> gcd = 1 の 場合は, 少なくとも, "setwise coprime" になる. int gcd = a[0]; repx(i, 1, N) gcd = __gcd(gcd, a[i]); if(gcd != 1) idx--; // 5. 出力. printf("%s\n", ans[idx].c_str()); 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 47 48 49 50 51 52 53 54 55 56 57 58 |
[入力例] 3 3 4 5 [出力例] pairwise coprime ※AtCoderテストケースより [入力例] 3 6 10 15 [出力例] setwise coprime ※AtCoderテストケースより [入力例] 3 6 10 16 [出力例] not coprime ※AtCoderテストケースより [入力例] 3 6 6 9 [出力例] not coprime [入力例] 3 2 3 6 [出力例] setwise coprime [入力例] 5 1 2 1 3 1 [出力例] pairwise coprime [入力例] 10 2 3 5 7 11 13 17 19 23 29 [出力例] pairwise coprime [入力例] 17 2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 57 [出力例] setwise coprime |
■参照サイト
AtCoder Beginner Contest 177