很经典的一道多个同余式求解的题目。
中国剩余定理只能解决每个mod的数两两互素的情况,但是我们通过方程合并可以很简单地处理不互素地情况。
比如有两个同余方程: x = a1modb1 x = a2modb2
那么我们可以把他们改写成这样的形式:x=a1+b1y1 x=a2+b2y2
合并为一个方程:a1+b1y1 = a2+b2y2
化简:b1y1 = a2-a1 +b2y2
那么我们有:b1y1 = (a2-a1) mod b2
这个可以通过逆元或者拓展欧几里得判断有无解,并且算出y1的值。
算出y1的值后,我们能形成一个新的方程: x = p mod lcm(b1,b2)。
这样通过一步一步合并我们就能得到最终的方程。
这题有一个陷阱,需要特判最后的答案是否是0。
#include "cstdio"
#define LL long long
LL gcd(LL a,LL b){
return b==0 ? a:gcd(b,a%b);
}
void extend_gcd(LL a,LL b,LL &d,LL &x,LL &y){
if(!b){d=a;x=1;y=0;}
else{extend_gcd(b,a%b,d,y,x);y-=x*(a/b);}
}
LL inv(LL a,LL n){
LL d,x,y;
extend_gcd(a, n, d, x, y);
return d == 1? (x+n)%n:-1;
}
int main(){
int n;
int i;
LL a1,m1,a2,m2;
LL d;
LL tem;
int flag;
LL K;
while(scanf("%d",&n)!=EOF){
flag=0;
scanf("%lld%lld",&m1,&a1);
for(i=1;i<n;i++){
scanf("%lld%lld",&m2,&a2);
d=gcd(m1,m2);
if((a2-a1)%d==0){
tem=inv(m1/d, m2/d);
K=(tem*(a2-a1)/d)%(m2/d);
a1=(K*m1+a1)%(m1*m2/d);
m1=m1*m2/d;
}else{
flag=1;
}
}
if(flag==1) printf("-1\n");
else{
a1=(a1+m1)%m1;
if(a1==0) a1=m1;
printf("%lld\n",a1);}
}
return 0;
}