【ZJOI2002】poj1061 青蛙的约会

我们把这两只青蛙分别叫做青蛙A和青蛙B，并且规定纬度线上东经0度处为原点，由东往西为正方向，单位长度1米，这样我们就得到了一条首尾相接的数轴。设青蛙A的出发点坐标是x，青蛙B的出发点坐标是y。青蛙A一次能跳m米，青蛙B一次能跳n米，两只青蛙跳一次所花费的时间相同。纬度线总长L米。现在要你求出它们跳了几次以后才会碰面。

x≠y < 2* 10^9，0 < m,n < 2* 10^9，0 <L < 2.1*10^9。


讲解：
容易得到：
    (n-m)*t=y-x (mod l)
所以有：
    (n-m)* t+l * ? =y-x
改变一下变量名称，有：
    (n-m)*p+ l *q =y-x
其中p=t 。

因为可以求出：
    (n-m)*p+l * q=gcd(n-m,l);
要求：
    (n-m)*p +l *q=x-y
可以写成：
           (n-m)*p *(x-y)/gcd+l *q *(x-y)/gcd=x-y;
因此这里面的变量等于：
            p'=p *(x-y)/gcd;
            q'=q *(x-y)/gcd;
已求出特解。

接下来求齐次方程：
            (n-m)*p'+l *q'=x-y;
            那么:
            p'=l/gcd(p',q') * v=l/gcd(p,q) * v;
            q'=(n-m)/gcd(p',q') * v= (n-m)/gcd(p,q) *v;

因此通解为：
            p'=p *(x-y)/gcd + l/gcd(p,q) * v.

简化一下，令：
            t=l/gcd,
            则p'是t的倍数+p *(x-y)/gcd；

因此：
            p'=((p*(x-y)/gcd%t)+t)%t;
            为解得最小整数解。

#include<cstdio>
#include<cmath>
#define L long long
L abs(L a)
{
    return a>0?a:-a;
}
void euc(L a,L b,L &d,L &x,L &y)
{
    if (b==0)
    {
        d=a;
        x=1;
        y=0;
        return;
    }
    euc(b,a%b,d,y,x);
    y-=x*(a/b);
}
int main()
{
    L x,y,m,n,p,q,l,d,t;
    scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l);
    euc(m-n,l,d,p,q);
    if ((y-x)%d)
    {
        printf("Impossible\n");
        return 0;
    }
    t=abs(l/d);
    printf("%lld\n",((p*(y-x)/d)%t+t)%t);
}