HDU2844 Coins【多重背包新姿势】

题目描述：

n种硬币，面值$A_i$，个数$C_i$，求可以凑出$[1,m]$中多少个数。
$n\le 100,m\le 100000,A_i\le 100000,C_i\le 1000$

题目分析：

题面显然是要求一个多重背包。

• 方法一：二进制优化。$O(nmlogC)$
  将$C_i$个物品拆成$2^0+2^1+2^2...+rest$形式的物品做01背包。
• 方法二：记录凑出这个数最少需要多少当前物品。$O(nm)$
  此题是求是否可以凑出，在添加一个物品时记录当前状态用过多少次这个物品记为$cnt[x]$，如果$x$在上一轮没有被凑出而在这一轮$x-A[i]$已经被凑出并且$cnt[x-A[i]]<C[i]$，那么$x$在这一轮就可以被凑出，同时令$cnt[x]=cnt[x-A[i]]+1$。
• 方法三：统计方案数。$O(nm)$（常数比方法二大）
  记这一轮的背包数组为$g$，上一轮为$f$，那么$g[x]=f[x]+g[x-A[i]]-f[x-(C[i]+1)\ast A[i]]$，相当于减去用了$C[i]+1$个当前物品的方案数。此题没有模数可以采用自然溢出。\
• 方法四：单调队列优化（$O(nm)$且是通用方法）
  注意到$f_t[x]=\sum _{k=1}^{C_i}{f}_{t-1}[x-k\ast A_i]$（$\sum$是或的意思，找不到好的符号了，$t$表示考虑前$t$个物品），对于$xmodk$余数相同的$x$可以看做连续的序列，那么对于$xmodk=y$的每个$y$，存当前最大的$f_{t-1}[pk+y]$的$p$，然后就可以$O(1)$检验了。

Code1：见这里
Code2：

#include<bits/stdc++.h>
#define maxn 115
#define maxm 100015
using namespace std;
int n,m,A[maxn],C[maxn],cnt[maxm];
bool f[maxm];
int main()
{
    while(~scanf("%d%d",&n,&m)&&n+m){
        for(int i=1;i<=n;i++) scanf("%d",&A[i]);
        for(int i=1;i<=n;i++) scanf("%d",&C[i]);
        memset(f,0,sizeof f),f[0]=1;
        int ans=0;
        for(int i=1;i<=n;i++){
            memset(cnt,0,(m+1)<<2);
            for(int j=A[i];j<=m;j++)
                if(!f[j]&&f[j-A[i]]&&cnt[j-A[i]]<C[i])
                    f[j]=1,cnt[j]=cnt[j-A[i]]+1,ans++;
        }
        printf("%d\n",ans);
    }
}

Code3：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 105
#define maxm 100005
using namespace std;
int n,m,A[maxn],C[maxn];
int a[2][maxm],*f=a[0],*g=a[1];
int main()
{
    while(~scanf("%d%d",&n,&m)&&n){
        for(int i=1;i<=n;i++) scanf("%d",&A[i]);
        for(int i=1;i<=n;i++) scanf("%d",&C[i]);
        memset(f,0,(m+1)<<2),f[0]=1;
        for(int i=1;i<=n;i++,swap(f,g)){
            memcpy(g,f,A[i]<<2);
            for(int j=A[i];j<=m;j++)
                g[j]=g[j-A[i]]+f[j]-(j>=(C[i]+1)*A[i]?f[j-(C[i]+1)*A[i]]:0);
        }
        int ans=0;
        for(int i=1;i<=m;i++) if(f[i]) ans++;
        printf("%d\n",ans);
    }
}