关于线性筛的简单应用

问题

Forsaken有一个有趣的数论函数。对于任意一个数x，f(x)会返回x的最小质因子。如果这个数没有最小质因子，那么就返回0。
现在给定任意一个n，Forsaken想知道 f(1) + f(2) + … + f(n) 等于多少。

输入

一个整数 n

输出

一个整数代表求和答案

算法思想

这道题是一道比较简单的线性筛的模板应用，

代码实现

# include<bits/stdc++.h>
using namespace std;
int flag[30003000]={0};
int zuixiaoyinzi[30003000]={0};
int zhiyinzi[30003000]={0};
int main()
{
    long long int n,sum=0;
    cin>>n;
    long long int i,j=0;
    for(i=2;i<=n;i++)
    {
        if(flag[i]==0)
        {
            zhiyinzi[j]=i;
            zuixiaoyinzi[i]=i; //增加了一个数组存放每个数的最小质因子
            j++;
        }
        for(long long int k=0;zhiyinzi[k]<=n/i;k++)
        {
            flag[zhiyinzi[k]*i]=1;
            zuixiaoyinzi[zhiyinzi[k]*i]=zhiyinzi[k];
            if(zhiyinzi[k]!=0&&i%zhiyinzi[k]==0) //为了使得每个数被最小质因子筛掉 
            break;
            /*如果去掉了这一行则会导致i*zhiyinzi[m](m>n,i%zhiyinzi[n]==0)
            的最小质因子是zhiyinzi[m],但实际上，i*zhiyinzi[m]/(zhiyinzi[n]*zhiyinzi[m])
            可能是比zhiyinzi[m]更小的质因子*/
        }
    }
    for(long long int i=0;i<=n;i++)
    sum+=zuixiaoyinzi[i];
    cout<<sum<<endl;
    return 0;
}