扩展欧几里德 POJ1061

#include<stdio.h>
typedef long long ll; //重新定义long long
ll x,y;
ll gcd(ll a,ll b)  //找ab的最大公约数
{
if(b==0)  return a;
else return gcd(b,a%b);
}
void extend_gcd(ll a,ll b)  //函数求解ax+by=d的一组解
{
if(b==0) //最底层  
{
x=1;y=0;
return;
}
extend_gcd(b,a%b);
ll temp=x;
x=y;  //x=y'
y=temp-a/b*y; //y=x'-a/b*y'
return;
}
int main()
{
 ll a,b,c,x0,y0,m,n,L;
 scanf("%lld%lld%lld%lld%lld",&x0,&y0,&m,&n,&L);
 a=m-n;b=-L;c=y0-x0; 
 ll d=gcd(a,b);
 if(c%d!=0) {
  printf("Impossible\n");
  return 0;
 }
 {
 extend_gcd(a,b);
 x=x*(c/d);
 y=y*(c/d);
 x=(x%(b/d)+(b/d))%(b/d); //求通解中的最小正整数x 
 printf("%lld\n",x);
 }
return 0; 
}