BIT1052 poj1061 青蛙的约会

拿了以前在pojAC的代码一交，WA了。。。

重写了、

根据题意有

p+m*k=q+n*k (mod L)

化简

p-q+(m-n)k=(k')*L

再化简

(m-n)*k+k'*(-L)=q-p

另a=m-n,b=-L

a*k+b*k'=q-p

这就可以用扩展欧几里德了

#include<iostream>
#include<algorithm>
using namespace std;
long long x,y;
long long extended_gcd(long long a,long long b)
{
long long t,gcd;
if(b==0)
{
x=1;
y=0;
return a;
}
gcd=extended_gcd(b,a%b);
t=x;
x=y;
y=t-a/b*y;
return gcd;
}
int main()
{
long long p,q,m,n,l;
while(~scanf("%lld %lld %lld %lld %lld",&p,&q,&m,&n,&l))
{
long long a=m-n;
long long b=-1*l;
long long d=extended_gcd(a,b);//求出d公约数，同时求出(a/d)*k+(b/d)*k'=1的一组解
if((q-p)%d!=0)
{
printf("Impossible\n");
continue;
}
long long k=x*(q-p)/d;//获得k的一个解
//k的通解为k+b/d*n   n为1，2，3，，
k%=(b/d);
if(b/d>0)
{
while(k<0)
{
k+=b/d;
}
}
if(b/d<0)
{
while(k<0)
{
k-=b/d;
}
}
printf("%lld\n",k);
}
return 0;
}