青蛙的约会与欧几里德定理

参考出处：http://hi.baidu.com/newmyl/blog/item/d4e7ac07227afcc87a894735.html/cmtid/b5e81360995e8449eaf8f823

青蛙的约会：

在网上看了很多关于不定方程方程求解的问题，可都没有说全，都只说了一部分，看了好多之后才真正弄清楚不定方程的求解全过程，步骤如下：

　　求a * x + b * y = n的整数解。

　　1、先计算Gcd(a,b)，若n不能被Gcd(a,b)整除，则方程无整数解；否则，在方程两边同时除以Gcd(a,b)，得到新的不定方程a' * x + b' * y = n'，此时Gcd(a',b')=1;

    2、利用上面所说的欧几里德算法求出方程a' * x + b' * y = 1的一组整数解x0,y0，则n' * x0,n' * y0是方程a' * x + b' * y = n'的一组整数解；

　　3、根据数论中的相关定理，可得方程a' * x + b' * y = n'的所有整数解为：

        x = n' * x0 + b' * t
y = n' * y0 - a' * t
(t为整数)

　　　 上面的解也就是a * x + b * y = n 的全部整数解。

#include<iostream> #include<cmath> using namespace std; int Gcd(int a, int b) { while(b) { int m = a; a = b; b = m%b; } return a; } void Extend_gcd(int a, int b, int &x, int &y) { int tmp; if(b==0) { x = 1; y = 0; return; } else { Extend_gcd(b,a%b,x,y); tmp = x; x = y; y = tmp - (a/b)*y; } } int main() { int x,y,m,n,L,step,loops,t; cin>>x>>y>>m>>n>>L; int a = n - m; int c = x - y; int gcd = Gcd(a,L); if(c%gcd) { cout<<"Impossible"<<endl; } else { a /= gcd; L /= gcd; c /= gcd; Extend_gcd(a,L,step,loops); t=c*step/L;//注1 step=c*step-t*L; if(step<0) { step += L; } cout<<step<<endl; } return 0; }

注1：（原注1做废）

              此时方程的所有解为：x=c*k1-b*t,x的最小的可能值是0，令x=0可求出当x最小时的t的取值，但由于x=0是可能的最小取值，实际上可能x根本取不到0，那么由计算机的取整除法可知：由 t=c*k1/b算出的t，代回x=c*k1-b*t中，求出的x可能会小于0，此时令t=t+1，求出的x必大于0；如果代回后x仍是大于等于0的，那么不需要再做修正。