洛谷-1403 [AHOI2005]约数研究

最新推荐文章于 2022-05-08 11:22:13 发布

mkopvec

最新推荐文章于 2022-05-08 11:22:13 发布

阅读量130

点赞数

分类专栏：算法题

本文链接：https://blog.csdn.net/mkopvec/article/details/90478708

版权

算法题专栏收录该内容

311 篇文章 3 订阅

订阅专栏

题目描述
科学家们在Samuel星球上的探险得到了丰富的能源储备，这使得空间站中大型计算机“Samuel II”的长时间运算成为了可能。由于在去年一年的辛苦工作取得了不错的成绩，小联被允许用“Samuel II”进行数学研究。
小联最近在研究和约数有关的问题，他统计每个正数N的约数的个数，并以f(N)来表示。例如12的约数有1、2、3、4、6、12。因此f(12)=6。下表给出了一些f(N)的取值：
f(n)表示n的约数个数，现在给出n，要求求出f(1)到f(n)的总和。
输入输出格式
输入格式：
输入一行，一个整数n
输出格式：
输出一个整数，表示总和

输入输出样例
输入样例#1：
3

输出样例#1：
5

说明
【数据范围】
20%N<=5000
100%N<=1000000

解释： $s u m [i$ ]为 $i$ 的因数个数，那么如果i是素数那就是 $s u m [i] = 2, 否则 s u m [i] = s u m [t e m p] * (1 + n u m), 设 p 是 i 最小的素因$ 子， $t e m p 就是去掉 p 的数， n u m 是 p 的指数个数$

#include<cstdio>
#include<iostream>
#include<cstring>
#define N 1000004
using namespace std;
int fact[N]={0};
bool Prime[N]={0};
int P[N]={0};
int sum[N]={0};
int t=0;
int n=0;
void make_prime(){
    memset(Prime,true,sizeof(Prime));
    Prime[0]=Prime[1]=false;
    for(int i=2;i<N;i++){
        if(Prime[i]) P[t++]=i;
        for(int j=0;j<t&&P[j]*i<N;j++){
            Prime[P[j]*i]=false;
            fact[P[j]*i]=P[j];
            if(i%P[j]==0) break;
        }
    }
}
int ret=1;
int main(){
    make_prime();
    cin>>n;
    sum[1]=1;
    for(int i=2;i<=n;i++){
        if(Prime[i]) sum[i]=2;
        else{
            int num=0;
            int temp=i;
            while(temp%fact[i]==0){
                temp/=fact[i];num++;
            }
            sum[i]=sum[temp]*(num+1);
        }
        ret+=sum[i];
    }
    cout<<ret<<endl;
    return 0;
}