欧几里得算法和扩展的欧几里得算法C++递归实现

欧几里得算法和扩展的欧几里得算法C++递归实现

关于欧几里得算法的流程不再赘述,不清楚的可以搜得到。本篇主要通过C++代码利用递归的思想实现,参考书籍是《密码编码与信息安全:C++实践》。

1、欧几里得算法实现

欧几里得算法比较简单,主要用于求两个数(多项式)的最大公因数(式),直接上代码。

#include <iostream>
using namespace std;
int Euclidean(int a, int b){
	if(b==0){
		return a;
	}else{
		return Euclidean(b, a%b);
	}
}

2、扩展的欧几里得算法实现

扩展的欧几里得算法主要用于求模逆运算。我第一次实现扩展的欧几里得算法是通过辗转相除,然后再回溯求出了 a a a b b b的系数。感觉可以像普通欧几里得算法一样可以递归编程,但是总是没想出来。后来借助参考书实现了。主要思想是要写出此时的递推关系式。

2.1递推关系式的说明

扩展欧几里得算法本质上是要求得
a × s + b × t = g c d a\times s+b\times t=gcd a×s+b×t=gcd
这个式子。按照递归的思想,我们应该这么考虑:要利用递归得到 a a a b b b的式子,结合辗转相除的原理,我们肯定是先得到 b b b a ( m o d b ) a\pmod b a(modb)的关系式。假设已经有了
s ′ × b + t ′ × ( a ( m o d b ) ) = g c d s'\times b+t'\times (a\pmod b)=gcd s×b+t×(a(modb))=gcd

如何得到我们需要的式子?

做一个简单的变形就出来了:我们知道 a ( m o d b ) = a − ⌊ a / b ⌋ ∗ b a\pmod b=a-\lfloor a/b\rfloor *b a(modb)=aa/bb,代进上式,有
s ′ × b + t ′ × ( a ( m o d b ) ) = s ′ × b + t ′ × ( a − ⌊ a / b ⌋ ∗ b ) = t ′ × a + ( s ′ − t ′ × ⌊ a / b ⌋ ) × b = g c d s'\times b+t'\times (a\pmod b)=s'\times b+t'\times (a-\lfloor a/b\rfloor *b)\\ =t'\times a+(s'-t'\times \lfloor a/b\rfloor)\times b =gcd s×b+t×(a(modb))=s×b+t×(aa/bb)=t×a+(st×a/b)×b=gcd
所以递归就是根据这个式子编出的。具体代码如下:

#include <iostream>
using namespace std;
struct ExtEuc
{
	int gcd;
	int s;
	int t;
};
//the most important is the "recursion formula",especially in Extends_Eclidean
ExtEuc Extends_Eclidean(int a, int b){
	ExtEuc obj, obj1;
	if(b == 0){
		obj1.gcd = a;
		obj1.s = 1;
		obj1.t = 0;
		return obj1;
	}
	obj1 = Extends_Eclidean(b, a%b);
	//attention!!
	obj.gcd = obj1.gcd;
	obj.s = obj1.t;
	obj.t = obj1.s - (a/b)*obj1.t;
	return obj;
}

3、整体代码如下

#include <iostream>
using namespace std;
struct ExtEuc
{
	int gcd;
	int s;
	int t;
};

//the most important is the "recursion formula",espeacially in Extends_Eclidean
int Euclidean(int a, int b){
	if(b==0){
		return a;
	}else{
		return Euclidean(b, a-(a/b)*b);
	}
}

ExtEuc Extends_Eclidean(int a, int b){
	ExtEuc obj, obj1;
	if(b == 0){
		obj1.gcd = a;
		obj1.s = 1;
		obj1.t = 0;
		return obj1;
	}

	obj1 = Extends_Eclidean(b, a%b);
	//attention!!
	obj.gcd = obj1.gcd;
	obj.s = obj1.t;
	obj.t = obj1.s - (a/b)*obj1.t;
	return obj;
	
}

int main(int argc, char const *argv[])
{
	/* code */
	ExtEuc f = Extends_Eclidean(2,3);
	cout<<f.s<<"*2+"<<f.t<<"*3="<<f.gcd<<"  "<<endl;
	return 0;
}

[1] 王静文, 吴晓艺. 密码编码与信息安全:C++实践[M]. 清华大学出版社, 2015.

已标记关键词 清除标记
相关推荐
【为什么还需要学习C++?】 你是否接触很多语言,但从来没有了解过编程语言的本质? 你是否想成为一名资深开发人员,想开发别人做不了的高性能程序? 你是否经常想要窥探大型企业级开发工程的思路,但苦于没有基础只能望洋兴叹?   那么C++就是你个人能力提升,职业之路进阶的不二之选。 【课程特色】 1.课程共19大章节,239课时内容,涵盖数据结构、函数、类、指针、标准库全部知识体系。 2.带你从知识与思想的层面从0构建C++知识框架,分析大型项目实践思路,为你打下坚实的基础。 3.李宁老师结合4大国外顶级C++著作的精华为大家推出的《征服C++11》课程。 【学完后我将达到什么水平?】 1.对C++的各个知识能够熟练配置、开发、部署; 2.吊打一切关于C++的笔试面试题; 3.面向物联网的“嵌入式”和面向大型化的“分布式”开发,掌握职业钥匙,把握行业先机。 【面向人群】 1.希望一站式快速入门的C++初学者; 2.希望快速学习 C++、掌握编程要义、修炼内功的开发者; 3.有志于挑战更高级的开发项目,成为资深开发的工程师。 【课程设计】 本课程包含3大模块 基础篇 本篇主要讲解c++的基础概念,包含数据类型、运算符等基本语法,数组、指针、字符串等基本词法,循环、函数、类等基本句法等。 进阶篇 本篇主要讲解编程中常用的一些技能,包含类的高级技术、类的继承、编译链接和命名空间等。 提升篇: 本篇可以帮助学员更加高效的进行c++开发,其中包含类型转换、文件操作、异常处理、代码重用等内容。
©️2020 CSDN 皮肤主题: 1024 设计师:白松林 返回首页