poj 1061 扩展欧几里得解同余方程（求最小非负整数解）

题目可以转化成求关于t的同余方程的最小非负数解：

x+m*t≡y+n*t (mod L)

该方程又可以转化成：

k*L+(n-m)*t=x-y

利用扩展欧几里得可以解决这个问题：

eg：对于方程ax+by=c

设tm=gcd(a,b)

若c%tm!=0，则该方程无整数解。

否则，列出方程：

a*x0+b*y0=tm

易用extend_gcd求出x0和y0

然后最终的解就是x=x0*(c/tm)，y=y0*(c/tm)

注意：若是要求最小非负整数解？

例如求y的最小非负整数解，

令r=a/tm，则y=(y%r+r)%r;

 #include <iostream>
 using namespace std;
 __int64 x,y,m,n,L;

 __int64 extend_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y){
     if (b==){
         x=;y=;
         return a;
     }
     else{
         __int64 r=extend_gcd(b,a%b,y,x);
         y=y-x*(a/b);
         return r;
     }
 }

 int main()
 {
     cin>>x>>y>>m>>n>>L;

     __int64 k,t;
     __int64 tm=extend_gcd(n-m,L,t,k);
     if ((x-y)%tm!=)
         cout<<"Impossible"<<endl;
     else
     {
         __int64 ans=t*((x-y)/tm);
         __int64 r=L/tm;
         ans=(ans%r+r)%r;    //求出最小非负整数解
         //while (ans<0)   ans+=k/tm;　　　　//这样做是错的= =
         cout<<ans<<endl;
     }

     return ;
 }

Reference：http://www.cnblogs.com/yueshuqiao/archive/2011/08/23/2150960.html