gcd与快速gcd

GCD与快速GCD

求最大公约数有两种方法：更相减损术和辗转相除法。

转转相除法

欧几里得算法又称辗转相除法，是指用于计算两个非负整数a，b的最大公约数。应用领域有数学和计算机两个方面。计算公式gcd(a,b) = gcd(b,a mod b)。

通过递归实现

int gcd(int a, int b)
{ 
    return b == 0 ? a : gcd( b, a % b);
} 

通过循环实现

int gcd(int x, int y)
{ 
    if (!x || !y) return x > y ? x : y; 
    for (int t; t = x % y; x = y, y = t); 
    return y; 
} 

通过逻辑运算实现

我们观察gcd的递归公式可以发现他的操作步骤：

1. a = a%b
2. 交换a b
3. 重复前两步，直到b等于0，此时a就是最大公约数
   我们可以通过逻辑运算实现交换a，b

b = a ^ b;
a = a ^ b;
b = a ^ b;

连续异或两个相同的数就相当于没有进行任何操作。
b = a ^ b;
a = a ^ b = a ^ (a ^ b) = b;
b = (b) ^ (a ^ b) = a;
所以我们通过一个while循环就能实现a % = b 然后交换a和b

int gcd(int a,int b)
{
	while(b ^= a ^= b ^= a %= b);
	return a;
}

更相减损术

可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。

这就是更相减损术的原文，我们在这个基础上优化后叫快速gcd。

正常情况下

第一步：任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简；若不是则执行第二步。
第二步：以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到所得的减数和差相等为止。
则第一步中约掉的若干个2的积与第二步中等数的乘积就是所求的最大公约数。
其中所说的“等数”，就是公约数。求“等数”的办法是“更相减损”法。

如何优化

按照这个思路我们来考虑两个数a和b
1.如果两个数都是偶数，按照文中所说两个数都除以二然后结果再乘以二就行了
2.如果两个数是一奇一偶，那么很明显2不是他们的公约数。偶数除以二并不会对最大公约数造成影响。
3.如果两个都是奇数，那么两数相减一定是偶数。这样按照更相减损术就又能得到一奇一偶，然后就可以继续进行除二操作。

代码

ll kgcd(ll a, ll b)
{ 
    if (a == 0) return b; 
    if (b == 0) return a; 
    if (!(a & 1) && !(b & 1)) 
        return kgcd(a>>1, b>>1)<<1;
    else if (!(b & 1))
        return kgcd(a, b>>1); 
    else if (!(a & 1)) 
        return kgcd(a>>1, b); 
    else 
        return kgcd(abs(a - b), min(a, b)); 
}