Acwing-基础算法课笔记之数学知识（约数）

一、试除法求约数

1、概念

$\bullet$利用试除法即可
如果我们要枚举较小的约数，则$d\le \frac{n}{d}$，即$d\le \sqrt{n}$，那么较大的约束可以直接算出来。

2、代码模板

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;
}

Acwing-试除法求约数

二、约数个数

1、概述

利用短除法分解质因数
设要求的数是$N$
利用短除法分解质因数得$N=p_1^{n_1}\cdot p_2^{n_2}\cdot p_3^{n_3}\cdot \cdot \cdot \cdot \cdot p_n^{n_{\alpha }}$
所以$N$的约数个数为$(n_1+1)\cdot (n_2+1)\cdot (n_3+1)\cdot \cdot \cdot \cdot \cdot (n_{\alpha }+1)$

2、过程模拟

要求360的约数个数，步骤如下：
$360=2^3\times 3^2\times 5^1$
所以$360$的约数个数是$(3+1)\times (2+1)\times (1+1)=24$

三、约数之和

1、概述

利用短除法分解质因数
设要求的数是$N$
利用短除法分解质因数得$N=p_1^{n_1}\cdot p_2^{n_2}\cdot p_3^{n_3}\cdot \cdot \cdot \cdot \cdot p_n^{n_{\alpha }}$
所以$N$的约数之和为$(p_1^0+p_1^1+p_1^2+\cdot \cdot \cdot +p_1^{n_1})\cdot (p_2^0+p_2^1+p_2^2+\cdot \cdot \cdot +p_2^{n_2})\cdot (p_3^0+p_3^1+p_3^2+\cdot \cdot \cdot +p_3^{n_3})\cdot \cdot \cdot \cdot \cdot (p_n^0+p_n^1+p_n^2+\cdot \cdot \cdot +p_n^{n_{\alpha }})$

2、过程模拟

要求360的约数个数，步骤如下：
$360=2^3\times 3^2\times 5^1$
所以$360$的约数个数是$(2^0+2^1+2^2+2^3)\times (3^0+3^1+3^2)\times (5^0+5^1)$

四、最大公约数（欧几里得算法即辗转相除法）

1、举例

$\bullet$欧几里得在几何原本当中举的例子（铺瓷砖）
例如，求$(6,16)$的最大公约数，步骤如下：

所以，$(6,16)$的最大公约数是$2$。

$\bullet$当两个数比较大时，用辗转相除法
例如，求$(12921,4234)$的最大公约数，步骤如下：

所以，$(12921,4234)$的最大公约数是$73$。

由此可知，设$A$为被除数，$B$为除数，$Q$为商数，$R$为余数
所以，$(A,B)=(B,R)$

2、证明欧几里得算法

$gcd(a,b)=gcd(b,amodb)$

不妨设$a>b$，令$a=k\cdot b+c$，其中$k$是商数，$c$是余数，$c=amodb$
不妨设$d=gcd(a,b)$，则设$\{\begin{array}{l}a=p_1\cdot d\\ b=p_2\cdot d\end{array}$，（$p_1$，$p_2$互质）

代入得，$p_1\cdot d=k\cdot p_2\cdot d+c$$\Rightarrow$$c=d\cdot (p_1-k\cdot p_2)$

结论：$gcd(p_1\cdot d,p_2\cdot d)=gcd(p_2\cdot d,d\cdot (p_1-k\cdot p_2))=d$

证明：$p_2$与$p_1-k\cdot p_2$互质即可（反证法）

若$p_2$与$p_1-k\cdot p_2$不互质，则设其最大公因数为$t$，设$\{\begin{array}{l}{p}_1=t\cdot m，(m,n互质)\\ p_1-k\cdot p_2=t\cdot n，（t\ge 2）\end{array}$

所以，$p_1-k\cdot t\cdot m=t\cdot n$，$\{\begin{array}{l}{p}_1=t\cdot (n+k\cdot m)\\ p_2=t\cdot m\end{array}$

所以，$p_1=t\cdot (n+k\cdot m)$，$\{\begin{array}{l}a=p_1\cdot d=t\cdot d\cdot (n+k\cdot m)\\ b=p_2\cdot d=t\cdot m\cdot d=t\cdot d\cdot m\end{array}$

所以，$gcd(a,b)=t\cdot d\Rightarrow x$。

3、代码模板

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