poj 1061 青蛙约会（扩展欧几里德）

题目链接：

　　http://poj.org/problem?id=1061

题目大意：

　　中文题目，题意一目了然，就是数据范围大的出奇。

解题思路：

　　假设两只青蛙都跳了T次，可以列出来不定方程：p*l + (n-m)*T == x - y。列出等式以后，利用扩展欧几里德计算不定方程的解。在求出整数最小解的地方卡了好久，好久。

想具体了解扩展欧几里德的用法和证明的话，可以看一下神牛的博文，我自认弱绞尽脑汁也写不来这么好，附上链接：http://www.cnblogs.com/frog112111/archive/2012/08/19/2646012.html

代码：

 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <algorithm>
 #include <iostream>
 #include <cmath>
 #include <queue>
 using namespace std;
 #define LL long long
 LL gcd (LL a, LL b, LL &x, LL &y)//扩展欧几里德模板
 {
     if (b == )
     {
         x = ;
         y = ;
         return a;
     }
     LL r = gcd (b, a%b, x, y), t;
     t = x;
     x = y;
     y = t - a / b * y;
     return r;
 }
 int main ()
 {
     LL x, y, m, n, l;
     while (scanf ("%lld %lld %lld %lld %lld", &x, &y, &m, &n, &l) != EOF)
     {
         LL u, v, a, b, c;
         a = n - m;
         b = l;
         c = gcd(a, b, u, v);//求出来的u是gcd(a, b) = au + bv中的。
         if ((x - y) % c != )//x*gcd(a, b) = a*u*x + b*v*x成立，才能求出题目中的解
         {
             printf ("Impossible\n");
             continue;
         }
         LL s = b / c;
         u = u * (x - y) / c;//这个时候才是(x-y) = (n-m)*t + p*l中t的解
         u = (u % s + s) % s;
         printf ("%lld\n", u);
     }
     return ;
 }