非原创：最大公约数问题



声明：算法思想来源于《编程之美》，在此仅作整理分享之用。

最大公约数问题在初等数学中是一个基础问题，当时我们用的方法是短除法。后来接触编程，常用的方法是转转相除法，此法由公元前大数学家欧几里得在《几何原本》给出解释。如果一个数能同时整除X和Y，则必能同时整除X%Y和Y，反之亦然，推出：X和Y的最大公约数与X%Y的最大公约数相等。

解法一：F（X,Y） = F(Y,X%Y)

<pre class="java" style="margin: 4px 0px; background-color: rgb(240, 240, 240);" name="code">public static int gcd1(int x,int y)
	{
		if(y==0)
			return x;
		else
			return gcd(y,x%y);
	}

优点：递推思路清晰，代码精简。

缺点：递归取模，开销昂贵。

为了避免使用取模，可以用减法替代取模运算。

解法二：F(X,Y) = F(X-Y,Y) (X>Y)

public static int gcd2(int x,int y)
	{
		if(x<y)
			return gcd2(y,x);
		else
		{
			if(y==0)
				return x;
			else
				return gcd2(x-y,y);
		}
	}

优点：避免了取模运算；

缺点：迭代次数很多，尤其是X和Y相差悬殊的时候。

看来为了减少迭代次数，除法运算不可避免。但为了减少开销，又要尽量避免除法运算，所以顺理成章的移位运算成了首选，也就是除以2^n.

对于X和Y奇偶情况，分为以下几种情况。

1，若X和Y均为偶数，F（X,Y） = 2*F(X>>1，Y>>1)

2，若X为偶数，Y为奇数，F（X,Y） = F(X>>2，Y)

3，若X为奇数，Y为偶数，F（X,Y） = F(X,Y>>2)

4，若X和Y均为奇数，F（X,Y） = F(X-Y,Y)

以上思路为解法一和解法二的综合版，但有一种特殊情况还是不能避免，如 F（5000,2） = 2*F（2500,1） = F（1250,1） =F（625,1） = F（625-1,1） = ......

所以在推出判断的时候加一个判断条件，判断是否有X或Y为1.

public static int gcd3(int x,int y)
	{
		if(x<y)
			return gcd3(y,x);
		if(y==0||x==1)
			return x;
		else
		{
			if(isEven(x)&&isEven(y))
				return (gcd3(x>>1,y>>1))<<1;
			else if(isEven(x))
				return gcd3(x>>1,y);
			else if(isEven(y))
				return gcd3(x,y>>1);
			else
				return gcd3(x-y,y);
				
		}
	}

总之，求最大公约数是一个很有趣的问题，尤其是当两个数比较大或者是要求的数比较多时，一个好的算法就显得尤为重要。