POJ 1061 青蛙的约会

问两只青蛙能否相遇，即问x+km=y+kn(mod l) 是否有解，如果有解，还需要输出最小解

移项后，变成一次同余式ax=b(mod l)

根据定理，当d=gcd(a,l)|b时，方程有解,且必然有一解为x'=u*b/d,其中au+lq=d,且共有l/d个解，X=x'+k*（l/d),k=0,1,2……(l/d-1)
由于X=x' mod （l/d)，最小解肯定在(0,l/d)中间

代码:

#include<iostream>
#include<memory.h>
#include<string>
#include<cstdio>
#include<algorithm>
#include<math.h>
#include<stack>
#include<queue>
#include<vector>
#include<map>
using namespace std;
long long gcd(long long a,long long b,long long &x,long long &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
long long d=gcd(b,a%b,x,y);
long long tmp=x;
x=y;
y=tmp-a/b*y;
return d;
}
int main()
{
long long n,m,a,b,l,x,y;
while(scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l)!=EOF)
{
a=((m-n)%l+l)%l;
b=((y-x)%l+l)%l;
long long d=gcd(a,l,x,y);
if(b%d!=0)
{
puts("Impossible");
continue;
}
long long ans=(b*x/d+l)%l;
//cout<<ans<<endl;
l/=d;
ans=(ans%l+l)%l;
printf("%lld\n",ans);
}
return 0;
}