poj 1061 exgcd解同余方程

点击打开链接http://poj.org/problem?id=1061

#include <iostream>
#include <algorithm>
using namespace std;
__int64 t,p,d;
void exgcd(__int64 a,__int64 b)
{
	if(b==0)
	{
		t=1;
		p=0;
		d=a;//gcd
	}
	else
	{
		exgcd(b,a%b);
		__int64 tmp;
		tmp=t;
		t=p;
		p=tmp-(a/b)*p;	
	}
}
int main()
{
	__int64 x,y,m,n,l;    // t秒后能相遇  (x+mt)modl = (y+nt)modl 
	cin>>x>>y>>m>>n>>l; //  即x+mt ≡y+nt (mod l) 
						//  pl=x+mt-y-nt ->  (n-m)t+pl=x-y; 
										//可以先求出 (n-m)t1+lp1=gcd(n-m,l) 
	exgcd(n-m,l);							//如果 d|(x-y) 则 t=t1*((x-y)/d) 
	if((x-y)%d)										//	t=(t%(l/d)+l/d)%(l/d);  (没看懂) 
	{
		cout<<"Impossible"<<endl;
	}
	else
	{
		__int64 k=(x-y)/d;
		
		t=t*k;
		t=(t%(l/d)+l/d)%(l/d); // 不懂... 
		cout<<t<<endl;
	}
	return 0;
}