最大公约数的最佳解决算法

 

首先从分析公约数的特点入手：

对于y和x来说，如果y=k * y₁，x=k * x₁。那么有f（y, x）= k * f（y₁, x₁）。

另外，如果x = p * x₁，假设p是素数，并且y % p ! = 0（即y不能被p整除），那么f（x, y）= f（p * x₁, y）= f（x₁, 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）

若x为奇数，y为偶数，f（x, y）= f（x, y/2）= f（x, y>>1）

若x, y均为奇数，f（x, y）= f（x, x - y），

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

因此，最坏情况下的时间复杂度是O（log₂（max（x, y））。

考虑如下的情况：

       f（42, 30）= f（101010₂, 11110₂）

                       = 2 * f（10101₂, 1111₂）

                       = 2 * f（1111₂, 110₂）

                       = 2 * f（1111₂, 11₂）

                       = 2 * f（1100₂, 11₂）

                       = 2 * f（11₂, 11₂）

                       = 2 * f（0₂, 11₂）

                       = 2 * 11₂

                       = 6

根据上面的规律，具体代码实现如下：

代码清单2-16

BigInt gcd(BigInt x, BigInt y)

{

    if(x < y)

        return gcd(y, x);

    if(y == 0)

        return x;

    else

    {

        if(IsEven(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);

        }

    }

}

BigInt见解法二中的解释，IsEven（BigInt x）函数检查x是否为偶数，如果x为偶数，则返回true，否则返回false。