luogu3188/bzoj1190 梦幻岛宝珠 (分层背包dp)

他都告诉你能拆了 那就拆呗。把每个重量拆成$a*2^b$的形式

然后对于每个不同的b，先分开做30个背包

再设f[i][j]表示b<=i的物品中 容量为$ j*2^i+W&((1<<(i-1))-1) $（就是这一位是j+W的前i-1位）的最大权值（这个容量没必要填满）

然后f[i][j]就可以从f[i-1][j*2+W的第i-1位]转移过来，再拿着这个去更新本层的其他容量

最后答案就是f[x][1]，x是W的最高位1的位数

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=110,maxs=1e3+10,maxp=33;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

int N,W,sum[maxs];ll f[maxp][maxs];

int main(){
    //freopen(".in","r",stdin);
    int i,j,k;
    while(1){
        N=rd(),W=rd();
        if(N==-1) break;
        CLR(f,0);CLR(sum,0);
        int x=log2(W);
        for(i=1;i<=N;i++){
            int a=rd(),b=0,v=rd();
            while(a&&a%2==0) b++,a>>=1;
            if(b>x) continue;
            for(j=min(sum[b],1000-a);j>=0;j--){
                f[b][j+a]=max(f[b][j+a-1],max(f[b][j+a],f[b][j]+v));
            }
            sum[b]+=a;
        }
        
        for(i=0;i<x;i++){
            for(j=1;j<=1000;j++) f[i][j]=max(f[i][j-1],f[i][j]);
            for(k=1000;k>=0;k--){
                if(k&&!f[i+1][k]) continue; 
                for(j=1000;j>=k;j--){
                    f[i+1][j]=max(f[i+1][j],f[i+1][k]+f[i][min(1000,2*(j-k)+((W>>i)&1))]);
                }
            }
        }
        printf("%lld
",f[x][1]);
    }
    return 0;
}