其他分享
首页 > 其他分享> > 数学知识 约数的和

数学知识 约数的和

作者:互联网

约数之和

给定 n 个正整数 ai,请你输出这些数的乘积的约数之和,答案对 109+7 取模

输入格式
第一行包含整数 n

接下来 n 行,每行包含一个整数 ai

输出格式
输出一个整数,表示所给正整数的乘积的约数之和,答案需对 109+7 取模

数据范围
1≤ n ≤100,
1≤ ai ≤2×109
输入样例:

3
2
6
8

输出样例:

252

约数和的定理:

n = p1a1+p2a2+p3a3+p4a4+…pkak

p1a1 的约数有:p10,p11,p12,p13…p1k 共(a1+1)个;同理p2a2的约数个数有(a2+1)个。
n的约数是在p1a1、p2a2、…、pkak每一个的约数中分别挑一个相乘得来。所以共有(a1+1)*(a2+1)*(a3+1)*…(ak+1)种挑法。
所以约数的和为:
(p10+p11+p12…p1a1)*(p20+p21+p22+…p2a2)*…(pk0+pk1+pk2…pkak

eg:

23 * 32 = 72
23 含有的约数为(20,21,22,23)即(1,2,4,8)
32含有的约数为(30,31,32)即(1,3,9)
而72的约数为:
(1,2,3,4,6,8,9,12,18,24,36,72)=(20,21,22,23)*(30,31,32

#include<bits/stdc++.h>
using namespace std;
const int M = 1e9+7;
typedef long long int LL;
int n;
LL res=1;
int main()
{
    cin>>n;
    unordered_map<int,int> p;
    while(n--)
    {
        int x;
        cin>>x;
        for(int i=2;i<=x/i;i++)
        {
            while(x%i==0)
            {
                p[i]++;
                x/=i;
            }
        }
        if(x>1) p[x]++;
    }
    for(auto t:p)
    {
        LL i=t.first,k=t.second;
        LL m=1;
        while(k--)
        {
            m=(m*i+1)%M;
        }
        res=res*m%M;
    }
    cout<<res;
    return 0;
}

求( i0 + i1 + i2… + ik

while (k -- )  m=(m*i+1) % mod;

m = m*i + 1
m = 1
m = i + 1
m = i2 + i + 1
m = i3 + i2 + i + 1
……
m = ik + ik-1 +…+ i + 1

标签:约数,p2a2,23,int,32,p1a1,数学知识
来源: https://blog.csdn.net/qq_45806470/article/details/123071651