中国剩余定理

定理内容

中国剩余定理最早记载于《孙子算经》的“物不知数”问题，所以又叫做孙子剩余定理。
设$m_1,m_2,m_3,\cdots ,m_k$是$k$个两两互素的正整数，则对于任意的$b_1,b_2,b_3,\cdots ,b_k$，同余式组：
$\begin{array}{r}\{\begin{array}{l}x\equiv b_1\mathrm{m}\mathrm{o}\mathrm{d}{m}_1\\ x\equiv b_2\mathrm{m}\mathrm{o}\mathrm{d}{m}_2\\ x\equiv b_3\mathrm{m}\mathrm{o}\mathrm{d}{m}_3\\ ⋮\\ x\equiv b_k\mathrm{m}\mathrm{o}\mathrm{d}{m}_k\end{array}\end{array}$

有且仅有一个解，为：
$$x^{\prime }\equiv b_1\cdot M_1^{\prime }\cdot M+b_2\cdot M_2^{\prime }\cdot M_2+\cdots +b_k\cdot M_k^{\prime }\cdot M_k\mathrm{m}\mathrm{o}\mathrm{d}m$$
其中$M_i=\prod _{j=1}^k{m}_j(j\ne i),M_i^{\prime }\cdot M_i\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}{m}_i$.
对于任意一个同余式$x\equiv b_i\mathrm{m}\mathrm{o}\mathrm{d}{m}_i$，我们知道，对于$M_j(j\ne i)$是$m_i$的倍数，所以，$b_j\cdot M_j^{\prime }\cdot M_j\equiv 0\mathrm{m}\mathrm{o}\mathrm{d}{m}_i.$
$\begin{array}{rl}{x}^{\prime }& \equiv b_1\cdot M_1^{\prime }\cdot M_1+b_2\cdot M_2^{\prime }\cdot M_2+\cdots +b_k\cdot M_k^{\prime }\cdot M_k\mathrm{m}\mathrm{o}\mathrm{d}m\\ & \equiv b_i\cdot M_i^{\prime }\cdot M_i\mathrm{m}\mathrm{o}\mathrm{d}{m}_i\\ & \equiv b_i\mathrm{m}\mathrm{o}\mathrm{d}{m}_i\end{array}$
这说明$x^{\prime }$的确是同余式组的解，显然，若$x^{\prime \prime }$是同于式组的解，那么也必然满足上式，说明$x^{\prime }\equiv x^{\prime \prime }\mathrm{m}\mathrm{o}\mathrm{d}m$，即同余式组的解式唯一的。

c++代码演示

#include<iostream>
#include<vector>
using namespace std;
//We need the extedned gcd algorithm to compute the inverse of M
vector<int> exgcd(int a,int b,int s1=1,int s2=0,int t1=0,int t2=1){
	int r=a%b;
	if(r==0){
		return {b,s2,t2};
	}
	int q=a/b;
	int s3=-q*s2+s1;
	int t3=-q*t2+t1;
	return exgcd(b,r,s2,s3,t2,t3);
}

//Warnning: the elements in vector ms must be coprime!
int crt(vector<int> ms,vector<int> bs){
	vector<int> Ms(ms.size());
	vector<int> Ms_prime(ms.size());
	int m=1;
	for(auto i:ms){
		m=m*i;
	}
	for(int i=0;i<ms.size();i++){
		Ms[i]=m/ms[i];
		auto temp=exgcd(Ms[i],ms[i]);
		Ms_prime[i]=temp[1];
		cout<<temp[1]<<endl;
	}
	int res=0;
	for(int i=0;i<ms.size();i++){
		res+=bs[i]*Ms[i]*Ms_prime[i]%m;
	}
	return res;
}

int main(){
	vector<int> ms={3,5,7};
	vector<int> bs={2,4,6};
	int res=crt(ms,bs);
	cout<<res<<endl;
}