POJ 2891 Strange Way to Express Integers[中国剩余定理(非互质)][扩展欧几里得]

题意

给出n个一元一次同余方程$X≡a_i(mod$ $r_i)$

求最小的$X$,模数不一定互质

题解

模数非互质不能直接上中国剩余定理

任意找两个式子出来:

$\left\{\begin{array}\\X≡a_1(mod \ r_1)\\{X≡a_2(mod\ r_2)}\end{array}\right.$

可以写成

$\left\{\begin{array}\\X=k_1*r_1+a_1\\{X=k_2*r_2+a_2}\end{array}\right.$

即 $k_1*r_1+a_1=k_2*a_2+a_2$

写成 $k_1*a_1≡r_2-r_1(mod\ a_2)$

用扩欧找出最小的非负$k_1$,同时判一下无解的情况.

代到原来的方程中,得$x_0=k_1*r_1+a_1$.

则$X=x_0(mod\ lcm(a_1,a_2))$这其实就相当于得到一个新的方程了,新的余数就是$x_0$.

然后往后一直做,最后得到的$x_0$就是答案.

#include<cstdio>
#include<algorithm>
#define LL long long
#define N 1005
using namespace std;

int n;
LL a[N],r[N];

inline LL readl(){
    LL a=0;char f=1,c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){a=a*10ll+c-'0';c=getchar();}
    return a*f;
}

inline LL exgcd(LL a,LL b,LL &x,LL &y){
    if(!b){
        x=1,y=0;return a;
    }
    LL r=exgcd(b,a%b,x,y),t;
    t=x,x=y,y=t-a/b*y;
    return r;
}

//inline LL lcm(LL A,LL B){LL x,y,t=exgcd(A,B,x,y);return A*B/t;}

inline LL work(){
    LL A=a[1],R=r[1],x,y,d;
    for(int i=2;i<=n;++i){
        d=exgcd(A,a[i],x,y);
        if((R-r[i])%d) return -1;
        x=(R-r[i])/d*x%a[i];
        R-=A*x;
        A=A*a[i]/d;
        R=(R+A)%A;
    }
    return R;
}

int main(){
    while(~scanf("%d",&n)){
        for(int i=1;i<=n;++i) a[i]=readl(),r[i]=readl();
        printf("%lld\n",work());
    }
    return 0;
}