题目链接:http://poj.org/problem?id=1061
本题是简单的扩展欧几里得算法。
基本全裸。
请读者自行推导。
即推出ax+by=z的形式,然后利用扩展欧几里得算法算出其中x,y的一组解,然后求其他满足条件的解。
下面给出AC代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
__int64 exgcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
__int64 d=a;
if(b!=0)
{
d=exgcd(b,a%b,y,x);
y-=(a/b)*x;
}
else
{
x=1;
y=0;
}
return d;
}
int main()
{
__int64 s1,s2,v1,v2,m,x,y;
while(~scanf("%I64d%I64d%I64d%I64d%I64d",&s1,&s2,&v1,&v2,&m))
{
__int64 gcd=exgcd(v2-v1,m,x,y);
if((s1-s2)%gcd)
{
printf("Impossible\n");
}
else
{
__int64 d=(s1-s2)/gcd;
x*=d;
__int64 r=m/gcd;
printf("%I64d\n",(x%r+r)%r);
}
}
return 0;
}