中国剩余定理

中国剩余定理

模数互质

定理：若m₁,m₂,…,m_n是两两互质的正整数，则对于任意的n个整数a₁,a₂,…,a_n,同余方程组x≡a_i(mod m_i)(i=1,2,3,…,n)有整数解，并且在模M下解唯一。构造方程组的解为：x = a₁$M_1$$x_1$+$a_2{M}_2{x}_2$+…+$a_n{M}_n{x}_n$,其中：M = $\prod m_i,M_i=M/m_i$,$x_i$是线性同余方程，$M_i{x}_i$ ≡ 1(mod $m_i$)的一个解(因为$M_i和m_i$)互质，必有解）。

int ex_gcd(int a,int b,int &x,int &y){
    if(b==0) { x=1,y=0; return a}
    int d = ex_gcd(b,a%b,x,y);
    int temp = x;
    x = y;
    y = temp - a/b*y;
    return d;
}
int chinese_remain(int n){
    int M=1,x,y,ans=0;
    for(int i=0;i<n;i++) M*=mi[i];
    for(int i=0;i<n;i++){
        int MI = M/mi[i];
        ex_gcd(MI,mi[i],x,y);
        ans = (ans+MI*x*ai[i])%M;
    }
    return (ans%M+M)%M;
}

模数不互质

$$\{\begin{array}{l}x=a_1(mod{n}_1)\\ x=a_2(mod{n}_2)\end{array}\to \{\begin{array}{l}x=k_1{n}_1+a_1\\ x=k_2{n}_2+a_2\end{array}\to k_1{n}_1=(a_2-a_1)(modn2)$$

根据线性同余方程相关定理可得：当且仅当 $gcd(n_1,n_2)|a2-a1$，方程有解。

令 d = $gcd(n1,n2)$ c = $a_2-a_1$

$k_1{n}_1=(a_2-a_1)(mod{n}_2)\to \frac{k_1{n}_1}{d}=\frac{a_2-a_1}{d}(mod\frac{n_2}{d})\to k_1=\frac{(a2-a1)}{d}\ast inv(\frac{n_1}{d})(mod\frac{n_2}{d})$

int chinese_remain(int n){
	int a1=ai[0],k1=mi[0],x,y,ans=0;
	for(int i=1;i<n;i++){
		int d = ex_gcd(k1,mi[i],x,y);
		int c = ai[i]-a1;
		if(c%d){// 当且仅当gcd(n1,n2)|a2-a1,方程有解
			return -1;
		}
		int temp = c/d*x%(mi[i]/d);  // 由费小马定理可得求出的x为n1/d的逆元
		ans += k1*temp;
		k1 = (k1*mi[i])/d;
		ans %= k1;
	}
	return (ans%k1+k1)%k1;
}