BITACM 1053 Han Xin Count Soldiers

很经典的一道多个同余式求解的题目。

中国剩余定理只能解决每个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;
}