欧几里得原理及扩展欧几里得原理(Euclidean Theory and Extended Euclidean Theory)学习笔记

题记：这是我第四次复习扩展欧几里得原理，因为不常用到，想要使用的时候又想不起细节，总是要查资料，于是索性这次整理一下欧几里得原理及其扩展原理存档至博客以备查用。

一、欧几里得原理

欧几里得原理(Euclidean Theory)是初等数论中求两正整数最大公约数(Greatest Common Divisor, GCD)的方法。欧几里得原理在中国古代又称“辗转相除法”，这一称法揭示了其求最大公约数的过程。

对于两个正整数a,b，记其最大公约数为gcd (a,b)。

那么我们有 gcd (a,b) = gcd (b,a%b)

即a与b的最大公约数也是a除以b的余数与b的最大公约数。

欧几里得原理可以通过递归实现，下面给出其实现的C代码：

int gcd(int a, int b)
{if(b==0) return a;return gcd(b,a%b);
}

函数返回参数a,b的最大公约数。注意该gcd函数调用时参数a,b不应该同时为0，否则该函数将不存在实际意义。

二、扩展欧几里得原理

欧几里得原理的扩展版本最经典的用法在于求二元一次不定方程的整数解。

（1）先来看一个最简单的二元一次不定方程：

a*x + b*y = 1 (*)，且 gcd (a,b) = 1，设它的一组整数解为 x = x1， y = y1。

对于另一个二元一次不定方程：

b*x + (a%b)*y = 1 (**)，设它的一组整数解为 x = x2， y = y2。

那么我们可以得到 a*x1 + b*y1 = b*x2 +(a%b)*y2 = 1，又因为 a%b = a - a/b*b，带入前式，

我们可以得到 a*x1 + b*y1 = a*y2 + b*(x2-a/b*y2)，由于a,b是任意的互质正整数，故可由该式得到方程(*)与方程(**)的解的关系如下：

x1 = y2，

y1 = x2-a/b*y2，即方程(*)的解可由方程(**)的解得到。

若我们想要求得a*x + b*y = 1的解，可以先求b*x + (a%b)*y = 1的解。我们可以发现这两个方程其未知量x,y 的系数之间是“辗转相除”的关系。那么最终我们只需要求gcd (a,b)*x + 0*y = 1的解（这里gcd (a,b) = 1，其解为 x = 1，y任意，一般取0），逆推之就可得到原方程的解。而逆推的过程就是扩展欧几里得原理的过程。

（2）我们现在已经解决了求二元一次不定方程 a*x + b*y = 1, gcd (a,b) = 1的解。下面考察二元一次不定方程的一般情况：

对于一般的二元一次方程 m*x + n*y = t，

若gcd (m,n) | t，亦即 t%gcd(m,n) == 0，方程的解将与方程 (m/gcd(m,n))*x + (n/gcd(m,n))*y = t/gcd(m,n) 的解相同。记新方程为 a*x + b*y = d，那么其解将是方程 a*x + b*y = 1的d倍。而方程 a*x + b*y = 1, gcd(a,b) = 1的解在（1）中已经求得。
否则方程无整数解

上述过程即用扩展欧几里得原理求一般二元一次不定方程整解的过程。

欧几里得原理及扩展欧几里得原理(Euclidean Theory and Extended Euclidean Theory)学习笔记相关推荐

欧几里得定理与扩展欧几里得定理
欧几里的定理(辗转相除法): gcd(a,b) = gcd(b,a%b) gcd(a,b)表示a,b的最大公约数证明: 设 a > b c = a%b a = k * b + c (k为某个整 ...
欧几里得定理与扩展欧几里得
3,欧几里德定理:(射影定理) 定理指出素数是无限的. a*b*c+1要么是素数要么其质因子就是素数. 扩展欧几里得: 扩展欧几里得算法是欧几里得(又叫辗转相除法)的扩展.已知整数a.b,扩展欧几里得 ...
扩展欧几里得及其应用
欧几里得算法欧几里德算法又称辗转相除法,用于计算两个整数a,b的最大公约数.其计算原理依赖于下面的定理:gcd(a,b)=gcd(a,a mod b) gcd函数就是用来求(a,b)的最大公约数的. ...
扩展欧几里得定理的证明和代码
1.欧几里得算法,gcd(a,b)为a b(a>b)的最大公约数,则gcd(a,b) = acd(b, a%b) 利用这个定理我们可以反复对ab模下去求得a和b的最大公约数代码如下 int G ...
Strange Optimization(扩展欧几里得)
Strange Optimization Accepted : 67 Submit : 289 Time Limit : 1000 MS Memory Limit : 65536 KB Str ...
欧几里得算法扩展欧几里得
本文包括: 1.几里德算法递归实现 2.扩展欧几里德算法递归实现 -实际上就是把别人总结的,我认为有助于自己理解的内容copy过来,再加上几句自己的理解. 欧几里得算法欧几里德算法又称辗转相除法 ...
【欧几里得扩展欧几里得】
欧几里得 LL gcd(LL a,LL b){return (b==0) ? a : gcd(b,a%b);} 扩展欧几里得 int ex_gcd(int a,int b,int &x,int ...
欧几里得算法（扩展欧几里得、欧拉定理、费马小定理）
问题描述: a 和 b的最大公约数是多少? 古代解法:辗转相除法迭代过程:例如: {a = 15 和 b = 12 } =>{ a = 12,b = 15 - (15/12)* 12 = ...
扩展欧几里得定理新手向证明及代码
知识储备 1 . 朴素欧几里得原理:gcd(a,b) == gcd(b,a % b) 2 . 负数取模:忽略符号返回绝对值就好了 3 . 模数原理:对于整数a,b必然存在整数k使得a % b == a ...

欧几里得原理及扩展欧几里得原理(Euclidean Theory and Extended Euclidean Theory)学习笔记

欧几里得原理及扩展欧几里得原理(Euclidean Theory and Extended Euclidean Theory)学习笔记相关推荐

最新文章

热门文章