AcWing 871.约数之和（模板题）

最新推荐文章于 2025-02-27 01:46:34 发布

原创最新推荐文章于 2025-02-27 01:46:34 发布 · 364 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

enroll 专栏收录该内容

35 篇文章

订阅专栏

博客介绍了如何利用算术基本定理和秦九韶算法解决约数之和的问题。首先，通过质因数分解将数分解为质数的乘积形式，然后应用约数之和公式计算每个质因数的指数之和的乘积，最终得到所有可能的组合和。

Issue：

在这里插入图片描述

Meaning：

输出约数之和

Thinking：

用算术基本定理分解质因数；
unordered_map<int, int>一个用来存 $P$ ,一个用来存对应的 $α\alpha$ ；
记住约数之和公式： $(P10+P11+⋯+P1α1)×⋯×(Pk0+Pk1+⋯+Pkαk)(P_1^0 + P_1^1 + \cdots + P_1^{\alpha_1}) \times \cdots \times (P_k^0 + P_k^1 + \cdots + P_k^{\alpha_k})$ ；
其中的每一小项用秦九韶算法来计算：while (a -- ) t = (t * p + 1) % mod;

Code：

//记住两个公式,约数个数，约数之和
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;

const int mod = 1e9+7;

int main()
{
    int n;
    cin >> n;
    long long res = 1;
    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]++;
    }
    for(auto prime : primes)
    {
        int p = prime.first, a = prime.second;
        long long t = 1;
        while(a--)
            t = (t*p+1) % mod;    //秦九韶算法
        res = res * t % mod;
    }
    cout << res << endl;
    return 0;
}