最大公约数与最小公倍数求解算法

最大公约数与最小公倍数求解算法

1.最大公约数

1）辗转相除法

使用辗转相除法求解，这条算法基于一个定理：

两个正整数a和b（a>b），它们的最大公约数等于a除以b的余数c和b之间的最大公约数。
比如10和25，25除以10商2余5,那么10和25的最大公约数，等同于10和5的最大公约数。

使用递归方法来实现此算法：

首先，我们先计算出a除以b的余数c，把问题转化成求出b和c的最大公约数；然后计算出b除以c的余数d，把问题转化成求出c和d的最大公约数；再然后计算出c除以d的余数e，把问题转化成求出d和e的最大公约数…
以此类推，逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算，直到两个数可以整除，或者其中一个数减小到1为止。

int grd(int a ,int b){
	//这里规定传入参数a>=b，不再对a,b大小进行比较; 
	//规定b必须大于0，否则会出现对零取余错误 
	if(a%b==0) return b;
	else
	{
		return grd(b,a%b);//这里最多会递归到b=1。 
	}
	//返回值即为最大公约数 
}

2）更相减损术

当俩个整数型较大时，做a%b取模运算的性能比较低，这时候我们可以采用更相减损术。
原理：

两个正整数a和b（a>b），它们的最大公约数等于a-b的差值c和较小数b的最大公约数。
比如10和25，25减去10的差是15,那么10和25的最大公约数，等同于10和15的最大公约数。

使用递归方法实现此算法：

首先，我们先计算出a和b的差值c（假设a>b），把问题转化成求出b和c的最大公约数；然后计算出c和b的差值d（假设c>b），把问题转化成求出b和d的最大公约数；再然后计算出b和d的差值e（假设b>d），把问题转化成求出d和e的最大公约数…
以此类推，逐渐把两个较大整数之间的运算简化成两个较小整数之间的运算，直到两个数可以相等为止，最大公约数就是最终相等的两个数。

int grd(int a ,int b){
	if(a==b) return b;
	else
	{
		if(a>b)
		{
			return grd(a-b,b);
		}
		else
		{
			return grd(b-a,a);
		}
	}
	//返回值即为最大公约数 
}

3）Stein算法（结合辗转相除法和更相减损法的优势以及移位运算）

当俩个数相差悬殊时，比如1和10000，就需要递归9999次，这时候就需要使用Stein算法。
众所周知，移位运算的性能非常快。对于给定的正整数a和b，不难得到如下的结论。其中gcb(a,b)的意思是求a,b的最大公约数的函数。

1. 当a和b均为偶数，gcb(a,b) = 2gcb(a/2, b/2) = 2gcb(a>>1, b>>1)
2. 当a为偶数，b为奇数，gcb(a,b) = gcb(a/2, b) = gcb(a>>1, b)
3. 当a为奇数，b为偶数，gcb(a,b) = gcb(a, b/2) = gcb(a, b>>1)
4. 当a和b均为奇数，利用更相减损术运算一次，gcb(a,b) = gcb(b, a-b)， 此时a-b的结果必然是偶数，又可以继续进行移位运算。

int gcd(int a,int b){
    if(a == 0) return b;
    if(b == 0) return a;
    if(a % 2 == 0 && b % 2 == 0) return 2 * gcd(a >> 1, b >> 1);
    else if(a % 2 == 0)  return gcd(a >> 1, b);
    else if(b % 2 == 0) return gcd(a, b >> 1);
    else return gcd(abs(a - b), min(a, b));
}

2.最小公倍数

先判断a,b俩个整数是否能够整除，如果不能求则a,b最大公约数c，再a*b/c即为最小公倍数。

参考文章：
https://zhuanlan.zhihu.com/p/31824895
https://zhuanlan.zhihu.com/p/338809271