我们把这两只青蛙分别叫做青蛙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);
}