【初等数论】中国剩余定理及其扩展

原创已于 2025-12-12 21:05:06 修改 · 992 阅读

27 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

标签

#算法

于 2025-12-07 10:51:51 首次发布

中国剩余定理

中国剩余定理（CRT）
扩展中国剩余定理（EXCRT）

中国剩余定理（CRT）

“今有物不知其数，三三数之剩二（除以3余2），五五数之剩三（除以5余3），七七数之剩二（除以7余2），问物几何？” ——《孙子算经》

定理内容

一次同余方程组 :

$\left\{ \begin{array}{lcl} x≡a_{1} (mod\space m_{1})\\ x≡a_{2} (mod\space m_{2}) \\ \cdots \\ x≡a_{n} (mod\space m_{n}) \end{array} \right.$

如果满足 $m_{i}$ 两两互质，则在 $mod\space M(M=m_{1}\cdot m_{2}\cdot m_{3}\cdots m_{n})$ 意义下，存在唯一解.

定理证明

构造解 $M=m_{1}\cdot m_{2}\cdot m_{3}\cdots m_{n}$ ， $M_{i}=M/m_{i}$ ， $t_{i}$ 为 $M_{i}$ 在 $mod\space m_{i}$ 意义下的乘法逆元， $gcd(M_{i},m_{i})=1$ .

则 $x=(\sum\limits_{i=1}^na_{i}t_{i}M_{i})$ $(mod\space M)$

带入验证，对于每一个一次同于方程 $x≡a_{j} (mod\space m_{j})$ .

当 $i \neq = j$ 时，一定有 $m_{j}\space |\space a_{i}t_{i}M_{i}$ ， $a_{i}t_{i}M_{i}\space mod\space m_{j}=0$ .

当 $i = j$ 时，一定有 $a_{i}t_{i}M_{i}=a_{j}$ ， $a_{i}t_{i}M_{i}\space mod\space m_{j}=a_{j}$ .

由此，同余方程组成立 .

典型例题

【模板】中国剩余定理（CRT）/ 曹冲养猪

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=15;
namespace zyx_{

ll read(){
	ll x=0;int f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
void print(ll x){
	if(x<0) putchar('-'),x=-x;
	if(x>9) print(x/10);
	putchar('0'+x%10);
	return ;
}
ll mul(ll a,ll b,ll mod){
	ll res=0;
	a%=mod;
	while(b){
		if(b&1) res=(res+a)%mod;
		a=(a+a)%mod;
		b>>=1;
	}
	return res;
}
ll exgcd(ll a,ll b,ll &x,ll &y){
	if(!b){
		x=1,y=0;
		return a;
	}
	ll d=exgcd(b,a%b,y,x);
	y-=(a/b)*x;
	return d;
}
ll inv(ll a,ll b){
	ll x,y;
	exgcd(a,b,x,y);
	return (x%b+b)%b;
}
ll n,M=1,ans,a[N],b[N];
int my_main(){
	n=read();
	for(int i=1;i<=n;i++)
		a[i]=read(),b[i]=read(),M*=a[i];
	for(int i=1;i<=n;i++){
		ll Mi=M/a[i];
		ll tmp=mul(b[i],Mi,M);
		tmp=mul(tmp,inv(Mi,a[i]),M);
		ans=(ans+tmp)%M;
	}
	print(ans);
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}

扩展中国剩余定理（EXCRT）

中国剩余定理的局限性

中国剩余定理(CRT)不能解决模数不互质情况的模线性同余方程组：

$\sum\limits_{i=1}^na_{i}\cdot M_{i}\cdot inv(M_{i},m_{i})\space (mod\space M)$

其中的 $inv(M_{i},m_{i})$ 在 $M_{i},m_{i}$ 不互质的情况下，根本不存在。而只要有任意两个模数不互质，就会破坏掉整个解。

解决方案

扩展中国剩余定理的核心思想就是将两个同余方程组两两合并：

对于两个同余方程：

$\left\{ \begin{array}{lcl} x≡a_{1} (mod\space m_{1})\\ x≡a_{2} (mod\space m_{2}) \end{array} \right.$

可以将其转化为不定方程组：

$\left\{ \begin{array}{lcl} x=k_{1}m_{1}+a_{1}\\ x=k_{2}m_{2}+a_{2} \end{array} \right.$

合并，移项后，得到新的不定方程：

$k_{1}m_{1}-k_{2}m_{2}=a_{2}-a_{1}$

由裴蜀定理得，当且仅当 $gcd(m_{1},m_{2})\space |\space (a_{2}-a_{1})$ 时，有解；

我们就能够根据扩展欧几里德算法计算出系数：

设 $d=gcd(m_{1},m_{2})$

将方程左右两边同时除以 $d$ ，得

$\frac{m_{1}}{d}k_{1}+\frac{m_{2}}{d}(-k_{2})=\frac{a_{2}-a_{1}}{d}$

为了简化方程，令：

$m_{1}^{'}=\frac{m_{1}}{d}k_{1},m_{2}^{'}=\frac{m_{2}}{d}k_{1},c=\frac{a_{2}-a_{1}}{d}$

方程变为：

$m_{1}^{'}k_{1}+m_{2}^{'}(-k_{2})=c$

其中 $gcd(m_{1}^{'},m_{2}^{'})=1$ ，根据裴蜀定理，一定有：

$m_{1}^{'}u+m_{2}^{'}v=1$

则通解：

$k_{1}=cu+tm_{2}^{'},k2=-(cv-tm_{1}^{'})\space (t\in Z)$

代回求 $x$ ：

$x=a_{1}+k_{1}m_{1}=a_{1}+(cu+tm_{2}^{'})m_{1}$

化简，得：

$x=a_{1}+cu\cdot m_{1}+t\cdot \frac{m_{1}m_{2}}{d}$

令 $x_{0}=a_{1}+cu\cdot m_{1}$ ，则合并结果为：

$x≡x_{0}\space (mod\space lcm(m_{1},m_{2}))$

典型例题

【模板】扩展中国剩余定理（EXCRT）

#include<bits/stdc++.h>
using namespace std;
typedef __int128 ll;
const int N=1e5+5;
namespace zyx_{

ll read(){
	ll x=0;int f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
void print(ll x){
	if(x<0) putchar('-'),x=-x;
	if(x>9) print(x/10);
	putchar(x%10+'0');
	return ;
}
ll n,A=1,B=0,x,y;
ll exgcd(ll a,ll b,ll &x,ll &y){
	if(!b){
		x=1,y=0;
		return a;
	}else{
		ll d=exgcd(b,a%b,y,x);
		y-=a/b*x;
		return d;
	}
}
int my_main(){
	n=read();
	for(int i=1;i<=n;i++){
		ll a,b;
		a=read(),b=read();
		ll d=exgcd(A,a,x,y);
		x=(B-b)/d*x;
		B=B-A*x;
		A=a/d*A;
		B=(B%A+A)%A;
	}
	ll ans=(B%A+A)%A;
	print(ans);
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}