OJ青蛙的约会

Description

两只青蛙在网上相识了，它们聊得很开心，于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直到碰面为止。可是它们出发之前忘记了一件很重要的事情，既没有问清楚对方的特征，也没有约定见面的具体位置。不过青蛙们都是很乐观的，它们觉得只要一直朝着某个方向跳下去，总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上，不然是永远都不可能碰面的。为了帮助这两只乐观的青蛙，你被要求写一个程序来判断这两只青蛙是否能够碰面，会在什么时候碰面。
我们把这两只青蛙分别叫做青蛙A和青蛙B，并且规定纬度线上东经0度处为原点，由东往西为正方向，单位长度1米，这样我们就得到了一条首尾相接的数轴。设青蛙A的出发点坐标是x，青蛙B的出发点坐标是y。青蛙A一次能跳m米，青蛙B一次能跳n米，两只青蛙跳一次所花费的时间相同。纬度线总长L米。现在要你求出它们跳了几次以后才会碰面。

Input

输入只包括一行5个整数x，y，m，n，L，其中x≠y < 2000000000，0 < m、n < 2000000000，0 < L < 2100000000。

Output

输出碰面所需要的跳跃次数，如果永远不可能碰面则输出一行"Impossible"

Sample Input

1 2 3 4 5

Sample Output

4

分析与解答：

此题其实就是扩展欧几里德算法－求解不定方程。

设经过s步后，A、B两只青蛙相遇，则可得出相应的方程：x+m*s-(y+m*s)=k*L,(其中k=0,1,2,3……)，稍微变一下形得：

(n-m)*s+k*l=x-y，则令n-m=a,l=b,x-y=c，该方程变为：a*s+b*k=c

即转化为

a*x+b*y=c  （1）

求最小x的模型了。

之所以要转化为这个方程，是为了便于建立数学模型，所以在OJ上遇到任何问题，首先要分析，并建立相应的数学模型。

所以只要上式存在整数解，则两青蛙能相遇，否则不能。

那么如何求上式方程的整数解x呢？

思路：用欧几里德算法和扩展欧几里德算法。

欧几里德算法又称辗转相除法，用于计算两个整数a,b的最大公约数。其计算原理依赖于下面的定理：

　　定理：gcd(a,b) = gcd(b,a%b)（证明略）

    首先，设d = gcd(a, b)，方程两边除以d得到a/d*x + b/d*y = c/d，很显然嘛，a是整除d的，b也是整除d的，而x、y都是整数解，所以要求c/d也是整数嘛。如果c不整除d，即没有整数解x，当然就是Impossible。

    然后用扩展欧几里德算法求出其中一组解x0 ,y0, 有方程a*x0+b*y0=d，并令d = gcd(a, b)，那么两边乘以c/d即a(c/d * x0) + b(c/d * y0) = c  （2）

    对方程（1）和（2），可知x = (c/d * x0)，y = (c/d * y0)。

    最后的问题就是如何用扩展欧几里得算法求出x0和y0呢？

    对于的非负整数a, b， gcd(a, b)表示a, b 的最大公约数。那么存在整数x0， y0使得a*x0+b*y0=gcd(a,b)

    1、当b = 0 时，gcd(a, b) = a , 此时 x0 = 1, y0 = 0

    2、b不为0时，有a*x0+b*y0=gcd(a,b)= gcd(b, a%b) = bx1 + (a%b)y1，而a%b又可以写成a-a/b*b，接着不停的递归就可以就出x0和y0。

由此可以得出扩展欧几里德的递归程序：

_int64 gcd(_int64 a, _int64 b, _int64 &x, _int64 &y)
{
   _int64 d, t;
   if (b == 0)

     {

       x = 1; y = 0;

return a;

      }

     else
        d = gcd(b, a % b, x, y);
        t = x - a/b*y;

        x = y;

        y = t;
   return d;
}

源程序：


#include

using namespace std;

_int64 gcd( _int64 a, _int64  b, _int64  &x, _int64 &y)

{

    _int64  d, t;

    if (b == 0)
{
    x = 1; y = 0; return a;
}

   else

    d = gcd(b, a % b, x, y);

    t = x - a / b * y;
     x = y; y = t;

    return d;

}

int main()

{

    _int64  x, y, m, n, L, X, Y, d, r;

    while (scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&L) )

    {

        d = extgcd(n - m, L, X, Y);
        r = L / d;

        if ((x - y) % d)

           printf("%I64d\n","impossible");

        else
            printf("%I64d\n",((x - y) / d * X % r + r) % r);
     }

    return 0;

}