中国剩余定理（孙子定理）（模板）

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/yz467796454/article/details/82014536

中国剩余定理是求解一次同余式组的方法

当a互质的时候：

#include<iostream>
#include<cstdio>
#include<climits>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m[105],a[105],lcm=1;
int exgcd(int a,int b,int &x,int &y){//扩欧
    if(!b){x=1,y=0;return a;}
    int re=exgcd(b,a%b,x,y),tmp=x;
    x=y,y=tmp-(a/b)*y;
    return re;
}
int work(){
    int i,j,d,x,y,re=0;
    for(i=1;i<=n;i++)lcm=lcm*m[i];//因为互质所以直接这么写了
    for(i=1;i<=n;i++){
        int kl=lcm/m[i];
        d=exgcd(kl,m[i],x,y);
        x=(x%m[i]+m[i])%m[i];
        re=(re+a[i]*x*kl)%lcm;
    }
    return re;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d%d",&m[i],&a[i]);;
    printf("%d",work());
    return 0;
}

当a不互质的时候：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+5;

ll exgcd(ll a,ll b,ll &x,ll &y){
    if(!b){x=1,y=0;return a;}
    ll re=exgcd(b,a%b,x,y),tmp=x;
    x=y,y=tmp-(a/b)*y;
    return re;
}
ll m[maxn],a[maxn];
int n;

ll work(){
    ll M=m[1],A=a[1],t,d,x,y;
	int i;
    for(i=2;i<=n;i++){
        d=exgcd(M,m[i],x,y);//解方程
        if((a[i]-A)%d)return -1;//无解
        x*=(a[i]-A)/d,t=m[i]/d,x=(x%t+t)%t;//求x
        A=M*x+A,M=M/d*m[i],A%=M;
    }
    A=(A%M+M)%M;
    return A;
}
int main(){
    while(scanf("%d",&n)!=EOF){
        for(int i=1;i<=n;i++){
			scanf("%lld%lld",&m[i],&a[i]);
        }
        printf("%lld\n",work());
    }
    return 0;
}