扩展欧几里得算法——pku1061

#include<stdio.h>

__int64 Ext_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
if(b==0) { x=1, y=0; return a; }
__int64 res= Ext_gcd(b,a%b,y,x);
y-= a/b*x;
return res;
}
__int64 Ax_b_mod_n(__int64 a,__int64 b,__int64 n)//同余方程 (m-n)*t = (y-x) (mod L) 的最小正整数解
{
__int64 res,x,y,t;
res= Ext_gcd(a,b,x,y);
if(n%res==0)
{
x= x*n/res;
t= b/res;
if(t<0) t=-t;
x= (x%t+t)%t;
return x;
}
return -1;
}
int main()
{
__int64 x,y,m,n,len,t;
scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&len);
t= Ax_b_mod_n(m-n,len,y-x);
if(t==-1) printf("Impossible\n");
else printf("%I64d\n",t);
}