【Luogu 4777】扩展中国剩余定理

假设已经求出前面个方程的通解  , 。

不难得出通解为 

现在加入第个方程，等价于求解

通过扩展欧几里得求出方程的解，如果方程无解，则整个方程无解，如果有解 

#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;
}