求两个正整数的最大公约数
大四毕业前夕,计算机学院的小灰又一次顶着炎炎烈日,
去某IT公司面试研发工程师岗位……
半小时后,公司会议室,面试开始……
小灰奋笔疾书,五分钟后……
小灰的思路十分简单。他使用暴力枚举的方法,试图寻找到一个合适的整数 i,看看这个整数能否被两个整型参数numberA和numberB同时整除。
这个整数 i 从2开始循环累加,一直累加到 numberA 和 numberB 中较小参数的一半为止。循环结束后,上一次寻找到的能够被两数整除的最大 i 值,就是两数的最大公约数。
事后,垂头丧气的小灰去请教同系的学霸大黄……
辗转相除法, 又名欧几里得算法(Euclidean algorithm),目的是求出两个正整数的最大公约数。它是已知最古老的算法, 其可追溯至公元前300年前。
这条算法基于一个定理:两个正整数a和b(a>b),它们的最大公约数等于a除以b的余数c和b之间的最大公约数。比如10和25,25除以10商2余5,那么10和25的最大公约数,等同于10和5的最大公约数。
有了这条定理,求出最大公约数就简单了。我们可以使用递归的方法来把问题逐步简化。
首先,我们先计算出a除以b的余数c,把问题转化成求出b和c的最大公约数;然后计算出b除以c的余数d,把问题转化成求出c和d的最大公约数;再然后计算出c除以d的余数e,把问题转化成求出d和e的最大公约数……
以此类推,逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算,直到两个数可以整除,或者其中一个数减小到1为止。
五分钟后,小灰改好了代码……
更相减损术, 出自于中国古代的《九章算术》,也是一种求最大公约数的算法。
他的原理更加简单:两个正整数a和b(a>b),它们的最大公约数等于a-b的差值c和较小数b的最大公约数。比如10和25,25减去10的差是15,那么10和25的最大公约数,等同于10和15的最大公约数。
由此,我们同样可以通过递归来简化问题。首先,我们先计算出a和b的差值c(假设a>b),把问题转化成求出b和c的最大公约数;然后计算出c和b的差值d(假设c>b),把问题转化成求出b和d的最大公约数;再然后计算出b和d的差值e(假设b>d),把问题转化成求出d和e的最大公约数……
以此类推,逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算,直到两个数可以相等为止,最大公约数就是最终相等的两个数。
五分钟后,小灰重写了代码……
众所周知,移位运算的性能非常快。对于给定的正整数a和b,不难得到如下的结论。其中gcb(a,b)的意思是a,b的最大公约数函数:
当a和b均为偶数,gcb(a,b) = 2*gcb(a/2, b/2) = 2*gcb(a>>1, b>>1)
当a为偶数,b为奇数,gcb(a,b) = gcb(a/2, b) = gcb(a>>1, b)
当a为奇数,b为偶数,gcb(a,b) = gcb(a, b/2) = gcb(a, b>>1)
当a和b均为奇数,利用更相减损术运算一次,gcb(a,b) = gcb(b, a-b), 此时a-b必然是偶数,又可以继续进行移位运算。
比如计算10和25的最大公约数的步骤如下:
- 整数10通过移位,可以转换成求5和25的最大公约数
- 利用更相减损法,计算出25-5=20,转换成求5和20的最大公约数
- 整数20通过移位,可以转换成求5和10的最大公约数
- 整数10通过移位,可以转换成求5和5的最大公约数
- 利用更相减损法,因为两数相等,所以最大公约数是5
在两数比较小的时候,暂时看不出计算次数的优势,当两数越大,计算次数的节省就越明显。
最后总结一下上述所有解法的时间复杂度:
1.暴力枚举法:时间复杂度是O(min(a, b)))
2.辗转相除法:时间复杂度不太好计算,可以近似为O(log(min(a, b))),但是取模运算性能较差。
3.更相减损术:避免了取模运算,但是算法性能不稳定,最坏时间复杂度为O(max(a, b)))
4.更相减损术与移位结合:不但避免了取模运算,而且算法性能稳定,时间复杂度为O(log(max(a, b)))
本文原本只写到辗转相除法就终告结束,后来网友们指出还有更优化的解法,看来自己还是才疏学浅,很感谢大家指出问题。另外,方法的参数默认必定是正整数,所以在代码中省去了合法性检查。
文中描述的更相减损术是简化了的方式。在九章算术原文中多了一步验证:如果两数都是偶数,计算差值之前会首先让两个数都折半,使得计算次数更少。这种方法做到了部分优化,但古人似乎没想到一奇一偶的情况也是可以优化的。
由于篇幅所限,本文省略了关于辗转相除法原和更相减损术的原理及证明。其实证明过程并不复杂,细心的同学们也可以自己尝试研究一下。谢谢大家的捧场!
求两个正整数的最大公约数相关推荐
- c语言求出两个最大素数,求两个正整数的最大公约数 思路:这是一个很基本的问题,最常见的就是两种方法,辗转相除法和辗转相减法。通式分别为 f(x, y) = f(y, x%y...
求两个正整数的最大公约数 思路:这是一个很基本的问题,最常见的就是两种方法,辗转相除法和辗转相减法.通式分别为 f(x, y) = f(y, x%y), f(x, y) = f(y, x - y) ( ...
- php怎么求最小公倍数,C++_详解C语言求两个数的最大公约数及最小公倍数的方法,求两个正整数的最大公约数nbs - phpStudy...
详解C语言求两个数的最大公约数及最小公倍数的方法 求两个正整数的最大公约数 思路:这是一个很基本的问题,最常见的就是两种方法,辗转相除法和辗转相减法.通式分别为 f(x, y) = f(y, x%y) ...
- C语言丨求两个正整数的最大公约数
两个正整数的最大公约数(Greatest Common Divisor, GCD)是能够整除这两个整数的最大整数.两个正整数的最大公约数的求法有多种解答,本文就三种方法做详细介绍:穷举法.欧几里得算法 ...
- 三种算法求两个正整数的最大公约数和最小公倍数;求三个数的最大公约数和最小公倍数
第二次作业 题目:求两个正整数的最大公约数和最小公倍数. 基本要求:1.程序风格良好(使用自定义注释模板),两种以上算法解决最大公约数问题,提供友好的输入输出. 提高要求:1.三种以上算法解决两个正整 ...
- 求两个正整数的最大公约数和最小公倍数
基本要求:1.程序风格良好(使用自定义注释模板),两种以上算法解决最大公约数问题,提供友好的输入输出. 提高要求:1.三种以上算法解决两个正整数最大公约数问题. 2.求3个正整数的最大公约数和最小公倍 ...
- 【算法】求两个正整数的最大公约数和最小公倍数
题目 输入两个正整数m和n,求其最大公约数和最小公倍数. public class Prog6{public static void main(String[] args){int m,n;try{m ...
- 求两个正整数的最大公约数 C语言
算法,采用辗转相除法求最大公约数. (1)输入两个整数数m和n (2) m%n余数为r,然后m=n;n=r; (3)当n=0,则m是最大公约数,算法结束,否则转至执行(2)重复上述过程,直至n=0为止 ...
- 欧几里得算法求两个正整数的最大公约数
此题可以使用暴力破解法解决,但是效率不高. 欧几里德算法又称辗转相除法,其计算原理依赖于下面的定理: gcd(a,b) = gcd(b,a mod b) (a>b 且a mod b 不为0) ...
- Python 求两个正整数的最大公约数和最小公倍数
最大公约数 a . b 的顺序并不影响最终结果,如果 a .b 是顺序的话,会多一次递归(相当于第一次递归会用来调换大小顺序). # 最大公约数 greatest common divisor # 辗 ...
- c语言两个正整数的最小公倍数,C语言求两个正整数的最小公倍数
这里将介绍求两个正整数的最小公倍数(Least Common Multiple,LCM)的方法.提供两种主要思路,一种是直接根据最小公倍数的定义设计算法,一种是由最大公约数计算得出.下面来介绍这两种方 ...
最新文章
- 【ACM】杭电OJ 2048 2049
- 我的Dll(动态链接库)学习笔记
- Expo大作战(三十)--expo sdk api之Permissions(权限管理模块),Pedometer(计步器api)
- linux filesystem_如何使用cgdb + qemu调试linux内核模块
- shell 下使用 echo 打印彩色字体及彩色背景
- c# 第四课 interfaces
- 视图编辑php,PhpStorm视图模式
- linux 的间隔定时器函数setitimer
- 01_Android应用开发环境_01_android发展史及系统架构
- 图像识别-百度AI开放平台
- python中︿是什么意思_Python learning notes-0003-注释、变量、简单类型、运算符,学习,笔记...
- DELPHI10.2的LINUX数据库开发环境配置
- MaxCompute命令行工具——odpscmd的操作使用
- 解决vm安装过程显示无人参与应答文件包含的产品密钥无效
- Snake模型综合介绍
- unity 导入STL格式模型(STL 文本ASCII码格式文件)
- 华为路由器忘记密码怎么恢复
- 深圳大学计算机系在哪个校区,2021年深圳大学有几个校区,大一新生在哪个校区...
- C# 中的委托和事件[转自http://www.cnblogs.com/jimmyzhang/archive/2007/09/23/903360.html]
- 解决tensorflow.python.framework.errors_impl.InvalidArgumentError: indices[0,32] = -1 is not in [0, 50)
热门文章
- 学习Linux内核内存管理要看的布局图
- 计算机应用基础演示文稿内容,计算机应用基础(PowerPoint 2003电子演示文稿系统)...
- 副高相当于副处吗_事单副高职称相当于副处待遇?别逗了,连公务员副科实职都不如...
- 【微信小程序】小程序是什么?怎样快速上手开发一个小程序?
- 佳能无线打印机服务器,佳能LBP6018W打印机WIFI无线打印心得分享
- ValueError: The name None occurs multiple times, use a level number
- 模式识别的常用英文总结
- 《图解TCP/IP》读书笔记
- 帮蔡徐坤刷出一亿转发,APP开发者一审获刑
- 流落在帝都的那些80后北漂IT人,你们的未来在哪儿