EXCRT/扩展中国剩余定理学习笔记

当a1,a2a_{1},a_{2}a1​,a2​互质的时候
我们可以使用中国剩余定理合并
即∑ibi∗(Mai∗Inv(Mai))\sum_{i}b_{i}*(\frac{M}{a_{i}}*Inv_{(\frac{M}{a_{i}})})∑i​bi​∗(ai​M​∗Inv(ai​M​)​)
其中M=∏iaiM=\prod_{i}a_{i}M=∏i​ai​
但是不同于别的定律
从CRT到EXCRT是完全不同的

观察下面$T$ $mod$ a1a_{1}a1​=b1b_{1}b1​ $T$ $mod$ a2a_{2}a2​=b2b_{2}b2​
列出以下式子 a1∗x1+b1=a2∗x2+b2a_{1}*x_{1}+b_{1}=a_{2}*x_{2}+b_{2}a1​∗x1​+b1​=a2​∗x2​+b2​
转化为a1∗x1−a2∗x2=b2−b1a_{1}*x{1}-a_{2}*x_{2}=b_{2}-b_{1}a1​∗x1−a2​∗x2​=b2​−b1​
所以可以EXGCD
解出最小正整数解x1x_{1}x1​后回带求出$T$
这时有着一个新的同余式
$T$ $mod$ LCM(a1,a2)LCM(a_{1},a_{2})LCM(a1​,a2​) = bSolb_{Sol}bSol​
递推计算$n-1$次即答案

#include<bits/stdc++.h>
using namespace std;
typedef int INT;
#define int long long
const int N=1e5+10;
inline void read(int &x){
    x=0;
    int f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    x*=f;
}
struct Node{
    int ret,mod;
}P[N];
int n;
int EXGCD(int &x,int &y,int a,int b){
    if(!b){
        x=1;
        y=0;
        return a;
    }
    int GCD=EXGCD(x,y,b,a%b);
    int tmp=x;
    x=y;
    y=tmp-(a/b)*y;
    return GCD;
}
int mul(int x,int y,int mod){
    int ret=0;
    while(y){
        if(y&1)ret=(ret+x)%mod;
        x=(x+x)%mod;
        y/=2;
    }
    return ret;
}
void EXCRT(int tot){
    int ret=P[1].ret;
    int mod=P[1].mod;
    for(int i=2;i<=tot;++i){
        int x,y;
        int GCD=EXGCD(x,y,mod,P[i].mod);
        int now=((P[i].ret-ret)%P[i].mod+P[i].mod)%P[i].mod;
        x=(x%P[i].mod+P[i].mod)%P[i].mod;
        x=mul(x,now/GCD,P[i].mod)%P[i].mod;
        ret=(ret+mul(x,mod,mod/GCD*P[i].mod));
        mod=mod/GCD*P[i].mod;
        ret=ret%mod;
    }
    cout<<ret;
}
INT main(){
//	freopen("test.in","r",stdin);
    read(n);
    for(int i=1;i<=n;++i){
        read(P[i].mod);
        read(P[i].ret);
    }
    EXCRT(n);
}

转载于:https://www.cnblogs.com/Leo-JAM/p/10079085.html