漫画算法:辗转相除法是什么鬼?

大四毕业前夕,计算机学院的小灰又一次顶着炎炎烈日,

去某IT公司面试研发工程师岗位……

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%b5%b6%e8%b7%af

半小时后,公司会议室,面试开始……

%e5%b0%8f%e4%bb%93%e9%bc%a0%e9%9d%a2%e8%af%95%e5%ae%981

%e5%b0%8f%e4%bb%93%e9%bc%a0%e9%9d%a2%e8%af%95%e5%ae%987

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%9b%9e%e7%ad%944

%e5%b0%8f%e4%bb%93%e9%bc%a0%e9%9d%a2%e8%af%95%e5%ae%985

%e9%9d%a2%e8%af%95%e5%ae%98%e5%82%b2%e6%85%a22

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%8b%a6%e6%80%9d

小灰奋笔疾书,五分钟后……

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%9b%9e%e7%ad%945

%e4%bb%a3%e7%a0%811

小灰的思路十分简单。他使用暴力枚举的方法,试图寻找到一个合适的整数 i,看看这个整数能否被两个整型参数numberA和numberB同时整除。

这个整数 i 从2开始循环累加,一直累加到 numberA 和 numberB 中较小参数的一半为止。循环结束后,上一次寻找到的能够被两数整除的最大 i 值,就是两数的最大公约数。

%e5%b0%8f%e4%bb%93%e9%bc%a0%e9%9d%a2%e8%af%95%e5%ae%988

%e9%9d%a2%e8%af%95%e5%ae%98%e5%82%b2%e6%85%a2

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%8b%a6%e6%80%9d

%e5%b0%8f%e4%bb%93%e9%bc%a0%e6%97%a0%e5%a5%88

%e5%b0%8f%e4%bb%93%e9%bc%a0%e9%9d%a2%e8%af%95%e5%ae%980

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%a4%b1%e6%9c%9b

事后,垂头丧气的小灰去请教同系的学霸大黄……

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a33

%e5%b0%8f%e4%bb%93%e9%bc%a0%e6%97%a0%e5%a5%882

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a34

辗转相除法, 又名欧几里得算法(Euclidean algorithm),目的是求出两个正整数的最大公约数。它是已知最古老的算法, 其可追溯至公元前300年前。

这条算法基于一个定理:两个正整数a和b(a>b),它们的最大公约数等于a除以b的余数c和b之间的最大公约数。比如10和25,25除以10商2余5,那么10和25的最大公约数,等同于10和5的最大公约数。

%e6%ac%a7%e5%87%a0%e9%87%8c%e5%be%97

有了这条定理,求出最大公约数就简单了。我们可以使用递归的方法来把问题逐步简化。

首先,我们先计算出a除以b的余数c,把问题转化成求出b和c的最大公约数;然后计算出b除以c的余数d,把问题转化成求出c和d的最大公约数;再然后计算出c除以d的余数e,把问题转化成求出d和e的最大公约数……

以此类推,逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算,直到两个数可以整除,或者其中一个数减小到1为止。

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a35b

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%9b%9e%e7%ad%946

 

五分钟后,小灰改好了代码……

%e4%bb%a3%e7%a0%812

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a37

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%8b%a6%e6%80%9d3

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a38

更相减损术, 出自于中国古代的《九章算术》,也是一种求最大公约数的算法。

他的原理更加简单:两个正整数a和b(a>b),它们的最大公约数等于a-b的差值c和较小数b的最大公约数比如10和25,25减去10的差是15,那么10和25的最大公约数,等同于10和15的最大公约数。

%e4%b9%9d%e7%ab%a0%e7%ae%97%e6%9c%af

由此,我们同样可以通过递归来简化问题。首先,我们先计算出a和b的差值c(假设a>b),把问题转化成求出b和c的最大公约数;然后计算出c和b的差值d(假设c>b),把问题转化成求出b和d的最大公约数;再然后计算出b和d的差值e(假设b>d),把问题转化成求出d和e的最大公约数……

以此类推,逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算,直到两个数可以相等为止,最大公约数就是最终相等的两个数。

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a39

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%9b%9e%e7%ad%946

五分钟后,小灰重写了代码……

