实现扩展欧几里德算法实验报告
实现扩展欧几里德算法实验报告
一、 扩展欧几里德算法问题简介
对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，
必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。
二、 扩展欧几里德算法设计方法简介
欧几里德算法又称辗转相除法，用于计算两个整数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)的公约数是一样的，其最大公约数也必然相等，得证
三、 程序代码
#include
#include
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;
cout<<"请输入a"<
cout<<"请输入b"<
if(a < b)
{
int temp=a;
a=b;
b=temp;
}
extend_Eulid(a,b);
printf("最大公约数=%d=(%d)*%d+(%d)*%d\n",q,x,a,y,b);
return 0;
}
四、 算法介绍
扩展欧几里德算法是用来在已知a, b求解一组x，y使得ax+by = Gcd(a, b)
=d(解一定存在，根据数论中的相关定理)。扩展欧几里德常用在求解模线性方程
及方程组中。
五、 实验数据
a=22
b=33
六、 实验结果
七、 实验体会
欧几里得算法应该算是在这两个实验中比较简单的一个程序，但是就是这样
的算法，我在实际操作中也遇到了很多的问题，有些是输入法问题，有些是运算
问题，所幸的是最后终于完成了。
扩展欧几里得算法
序号 i 商 q i −1余数 r i s i t i 0240101 4601扩展欧⼏⾥得算法扩展欧⼏⾥得算法 扩展欧⼏⾥得算法（英语：Extended Euclidean algorithm ）是欧⼏⾥得算法（⼜叫辗转相除法）的扩展。
已知整数a 、b ，扩展欧⼏⾥得算法可以在求得a 、b 的最⼤公约数的同时，能找到整数x 、y （其中⼀个很可能是负数），使它们满⾜贝祖等式： ax + by = gcd(a, b)如果a 是负数，可以把问题转化成: |a|(-x) + by = gcd(|a|, x) （|a|为a 的绝对值），然后令x ' = (-x)。
通常谈到最⼤公约数时，我们都会提到⼀个⾮常基本的事实（贝祖等式给定）：给定⼆个整数a 、b ，必存在整数x 、y 使得ax + by =gcd(a,b)。
众所周知，已知两个数a 和b ，对它们进⾏辗转相除，可得它们的最⼤公约数。
不过，在欧⼏⾥得算法中，我们仅仅利⽤了每步带余除法所得的余数。
扩展欧⼏⾥得算法还利⽤了带余除法所得的商，在辗转相除的同时也能得到贝祖等式中的x 、y 两个系数。
以扩展欧⼏⾥得算法求得的系数是满⾜裴蜀等式的最简系数。
另外，扩展欧⼏⾥得算法是⼀种⾃验证算法，最后⼀步得到的s i+1和t i+1(s i+1和t i+1的含义如下)乘以gcd(a,b)之后恰好是a 和b,可以⽤来验证计算结果是否正确。
扩展欧⼏⾥得算法可以⽤来计算模反元素(也叫模逆元)，求出模反元素是RSA 加密算法中获得所需公钥、私钥的必要步骤。
算法和举例在标准的欧⼏⾥得算法中，我们记欲求最⼤公约数的两个数为a 和b ，第i 步带余除法得到的商为q i ，余数为r i+1，则欧⼏⾥得算法可以写成如下形式： r 0 = a r 1 = b ....... r i+1 = r i-1 - q i r i 且 0 ≤ r i+1 < |r i | ......当某步得到的 r i+1 = 0 时，计算结束。
扩展欧几里得算法与线性同余方程的求解
扩展欧几里得算法与线性同余方程的求解扩展欧几里得算法（Extended Euclidean Algorithm）是欧几里得算法（Euclidean Algorithm）的一个扩展版本，它不仅用于计算两个整数的最大公约数（GCD），还能同时找到满足贝祖等式（Bézout's Identity）的整数解。
这个等式是，其中a和b是任意两个整数，x和y是待求的整数。
当我们要用这个算法来求解线性同余方程ax≡c(mod b)时，我们首先需要找到一个整数x使得axmod b=c。
这可以转化为求解ax−by=c，其中y是另一个待求的整数。
如果不能整除c，则方程ax≡c(mod b)无解。
否则，我们可以先使用扩展欧几里得算法找到x′和y′满足，然后通过调整x′来找到x，使得ax≡c(mod b)成立。
扩展欧几里得算法的基本步骤1.初始化：设x2=1,y2=0,x1=0,y1=1，并且r1=a,r2=b（其中a是被除数，b是除数）。
2.迭代过程：当r2≠0时，执行以下步骤：3.计算商（向下取整）。
4.计算余数r=r1mod r2。
5.更新x和y：x=x1−q∙x2,y=y1−q∙y2。
6.更新r1,r2,x1,x2,y1,y2为r2,r,x2,x,y2,y。
7.终止条件：当r2=0时，算法结束。
此时r1就是a和b的最大公约数，而x1和y1就是满足的整数解。
求解线性同余方程一旦你有了满足的x1和y1，你可以通过以下方式找到x：注意，这里使用了模b来确保x是一个在模b意义下的有效解。
这是因为如果，那么乘以后，等式的左边就变成了ax（其中），而等式的右边则变成了c。
然而，由于我们只关心模b的结果，因此需要对x取模b。
但是，更常见的做法是使用模的逆元来避免直接的除法运算（这在b不是质数时可能不可用）。
如果b和互质，则可以使用扩展欧几里得算法找到关于模b的逆元c′，然后计算x=x1∙c′mod b。
然而，在这个特定的情况下，由于我们已经知道能整除c，我们可以简单地通过乘以并取模b来找到解，而不需要显式地计算逆元。
辗转相除法实验报告
一、实验目的1. 理解辗转相除法的原理。
2. 掌握使用辗转相除法求解两个正整数最大公约数的方法。
3. 通过实验加深对数学理论知识的理解。
二、实验原理辗转相除法，也称欧几里得算法，是一种用于求两个正整数a和b（a>b）的最大公约数（GCD）的算法。
其基本原理是：用较大数a除以较小数b，得到余数r1，再用b除以r1，得到余数r2，如此重复，直到余数为0，此时的除数即为最大公约数。
三、实验器材1. 计算机2. 编程软件（如Python、C++等）3. 笔记本四、实验步骤1. 输入两个正整数a和b：确保a>b，否则交换a和b的值。
2. 初始化余数r为a。
3. 循环执行以下步骤：a. 计算r除以b的余数，赋值给r。
b. 如果r等于0，则b即为最大公约数，结束循环。
c. 如果r不等于0，则将b赋值给r，b赋值给新的临时变量。
4. 输出最大公约数。
五、实验代码（Python）```pythondef gcd(a, b):while b != 0:r = a % ba = bb = rreturn a# 输入两个正整数a = int(input("请输入第一个正整数a："))b = int(input("请输入第二个正整数b（a>b）："))# 调用gcd函数计算最大公约数result = gcd(a, b)# 输出结果print("最大公约数为：", result)```六、实验结果与分析1. 当输入a=56，b=48时，程序输出最大公约数为8。
2. 当输入a=100，b=50时，程序输出最大公约数为50。
3. 当输入a=101，b=102时，程序输出最大公约数为1。
通过实验，我们可以看到辗转相除法能够快速且准确地计算出两个正整数的最大公约数。
该算法在数学领域有着广泛的应用，如密码学、计算机科学等。
七、实验总结本次实验让我们深入理解了辗转相除法的原理，并掌握了使用Python编程实现该算法的方法。
python扩展欧几里得算法求逆元
python扩展欧几里得算法求逆元以下是使用Python实现扩展欧几里得算法求逆元的代码：```pythondef extended_gcd(a, b):if b == 0:return a, 1, 0else:gcd, x, y = extended_gcd(b, a % b)return gcd, y, x - (a // b) * ydef modular_inverse(a, m):gcd, x, _ = extended_gcd(a, m)if gcd != 1:# 如果a和m不互质，则不存在逆元return Noneelse:return (x % m + m) % m# 示例a = 7m = 26inverse = modular_inverse(a, m)if inverse is not None:print(f"The modular inverse of {a} (mod {m}) is {inverse}.") else:print(f"The modular inverse of {a} (mod {m}) does not exist.") ```上述代码定义了两个函数：`extended_gcd`用于计算扩展欧几里得算法，`modular_inverse`用于求取逆元。
在示例中，我们求取了数字7在模26下的逆元。
如果逆元存在，则打印出逆元的值；如果逆元不存在（即a和m不互质），则打印出不存在的提示。
请注意，逆元的计算要求a和m是互质的（即它们的最大公约数为1）。
如果a和m不互质，则不存在逆元。
拓展的欧几里得算法求逆元c++
拓展的欧几里得算法求逆元c++拓展的欧几里得算法用于求解两个数a和b的最大公约数，同时可以求解出满足a * x + b * y = gcd(a, b)的整数x和y，其中a和b互质时，x即为b的乘法逆元。
以下是一个用C++实现的拓展的欧几里得算法，求解逆元的函数为`modular_inverse`：```cpp#include <iostream>using namespace std;// 拓展的欧几里得算法int extended_euclidean_algorithm(int a, int b, int &x, int &y) {if (b == 0) {x = 1;y = 0;return a;}int x1, y1;int gcd = extended_euclidean_algorithm(b, a % b, x1, y1);x = y1;y = x1 - (a / b) * y1;return gcd;}// 求解逆元int modular_inverse(int a, int m) {int x, y;int gcd = extended_euclidean_algorithm(a, m, x, y);if (gcd != 1) {throw invalid_argument("Inverse doesn't exist");}return (x % m + m) % m;}int main() {int a, m;cout << "Enter a: ";cin >> a;cout << "Enter m: ";cin >> m;try {int inverse = modular_inverse(a, m);cout << "Inverse of " << a << " (mod " << m << ") is: " << inverse << endl;} catch (const invalid_argument& e) {cout << e.what() << endl;}return 0;}```在这个代码示例中，用户需要输入一个整数a和模数m。
欧几里德算法及其扩展
欧几里德算法欧几里德算法又称辗转相除法，用于计算两个整数a,b 的最大公约数。
基本算法：设a=qb+r,其中a, b, q, r都是整数，则gcd(a,b)=gcd(b,r)，即gcd(a,b)=gcd(b,a%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)的公约数，贝Ud | b , d |r ,但是a = kb +r因此d也是(a,b)的公约数因此(a,b)和(b,a modb)的公约数是一样的，其最大公约数也必然相等，得证第二种证明：要证欧几里德算法成立,即证: gcd(a,b)=gcd(b,r), 其中gcd 是取最大公约数的意思, r=a mod b下面证gcd (a, b)=gcd(b, r)设 c 是a, b的最大公约数，即c=gcd (a, b),贝U有a=mc, b=nc,其中m, n 为正整数，且m，n 互为质数由r= a mod b 可知，r= a- qb 其中，q 是正整数，则r=a-qb=mc-qnc= ( m-qn) cb=nc,r=(m-qn)c ,且n, (m-qn)互质(假设n, m-qn不互质，贝U n=xd, m-qn=yd其中x,y,d 都是正整数，且d>1贝a=mc=(qx+y)dc, b=xdc,这时a,b的最大公约数变成dc,与前提矛盾，所以n , m-qn —定互质)贝gcd( b,r ) =c=gcd(a,b )4 {5 int r = b; 得证。
算法的实现：最简单的方法就是应用递归算法，代码如下：1 int gcd( int a, int b)2 {3if (b==0) 4return a; 5return 6gcd(b,a%b);7 } 代码可优化如下：1 int gcd( int a, int b)2 {3return b ? gcd(b,a%b) : a;4 }当然你也可以用迭代形式： 1 int Gcd( int a, int b)2 {3while (b != 0) 6b = a % b; 7a = r;8 }a;9 return10 }扩展欧几里德算法基本算法：对于不完全为0的非负整数a , b, gcd (a, b)表示a , b的最大公约数，必然存在整数对x ，y ，使得gcd ( a，b) =ax+by。
扩展欧几里得算法详解
扩展欧⼏⾥得算法详解本篇将附上扩展欧⼏⾥得算法的思想与推导;对于⼀个⽅程a∗x+b∗y=gcd(a,b)来说,我们可以做如下的推导:设有a∗x1+b∗y1=gcd(a,b);同时我们有b∗x2+(a%b)∗y2=gcd(b,a%b);对于这个⽅程组,我们希望知道的是x1,x2,y1,y2之间的关系,这样我们才可以递归解决这个问题我们观察b∗x2+(a%b)∗y2这个式⼦,我们可以将(a%b)写作(a−⌊ab⌋∗b),将括号打开常数a,b合并,合并之后的结果为a∗y2+b∗(x2−⌊ab⌋∗y2))由于欧⼏⾥得算法的原理gcd(a,b)==gcd(b,a%b),我们将两式⼦联⽴,对⽐系数即可得到x1=y2,y1=x2−⌊ab⌋∗y2这个递归的边界是什么呢?我们知道,当朴素欧⼏⾥得到达边界时,return gcd(a,0)=a,那么边界条件就是对a∗x0+b∗y0=a求解,很显然,此时x0=1,y0=0当我们递归求出了⼀个⽅程的特解时,如何求出这个⽅程的通解呢?⽅程a∗x+b∗y=gcd(a,b)中,如果将x加上⼀个常数k1,y减去⼀个常数k2,仍然保持原⽅程成⽴,那么x+k1,y−k2就是⽅程的⼀个新解,这个k应该如何选择呢?实际上很简单,a∗(x+k1)+b∗(y+k2)=gcd(a,b),打开括号,a∗x+a∗k1+b∗y−b∗k2=gcd(a,b);我们保证原⽅程成⽴,就需要a∗k1==b∗k2,那么显然k1=b,k2=a是⼀种合理的情况,但是这样是⽆法包含所有整数解的,因为我们加上的这个值并⾮是最⼩值那我们应该加上什么值才⾏呢?我们发现当a∗k1==b∗k2=t∗lcm(a,b)可以保证得到所有解,于是每次寻找解就可以分别在x加上bgcd(a,b),在y减去agcd(a,b)就可以了对于⽅程a∗x+b∗y=c我们⼜该如何求解?我们发现如果(c%gcd(a,b)!=0)那么这个⽅程是⽆解的,⽽如果gcd(a,b)∗t==c,我们就可以按上⾯的⽅法求解之后对我们的解乘上⼀个t(t=c gcd(a,b))code:int exgcd(int a,int b,int &x1,int &y1) {if(!b){x1=1,y1=0;return a;}int x2,y2;int d=exgcd(b,a%b,x2,y2);x1=y2,y1=x2-(a/b)*y2;return d;}Processing math: 100%。
扩展欧几里得算法求负数
扩展欧几里得算法求负数扩展欧几里得算法（Extended Euclidean Algorithm）是一个用于求解两个整数的最大公约数（GCD）以及找到使得ax + by = GCD(a, b)成立的整数x和y 的算法。
在正整数情况下，欧几里得算法非常方便和易于理解，但是当我们需要在负数集合上运算时，情况会变得复杂一些。
在这篇文章中，我们将介绍如何使用扩展欧几里得算法来处理负数的情况，并且逐步解释算法的具体执行过程。
首先，我们需要了解欧几里得算法的原理以及如何处理正整数情况。
欧几里得算法基于一个简单的观察结果：两个整数a和b的最大公约数等于b和a mod b 的最大公约数。
换句话说，如果a能够被b整除，那么最大公约数就是b；否则，最大公约数是b和a mod b的最大公约数。
现在，我们将这个观察结果扩展到负数集合上。
假设我们想要求解两个负数a 和b的最大公约数gcd(a, b)。
首先，我们需要注意到，最大公约数不受负号的影响。
也就是说，gcd(a, b) = gcd( a , b )。
因此，我们可以将负数转化为其绝对值来简化计算。
接下来，我们需要找到使得ax + by = gcd(a, b)成立的整数x和y。
这个方程被称为贝祖等式（Bézout's identity）。
为了使用扩展欧几里得算法，我们需要事先计算a和b的最大公约数和贝祖等式的解。
首先，我们使用欧几里得算法来计算gcd( a , b )。
这一步和正整数情况没有任何区别。
我们重复使用Euclid's Algorithm，直到余数为0为止。
也就是说，我们不断用较小的数去除较大的数，并用较大的数对较小的数取余，直到余数为0。
此时，较大的数就是最大公约数。
假设我们得到的最大公约数是d。
接下来，我们要找到整数x和y，使得ax + by = d成立。
我们可以通过设置初始的x和y的值，然后在每一步通过更新x和y 的值来求解。
欧几里得算法与扩展欧几里得算法_C++
欧⼏⾥得算法与扩展欧⼏⾥得算法_C++注：以下讨论的数均为整数⼀、欧⼏⾥得算法（重点是证明，对后续知识有⽤） 欧⼏⾥得算法，也叫辗转相除，简称 gcd，⽤于计算两个整数的最⼤公约数 定义 gcd(a,b) 为整数 a 与 b 的最⼤公约数 引理：gcd(a,b)=gcd(b,a%b) 证明： 设 r=a%b , c=gcd(a,b) 则 a=xc , b=yc , 其中x , y互质 r=a%b=a-pb=xc-pyc=(x-py)c ⽽b=yc 可知：y 与 x-py 互质 证明： 假设 y 与 x-py 不互质 设 y=nk , x-py=mk , 且 k>1 （因为互质） 将 y 带⼊可得 x-pnk=mk x=(pn+m)k 则 a=xc=(pn+m)kc , b=yc=nkc 那么此时 a 与 b 的最⼤公约数为 kc 不为 k 与原命题⽭盾，则 y 与 x-py 互质 因为 y 与 x-py 互质，所以 r 与 b 的最⼤公约数为 c 即 gcd(b,r)=c=gcd(a,b) 得证 当a%b=0时，gcd(a,b)=b 这样我们可以写成递归形式1 inline int gcd(int a,int b)2 {3return b?gcd(b,a%b):a;4 }⼆、扩展欧⼏⾥得算法 扩展欧⼏⾥得算法，简称 exgcd，⼀般⽤来求解不定⽅程，求解线性同余⽅程，求解模的逆元等 引理：存在 x , y 使得 gcd(a,b)=ax+by 证明： 当 b=0 时，gcd(a,b)=a，此时 x=1 , y=0 当 b!=0 时， 设 ax1+by1=gcd(a,b)=gcd(b,a%b)=bx2+(a%b)y2 ⼜因 a%b=a-a/b*b 则 ax1+by1=bx2+(a-a/b*b)y2 ax1+by1=bx2+ay2-a/b*by2 ax1+by1=ay2+bx2-b*a/b*y2 ax1+by1=ay2+b(x2-a/b*y2) 解得 x1=y2 , y1=x2-a/b*y2 因为当 b=0 时存在 x , y 为最后⼀组解 ⽽每⼀组的解可根据后⼀组得到 所以第⼀组的解 x , y 必然存在 得证 根据上⾯的证明，在实现的时候采⽤递归做法 先递归进⼊下⼀层，等到到达最后⼀层即 b=0 时就返回x=1 , y=0 再根据 x=y’ , y=x’-a/b/y’ ( x’ 与 y’ 为下⼀层的 x 与 y ) 得到当层的解 不断算出当层的解并返回，最终返回⾄第⼀层，得到原解1 inline void exgcd(int a,int b)2 {3if (b)4 {5 exgcd(b,a%b);6int k=x;7 x=y;8 y=k-a/b*y;9 }10else y=(x=1)-1;11 }三、exgcd 解不定⽅程（使⽤不将a与b转为互质的⽅法） 对于 ax+by=c 的不定⽅程，设 r=gcd(a,b) 当 c%r!=0 时⽆整数解 当 c%r=0 时，将⽅程右边 *r/c 后转换为 ax+by=r 的形式 可以根据扩展欧⼏⾥得算法求得⼀组整数解 x0 , y0 ⽽这只是转换后的⽅程的解，原⽅程的⼀组解应再 *c/r 转变回去 （如 2x+4y=4 转换为 2x+4y=2 后应再将解得的 x , y 乘上2） 则原⽅程解为 x1=x0*c/r , y1=y0*c/r 通解 x=x1+b/r*t , y=y1-a/r*t ，其中 t 为整数 证明： 将 x , y 带⼊⽅程得 ax+ab/r*t+by-ab/r*t=c ax+by=c 此等式恒成⽴ 得证 这⾥ b/r 与 a/r 为最⼩的系数，所以求得的解是最多最全⾯的 证明： 为了推出证明中的 ax+by=c ，且想达到更⼩的系数，只能将 b/r 与 a/r 同除以⼀个数 s ⽽ b/r 与 a/r 互质，且 s 为整数，则 s=1 ，不影响通解 那么 b/r 与 a/r 就为最⼩的系数 得证四、exgcd 解线性同余⽅程 关于 x 的模⽅程 ax%b=c 的解 ⽅程转换为 ax+by=c 其中 y ⼀般为⾮正整数 则问题变为⽤ exgcd 解不定⽅程 解得 x1=x0*c/r 通解为 x=x1+b/r*t 设 s=b/r （已证明 b/r 为通解的最⼩间隔） 则 x 的最⼩正整数解为 (x1%s+s)%s 证明： 若 x1>0，则 (x1%s+s)%s=x1%s%s+s%s=x1%s=x1-ts (t∈N) 若 x1<0，因在 C++ ⾥ a%b=-(-a%b)<0 (a<0 , b>0) 如 -10%4=-2 则 (x1%s+s)%s=(-(-x1%s)+s)%s=(-(ts-x1)+s)%s=ts-x1 (t∈N) 即为 x1 通过加或减上若⼲个 s 后得到的最⼩正整数解 得证 亦可伪证 x1<0 的情况：设 x1=-5 , s=2 则 (x1%s+s)%s=(-5%2+2)%2=(-1+2)%2=3%2=1 即为 x1 加上 3 个 s 后的到的最⼩正整数解版权所有，转载请联系作者，违者必究。
拓展的欧几里得算法
拓展的欧几里得算法
哎呀呀，同学们，你们知道什么是拓展的欧几里得算法吗？反正一开始我是完全摸不着头脑呢！
就好像有一天，我走进了一个神秘的数学城堡。
这个城堡里到处都是数字和符号，让我眼花缭乱。
然后我就碰到了这个叫拓展的欧几里得算法的家伙。
我就问老师：“老师，这拓展的欧几里得算法到底是个啥呀？”老师笑着说：“孩子，它就像是一把神奇的钥匙，可以帮我们解决很多难题。
”我心里直犯嘀咕，能有多神奇？
后来老师给我们举了个例子。
比如说，我们要找两个数的最大公约数，普通方法可能很麻烦，但是用拓展的欧几里得算法，就能很快找到答案。
这就好比我们在一堆乱七八糟的玩具里找一个特定的玩具，没有好方法的时候，那简直是大海捞针，有了这个算法，就像是有了指南针，一下子就能找到目标啦！
有一次，我和同桌一起研究这个算法。
我着急地说：“哎呀，我怎么还是搞不懂啊！”同桌拍着我的肩膀说：“别着急，咱们一起慢慢琢磨。
”于是我们一步一步地算，算着算着，突然就好像有一道光照进了黑暗的角落，我们一下子就明白了！
你说，这拓展的欧几里得算法是不是很有趣？它就像是一个藏在数学世界里的宝藏，等待着我们去发现。
我觉得呀，虽然一开始它让我觉得头疼，但是当我真正理解它的时候，那种感觉太棒啦！就好像我征服了一座高山，看到了别人看不到的美丽风景。
所以说，数学的世界真是充满了惊喜和挑战，只要我们勇敢去探索，就能收获满满的快乐和知识！。