题意:给出k个模方程组:x mod ai = ri。求x的最小正值。如果不存在这样的x,那么输出-1.
题解:
由于这道题目里面的ai、ri之间不满足两两互质的性质,所以不能用中国剩余定理直接求解。
不过,我们可以模仿中国剩余定理的做法来解决这个问题。
如果只有一个方程:x mod a0 = r0。那么,显然x的最小正值为a0+r0。
根据模的性质,我们容易得知,x+a0*k均为该方程的解。(k为正整数)
如果多了一个方程:x mod a1 = r1。那么,我们为了使之间求得的解x0=a0+r0能够同时满足这两个方程,只好令
x0=r0+a0*k,显然这样做x0仍然满足第一个方程。这时候我 们相当于要求解这样一个模方程:(r0+a0*k) mod a1 =
r1。这个方程我们可以用拓展欧几里得算法求得k的值。这样,只要令x0变成x0+a0*k,就能同时满足这两个方程了。
推而广之,对于方程x mod ai = ri,假如我们之前求得的解为X,那么我们要令X变成X+k*LCM(a0,a1,a2...ai-1),使得
它满足这个方程。k我们可以用拓展欧几里得 算法求解,LCM可以在每一次更新,这样就能在接近O(klogk)的时间复
杂度内解决这个问题了。
无解的判断:若某个(X+k*LCM) mod ai = ri无整数解,那么原方程组无解。
在<<ACM-ICPC程序设计系列--数论及应用>>的54页有解释。
见黑书P230页。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define LL long long
using namespace std;
LL extend_gcd(LL a,LL b,LL &x,LL &y){
if(a==0&&b==0) return -1;
if(b==0){
x=1,y=0;return a;
}
LL d=extend_gcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
int main(){
LL i,n,a1,r1,a2,r2,ans,d,x,y;
while(~scanf("%lld",&n)){
bool have=1;
scanf("%lld%lld",&a1,&r1);
for(i=1;i<n;i++){
scanf("%lld%lld",&a2,&r2);
LL a,b,c;
a=a1,b=a2,c=r2-r1;
d=extend_gcd(a,b,x,y);
if(!have) continue;
if(c%d!=0){
have=0;
}
int t=b/d;
x=(x*(c/d)%t+t)%t;
r1=a1*x+r1;
a1=a1*(a2/d);
}
if(!have){
puts("-1");
}
else printf("%lld\n",r1);
}
return 0;
}