数论--约数相关知识及代码实现

食用这篇文章前建议有素数相关的知识，若不了解的可以点开链接快速学习 数论 – 质数判定及其筛法求解

基本概念

• 设 a 和 b 是两个整数，如果 d | a 且 d | b，则称 d 为 a 与 b 的 $公约数$。

• 除 0 之外 ，任何整数只有有限个因子。因而两个不全为 0 的整数只有有限个公因子，其中最大的称作$最大公约数$

试除法求约数

• 给定一个数 N ，求 [1 ，N ]的所有约数

$注意$

$n为完全平方数的时候sqrt(n)|n，但此时sqrt(n)和n/sqrt(n)只能算作一个约数，因为他俩相等。此算法时间复杂度为固定值O(\surd n)$

代码模板

vector<int> get_dividors(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()); // 排序
}

约数个数

$算术基本定理$

任何一个大于 1 的自然数 N， 如果 N 不为 质数，那么 N 可以 唯一分解成有限个质数的乘积

设 a > 1 ，则
$$a=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdot \cdot \cdot p_k^{{\alpha }_k}$$
其中

​ $p_1,p_2,\cdot \cdot \cdot \cdot ,p_k$是不同的素数

​ ${\alpha }_1,{\alpha }_2,\cdot \cdot \cdot \cdot ,{\alpha }_k$是正整数

并且在不记顺序的情况下，该表示是唯一的

$约数个数定理$

对于 N分解质因数 ：
$$N=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdot \cdot \cdot p_k^{{\alpha }_k}$$
由 约数定义 可知 ，对于 $ p_1^{α_1}$ 来说 ，它的约数 为
$$p_1^0，p_1^1\cdot \cdot \cdot p_1^{{\alpha }_1}$$
共 $（{\alpha }_1+1）$个

同理 对于 N 其他的质因数来说 ，它们的约数 为 $...({\alpha }_n+1)$ 个

故根据乘法原理可知：

$$N的约数的个数为({\alpha }_1+1)({\alpha }_2+1)\cdot \cdot \cdot ({\alpha }_n+1)$$

代码模板

LL res ;
const long long mod = // 很大的一个数，用来保证答案在此范围
    
void get_primes(int x)
{
    for(int i = 2; i <= x/i ; i ++ )
        if(x % i == 0)
        {
            int s = 0;
            
            while(x % i == 0) s ++, x /= i; // 计算 质因子 的个数 即 某个 α
            
            res *= (LL) (s + 1) % mod; // N 的 约数的个数为 (α1 + 1)(α2 + 1)···(αn + 1)
        }
    
    if(x > 1) res *= (LL) (1 + 1) % mod;// x 若 最后为 大于 1 的质数，约数的个数为 1 
    
}

约数之和

$约数和定理$

对于 N 分解质因数 ：
$$N=p_1^{{\alpha }_1}\ast p_2^{{\alpha }_2}\ast \cdot \cdot \cdot p_k^{{\alpha }_k}$$
由约数个数定理可知 N 的 正约数个数有
$$({\alpha }_1+1)({\alpha }_2+1)\cdot \cdot \cdot ({\alpha }_n+1)$$
那么 N 的 d（n）个正约数的和为
$$\delta (n)=(p_1^0+p_2^1+\cdot \cdot \cdot p_1^{\alpha 1})\ast \cdot \cdot \cdot \ast (p_n^0+p_n^1+\cdot \cdot \cdot p_n^{{\alpha }_n})$$

$秦九韶算法$

计算一元n次多项式：
$$\begin{gathered} f(x)=a_n{x}^n+a_{n-1}{x}^{n-1}+a_{n-2}{x}^{n-2}+...a_1{x}^1+c \\ =(...(a_{n}x+a_{n-1}x+a_{n-2})x...+a_1)x+c \end{gathered}$$

秦九韶算法实现

秦九韶算法的实现只是一个模拟算法的过程。假设 $a[i]$ 数组输入了从多项式在次数为$\ i\ $那项的系数，那么秦九韶算法的核心代码块就是：

LL a[i];
void  qjs()
{
    for(int i = n - 1; i >= 1; i -- )
        ans *= x, ans += a[i]; // 也可以写成 ans =(ans * x + a[i])
    
}

代码模板

LL get_sum(int a, int b) // a 为 p ， b 为各种 α 
{
    
    LL t = 1 ;
    
    while(b --)
    {
        t = (t * a + 1) % mod;// 这里 x 为 p ， a【i】 恒为 1 
    }
    
}

---