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

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

题目大意：中文题，不解释

方法：设t为A青蛙和B青蛙的跳的次数，k为绕地球绕的圈数

则得出公式：（x+m*t）- （y+n*t）= k*L

化简得：（x-y）+（m-n）*t = k*L

    化简：（m-n）*t - k*L = y-x（由扩展欧几里得定理可知:存在x0，y0使得a*x0+b*y0=gcd（a,b)）

最终t = （x0+k(L/d)）

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define LL long long

using namespace std;

LL exgcd(LL a,LL b,LL &x,LL &y)
{
    LL r,t;
    if(b==0)
    {
        x=1, y=0;
        return a;
    }
    r=exgcd(b, a%b, x, y);
    t=x;
    x=y;
    y=t-a/b*y;
    return r;
}
int main()
{
    LL n, L, m, x, y, t;
    LL xx, yy;
    cin>>x>>y>>m>>n>>L;
    LL d = exgcd(n-m, L, xx, yy);
    if((x-y) % d != 0)
        cout<<"Impossible"<<endl;
    else
        cout<<((xx*(x-y)/d)%(L/d)+(L/d))%(L/d)<<endl;
    return 0;
}