质因子个数统计

#include<cstdio>
#include<cstring>

#define Max 1000000

int prime[Max+5];
int dp[Max+5];
void prime_init()
{
    memset(prime,0,sizeof(prime));
    prime[0]=prime[1]=1;
    int i,j;
    for(i=2;i<=Max;i++)
    {
        if(prime[i]!=0) continue;//判断是否为素数
        if(i>Max/i) continue;//防止越界
        for(j=i*i;j<=Max;j+=i) prime[j]=i;//保存i的其中一个因子，并标记为非素数
    }
}
void dp_init()//dp[i]保存i的质因子个数
{
    int i,j;
    memset(dp,0,sizeof(dp));
    for(i=2;i<=Max;i++)
    {
        if(prime[i]==0) dp[i]=1;
        else dp[i]=dp[prime[i]]+dp[i/prime[i]];//状态转移
    }
}
int main()
{
    prime_init();
    dp_init();
    int i;
    for(i=2;i<=Max;i++) printf("dp[%d]=%d\n",i,dp[i]);
    return 0;
}//YSF友情赞助!!!!