excrt

求一个形如\(x\equiv a_i\ (mod\ \ p_i)\) 的线性方程组的解
 
我们考虑合并两个方程：

先考虑将\(x\)设为\(a_1\)，再不断的加\(p_1\)的倍数使得\(x\)在模\(p_2\)意义下为\(a_2\)

我们可以用exgcd求得一组\(mp_1+np_2= gcd(p_1,p_2)\)的解，即\(mp_1\equiv gcd(p1,p2)\ (mod\ \ p2)\)

也就是说我们每次增加\(mp_1\)即可让\(x\)在模\(p_1\)意义下增加\(gcd(p_1,p_2)\)

令\(d=gcd(p_1,p_2)\)，那么我们应该增加\((a_2-a_1)/d\) 次，即\(x\)应该加上\((a_2-a_1)/d*mp_1\)，然后我们知道方程合并后的模数应该是\(P=p_1*p_2/gcd(p_1,p_2)\)，所以再将结果对\(P\)取模即可。
 
代码：

LL exgcd(LL x,LL y,LL&a,LL&b){
	LL z;
	return y?(z=exgcd(y,x%y,b,a),b-=x/y*a,z):(a=1,b=0,x);
}
LL mult(LL a,LL b,LL p){
	LL ans=0;
	while(b){
		if(b&1)ans=(ans+a)%p;
		b>>=1,a=(a<<1)%p;
	}
	return ans;
}
struct CRT{
	LL a[N],p[N],P,x;
	bool excrt(LL a1,LL p1,LL a2,LL p2){
		LL m,n,d=exgcd(p1,p2,m,n); P=p1/d*p2;
		LL a=(a2-a1%p2+p2)%p2;
		if(a%d!=0)return false;
		x=(x+mult(a/d*p1,m,P))%P; //注意这里一般先计算a/d*p1，然后再乘m，因为题目一般保证了所有p_i的最小公倍数小于某个值，这样前面的运算就不会溢出，而后面的乘m可能会溢出，所以还要用快速乘。
		return true;
	}
	bool run(){
		x=a[1],P=p[1];
		for(int i=2;i<=n;i++)
			if(!excrt(x,P,a[i],p[i]))
				return false;
		return true;
	}
}crt;

原文地址：https://www.cnblogs.com/lishuyu2003/p/13287249.html