【模板】中国剩余定理（CRT）

【模板】中国剩余定理

Part 1. 中国剩余定理是干什么的

用于求解如下问题：

给定 $n$ 个线性同余方程，如：
$$\{\begin{array}{l}x\equiv a_1(\mathrm{mod}{m}_1)\\ x\equiv a_2(\mathrm{mod}{m}_2)\\ ⋮\\ x\equiv a_n(\mathrm{mod}{m}_n)\end{array}$$
保证 $\forall i,j$ 满足 $\gcd (m_i,m_j)=1$，求出 $x$ 的最小非负整数解。

Part 2. 中国剩余定理的求解

令 $M=\prod _{i=1}^n{m}_i,M_i=\frac{M}{m_i},M_i{t}_i\equiv 1(\mathrm{mod}{m}_i)$。
则任意解 $x=Mk+\sum _{i=1}^n{a}_i{M}_i{t}_i$，最小非负整数解 $x_{\min }=\sum _{i=1}^n{a}_i{M}_i{t}_i$。

算法流程：

1. 计算出 $M=\prod _{i=1}^n{m}_i$。
2. 遍历每一个方程，求出 $M_i=\frac{M}{m_i}$，用扩展欧几里得算法求出 $t_i$。
3. 得到 $\sum _{i=1}^n{a}_i{M}_i{t}_i$。

总时间复杂度 $\mathcal{O}(n\log n)$。

Part 3. 证明

考虑第 $i$ 个方程，对于 $\forall j\ne i\in [1,n]$，显然有 $a_j{M}_j{t}_j\equiv 0(\mathrm{mod}{m}_i)$。

所以 $\sum _{j=1}^n{a}_j{M}_j{t}_j\equiv a_i{M}_i{t}_i\equiv a_i(\mathrm{mod}{m}_i)$。

Part 4. 代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e5+5;
inline int read(){
    char op=getchar();
    int w=0,s=1;
    while(op<'0'||op>'9'){
        if(op=='-') s=-1;
        op=getchar();
    }
    while(op>='0'&&op<='9'){
        w=(w<<1)+(w<<3)+op-'0';
        op=getchar();
    }
    return w*s;
}
int Pow(int a,int k,int mod){
    int ans=1;
    while(k){
        if(k&1) ans=(ans%mod*a%mod)%mod;
        a=(a%mod*a%mod)%mod;
        k>>=1;
    }
    return ans;
}
void exgcd(int a,int b,int &x,int &y){
    if(b==0){
        x=1,y=0;
        return;
    }
    exgcd(b,a%b,x,y);
    int t=x;
    x=y;
    y=t-a/b*y;
}
int inv(int x,int mod){return Pow(x,mod-2,mod);}
int M=1ll,a[N],b[N],m[N],t[N],ans=0;
signed main(){
    int n=read();
    for(register int i=1;i<=n;i++) a[i]=read(),b[i]=read();
    for(register int i=1;i<=n;i++) M=M*a[i];
    for(register int i=1;i<=n;i++){
        m[i]=M/a[i];
        int y=0;
        exgcd(m[i],a[i],t[i],y);
        if(t[i]<0) t[i]+=a[i];
    }
    for(register int i=1;i<=n;i++) ans=(ans+(b[i]%M*m[i]%M)%M*t[i]%M)%M;
    printf("%lld",(ans+M)%M);
}