CRT exCRT学习笔记

CRT

中国剩余定理用于求解这样一组同余式的最小正整数解：

$\{\begin{array}{l}x\equiv a_1(\mathrm{mod}{m}_1)\\ x\equiv a_2(\mathrm{mod}{m}_2)\\ \dots .\dots \\ x\equiv a_n(\mathrm{mod}{m}_n)\end{array}$

其中 $a_1,a_2,\dots \dots a_n$ 两两互质。

设 $M={\prod }_{i=1}^n$ $m_i$，$M_i=M÷m_i$ （原谅我的变量重名但只有这样看着顺眼qwq

构造逆元 $t_i$ 使得 $M_i{t}_i\equiv 1(\mathrm{mod}{m}_i)$，那么当前方程显然就有解 $M_i{t}_i{a}_i$，所以方程组就有解 $x={\sum }_{i=1}^n$ $M_i{t}_i{a}_i$，这是由于 $M_i$ 的存在使得第 $i$ 个方程的解对其余方程均无影响，所以加起来一定正确。

能求出一个解，则其它解均为 $x+k\times M$的形式，最小正整数解就也可求了。

代码实现的话，$t_i$ 用扩展欧几里得求解，注意乘法可能爆long long，要龟速乘就好了。

洛谷模板

#include<bits/stdc++.h>
#define int long long
#define ff(i,s,e) for(int i=(s);i<=(e);++i)
using namespace std;
inline int read(){
	int x=0,f=1;
	char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return x*f;
}
const int N=15;
int n,a[N],b[N],m[N];
inline void exgcd(int a,int b,int &x,int &y){
	if(!b) return x=1,y=0,void();
	exgcd(b,a%b,x,y);
	int t=x;
	x=y,y=t-a/b*y;
}
inline int mul(int x,int y,int mod){
	int res=0;
	for(;y;y>>=1){
		if(y&1) res=(res+x)%mod;
		x=(x<<1)%mod;
	}
	return res;
}
signed main(){
	n=read();
	int M=1;
	ff(i,1,n) a[i]=read(),b[i]=read(),M*=a[i];
	ff(i,1,n) m[i]=M/a[i];
	int ans=0;
	ff(i,1,n){
		int x,y;
		exgcd(m[i],a[i],x,y);
		x=(x%M+M)%M;
		ans=(ans%M+mul(m[i],mul(x,b[i],M),M)+M)%M;
	}
	printf("%lld",ans);
}

exCRT

中国剩余定理的求解基于求逆元，而当任意两个 $m_i$ 不互质的时候，$M_i$ 在 $m_i$ 的剩余系下的逆元不存在，所以 $m_i$ 不两两互质的情况下无法用普通CRT求解。

于是考虑另一种思路：假设已经求出前 $i-1$ 组方程的解 $x$，${\mathrm{lcm}}_{j=1}^{i-1}$$m_i=t$，则通解为 $x+k\times t$，那么加入第 $i$ 个方程就是要求解 $k$ 使得 $x+k\times t\equiv a_i(\mathrm{mod}{m}_i)$，即 $k\times t\equiv a_i-x(\mathrm{mod}{m}_i)$，扩欧求解即可，若有任意一个同余式无解，则整个方程组无解。

洛谷模板++

一定要注意好当前是在谁的剩余系里啊/kk

#include<bits/stdc++.h>
#define int long long
#define ff(i,s,e) for(int i=(s);i<=(e);++i)
using namespace std;
inline int read(){
	int x=0,f=1;
	char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return x*f;
}
const int N=1e5+5;
int n,a[N],b[N];
inline int lcm(int x,int y){
	return x/__gcd(x,y)*y;
}
inline void exgcd(int a,int b,int &x,int &y){
	if(!b) return x=1,y=0,void();
	exgcd(b,a%b,x,y);
	int t=x;
	x=y,y=t-a/b*y;
}
inline int mul(int x,int y,int mod){
	int res=0;
	for(;y;y>>=1){
		if(y&1) res=(res+x)%mod;
		x=(x+x)%mod;
	}
	return res;
}
signed main(){
	n=read();
	ff(i,1,n) a[i]=read(),b[i]=read();
	int ans=b[1],m=a[1];
	ff(i,2,n){
		int x,y;
		int g=__gcd(m,a[i]),c=(b[i]+a[i]-ans%a[i])%a[i];
		exgcd(m,a[i],x,y);
		x=mul(x,c/g,a[i]/g);
		ans+=x*m;
		m*=a[i]/g;
		ans=(ans%m+m)%m;
	}
	printf("%lld",ans);
	return 0;
}