poj1061

// 由题意可得：求解公式
//1.(x+k*m)-(y+k*n)=h*L(k,h=0,1,2,3....)
//2. 得：(n-m)*k+L*h=x-y;
//3.求解线性方程：扩展欧几里得
//4.令n-m=a;b=L;c=x-y;即：a*k+b*h=c;
//5.线性方程有解得条件是：c%gcd(a,b)==0;
//6.求解方程，扩展欧几里得a*x1+b*y1=b*x2+(a mod b)*y2=b*x2+(a-a/b*b)*y2;
//7.展开得：x1=x2;y1=x2-a/b*y2;
//8.根据数论中的相关定理，可得方程a' * x + b' * y = n'的所有整数解为：
//x = n' * x0 + b' * t
//y = n' * y0 - a' * t
//因此：在找最小值时：
//此时方程的所有解为：x=c*x'-b*t,x的最小的可能值是0，
//令x=0可求出当x最小时的t的取值，但由于x=0是可能的最小取值，实际上可能x根本取不到0，
//那么由计算机的取整除法可知：由 t=c*k1/b算出的t，代回x=c*k1-b*t中，求出的x可能会小于0，此时令t=t+1，
//求出的x必大于0；如果代回后x仍是大于等于0的，那么不需要再做修正。
#include<stdio.h>
__int64 x,y,m,n,L;
__int64 gcd(__int64 a,__int64 b)
{
    return b?gcd(b,a%b):a;
}
void  extend_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    if(b==0)
    {
     x=1;
     y=0;
     return ;
    }
    extend_gcd(b,a%b,x,y);
    __int64 t;
    t=x;
    x=y;
    y=t-a/b*x;
}
int main()
{
    __int64 a,b,c,Y,t;
    while(scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&L)!=EOF)
    {
        a=n-m;
        b=L;
        c=x-y;
        Y=gcd(a,b);
        if(c%Y)
        {
            printf("Impossible\n");
            continue;
        }
        a/=Y;
        b/=Y;
        c/=Y;
        extend_gcd(a,b,x,y);
    t=c*x/b;
    x=c*x-b*t;
   if(x<0)
     x+=b;
       printf("%I64d\n",x);
    }
    return 0;
}