能化出来线性同余方程:(m-n)t=(y-x) (mod L)
之后带公式进行求解就行,要求得最小正数解
由于网上都是n-m, x-y的方式来计算,求出来的gcd都是正的,但我也没从题里读出来大小关系啊…所以一开始用m-n,y-x算出来一直是负的,带网上的代码出来也是负的,所以一开始没有特别理解如何解线性同余方程还是比较麻烦的。
因此这里专门放个m-n和y-x的代码:(具体解的原理看其他博客或者看书吧)
AC代码:
#include <iostream>
using namespace std;
typedef long long LL;
//扩展欧几里得
//输入a,b,任意x,y
//返回ax+by=gcd(a,b)的一组特解x0,y0
LL exgcd(LL a,LL b,LL &x,LL &y)
{
if(b==0){
x=1,y=0;
return a;
}
LL d=exgcd(b,a%b,x,y);
LL z=x;x=y;y=z-y*(a/b);
return d;
}
int main()
{
LL x,y,m,n,L;
while(cin>>x>>y>>m>>n>>L){
LL X,Y;
LL d=exgcd(m-n,L,X,Y);
if((y-x)%d!=0){
cout<<"Impossible"<<endl;
return 0;
}
LL ans;
if(L/d>0) ans=(((y-x)/d*X)%(L/d)+L/d)%(L/d);
else ans=(((y-x)/d*X)%(L/d)-L/d)%(L/d);
cout<<ans<<endl;
}
}