中国剩余定理是求解一次同余式组的方法
当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;
}