约数个数和约数之和算法总结

ykycode

已于 2024-01-07 18:38:02 修改

阅读量648

点赞数 12

分类专栏：经典算法总结数论文章标签：算法数学试除法约数个数约数之和数论约数

于 2024-01-07 14:00:52 首次发布

本文链接：https://blog.csdn.net/u012181348/article/details/135438918

版权

经典算法总结同时被 2 个专栏收录

44 篇文章 0 订阅

订阅专栏

数论

17 篇文章 0 订阅

订阅专栏

知识概览

约数个数

基于算数基本定理，假设N分解质因数的结果为

$N = p_1^{\alpha_{1}}p_2^{\alpha_{2}} \cdots p_k^{\alpha_{k}}$

可得对于N的任何一个约数d，有

$d = p_1^{\beta_{1}}p_2^{\beta_{2}} \cdots p_k^{\beta_{k}}, 0\leqslant \beta_{i}\leqslant \alpha_{i}$

因为N的每一个约数和 $\beta_{1}$ ~ $\beta_{k}$ 的一种选法是一一对应的，根据乘法原理可得，

一个数的约数个数为

$(\alpha_{1} + 1)(\alpha_{2} + 1) \cdots (\alpha_{k} + 1)$

约数之和

一个数的约数之和公式为

$(p_1^{0} + p_1^{1} + \cdots + p_1^{\alpha_1})\cdots (p_k^{0} + p_k^{1} + \cdots + p_k^{\alpha_k})$

多项式乘积的每一项为

$p_1^{\beta_{1}}p_2^{\beta_{2}} \cdots p_k^{\beta_{k}}$

正好对应的是一个数的每一个约数。

例题展示

约数个数

题目链接

活动 - AcWing系统讲解常用算法与数据结构，给出相应代码模板，并会布置、讲解相应的基础算法题目。https://www.acwing.com/problem/content/872/

代码

#include <iostream>
#include <algorithm>
#include <unordered_map>

using namespace std;

typedef long long LL;

const int mod = 1e9 + 7;

int main()
{
    int n;
    cin >> n;
    
    unordered_map<int, int> primes;
    while (n--)
    {
        int x;
        cin >> x;
        
        for (int i = 2; i <= x / i; i++)
            while (x % i == 0)
            {
                x /= i;
                primes[i]++;
            }
            
        if (x > 1) primes[x]++;
    }
    
    LL res = 1;
    for (auto prime : primes) res = res * (prime.second + 1) % mod;
    
    cout << res << endl;
    
    return 0;
}

约数之和

题目链接

活动 - AcWing系统讲解常用算法与数据结构，给出相应代码模板，并会布置、讲解相应的基础算法题目。https://www.acwing.com/problem/content/873/

代码

#include <iostream>
#include <algorithm>
#include <unordered_map>

using namespace std;

typedef long long LL;

const int mod = 1e9 + 7;

int main()
{
    int n;
    cin >> n;
    
    unordered_map<int, int> primes;
    while (n--)
    {
        int x;
        cin >> x;
        
        for (int i = 2; i <= x / i; i++)
            while (x % i == 0)
            {
                x /= i;
                primes[i]++;
            }
            
        if (x > 1) primes[x]++;
    }
    
    LL res = 1;
    for (auto prime : primes)
    {
        int p = prime.first, a = prime.second;
        LL t = 1;
        while (a--) t = (t * p + 1) % mod;
        res = res * t % mod;
    }
    
    cout << res << endl;
    
    return 0;
}