题目梗概
你有一个大小为n的背包,你有n种物品,第i种物品的大小为i,且有i个,求装满这个背包的方案数有多少
解题思路
分块考虑这个问题,因为大于 n√ 的数可以无限使用。
小于 n√ 的背包用前缀和维护。
大于 n√ 的背包变换DP的定义。
f[i][j] 表示放入了i个数,和为j,那么只有两种情况:
f[i][j−i]−>f[i][j] 表示i个数都加1
f[i−1][j−n√−1]−>f[i][j] 表示新加入一个数,因为枚举的数都大于 n√ 所以新数为 n√+1
然后运用乘法原理合并答案就可以了
#include<cstdio>
#include<cmath>
#define maxn 100005
#define maxm 320
#define tt 23333333
#define nxt 1-(i&1)
#define now i&1
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,g[maxm][maxn],f[2][maxn],sum[maxn],ans;
int main(){
freopen("exam.in","r",stdin);
freopen("exam.out","w",stdout);
scanf("%d",&n);m=sqrt(n);
f[0][0]=1;
for(int i=1;i<=m;i++){
memset(sum,0,sizeof(sum));
for(int j=0;j<=n;j++){
f[now][j]=(sum[j%i]+f[nxt][j])%tt;
sum[j%i]=(sum[j%i]+f[nxt][j])%tt;
if (j>=i*i) sum[j%i]=(sum[j%i]-f[nxt][j-i*i]+tt)%tt;
}
}
g[0][0]=1;
for (int i=1;i<=m;i++)
for (int j=0;j<=n;j++){
if (j>=i) (g[i][j]+=g[i][j-i])%=tt;
if (j>=m+1) (g[i][j]+=g[i-1][j-m-1])%=tt;
}
for (int i=0;i<=m;i++)
for (int j=0;j<=n;j++)
ans=(ans+(long long)f[m&1][j]*g[i][n-j]%tt)%tt;
printf("%d\n",ans);
}