很妙的题。考虑到大于
n−−√
的物品,取不到数量上限就一定超出背包大小了,就可以看成无限背包的问题。且取的物品总数一定小于
n−−√
。所以就可以
nn−−√
DP搞了。
而
n−−√
之内的物品,直接看成一般的多重背包问题,用按模分类的方法优化,可以去掉个数的复杂度,为
nn−−√
.
两半求好后统计一下即可。
所以总复杂度
O(nn−−√)
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int maxn=100005,MOD=23333333;
int n,m,f1[maxn],sum[maxn],ans,f[330][maxn],f2[maxn];
int main(){
freopen("51nod1597.in","r",stdin);
freopen("51nod1597.out","w",stdout);
scanf("%d",&n); m=sqrt(n);
f1[0]=1;
for(int i=1;i<=m;i++)
for(int j=0;j<=i-1;j++){
sum[0]=f1[j];
for(int k=1;j+k*i<=n;k++) sum[k]=(sum[k-1]+f1[j+k*i])%MOD, (f1[j+k*i]+=sum[k-1]-(k-i-1>=0?sum[k-i-1]:0))%=MOD;
}
f[0][0]=1;
for(int i=1;i<=m;i++)
for(int j=i*(m+1);j<=n;j++)
(f[i][j]+=f[i][j-i]+f[i-1][j-(m+1)])%=MOD;
f2[0]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++) (f2[i]+=f[j][i])%=MOD;
for(int i=0;i<=n;i++) (ans+=(LL)f1[i]*f2[n-i]%MOD)%=MOD;
printf("%d\n",(ans+MOD)%MOD);
return 0;
}