bzoj1025 [SCOI2009]游戏——因数DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1025

这篇博客写得真好呢：https://www.cnblogs.com/phile/p/4473192.html

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,pri[],cnt;
long long f[][],ans;
bool vis[];
void init()
{
    for(int i=;i<=n;i++)
    {
        if(!vis[i])vis[i]=,pri[++cnt]=i;
        for(int j=;j<=cnt&&i*pri[j]<=n;j++)
        {
            vis[i*pri[j]]=;
            if(i%pri[j]==)break;
        }
    }
}
int main()
{
    scanf("%d",&n);
    init();
    f[][]=;
    for(int i=;i<=cnt;i++)
    {
        for(int j=;j<=n;j++)f[i][j]+=f[i-][j];
        for(int j=;j<=n;j++)
            for(int k=pri[i];j-k>=;k*=pri[i])
                f[i][j]+=f[i-][j-k];
    }
    for(int i=;i<=n;i++)
        ans+=f[cnt][i];
    printf("%lld",ans);
    return ;
}