結果

問題 No.186 中華風 (Easy)
ユーザー dyktr_06
提出日時 2022-10-19 13:45:24
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,391 bytes
コンパイル時間 1,616 ms
コンパイル使用メモリ 170,476 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-29 16:43:23
合計ジャッジ時間 2,533 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

// CRT::extGCD(a, b, x, y) : ax + by = gcd(a, b) となる x, y を求めます : O(log a)
// CRT::ChineseRem(b, m) : 中国剰余定理、返り値を (r, m) とすると解は x ≡ r (mod. m)、解なしの場合は (0, -1) を返します
namespace CRT{
    long long extGCD(long long a, long long b, long long &x, long long &y){
        if(b == 0){
            x = 1;
            y = 0;
            return a;
        }
        long long d = extGCD(b, a % b, y, x);
        y -= a / b * x;
        return d;
    }

    pair<long long, long long> ChineseRem(const vector<long long> &b, const vector<long long> &m) {
        long long r = 0, M = 1;
        for(int i = 0; i < (int) b.size(); i++){
            long long p, q;
            long long d = extGCD(M, m[i], p, q);
            if((b[i] - r) % d != 0) return {0, -1};
            long long tmp = (b[i] - r) / d * p % (m[i] / d);
            r += M * tmp;
            M *= m[i] / d;
        }
        r %= M;
        if(r < 0) r += M;
        return {r, M};
    }
}

// example (中華風 (Easy) (yukicoder))
int main(){
    vector<long long> b(3), m(3);
    for(int i = 0; i < 3; i++){
        cin >> b[i] >> m[i];
    }
    pair<long long, long long> p = CRT::ChineseRem(b, m);
    if(p.second == -1){
        cout << -1 << endl;
    }else{
        cout << p.first << endl;
    }
}
0