POJ 2891 中国剩余定理扩展 模版

POJ 2891   扩展欧几里得 扩展中国剩余定理 两篇写的很好

x=a1∗x1+b1
x=a2∗x2+b2
a1,a2是模数，b1,b2是余数。 求x, 模数不互质。

#include<iostream>
using namespace std;
#define ll long long
const int maxn=1e5+5;
int n;
ll exgcd(ll a,ll b,ll &x,ll &y) {
	if(!b) {//解方程求 x,y a*x+b*y = gcd(a,b)
		x=1,y=0;
		return a;
	}
	ll re=exgcd(b,a%b,x,y),tmp=x;
	x=y,y=tmp-(a/b)*y;
	return re;
}
ll m[maxn],a[maxn];
ll work() {
	ll M=m[1],A=a[1],t,d,x,y;
	int i;
	for(i=2; i<=n; i++) {
		d=exgcd(M,m[i],x,y);//解方程
		if((a[i]-A)%d)return -1;//无解
		x*=(a[i]-A)/d,t=m[i]/d,x=(x%t+t)%t;//求x
		A=M*x+A,M=M/d*m[i],A%=M;//日常膜一膜 防爆
	}
	A=(A%M+M)%M;
	return A;
}
int main() {
    int i,j;
    while(scanf("%d",&n)!=EOF){
        for(i=1;i<=n;i++)scanf("%lld%lld",&m[i],&a[i]);// %m[i] 余 a[i]
        printf("%lld\n",work());
    }

	return 0;
}