社区讨论

84pt #15炸了 求调

P4777【模板】扩展中国剩余定理(EXCRT)参与者 1已保存回复 0

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
0 条
当前快照
1 份
快照标识符
@lo3515lx
此快照首次捕获于
2023/10/24 00:54
2 年前
此快照最后确认于
2023/10/24 00:54
2 年前
查看原帖
开了__int128还是过不去
CPP
#include <iostream>
using namespace std;

#define int long long
int exgcd(int a, int b, int &x, int &y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    int r = exgcd(b, a % b, x, y);
    int t = x;
    x = y;
    y = t - a / b * y;
    return r;
}
int a[100005], m[100005];

int exCRT(int a[], int m[], int n) {
    int M = m[1], A = a[1], x, y, d;
    for (int i = 2; i <= n; i++ ) {
        d = exgcd(M, m[i], x, y);
        if ((A - a[i]) % d != 0) {
            return -1;
        }
        x = (__int128)((A - a[i]) / d * x % (m[i] / d) + m[i] / d) % (m[i] / d);
         A -= (__int128)x * M;
        M = (__int128)M / d * m[i];
        A = (A % M + M) % M;
    }
    return A;
}
signed main() {
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++) {
        cin >> m[i] >> a[i];
    }
    cout << exCRT(a, m, n) << endl;
    return 0;
}

回复

0 条回复,欢迎继续交流。

正在加载回复...