[DP]51 Nod 1597——有限背包计数问题

题目梗概

你有一个大小为n的背包，你有n种物品，第i种物品的大小为i，且有i个，求装满这个背包的方案数有多少

解题思路

分块考虑这个问题，因为大于$\sqrt n$的数可以无限使用。

小于$\sqrt n$的背包用前缀和维护。

大于$\sqrt n$的背包变换DP的定义。

$f[i][j]$表示放入了i个数，和为j，那么只有两种情况:

$f[i][j-i]->f[i][j]$表示i个数都加1

$f[i-1][j-\sqrt n-1]->f[i][j]$表示新加入一个数，因为枚举的数都大于$\sqrt n$所以新数为$\sqrt 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);
}