結果
問題 |
No.3020 ユークリッドの互除法・改
|
ユーザー |
|
提出日時 | 2025-02-14 22:28:43 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,633 bytes |
コンパイル時間 | 3,998 ms |
コンパイル使用メモリ | 283,740 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2025-02-14 22:28:48 |
合計ジャッジ時間 | 4,233 ms |
ジャッジサーバーID (参考情報) |
judge7 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 16 WA * 5 |
ソースコード
#include<bits/stdc++.h> using namespace std; int h = 2, w = 2; vector<vector<long long>> gcd(vector<vector<long long>> a){ if(a[0][1] == 0 || a[0][0] == 0) return a; if(abs(a[0][0]) > abs(a[0][1])){ swap(a[0][0], a[0][1]); swap(a[1][0], a[1][1]); } if(a[0][0] < 0){ a[0][0] *= -1; a[0][1] *= -1; } long long k = a[0][1]/a[0][0]; a[0][1] -= k*a[0][0]; if(a[0][1] < 0){ k--; a[0][1] += a[0][0]; } a[1][1] -= k*a[1][0]; return gcd(a); } vector<vector<long long>> rotate(vector<vector<long long>> a){ auto b = a; for(int i=0; i<h; i++){ for(int j=0; j<w; j++){ b[i][j] = a[j][w-i-1]; } } return b; } int main(){ vector<vector<long long>> a(h, vector<long long>(w)); for(int i=0; i<h; i++){ for(int j=0; j<w; j++){ cin >> a[i][j]; } } for(int i=0; i<4; i++){ a = gcd(a); a = rotate(a); } for(int i=0; i<4; i++){ if(a[0][0] == 0) continue; if(a[0][1]%a[0][0] == 0) a[0][1] = 0; if(a[1][0]%a[0][0] == 0) a[1][0] = 0; } long long mi = -1; for(int i=0; i<h; i++){ for(int j=0; j<w; j++){ if(a[i][j] != 0){ if(mi == -1) mi = abs(a[i][j]); if(abs(a[i][j]) < mi) mi = abs(a[i][j]); } } } if(mi == -1){ cout << "0 0" << endl; return 0; } while(abs(a[0][0]) != mi) a = rotate(a); { long long k = a[0][1]/a[0][0]; a[1][1] -= a[1][0]*k; a[0][1] = a[1][0] = 0; } long long ans1 = abs(a[0][0]); long long ans2 = abs(a[0][0]); for(int i=0; i<h; i++){ for(int j=0; j<w; j++){ if(a[i][j] != 0){ ans1 = __gcd(ans1, abs(a[i][j])); ans2 /= __gcd(ans2, abs(a[i][j])); ans2 *= abs(a[i][j]); } } } cout << ans1 << " " << ans2 << endl; }