题意可以理解为给定N 将其拆为
pα11∗pα22∗pα33⋅⋅⋅pαmm
的形式 其中
p
为素数且
求
lcm(pα11∗pα22∗pα33⋅⋅⋅pαmm)
的不同个数 (仔细想想很容易理解)
转化问题后就可以用DP来解决啦 令
f[i][j]
表示当前用到第
i
个素数总和为
直接贴代码
#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;
}