2.7 编程之美--最大公约数的3种解法[efficient method to solve gcd problem]

【本文链接】

http://www.cnblogs.com/hellogiser/p/efficient-method-to-solve-gcd-problem.html

【题目】

　　求两个正整数的最大公约数Greatest Common Divisor (GCD)。如果两个正整数都很大，有什么简单的算法吗？例如，给定两个数1 100 100 210 001, 120 200 021,求出其最大公约数。

【解法】

---

 【1. 辗转相除法】

辗转相除法：f(x,y) = f(y , x % y)(x>y)

f(42,30) = f(30,12) = f(12,6) = f(6,0) = 6

【代码】

 C++ Code 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

/*     version: 1.0     author: hellogiser     blog: http://www.cnblogs.com/hellogiser     date: 2014/7/8 */ int gcd(int x, int y) {     if(x < y)         return gcd(y, x);     if(y == 0)         return x;     else         return gcd(y, x % y); }

此方法中用到了取模运算，对于大整数而言，取模运算（其中用到了除法）开销是非常昂贵的，将成为整个算法的瓶颈。

---

【2. 辗转相减法】

辗转相减法：f(x,y) = f(y, x-y) (x>y)

f(42,30) = f(30,12) = f(12,18) = f(18,12) = f(12,6)=f(6,6)=f(6,0)=6

【代码】

 C++ Code 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

/*     version: 1.0     author: hellogiser     blog: http://www.cnblogs.com/hellogiser     date: 2014/7/8 */ int gcd(int x, int y) {     if(x < y)         return gcd(y, x);     if(y == 0)         return x;     else         return gcd(y, x - y); }

这个算法免去了大整数除法的繁琐，但同样也有不足之处。最大的瓶颈是迭代的次数过多，如果出现（1 000 000 000,1）这类情况，那就相当让人郁闷了。

---

【3. 奇偶法】

奇偶法：

此种方法是将解法1)和解法2)结合起来，降低计算复杂度的同时也降低迭代次数。

1:若 x, y均为偶数，f (x,y) = 2 * f(x / 2, y / 2) = 2 * f(x >> 1, y >> 1)

2:若x为偶，而y为奇，f (x , y )  = f (x / 2, y) = f ( x >> 1, y)

3:若x为奇，y为偶，f ( x, y)  = f (x , y / 2) = f(x , y >> 1)

4:若x,y均为奇，f ( x, y ) = f (y , x - y)

在f(x, y) = f(y, x - y)之后，(x - y)是一个偶数，下一步一定会有除以2的操作。

因此最坏情况下时间复杂度为O(log2 (max(x,y)))。

f (42 , 30 ) = 2 * f (21,15)      

                   = 2 * f (15,6)    

                   = 2 * f (15,3)      

                   = 2 * f (3,12)  =2 * f (12,3)  

                   = 2 * f (6,3)         

                   = 2 * f (3,3)       

                   = 2 * f (3,0)         

                   = 2 * 3

                   = 6

【代码】

 C++ Code 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36

/*     version: 1.0     author: hellogiser     blog: http://www.cnblogs.com/hellogiser     date: 2014/7/8 */ bool IsEven(int x) {     return (x & 0x1) == 0; } int gcd(int x, int y) {     if(x < y)         return gcd(y, x);     if(y == 0)         return x;     else     {         if(IsEven(x))         {             if(IsEven(y)) //case 1,x,y均为偶数                 return 2 * gcd(x >> 1, y >> 1);             else  //case 2，x为偶，y为奇                 return gcd(x >> 1, y);         }         else         {             if(IsEven(y)) //case 3，x为奇，y为偶                 return gcd(x, y >> 1);             else  //case 4，x,y均为奇                 return gcd(y, x - y);         }     } }

【参考】

http://blog.csdn.net/ajioy/article/details/7478008

http://blog.csdn.net/rein07/article/details/6739688

【本文链接】

http://www.cnblogs.com/hellogiser/p/efficient-method-to-solve-gcd-problem.html