编程之美--最大公约数问题

思路：

最大公约数问题也是一个非常典型的递归算法的应用。每次递归使得原来求两个大数之间的公约数转变成求两个稍微小点的数之间的公约数，要求转换的过程要保证不会改变公约数的值。这就要看其中转换的原理了。

原理从《几何原本》中得出--辗转相除。假设f(x, y) 表示x,y的最大公约数是g，而k = x/y,b= x%y，则g必能整出b。因为x = ky + b，b = x - ky，b/g = (x-ky)/g一定为整数，所以必有g整除b。

如下所示：

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

代码如下：

int gcd(int x , int y)
{
return (y == 0 )?x :gcd(y , x % y) ;

}

书中还引伸出辗转相减法，原理跟上面所述差不多。

代码如下：

int gcd(int x , int y)
{
if(x < y)
return  gcd(y , x) ;
else if(y == 0)
return x ;
else
return gcd(x - y , y) ;

}

最后书中提到一种相除和相减相结合的方法，保证了不做过多冗余的步骤。这就是把2当作每次转换的数量级。

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

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

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

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

此解法结合了上述两种方法：第一种方法递归次数相当少但每次取模挺耗时（到底耗时情况如何，我没有具体了解过，但肯定比加减法耗时些），而第二种递归可能会相当多，但每次运算都是减法运算， 很快。综合两者，用2为基数，可以保证递归次数不那么多，同时运算也快。

代码如下：

int gcd(int x , int y)
{
if(x < y)
return gcd(y , x) ;
if(y == 0)
return x ;
if(isEven(x)) //x为奇
{
if(isEven(y)) //y为奇
return gcd(x - y, y) ;
else         //y为偶
return gcd(x , y >>1) ;
}
else         //x为偶
{
if(isEven(y)) //y为奇数
return gcd(x >> 1, y) ;
else         //y为偶
return 2 * gcd(x >> 1, y >> 1) ;

}

}