[HNOI2007]梦幻岛宝珠

LCVI.[HNOI2007]梦幻岛宝珠

好题。

明显它是01背包的模型，但值域过大。咋办呢？

我们考虑令 $f_{i,j}$ 表示只考虑 $a\times 2^i$ 类型的物品，关于 $a$ 做的一个背包。显然，暴力求出这个东西的时空复杂度都是可接受的。

我们再考虑 $g_{i,j}$ 表示有 $j\times2^i+w\operatorname{and}(2^i-1)$ 这么多的背包容量时的答案，即只考虑 $w$ 的下 $i$ 位，且第 $i$ 位上选了 $j$ 单位的物品。我们考虑由 $g_{i,j}$ 更新 $g_{i+1}$。

显然，如果 $w$ 的第 $i$ 位有一个 $1$，在第 $i+1$ 位上，$g_{i,j}$ 就与一个大小为 $\left\lfloor\dfrac{j}{2}\right\rfloor$ 的物品等价；否则，即 $w$ 的第 $i$ 位没有 $1$，它与 $\left\lceil\dfrac{j}{2}\right\rceil$ 等价。

于是我们就用 $g_i$ 和 $f_{i+1}$ 即可拼凑出 $g_{i+1}$。

代码：

#include<bits/stdc++.h>
using namespace std;
int n,m;
vector<pair<int,int> >v[40];
int f[40][2000],lim[40],g[40][2000];
void chmx(int &x,int y){if(x<y)x=y;}
int main(){
	while(true){
		scanf("%d%d",&n,&m);
		if(n==-1&&m==-1)break;
		memset(f,0,sizeof(f)),memset(g,0,sizeof(g)),memset(lim,0,sizeof(lim));
		for(int i=0;i<=30;i++)v[i].clear();
		for(int i=1,a,b,c;i<=n;i++){
			scanf("%d%d",&a,&b),c=0;
			while(!(a&1))a>>=1,c++;
			v[c].push_back(make_pair(a,b));
		}
		for(int i=0;i<=30;i++){
			for(auto k:v[i]){
				for(int j=lim[i];j>=0;j--)chmx(f[i][j+k.first],f[i][j]+k.second);
				lim[i]+=k.first;
//				printf("(%d,%d)\n",k.first,k.second);
			}
//			for(int j=0;j<=lim[i];j++)printf("%d ",f[i][j]);puts("");
		}
		for(int i=0;i<=lim[0];i++)g[0][i]=f[0][i];
		for(int i=0;i<=30;i++){
			for(int j=0;j<=lim[i];j++)for(int k=lim[i+1];k>=0;k--)chmx(g[i+1][k+((j+!((m>>i)&1))>>1)],f[i+1][k]+g[i][j]);
			lim[i+1]+=(lim[i]+!((m>>i)&1))>>1;
		}
		printf("%d\n",g[31][0]);
	}
	return 0;
}