%e4%bb%a3%e7%a0%813

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a310

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%8b%a6%e6%80%9d4

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a311

%e5%b0%8f%e4%bb%93%e9%bc%a0%e6%b2%89%e6%80%9d4

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a312

众所周知,移位运算的性能非常快。对于给定的正整数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的最大公约数的步骤如下:

  1. 整数10通过移位,可以转换成求5和25的最大公约数
  2. 利用更相减损法,计算出25-5=20,转换成求5和20的最大公约数
  3. 整数20通过移位,可以转换成求5和10的最大公约数
  4. 整数10通过移位,可以转换成求5和5的最大公约数
  5. 利用更相减损法,因为两数相等,所以最大公约数是5

在两数比较小的时候,暂时看不出计算次数的优势,当两数越大,计算次数的节省就越明显。

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a313

%e4%bb%a3%e7%a0%814

最后总结一下上述所有解法的时间复杂度:

1.暴力枚举法:时间复杂度是O(min(a, b)))

2.辗转相除法:时间复杂度不太好计算,可以近似为O(log(min(a, b))),但是取模运算性能较差。

3.更相减损术:避免了取模运算,但是算法性能不稳定,最坏时间复杂度为O(max(a, b)))

4.更相减损术与移位结合:不但避免了取模运算,而且算法性能稳定,时间复杂度为O(log(max(a, b)))

%e5%b0%8f%e4%bb%93%e9%bc%a0%e5%9b%9e%e7%ad%947

%e5%b0%8f%e4%bb%93%e9%bc%a0%e8%ae%b2%e8%a7%a314

本文原本只写到辗转相除法就终告结束,后来网友们指出还有更优化的解法,看来自己还是才疏学浅,很感谢大家指出问题。另外,方法的参数默认必定是正整数,所以在代码中省去了合法性检查。

文中描述的更相减损术是简化了的方式。在九章算术原文中多了一步验证:如果两数都是偶数,计算差值之前会首先让两个数都折半,使得计算次数更少。这种方法做到了部分优化,但古人似乎没想到一奇一偶的情况也是可以优化的。

由于篇幅所限,本文省略了关于辗转相除法原和更相减损术的原理及证明。其实证明过程并不复杂,细心的同学们也可以自己尝试研究一下。谢谢大家的捧场!

 

本人微信号:bjweimengshu

欢迎朋友们一起交流讨论,加好友请注明伯乐在线 :)

打赏支持我写出更多好文章,谢谢!

打赏作者

打赏支持我写出更多好文章,谢谢!

5 26 收藏 32 评论

关于作者:玻璃猫

互联网公司的码农一枚,喜欢算法和面向对象设计。个人微信号:13522239721 个人订阅号:dreamsee321欢迎一起交流讨论! 个人主页 · 我的文章 · 68 ·  

相关文章

可能感兴趣的话题



直接登录
最新评论
  • 王念一 初三学生 09/30

    位移啊!

    位移!!!

    • 玻璃猫 程序员 09/30

      移位怎么解决,求指教

      • meta python,前端 09/30

        编程之美2.7,求大数公约数,避开大整数相除的方法。

      • 王念一 初三学生 10/01

        你真的是程序员吗…………

  • 玻璃猫 程序员 10/04

    文章进行了扩充,增加了更相减损法和移位方法。感谢meta和王念一同学及时指出!

  • 小吒 产品设计 10/08

    不同的思路和方法,带来最优的解法

  • topillar   10/08

    博主,最后的

    return gcd(numberA, numberA - numberB)

    应该是

    return gcd(numberB, numberA - numberB)

     

    除此,真正的好文章。

  • 清风   10/11

    小数如果是素数,最大公约数一定是1,如果先判断小数是否是素数,是不是会更快一些?

  • Hobit sorry,coder 10/17

    太帅了,主要是这个漫画表情太帅了

  • Louislip   10/19

    在最后位移和更相减法的算法中,只有a=b=1是才会结束递归,若a=1000000,b=2,在递归第一次后其实已经知道结果,由于a非常大,所以会继续递归;加上判断 if(a>>1==0 || b>>1==0) return 1 结束递归?

跳到底部
返回顶部