POJ1061 青蛙的约会

题目：
http://poj.org/problem?id=1061
分析：
我们可以知道x0+m*t≡y0+n*t (mod L)
也就是x0+m*t1+t2*L=y0+n*t1
化简得(m-n)*t1+L*t2=y0-x0
这就是ax+by=c的形式，其中a=(m-n)，b=L，c=y0-x0
用扩展欧几里得即可求解t1
这里需要注意的是，m-n有可能为负数，求得的最大公因数可能有问题，此时需要令a=-a，c=-c，因为t2是未知数，所以方程两边同乘以-1没有问题，t2把-1吃掉了
用c%gcd(a,b)==0判断有解
用扩展欧几里得求得方程ax+by=gcd(a,b)的解x，然后x*=c/gcd(a,b)，将x变换为ax+by=c的解
x的周期为t=b/gcd
所以x>0时最小正整数解为x%t，x<0时最小正整数解为x%t+t
代码：

#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
long long x0,y0,m,n,L;
long long exgcd(long long a,long long b,long long& x,long long& y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    long long r=exgcd(b,a%b,y,x);
    y-=x*(a/b);
    return r;
}
int main()
{
    long long a,b,c,gcd,x,y,t;
    scanf("%lld%lld%lld%lld%lld",&x0,&y0,&m,&n,&L);
    a=m-n;
    b=L;
    c=y0-x0;
    if(a<0)
    {
        a=-a;
        c=-c;
    }
    gcd=exgcd(a,b,x,y);
    if(c%gcd!=0) printf("Impossible");
    else
    {
        x*=c/gcd;
        t=b/gcd;
        if(x>0) x%=t;
        else x=x%t+t;
        printf("%lld",x);
    }
    return 0;
}