扩展中国剩余定理

问题描述

 

求解同余方程组

解释思路

假设已经求出前k-1个方程组成的同余方程组的一个解为
M为前k-1个方程的m的最小公倍数 
则前k-1个方程的方程组通解为+iM(i∈Z)
那么对于加入第k个方程后的方程组
我们就是要求一个正整数t，使得
+tM == a[k] (mod m[k]) 
变形得Mt==(a[k]-)(mod m[k]) 
即求解不定方程 Mt + m[k]y = a[k] - 的t 
= tM+

 代码

ll m[N];
ll a[N];
int  n;
ll exgcd(ll a,ll b,ll &x,ll &y){//这一层为x1,y1 
	if(b == 0){
		x = 1; y = 0;
		return a;
	}
	ll d = exgcd(b, a%b, x, y);//这一层为x2,y2,下面赋值语句后面的都是修改过的值 
	ll xt= x;//先存下x2的值，后面会修改 
	x = y;//x1=y2
	y = xt - a / b * y;//y1=x2-a/b*y2
	//此时修改后得到本层的x,y的值 
	return d;
}
ll excrt(){
	ll x,y,t;
	ll M=m[1],ans=a[1];
	for(int k=2;k<=n;k++){//求前k个方程的解ans
		ll gcd=exgcd(M,m[k],x,y);
		if((a[k]-ans)%gcd)return -1;//无解
		t=(a[k]-ans)/gcd*x;
		ans=ans+t*M; 
		M=M*m[k]/gcd;//前k个m的最小公倍数 
		ans=(ans%M+M)%M;//防止溢出 
	}
	return ans;
}