题意
给出n个一元一次同余方程 X≡ai(mod ri)
求最小的 X ,模数不一定互质
题解
模数非互质不能直接上中国剩余定理
任意找两个式子出来:
可以写成
{X=k1∗r1+a1X=k2∗r2+a2
即 k1∗r1+a1=k2∗a2+a2
写成 k1∗a1≡r2−r1(mod a2)
用扩欧找出最小的非负 k1 ,同时判一下无解的情况.
代到原来的方程中,得 x0=k1∗r1+a1 .
则 X=x0(mod lcm(a1,a2)) 这其实就相当于得到一个新的方程了,新的余数就是 x0 .
然后往后一直做,最后得到的 x0 就是答案.
#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;
}