最大公约数和快速gcd

最大公约数有两种基本的求法:(1)辗转相除法(2)更相减损法，首先来证明一下:

$(1)证明gcd(a,b) = gcd(b,a%b)$

$假设a = k_1m,b = k_2m,a = sb+r$

$可以知道r = a-sb = (k_1-sk_2)m,所以m是(b,r)的公约数$

$如果存在n>m,使得n也是(b,r)的公约数,b=k_3n,r=k_4n那么a = sk_3n+k_4n$

$那么n也是(a,b)的公约数,这与m是(a,b)最大公约数矛盾。$

$(2)证明gcd(a,b) = gcd(a-b,b)$

$假设a = k_1m,b = k_2m,a - b= (k_1-k_2)m$

$可以知m是(a-b,b)的公约数,$

$如果存在n>m,使得n也是(a-b,b)的公约数,那么a = (k_3+k_4)n,则n也是(a,b)的最大公约数,这同样与m=gcd(a,b)矛盾$

$(3)那么相对于辗转相除法,显然更相减损法要慢的多$

然而当我们进行高精度计算的时候(比如几千位的运算)，取余操作时非常慢的。但是我们可以在更相减损的基础上加以改进，加快算法的实现,避免求余:起始所谓的快速gcd很简单,只要把a,b的奇数与偶数的性质分类讨论一下,就很容易优化。比如a,b同时为偶数，那么可以同时让a,b右移一位，同时记得将最后的结果乘2。如果一奇一偶,很容易就可以想到:2一定不是约数,所以可以将偶数右移。如果两个奇数,那么就可以相减，这样就会转换到一奇数一偶数的状态。c++的实现如下:

long long int gcd1(const long long int _a,const long long int _b){
    auto a = _a, b = _b;
    if (a < b)
        return gcd1(b, a);
    long long int t;
    while (b>0){
        t = a;
        a = b;
        b = t%b;
    }
    return a;
}
long long int gcd3(const long long int _a,const long long int _b){
    auto a = _a, b = _b;
    if (a < b){ gcd3(b, a); }

    int factor = 1;
    int cond = 0;
    int s;
    while (a != b){
        if (a < b) {
            long long int t = a;
            a = b;
            b = t;
        }
        cond = ((a & 1) << 1) + (b & 1);
        if (cond == 0) { factor <<= 1; a >>= 1; b >>= 1; }//even, even
        else if (cond == 1){ a >>= 1; }//a = even, b = odd
        else if (cond == 2){ b >>= 1; }//a = odd, b = even
        else {
            long long int t = a;
            a = b;
            b = t - b;       //#a = odd, b = odd
        }
    }
    return a*factor;
}

python实现

def gcd1(a,b):
    if a<b:return gcd1(b,a)

    while b:
        a,b = b,a%b
    return a
def gcd2(a,b):
    if a<b:return gcd2(b,a)

    while a!=b:
        a,b = b,a-b
        if a< b:a,b = b,a
    return a
def gcd3(a,b):
    if a<b:gcd3(b,a)
    factor = 0

    while a!=b:
        if a<b: a,b = b,a
        cond = ((a&1)<<1)+(b&1)
        if cond==0:a=1;b>>=1;factor += 1#even,even
        elif cond==1:a,b = a>>1,b #a = even,b = odd
        elif cond==2:a,b = a,b>>1 #a = odd,b = even
        else:a,b = b,a-b          #a = odd,b = odd
    return a*(1<<factor)