2013寒假练习 1052 青蛙的约会

地址：http://acm.bit.edu.cn/mod/programming/view.php?a=538

扩展欧几里得算法解ax+by=c，经典题。注意最后要输出x的最小的正整数解，故求出通解应该先%(b/gcd),再把负的加成正的。直接加减会超时。

#include<iostream>
#define ll __int64
using namespace std;
ll x0,y0;
ll exgcd(ll a,ll b)
{
	if(b==0)
	{
		x0=1,y0=0;
		return a;
	}
	ll gcd=exgcd(b,a%b);
	ll temp=x0;
	x0=y0,y0=temp-a/b*y0;
	return gcd;
}
int main()
{
	ll x,y,m,n,l,a,b,c;
	while(~scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&l))
	{
		a=m-n,b=(-1)*l,c=y-x;
		ll gcd=exgcd(a,b);
		if(c%gcd!=0) printf("Impossible\n");
		else
		{
			x0*=c/gcd;
			x0%=b/gcd;  //不然TLE
			if(b/gcd>0)
			while(x0<0)  x0+=b/gcd;
			printf("%I64d\n",x0);
		}
	}
	return 0;
}