1516 青蛙的约会

1516 青蛙的约会

我是来练习同余的…
两个小青蛙谈恋爱了，他们要奔现了，他们各自在数轴上跳，给定他们一次跳的距离m和n，还有他们所在的位置x和y，以及数轴的长度L，求他们最少几次能相遇
如果不能相遇就输出“Impossible”
也就是计算周期相遇问题
两个恋爱的青蛙只有跳到了同一个点才算相遇了，所以我们可以推出式子
(x+mt)-(y+nt)=pL
其中t是跳的次数，p是两人的跳的圈数之差，公式就是路程差等于纬线周长的整数倍
然后转化一下这个式子
(n-m)t+Lp=x-y
令a=n-m，b=L，c=gcd(a,b)，d=x-y
有：at+bp=d
求：t的最小整数解决
用exgcd求其中一个整数解t0，p0
并令c=gcd(a,b)
有：at0+bp0=c
然后将这个式子乘上(d/c)，得到at0d/c+bp0d/c=d
所以t0d/c是最小的解，但有可能是负数
所以解为(t0d/c%b+b)%b
那么如何求t0和p0呢
这里就可以用exgcd来求at0+bp0=c的解，t0和p0
ax+by=bx0+(a%b)y0=bx0+(a-a/bb)y0=ay0+(x0-a/by0)b
所以x=y0，y=x0-a/by0
不过要开 longlong
不过开了还是只有70分
负数要进行特判，真好

#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
using namespace std;
long long typedef ll;
ll x,y,m,n,L;
ll gcd; 
ll X,Y;
ll exgcd(ll a,ll b,ll &x,ll &y)//计算t0 p0 
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	int gcd=exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-a/b*y;
	return gcd;
}
int main()
{
	cin>>x>>y>>m>>n>>L;
	int a=x-y,b=n-m;
	if(b<0)
	{
		b=-b;
		a=-a;
	}
	if((a)%(gcd=exgcd(b,L,X,Y))!=0)
	{
		cout<<"Impossible"<<endl;
	}
	else
	{
		cout<<(a/gcd*X%(L/gcd)+(L/gcd))%(L/gcd)<<endl;
	}
	return 0;
}