辗转相除法
又名欧几里得算法(Euclidean algorithm)
- 该算法的目的是求出两个 正整数的最大公约数。它是已知最古老的算法, 其产生时间可追溯至公元前300年前。
- 这条算法基于一个定理:两个正整数a和b(a>b),它们的最大公约数等于a除以b 的余数c和b之间的最大公约数。
- 例如10和25,25除以10商2余5,那么10和25的最大公约数,等同于10和5的最大公约 数
更相减损术
- 更相减损术,出自中国古代的《九章算术》,也是一种求最大公约数的算法。
- 古希腊 人很聪明,可是我们炎黄子孙也不差。
- 它的原理更加简单:两个正整数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,不难得到如下的结论。
-
(从下文开始,获得最大公约数的方法getGreatestCommonDivisor被简写为gcd。)
-
当a和b均为偶数时,gcd(a,b) = 2×gcd(a/2, b/2) = 2×gcd(a>>1,b>>1)。
-
当a为偶数,b为奇数时,gcd(a,b) = gcd(a/2,b) = gcd(a>>1,b)。
-
当a为奇数,b为偶数时,gcd(a,b) = gcd(a,b/2) = gcd(a,b>>1)。
-
当a和b均为奇数时,先利用更相减损术运算一次,gcd(a,b) = gcd(b,a-b),此时a-b必然 是偶数,然后又可以继续