写一个程序,求两个正整数的最大公约数。如果两个正整数都很大,有什么简单算法?
例如,给定两个数 1100100210001,,10200021,求出其最大公约数。
分析
辗转相除法的原理,假设用f(x,y)表示x,y的最大公约数,取k=x/y, b=x%y,则 x=ky + b,如果一个数能够同时整除x和y,则必能同时整除b和y;而能够同时整除b和y的数也必能同时整除x和y,即x和y的公约数与 b和y的公约数是相同的,其最大公约数也是相同的,则有f(x,y)=f(y,x%y) (x>=y>0),如此便可把原问题转化为求两个更小数的最大公约数,直到其中一个数为0,剩下另一个数就是两者最大的公约数。
示例如下:
f(42,30)=f(30,12)=f(12,6)=f(6,0)=6
解法一:最简单的方法,直接用代码来实现辗转相除法。利用递归就能够解决这个问题。
int gcd(int x,int y) {
return (!y)?x:gcd(y, x%y);
}
解法二:
解法一中,用到了取模运算。但对于大整数而言,取模运算(其中用到除法)是非常昂贵开销,将成为整个算法的瓶颈。是否有办法不用取模运算?
采用类似前面辗转相除法的分析,如果一个数能够同时整除x和y,则必能同时整除x-y和y;而能够同时整除x-y 和 y的数也必能同时整除x和y,即x和y的公约数与x-y和y的公约数是相同的,其最大公约数也是相同的,即f(x,y)=f(x-y,y),那么就可以不再需要进行大整数的取模运算,而转换成简单得多的大整数的减法。
实际操作中,如果x<y,可以先交换(x,y),从而避免求一个正数和一个负数的最大公约数情况的出现。一直迭代下去,直到妻子一个数位0。
示例如下: f(42,30)=f(30,12)=f(12,18)=f(18,12)=f(12,6)=f(6,6)=f(6,0)=6
BigInt gcd(BigInt x, BigInt y) {
if(x<y)
return gcd(y,x);
if(y==0)
return x;
else
return gcd(x-y,y);
}
此算法,免去了大整数除法的繁琐,但是同样有不足之处。最大的瓶颈就是迭代的次数比之前的算法多了不少,如果遇到 (10 000 000 000 000,1)这类清空,就会很烦。
解法三:
结合上述的解法一和解法二从而使其成为一个最佳的算法。
从公约数的特点入手。对于y和x来说,如果y=k*y1,x=k*x1.那么有f(y,x)=k*f(y1,x1)。另外,如果x=p*x1,假设p是素数,并且y%p!=0(即y不能被p整除),那么f(x,y)=f(p*x1,y)=f(x1,y)。
注意到以上两点之后,可以利用这两点对算法进行改进。
2是一个素数,同时对于二进制表示的大整数而言,可以很容易地将除以2和乘以2的运算转换成移位运算,从而避免大整数除法,由此就可以利用2这个数字来进行分析。
取p=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)
若y为奇数,y为偶数,f(x,y)=f(x,y/2)=f(x,y>>1)
若x,y均为奇数,f(x,y)=f(y,x-y),
那么在f(x,y)=f(y,x-y)之后,(x-y)是一个偶数,下一步一定会有除以2的操作。
因此,最坏情况下的时间复杂度是O(log(max(x,y)))。
BigInt gcd(BigInt x, BigInt y) {
if(x<y)
return gcd(y,x);
if(y==0)
return x;
else {
if(IsEven(x)) //检查x是否为偶数,
{
if(IsEven(y))
return (gcd(x>>1,y>>1) <<1);
else
return gcd(x>>1,y);
}
else
{
if(IsEven(y))
return gcd(x,y>>1);
else
return gcd(y,x-y);
}
}
}
巧妙利用了移位运算和减法运算,避开了大整数除法,提高了算法的效率。