BZOJ1025 [SCOI2009]游戏 DP

题意可以理解为给定N 将其拆为$p_1^{α_1}*p_2^{α_2}*p_3^{α_3}···p_m^{α_m}$的形式 其中$p$为素数且 $\sum p_i^{α_i}\leq n$
求$lcm(p_1^{α_1}*p_2^{α_2}*p_3^{α_3}···p_m^{α_m})$的不同个数 （仔细想想很容易理解）

转化问题后就可以用DP来解决啦 令$f[i][j]$表示当前用到第$i$个素数总和为$j$我们枚举当前使用的$i$的指数$α$来转移 方程为$f[i][j]=sum\{ f[i-1][j-prime[i]^α]\}(prime[i]^α\leq j)$

直接贴代码

#include<bits/stdc++.h>
#define bug(x) cout<<(#x)<<" "<<(x)<<endl
#define ll long long
/*
char *TT,*mo,but[(1<<15)+2];
#define getchar() ((TT==mo&&(mo=(TT=but)+fread(but,1,1<<15,stdin),TT==mo))?-1:*TT++)//*/
using namespace std;
const int N=1e5+5;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
ll ans;
int n,cnt;
ll f[1005][1005];
int p[1005];
bool vis[1005];
void pre(){
    for(int i=2;i<=n;i++){
        if(!vis[i]) p[++cnt]=i;
        for(int j=1;j<=cnt&&i*p[j]<=n;j++){
            vis[i*p[j]]=1;
            if(i%p[j]==0) break;
        }
    } 
}
int main(){
#ifdef Devil_Gary
    freopen("in.txt","r",stdin);
#endif
    n=read(),pre(),f[0][0]=1;
//  bug(cnt);
    for(int i=1;i<=cnt;i++) for(int j=0;j<=n;j++){
        f[i][j]+=f[i-1][j];
        for(int k=p[i];k<=j;k*=p[i])
            f[i][j]+=f[i-1][j-k];
        }
    for(int i=0;i<=n;i++) ans+=f[cnt][i];
    cout<<ans<<endl;
}