ACM数论基础之扩展欧几里德详细证明.docxVIP

ACM数论基础之扩展欧几里德算法欧几里德算法概述：　欧几里德算法又称辗转相除法，用于计算两个整数a,b的最大公约数。其计算原理依赖于下面的定理：　gcd函数就是用来求(a,b)的最大公约数的。　gcd函数的基本性质：　gcd(a,b)=gcd(b,a)=gcd(-a,b)=gcd(|a|,|b|)欧几里得算法的公式表述　gcd(a,b)=gcd(b,a mod b)　证明：a可以表示成a = kb + r，则r = a mod b　假设d是a,b的一个公约数，则有　d|a, d|b，而r = a - kb，因此d|r　因此d是(b,a mod b)的公约数　假设d 是(b,a mod b)的公约数，则　d | b , d |r ，但是a = kb +r　因此d也是(a,b)的公约数　因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必然相等，得证欧几里德算法的C++语言描述　int Gcd(int a, int b)　{　if(b == 0)　return a;　return Gcd(b, a % b);　}　当然你也可以写成迭代形式：　int Gcd(int a, int b)　{　while(b != 0)　{　int r = b;　b = a % b;　a = r;　}　return a;　}扩展欧几里德定理　对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。c++语言实现　#include iostream　using namespace std;　int x,y,q;　void extend_Eulid(int a,int b)　{　if(b == 0)　{　x = 1;y = 0;q = a;　}　else　{　extend_Eulid(b,a%b);　int temp = x;　x = y;　y = temp - a/b*y;　}　}　int main()　{　int a,b;　cinab;　if(a b)　swap(a,b);　extend_Eulid(a,b);　printf(%d=(%d)*%d+(%d)*%d\n,q,x,a,y,b);　return 0;　}求解 x，y的方法的理解　设 ab。　1，显然当 b=0，gcd（a，b）=a。此时 x=1，y=0；　2，ab0 时　设 ax1+by1=gcd(a,b);　bx2+(a mod b)y2=gcd(b,a mod b);　根据朴素的欧几里德原理有 gcd(a,b)=gcd(b,a mod b);　则:ax1+by1=bx2+(a mod b)y2;　即:ax1+by1=bx2+(a-(a/b)*b)y2=ay2+bx2-(a/b)*by2;　根据恒等定理得：x1=y2; y1=x2-(a/b)*y2;　这样我们就得到了求解 x1,y1 的方法：x1，y1 的值基于 x2，y2.　上面的思想是以递归定义的，因为 gcd 不断的递归求解一定会有个时候 b=0，所以递归可以　结束。扩展欧几里德算法　扩展欧几里德算法是用来在已知a, b求解一组x，y使得ax+by = Gcd(a, b) =d(解一定存在，根据数论中的相关定理)。扩展欧几里德常用在求解模线性方程及方程组中。下面是一个使用C++的实现：　　int exGcd(int a, int b, int x, int y)　{　if(b == 0)　{　x = 1;　y = 0;　return a; ---很难找出一个这么实现的价值，因为扩展欧几里得还有更大的用途；个人认为定义全局数组更好,不用return r。　}　int r = exGcd(b, a % b, x, y);　int t = x;　x = y;　y = t - a / b * y;　return r;　}　把这个实现和Gcd的递归实现相比，发现多了下面的x,y赋值过程，这就是扩展欧几里德算法的精髓。　可以这样思考:　对于a = b, b = a % b 而言，我们求得 x, y使得 ax + by = Gcd(a, b)　由于b = a % b = a - a / b * b (注：这里的/是程序设计语言中的除法)　那么可以得到:　ax + by = Gcd(a, b) ===　bx + (a - a / b * b)y = Gcd(a, b) = Gcd(a, b) ===　ay +b(x - a / b*y) = Gcd(a, b)　因此对于a和b而言，他们的相对应的p，q分别是 y和(x-a/b*y)使用扩展欧几里德算法解决不定方程的办法　对于不定整数方程pa+qb=c，若 c mod Gcd(a, b)=0,则