假设已经求出前面个方程的通解
,
。
不难得出通解为
现在加入第个方程,等价于求解
通过扩展欧几里得求出方程的解,如果方程无解,则整个方程无解,如果有解
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 7;
typedef long long ll;
ll a[maxn], b[maxn];
int n;
ll exgcd(ll a, ll b, ll &x, ll &y) // 扩展欧几里得
{
if(b == 0) {
x = 1;
y = 0;
return a;
}
ll d = exgcd(b, a % b, x, y);
ll t = x;
x = y;
y = t - a / b * y;
return d;
}
ll quickmul(ll a, ll b, ll mod) // 快速乘
{
a = (a % mod + mod) % mod;
b = (b % mod + mod) % mod;
ll ans = 0;
while(b) {
if(b & 1) ans = (ans + a) % mod;
a = (a + a) % mod;
b >>= 1;
}
return ans;
}
ll EXCRT() // 扩展中国剩余定理
{
ll M, x, y, ans;
M = b[0], ans = a[0];
for(int i = 1; i < n; i++) { // n - 1 次运算
// t * M = ak - x mod bk
// x * aa = cc mod bb
ll aa = M, bb = b[i], cc = (a[i] - ans % bb + bb) % bb; // aax = cc mod bb
ll d = exgcd(aa, bb, x, y), g = bb / d; //
if(cc % d != 0) return - 1; // 无解
x = quickmul(x, cc / d, g);
ans = ans + x * M; // x = x + t * M, M 前 k - 1个 b[i] 的 lcm
M *= g; // 前 k 个 b[i] 的 lcm
ans = (ans % M + M) % M;
}
return ans;
}
int main()
{
while(scanf("%d", &n) == 1) {
for(int i = 0; i < n; i++) scanf("%lld %lld", &b[i], &a[i]);
ll ans = EXCRT();
printf("%lld\n", ans);
}
return 0;
}