poj 2891 Strange Way to Express Integers (扩展gcd)

题目链接

题意：给k对数，每对ai, ri。求一个最小的m值，令m%ai = ri;

分析：由于ai并不是两两互质的， 所以不能用中国剩余定理。

只能两个两个的求。

a1*x+r1=m=a2*y+r2
联立得：a1*x-a2*y=r2-r1;
设r=r2-r2;

互质的模线性方程组m=r[i](mod a[i])。两个方程可以合并为一个，新的a1为lcm(a1,a2)，

新的r为关于当前两个方程的解m，然后再和下一个方程合并……。(r2-r1)不能被gcd(a1,a2)整除时无解。   怎么推出的看了好多博客也没有介绍。

下面求最小的解的时候有一个定理，上一篇博客青蛙 也用到了：   对方程  ax  ≡  b (mod) n

d = gcd(a,n) 若 d | b 则 有 d 个解 ,最小的解为x0 = (x*(b/d) mod n + n )mod(n/d)

则所有满足方程 x = x0 (mod) (n/d) 的 x 都是原方程的解。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #define LL long long
 using namespace std;

 void exgcd(LL a, LL b, LL &d, LL &x, LL &y)
 {
     if(!b) {d = a; x = ; y = ;}
     else{ exgcd(b, a%b, d, y, x); y -= x*(a/b); }
 }
 int main()
 {
     int flag;
     LL k, a1, a2, r1, r2;
     LL c, d, x, y, q;
     while(~scanf("%lld", &k))
     {
         flag = ;
         scanf("%lld%lld", &a1, &r1);
         k--;
         while(k--)
         {
             scanf("%lld%lld", &a2, &r2);
             if(flag)
             continue;
             c = r2-r1;
             exgcd(a1, a2, d, x, y);
             if(c%d!=)
             flag = ;
             q = a2/d;
             x = (x*(c/d)%q + q)%q; //求最小的解
             r1 = x*a1 + r1;  //令r1为所求值，x;
             a1 = a1*a2/d;  //令a1为a1a2最小公倍数
         }
         if(flag)
         printf("-1\n");
         else
         printf("%lld\n", r1);
     }
     return ;
 }