标签：style blog http io color os ar sp strong
题目可以转化成求关于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;
1 #include <iostream> 2 using namespace std; 3 __int64 x,y,m,n,L; 4 5 __int64 extend_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y){ 6 if (b==0){ 7 x=1;y=0; 8 return a; 9 } 10 else{ 11 __int64 r=extend_gcd(b,a%b,y,x); 12 y=y-x*(a/b); 13 return r; 14 } 15 } 16 17 int main() 18 { 19 cin>>x>>y>>m>>n>>L; 20 21 __int64 k,t; 22 __int64 tm=extend_gcd(n-m,L,t,k); 23 if ((x-y)%tm!=0) 24 cout<<"Impossible"<<endl; 25 else 26 { 27 __int64 ans=t*((x-y)/tm); 28 __int64 r=L/tm; 29 ans=(ans%r+r)%r; //求出最小非负整数解 30 //while (ans<0) ans+=k/tm; //这样做是错的= = 31 cout<<ans<<endl; 32 } 33 34 return 0; 35 }
Reference：http://www.cnblogs.com/yueshuqiao/archive/2011/08/23/2150960.html
poj 1061 同余方程+扩展欧几里得+（求最小非负整数解）
标签：style blog http io color os ar sp strong
原文地址：http://www.cnblogs.com/pdev/p/4067571.html
