P3188 [HNOI2007]梦幻岛宝珠

传送门

注意到 $a,b$ 不大

考虑对每一个 $a*2^b$ 的 $b$ 分别背包

设 $f[i][j]$ 表示只考虑 $b=i$ 的物品时，容量为 $j= \sum a$ 的最大价值

这个就是普通的 $01$ 背包

考虑把 $f[i][j]$ 之间合并起来，为了得到容量为 $W$ 时的答案，我们要把 $f$ 的含义稍微变化一下

变成 $f[i][j]$ 表示当前考虑了 $b=2^1$ 到 $b=2^i$ 时的物品，容量为 $j*2^i$ 加上 $W$ 二进制下前 $i-1$ 位的值，此时的最大价值

考虑用 $f[i-1][]$ 更新 $f[i][j]$，枚举总体积为 $(j-k)*2^i$ 的 $b=2^i$ 的物品的最大价值（$f[i][j-k]$）

加上总体积为 $2k*2^{i-1}$ 的 $b<2^i$ 的物品的最大价值 （$f[i-1][k*2]$），注意我们还要考虑 $W$ 的容积，所以设 $W$ 第 $i-1$ 位为 $p$

那么转移为 $f[i][j]=f[i][j-k]+f[i-1][ min(sw[i-1],k*2+p) ]$，此时 $sw[i]$ 表示第 $b<=2^i$ 时物品的体积和上取整为 $2^{sw[i]}$

注意上面枚举 $j$ 的时候要从大到小枚举

转移的细节挺多的...，最终答案即为 $f[m][1]$， $m$ 表示 $W$ 的最高位

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=,M=;
int n,W,m,sw[];
struct Orb{
    int val,w;
};
vector <Orb> V[];
ll f[][M];
int main()
{
    while()
    {
        n=read(),W=read(); int w,v;
        if(n==-&&W==-) break;
        for(int i=;i<=;i++) V[i].clear(),sw[i]=;
        memset(f,,sizeof(f));
        for(int i=;i<=n;i++)
        {
            w=read(),v=read(); int cnt=;
            while(!(w&)) w>>=,cnt++;
            V[cnt].push_back((Orb){v,w}); sw[cnt]+=w;
        }
        int t=,cnt=; while(t<=W) t<<=,cnt++;
        m=cnt-;
        for(int i=;i<=m;i++)
        {
            int len=V[i].size();
            for(int j=;j<len;j++)
                for(int k=sw[i];k>=V[i][j].w;k--)
                    f[i][k]=max(f[i][k],f[i][k-V[i][j].w]+V[i][j].val);
        }
        // f[i][j]=f[i][j-k]+f[i-1][ (k<<1) | ( (W>>(i-1)) &1) ]
        for(int i=;i<=m;i++)
        {
            sw[i]+=(sw[i-]+)>>; int p=(W>>(i-))&;
            for(int j=sw[i];j>=;j--)
                for(int k=;k<=j;k++)
                    f[i][j]=max(f[i][j],f[i][j-k]+ f[i-][min(sw[i-],(k<<)|p)] );
        }
        printf("%lld\n",f[m][]);
    }
    return ;
}