约数、最大公约数

1. 约数

约数，又称因数。整数 $a$ 除以整数 $b(b\ne 0)$ 除得的商正好是整数而没有余数，我们就说 $a$ 能被 $b$ 整除，或 $b$ 能整除$a$。$a$ 称为 $b$ 的倍数，$b$ 称为 $a$ 的约数。

2. 试除法求约数

试除法求约数的原理和试除法判断质数的原理基本一致，但是注意求约数的时候要从 $1$ 开始，时间复杂度是 $O(sqrt(n))$。

#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;

// 求 x 的所有约数
vector<int> get_divisors(int x)
{
	vector<int> res;
	for (int i = 1; i <= x / i; i ++ )
		if (x % i == 0)
		{
			res.push_back(i);
			if (i != x / i) res.push_back(x / i);
		}
	sort(res.begin(), res.end());
	return res;
}

3. 最大公约数—欧几里得算法（辗转相除法）

首先，如果 $d|a$，并且 $d|b$，那么 $d|xa+yb$

我们先直接上结论，记 $x,y$ 的最大公约数为 $gcd(x,y)$，则 $gcd(a,b)=gcd(b,a\%b)$。

证明：
因为 $a\%b=a-\lfloor \frac{a}{b}\rfloor b$，因为 $\lfloor \frac{a}{b}\rfloor$ 是一个整数，所以我们可以直接将其记为 $c$。因此，如果 $gcd(a,b)=gcd(b,a\%b)$ 成立，则有 $gcd(a,b)=gcd(b,a-c\cdot b)$ 成立，所以我们只要证明 $gcd(a,b)=gcd(b,a-c\cdot b)$ 成立即可。

对于 $(a,b)$ 中的任何一个约数 $d$，都有 $d|a$ 和 $d|b$，而 $d|a=d|a-c\cdot b+c\cdot b$，又因为 $d|c\cdot b$ 所以 $d|a-c\cdot b$，因此 $(a,b)$ 与 $(b,a\%b)$ 的公约数集合是相同的，所以说 $gcd(a,b)=gcd(b,a-c\cdot b)$。

int gcd(int a, int b)
{
	// gcd(a, 0) = a，因为 0 可以整除任何一个数
	return b ? gcd(b, a % b) : a;
}