社区讨论
84pt #15炸了 求调
P4777【模板】扩展中国剩余定理(EXCRT)参与者 1已保存回复 0
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 0 条
- 当前快照
- 1 份
- 快照标识符
- @lo3515lx
- 此快照首次捕获于
- 2023/10/24 00:54 2 年前
- 此快照最后确认于
- 2023/10/24 00:54 2 年前
开了
CPP__int128还是过不去#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 条回复,欢迎继续交流。
正在加载回复...