> __int64 MLE(__int64 a,__int64 b,__int64 n)
> {
> 	__int64 x,y,x0;
> 	__int64 d=Extended_Euclid(a,n,x,y);
> 	if(b%d)
> 		return -1;
> 	
> 	x0=x*(b/d)%n+n;
> 	return x0%(n/d);
> 	
> 	
> }
> 
> 
> x0=x*(b/d)%n+n;
> 	return x0%(n/d);
> 就是这一步，不能严格证明，有点晕

算法导论上有详细证明。。。