中国剩余定理

含义：给定两两互质的正整数n1,n2,...,nk,要求找到最小的正整

数x,满足方程组x≡ai(mod ni) (i=1,2...k)

求法：

– 令n=n1，n2，...，nk, mi = n/ni
– 显然gcd(mi,ni)=1,利用扩展欧几里德算法计算出xi满
足mi*xi≡1(mod ni)
– x = (a1*x1*m1+a2*x2*m2+...+ak*xk*mk) mod n
– 此方程组任意两个解模 n 同余，因此x就是最小的

代码如下：

#include<iostream>
using namespace std;
int gcdEx(int a,int b,int& x,int& y)
{   //求ax+by=gcd(a,b) 的整数解 返回gcd(a,b); 
	if(b==0) {
		x=1; y=0; return a;
	}
	int x1,y1;
	int gcd=gcdEx(b,a%b,x1,y1);
	x=y1;
	y=x1-a/b*y1;
	return gcd;
}
int CRT(int* a,int* n,int k) 
{   //a[i]表示第i个余数，n[i]表示第i个除数，一共有k个数，返回x值 
	int n0=1;
	for(int i=0;i<k;i++) n0*=n[i]; //n0=n1*n2....*nk;
	int ans=0,x,y,m;
	for(int i=0;i<k;i++)
	{
		m=n0/n[i];
		gcdEx(m,n[i],x,y);
		ans=(ans+a[i]*x*m)%n0;
	}
	if(ans<=0) ans+=n0;
	return ans;
}