poj 1061 青蛙的约会

//poj 1061 扩展的欧几里德，数据要求long long ...
#include <iostream>
#include <algorithm>

using namespace std;

const __int64 inf = ((__int64)1)<<62;

__int64 ext_euclid(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    if (b==0){ x=1; y=0; return a;}
    __int64 d=ext_euclid(b,a%b,x,y);
    __int64 t=x;
    x=y;
    y=t-a/b*y;
    return d;
}

__int64 solve(__int64  a,__int64  b,__int64  n)
{
     __int64  x,y;
     __int64 d=ext_euclid(a,n,x,y);
    
     if (b%d) return inf;
  n/=d;
     __int64  e=(x*(b/d))%n;
  return (e+n)%n;
}

int main()
{
    __int64 n,m,x,y,l;
    while (scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&l)!=EOF)
 {
          __int64 a=m-n,b=y-x;
         if (a<0)
  {
      a=-a,b=-b;
      if (b<0) b+=l;
   }
         __int64 ans=solve(a,b,l);
          if (ans==inf) printf("Impossible/n");
          else printf("%I64d/n",ans);
    }
    system("pause");
    return 0;
}