poj 1061 青蛙的约会 (扩展gcd， 模线性方程)

好长时间没贴代码了，这段时间看了些数论，水了几道题，过一段时间把以前的题好好整理下。唉，算法学不成器，技术又不过关，菜鸟都要哭死了

这道题使用模线性方程做的，注意应该使用 long long 或 __int64数据类型，具体的exgcd 函数 和 解模线性方程 可以参照 《算法导论》。

//
#include <iostream>
#include <math.h>
using namespace std;
__int64 dis,step;

__int64 exgcd(__int64 a,__int64 b,__int64 & x,__int64 & y)// a > b
{
	if(0==b)
	{
		x = 1;
		y = 0;
		return a;
	}
	else
	{
		__int64 xt,yt,ret = 0;
		ret = exgcd(b,a % b,xt,yt);
		x = yt;
		y = xt - (a / b ) * yt;
		return ret;
	}
}

void MLES(__int64 a,__int64 b,__int64 n)// ax = b(mod n)
{
	__int64 x, y ,d;
	d = exgcd(a,n,x,y);
	if(0== b % d)
	{
		__int64 xt = x * (b / d) % n;
		while(xt < 0)
			xt += (n / abs((double)d));
		cout<<xt<<endl;
	}
	else
		cout<<"Impossible"<<endl;
}

int main()
{
	__int64 px,py,m,n,L;
	while(cin>>px>>py>>m>>n>>L)
	{
		 dis = (py - px) % L;//第二只 相对于第一只的 距离
		 step = (m - n) % L;//第一只 的步进速度
		if(dis==0)
			cout<<0<<endl;
		else if(0==step)
			cout<<"Impossible"<<endl;
		else
			MLES(step,dis,L);
	}
	return 0;
}