CRT
CRT 极简记录
CRT
CRT 原名为中国剩余定理,即 China Remainder Theorem。能够解决同于方程组的问题。
问题原型
解同余方程组:
$$ \left \{ \begin{aligned} x\equiv a_1 \bmod m_1\\ x\equiv a_2 \bmod m_2\\ x\equiv a_3 \bmod m_3\\ ...\\ x\equiv a_n \bmod m_n \end{aligned} \right. $$
其中满足所有 $m_i$ 互质。
设 $M=\prod m_i,M_i=\frac{M}{m_i},t_iM_i\equiv 1\bmod m_i$
则最小解 $x=\sum\limits_{i=1}^{n}M_it_ia_i\bmod M$
但注意普通的 CRT 必须满足 $m_i$ 互质,否则不成立,因为不互质不一定能够构造出来 $t_iM_i\equiv 1\bmod m_i$。
证明
证明是好理解的,首先,$\forall j\in [1,n],j\neq i$ 有 $M_jt_ja_j\equiv 0 \bmod m_i$,因为 $M_i$ 中有因子 $m_i$。之后又有 $M_it_ia_i\equiv a_i \bmod m_i$,因为 $t_iM_i\equiv 1\bmod m_i$。
这样加起来之后的数,一定能够保证其正确性。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=2e5+10;
void exgcd(ll a,ll b,ll &x,ll &y){
if(b==0){x=1,y=0;}
else exgcd(b,a%b,y,x),y-=a/b*x;
}
ll inv(ll x,ll mod){
ll q,p;
exgcd(x,mod,q,p);
return (q%mod+mod)%mod;
}
int n,a[maxn],m[maxn];
ll M=1;
int main(){
cin.tie(0);
ios::sync_with_stdio(false);
cin>>n;
for(int i=1;i<=n;i++)cin>>m[i]>>a[i],M*=1ll*m[i];
ll ans=0;
for(int i=1;i<=n;i++){
ll Mi=M/(1ll*m[i]);
ll t=inv(Mi,1ll*m[i]);
ans=(ans+Mi*t*1ll*a[i])%M;
}
cout<<ans<<"\n";
return 0;
}
EXCRT
那么有没有一个真正的同解,即没有 $m_i$ 限制的条件的解法呢?肯定有的,但是其逻辑和原来的 CRT 是完全不一样的,因为 CRT 其实只是简单地构造出来了一些等于 $0$ 和等于 $1$,的解,但是想要再进一步很困难,我们尝试回到两个同余方程组的本质。
$$ \left \{ \begin{aligned} x\equiv a_1 \bmod m_1\\ x\equiv a_2 \bmod m_2\\ \end{aligned} \right. $$
尝试将其合并,可以得到$a_1+k_1m_1=a_2+k_2m_2$,之后移项有 $k_1m_1-k_2m_2=a_2-a_1$。对于这个方程,可以简单用裴蜀定理来判断是否有解,即是否满足 $\gcd(m_1,m_2)|a_2-a_1$。
接着考虑如何解一个解出来,这里设 $\gcd(m_1,m_2)=d$ 则将两侧同时除以 $d$ 有 $k_1\dfrac{m_1}{d}-k_1\dfrac{m_2}{d}=\dfrac{a_2-a_1}{d}$。
左边两个系数是互质的,因此我们解出来 $X\dfrac{m_1}{d}+Y\dfrac{m_2}{d}=1$ 之后可以得到
$$ \left \{ \begin{aligned} k_1=X\dfrac{a_2-a_1}{d}\\ k_2=-Y\dfrac{a_2-a_1}{d}\\ \end{aligned} \right. $$
则有 $x=a_1+m_1X\dfrac{a_2-a_1}{d}$。
这样,将两个方程合并之后,我们有
$$ x \equiv a_1+m_1X\dfrac{a_2-a_1}{d} \bmod \text{lcm}(m_1,m_2) $$
至于证明唯一性,可以简单的设在 $\text{lcm}(m_1,m_2)$ 范围内有两个解,发现矛盾。
之后通过不断合并两个方程,我们就可以得到最后的解。
#include<bits/stdc++.h>
using namespace std;
#define ll __int128
inline ll read(){
ll res=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch)){
res=(res<<1)+(res<<3)+(ch-'0');
ch=getchar();
}
return res;
}
inline void pr(ll x){
if(x>=10)pr(x/10);
putchar(x%10+'0');
}
const int maxn=2e5+10;
int n;
ll a[maxn],b[maxn];
ll gcd(ll x,ll y){
if(y==0)return x;
return gcd(y,x%y);
}
void exgcd(ll a,ll b,ll &x,ll &y){
if(!b){x=1,y=0;}
else exgcd(b,a%b,y,x),y-=a/b*x;
}
ll lcm(ll x,ll y){
return x/gcd(x,y)*y;
}
int main(){
cin>>n;
for(int i=1;i<=n;i++)a[i]=read(),b[i]=read();
ll A=a[1],B=b[1];
for(int i=2;i<=n;i++){
ll d=gcd(A,a[i]),X=0,Y=0;
exgcd(A/d,a[i]/d,X,Y);
if(X<0)X+=a[i]/d;
else X%=a[i]/d;
ll LCM=lcm(A,a[i]);
B+=(b[i]-B)/d%LCM*A%LCM*X%LCM,A=LCM;
if(B<0)B+=A;
}
pr(B%A);
return 0;
}