中国剩余定理

中国剩余定理：若各方程的模数两两互质，设
M =
k ∏ i=1
mi,Mi = M mi,MiM′ i ≡ 1 (mod mi) 则同余方程组在模 M 下的解唯一且解为 x≡a1M1M′ 1 +a2M2M′ 2 + ... +akMkM′ k (mod M)

然后我们需要用到扩欧求Mi'的逆元就好了扩欧（求逆元的数，mod数，x，y）

#include<iostream>
#include<cstdio>
using namespace std;
int d,M=23*28*33,ans,t;
int m[4]={0,23,28,33};//存mod数 
int a[4];
void exgcd(int a,int b,int &d,int &x,int &y)
{
    if(!b)
    {
        x=1;
        y=0;
        d=a;
    }
    else
    {
        exgcd(b,a%b,d,x,y);
        int t=x;
        x=y;
        y=t-a/b*y;
    }
}
int intchina(int n)//中国剩余定理 
{
    ans=0;//注意初始化 
    int Mi,x,y,g;
    for(int i=1;i<=n;i++)
    {
        Mi=M/m[i];
        exgcd(Mi,m[i],g,x,y);
        ans=(ans+Mi*x*a[i])%M;
    }
    if(ans<=d) return 21252-(d-ans);
    else return ans-d;//ans的分类讨论 
}
int main()
{
    while(1)
    {
        t++;
        scanf("%d%d%d%d",&a[1],&a[2],&a[3],&d);
        if(a[1]==-1&&a[2]==-1&&a[3]==-1&&d==-1) break;
        a[1]%=23,a[2]%=28,a[3]%=33;
        ans=intchina(3);
        printf("Case %d: the next triple peak occurs in %d days.
", t, ans);
    }
}

https://www.luogu.org/problemnew/solution/P4777

把两个同余方程 x≡a1 (mod m1),x≡a2 (mod m2)

改写为 x = a1 +pm1,x = a2 +qm2

联立得 pm1 −qm2 = a2 −a1

通过拓展欧几里得算法得到 p,q 的一组解，

x′ = a1 +pm1 为这两个方程的一个解

于是得到新的同余方程 x≡x′ (mod lcm(m1,m2))

//niiick
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long lt;
const int maxn=100010;
int n;
lt ai[maxn],bi[maxn];

lt mul(lt a,lt b,lt mod)
{
    lt res=0;
    while(b>0)
    {
        if(b&1) res=(res+a)%mod;
        a=(a+a)%mod;
        b>>=1;
    }
    return res;
}

lt exgcd(lt a,lt b,lt &x,lt &y)
{
    if(b==0){x=1;y=0;return a;}
    lt gcd=exgcd(b,a%b,x,y);
    lt tp=x;
    x=y; y=tp-a/b*y;
    return gcd;
}

lt excrt()
{
    lt x,y,k;
    lt M=bi[1],ans=ai[1];
    for(int i=2;i<=n;i++)
    {
        lt a=M,b=bi[i],c=(ai[i]-ans%b+b)%b;
        lt gcd=exgcd(a,b,x,y),bg=b/gcd;
        if(c%gcd!=0) return -1;
        x=mul(x,c/gcd,bg);
        ans+=x*M;
        M*=bg;
        ans=(ans%M+M)%M;
    }
    return (ans%M+M)%M;
}

int main()
{
    cin>>n;
    for(int i=1;i<=n;++i)
    cin>>bi[i]>>ai[i];
    printf("%lld",excrt());
    return 0;
}