結果

問題 No.186 中華風 (Easy)
ユーザー とばり
提出日時 2018-09-21 06:48:14
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,440 bytes
コンパイル時間 1,836 ms
コンパイル使用メモリ 177,436 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 18:32:22
合計ジャッジ時間 2,377 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using int64 = long long;
using uint64 = unsigned long long;

inline int64 mod(int64 a, int64 m) { return (a % m + m) % m; }

std::tuple<int64, int64, int64> extGCD(int64 a, int64 b) {
    if (b == 0) return std::make_tuple((int64)1, (int64)0, a);

    int64 s, t, g;
    std::tie(s, t, g) = extGCD(b, a % b);
    return std::make_tuple(t, s - (a / b) * t, g);
}

std::tuple<int64, bool> CRT(std::vector<int64> b, std::vector<int64> m) {
    int n = b.size();

    int64 B = 0, M = 1;
    for (int i = 0; i < n; i++) {
        int64 p, q, d;
        std::tie(p, q, d) = extGCD(M, m[i]);

        if ((b[i] - B) % d != 0) return std::make_tuple(0, false);

        int64 tmp = mod((b[i] - B) / d * p, m[i] / d);
        B += M * tmp;
        M *= m[i] / d;
        B = mod(B, M);
    }
    return std::make_tuple(B, true);
}

int64 GCD(int64 a, int64 b) {
    if (b == 0) return a;
    else return GCD(b, a % b);
}

int64 LCM(int64 a, int64 b) {
    return (a / GCD(a, b)) * b;
}

int main() {
    std::cin.tie(nullptr);
    std::ios::sync_with_stdio(false);

    std::vector<int64> X(3), Y(3);
    for (int i = 0; i < 3; i++) {
        int64 x, y;
        std::cin >> x >> y;

        X[i] = x; Y[i] = y;
    }

    int64 x;
    bool ok;
    std::tie(x, ok) = CRT(X, Y);

    if (x == 0) {
        x = LCM(Y[0], LCM(Y[1], Y[2]));
    }

    if (!ok) x = -1;
    std::cout << (int64)x << std::endl;

    return 0;
}
0