POJ 1061 青蛙的约会（扩展欧几里得算法）

题目链接：http://poj.org/problem?id=1061

本题是简单的扩展欧几里得算法。
基本全裸。

请读者自行推导。
即推出ax+by=z的形式，然后利用扩展欧几里得算法算出其中x,y的一组解，然后求其他满足条件的解。
下面给出AC代码：

#include
#include
#include
using namespace std;

__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    __int64 d=a;
    if(b!=0)
    {
        d=exgcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else
    {
        x=1;
        y=0;
    }
    return d;
}

int main()
{
    __int64 s1,s2,v1,v2,m,x,y;
    while(~scanf("%I64d%I64d%I64d%I64d%I64d",&s1,&s2,&v1,&v2,&m))
    {
        __int64 gcd=exgcd(v2-v1,m,x,y);
        if((s1-s2)%gcd)
        {
            printf("Impossible\n");
        }
        else
        {
            __int64 d=(s1-s2)/gcd;
            x*=d;
            __int64 r=m/gcd;
            printf("%I64d\n",(x%r+r)%r);
        }
    }
    return 0;
}