最大公约数GCD&最小公倍数LCM

0x32约数

最大公约数

  若自然数 $d$ 同时是自然数 $a$ 和 $b$ 的约数，则称 $d$ 是 $a$ 和 $b$ 的公约数。在所有 $a$ 和 $b$ 的公约数中最大的一个，称为 $a$ 和 $b$ 的最大公约数，记为 $gcd(a,b)$。
  同理，我们也可以定义三个数以及多个数的最大公约数。

GCD(欧几里德算法)

欧几里德算法又称辗转相除法，是指用于计算两个非负整数 $a,b$ 的最大公约数。

辗转相除法的关键在于如下恒等式：$gcd(a,b)=gcd(b,amodb)$。它的边界条件是 $gcd(a,0)=a$;
我们可以得到如下程序

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

我们也可以直接使用C++的内置函数。
在头文件 $algorithm$ 中

#include <algorithm>
int gcd(int a,int b){
	return __gcd(a,b);
}

gcd函数的递归层数不超过 4.7851gN + 1.6723，其中 N = max(a,b)。所以不会栈溢出。

最小公倍数

若自然数 $d$ 同时是自然数 $a$ 和 $b$ 的倍数，则称 $d$ 是 $a$ 和 $b$ 的公倍数。在所有 $a$ 和 $b$ 的公倍数中最小的一个，称为 $a$ 和 $b$ 的最小公倍数，记为 $lcm(a,b)$。
  同理，我们也可以定义三个数以及多个数的最小公倍数。

LCM

利用 gcd 还可以求出两个整数的最小公倍数 $lcm(a,b)$。
由一分定理得：

• $gcd(a,b)\ast lcm(a,b)=a\ast b$

我们可以得到如下程序

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

注意：在这里不能写成 $a\ast b/gcd(a,b)$ 因为 $a\ast b$ 可能会产生溢出。而先除后乘可以避免这种情况发生。