poj1061 青蛙的约会 扩展欧几里得

中文题

思路: 扩展欧几里得

扩展欧几里得

• 给定整数a 和 b, 且满足 a*x1 + b*y1 = gcd(a, b), 求解x, y.
• 当 b == 0 的时候, gcd(a, b) = a. 此时x = 1, y = 0.
• 当a*b != 0的时候, 推理: 根据欧几里得可知, gcd(a, b) = gcd(b, a%b). 那么 a*x1 + b*y1 = b*x2 + a%b*y2

　　　化解上式得到 b*x2 + (a - (a/b)*b)*y2 = a*y2 + b*x2 - (a/b)*b*y2.

　　　而根据恒等定理可知: x1 = y2, y1 = x2 - a/b*y2 ,由此我们可以知道求得x1, x2 的结果是基于x2, y2 而得到的, 因此我们只需要不断的递归下去,即:

　　　gcd(a, b) =gcd(b, a%b)=gcd(a%b, b%(a%b)) 直到b == 0 返回x2, y2 的结果, 通过恒等定理求得x1, y1.

• 对于方程a*x1 + b*y1 = c , 我们可以根据上面求得的方程 a*x1 + b*y1 = gcd(a, b)转换得到(设d = gcd(a, b)).
• a*x1 + b*y1 = d 转换得到: a(x1*c/d) + b(y1*c/d) = c.
• 因此方程的解为x = x1*c/d, y = y1*c/d.

构建方程:

 (x + m* t) % L = (y + n*t) % L

转换得到:

x + m*t = y + n*t + w*L;
 (n - m)*t + w*L = x - y;
 即扩展欧几里得:a*x + b*t = c

 1 #include <iostream>

 2 using namespace std;

 3 

 4 long long extended_gcd(long long a, long long b, long long  &x, long long  &y)

 5 {

 6     long long ret , tmp;

 7     if (b == 0)

 8     {

 9         x = 1, y = 0;

10         return a;

11     }

12     ret = extended_gcd(b, a%b, x, y);

13     tmp = x;

14     x = y;

15     y = tmp - a /b *y;

16     return ret;

17 }

18 

19 int main()

20 {

21     long long x, y, d, X, Y, n, m, L;

22     /*

23     x + m*t = y + n*t + w*L;

24     (n - m)*t + w*L = x - y; 

25     扩展欧几里得:a*x + b*t = b

26     */

27     while (cin>>x>>y>>m>>n>>L)

28     {

29         d = extended_gcd(n-m, L, X, Y);

30         if ((x - y) % d != 0)

31             cout<<"Impossible"<<endl;

32         else

33             cout<<((x - y) / d * X % L + L) % L<<endl;

34     }

35 

36     

37 }

View Code