对拓展欧几里得算法的一点理解

首先需要明确的一点是:这是一种算法,而非一个证明题。

算法的需求与数学证明题是不一样的,数学证明题要求严谨完整,而算法只需要证明我用到的某个的性质成立即可,相当于是“恰好发现了这一点”。

于是对于拓展欧几里得,我们是从欧几里得算法中发现了一个递推的性质,从而受到启发,产生猜想:可不可以利用递推求出二元一次方程的解?

我们把猜想建立在欧几里得算法之上,利用该算法的递推过程,贯穿该过程来得到想要解决问题的答案。

也就是说,我们仅仅需要证明我们得到的答案是正确的,而并非深刻挖掘欧几里得算法的过程。

拓展欧几里得算法及证明如下:

ax+by=gcd(a,b)

bx`+{a/b}by`=gcd(b,{a/b}b)

由于gcd(a,b)=gcd(b,{a/b}b),

得到ax+by=bx`+{a/b}by`

∵{a/b}=a/b-[a/b]

∴ax+by=bx`+(a/b-[a/b])by`

∴ax+by=bx`+(a-[a/b]b)y`

∴a(x-y`)=b(x`-y-[a/b]y`)

由于当x=y`,y=[a/b]y`-x`时该等式一定成立

//注意这里并不代表x=y`,y=[a/b]y`-x`是该等式的唯一解,而是我们取出了在实际情况中所需要的解

所以我们主观上取这一组解做我们的递归项

//因为这样可以解决我们的问题,可以求解方程

于是我们得到了这个递归公式

★那么说到底我们为什么会想到利用这个方法?

因为当b=0时这个方程的根是已知的x=c,y=任意值.

我们要通过辗转相除法求出与原方程有关联的一个满足b=0的方程。

这里的关联就是递推式,先递到底得到满足b=0的方程,

再把这个b=0的方程的根带入递推式回归到原方程,得到原方程的根。

原文地址:https://www.cnblogs.com/AppleC/p/11332872.html

时间: 2024-10-06 00:07:07

对拓展欧几里得算法的一点理解的相关文章

HDU - 1356 The Balance(拓展欧几里得算法的解空间结构)

题目: Description Ms. Iyo Kiffa-Australis has a balance and only two kinds of weights to measure a dose of medicine. For example, to measure 200mg of aspirin using 300mg weights and 700mg weights, she can put one 700mg weight on the side of the medicin

拓展欧几里得算法及代码实现

扩展欧几里得算法就是求: ax + by = gcd(a, b) 的一组整数解(x, y) 一.非递归的实现: 首先看a = 60, b = 22的情况: 表格左边是欧几里得算法,右边等式计算ax + by = gcd(a, b)的解 a = 2 × b + 16 16 = a - 2b b = 1 × 16 + 6 6 = b - 1 × 16 = b - 1 × (a - 2b) = -a + 3b 16 = 2 × 6 + 4 4 = 16 - 2 × 6 = (a - 2b) - 2 ×

数学问题——拓展欧几里得算法

一.拓展欧几里得算法 该算法用来解决这样一个问题:给定两个非零整数 a 和 b,求一组整数解 (x,y) ,使得 ax + by = gcd(a,b) 成立,其中 gcd(a,b) 表示 a 和 b 的最大公约数. 递归边界:当 b 为 0 时,此时的 a 就等于 gcd,显然有 a*1+b*0=gcd 成立,此时 x=1,y=0: 递推公式:设当计算 gcd(a,b) 时,有 ax1 + by1 = gcd 成立:而在下一步计算gcd(b,a%b) 时,又有 bx2 + (a%b)y2 = g

如何使用拓展欧几里得算法求解模线性方程组(详解)

式子a≡b(mod n)称为a和b关于模n同余,它的充要条件是a-b是n的整数倍,即a-b=zn(其中z取整数). 而模线性方程组ax≡b(mod n)可以写成ax-b=zn(其中z取整数),移项可得 ax-zn=b,也即二元一次方程ax+by=c的形式,利用拓展欧几里得算法(extgcd)可以求解该方程是否有解及其一组解,并可根据该组解写出解系,进而求出一个特解,比如最小正整数解. 下面给出拓展欧几里得算法的程序. 1 typedef long long LL; 2 void extgcd(L

Sumdiv|同余|约数|拓展欧几里得算法

目录 Sumdiv|同余|约数|拓展欧几里得算法 Problem 分析 约数个数定理部分 约数和定理部分 等比数列部分 题目分析 扩展欧几里得算法部分 Code 呕,我吐了. Sumdiv|同余|约数|拓展欧几里得算法 Problem \[ 求A^{B}的所有约数之和 \ mod \ 9901\left(1\leqslant A,B \leqslant 5*10^{7}\right) \] 分析 约数个数定理部分 定理内容: 对于一个大于1的正整数n可以分解质因数: 则n的正约数个数为: 定理证

[数论]拓展欧几里得算法

欧几里得算法(辗转相除法) 用来求解最大公约数 1 int gcd(int a,int b){ 2 return b ? gcd(b,a%b) : a; 3 } 在 #include<algorithm> 中也可以直接调用 __gcd(a,b) 拓展欧几里得算法 求解不定方程: 引理:存在 x , y 使得 ax+by=gcd(a,b) 设a,b,c为任意整数,若方程ax+by=c的一组解是(x0,y0),则它的任意整数解都可以写成(x0+k*b/gcd(a,b),y0-k*a/gcd(a,b

欧几里得 &amp; 拓展欧几里得算法 讲解 (Euclid &amp; Extend- Euclid Algorithm)

欧几里得& 拓展欧几里得(Euclid & Extend-Euclid) 欧几里得算法(Euclid) 背景: 欧几里德算法又称辗转相除法,用于计算两个正整数a,b的最大公约数.                                                                                                                                               --百度百科 代码: 递推的代

欧几里得 &amp;amp; 拓展欧几里得算法 解说 (Euclid &amp;amp; Extend- Euclid Algorithm)

欧几里得& 拓展欧几里得(Euclid & Extend-Euclid) 欧几里得算法(Euclid) 背景: 欧几里德算法又称辗转相除法.用于计算两个正整数a.b的最大公约数. --百度百科 代码: 递推的代码是相当的简洁: int gcd(int a,int b) { return b == 0 ? a : gcd(b, a % b); } 分析: 方法说了是辗转相除法,自然没有什么好介绍的了. . Fresh肯定会认为这样递归下去会不会爆栈?实际上在这里是不会爆栈的,由于递归的层数是

拓展欧几里得算法求不定方程

对于 ax+by=gcd(a,b) 这样的方程,可以用扩展欧几里得算法exgcd求出一组通解. 根据欧几里得求gcd: gcd(a,b)=gcd(b,a%b) 可得 bx+(a%b)y=gcd(b,a%b) 根据 a%b=a−(a/b)∗b 可得 bx+ay−(a/b)b∗y=gcd(b,a%b) 化简得 ay+b(x−(a/b)y)=gcd(b,a%b) x′=y,y′=(x−(a/b)y) ax′+by′=gcd(b,a%b)<=>ax+by=gcd(a,b) 根据 gcd(a,0)=a