题解 P3188 【[HNOI2007]梦幻岛宝珠】 DP

~By Bartholomew~

其实,这道题目,评价:

1.难度: @@@@@

2.思维难度: @@@

3.细节难度: @@@@@

对于这道题目,我们发现其实就是普通的 01 背包,但是数据 十分的大!

那么我们可以对于 v = a×2b;的 b 相同的物品做一次 dp

那么就是* f[i][j] * 表示体积为 (j×2i) 的体积 只收纳 k×2i 的物品的最大值

那么最后就是 怎么合并是关键部分!

如果我们假设 f[i][j] 现在表示

表示20 ~~ 2i组内， 容量为 j×2i+(w&((1<< i)−1)) 时的最大价值

也就是f[i,j]表示体积为j× 2i再加上W**二进制第i位以下的体积**最多可以获得多少价值

那么 我们的 答案其实就是 dp[len][1] , len指的是 最大的i 是 2i 小于等于 m 的值!

那么 我们就是考虑转移 从 f[i−1][...] 开始转移!

转移方程:

f[i,j]=max(f[i,j],f[i,j−k]+f[i−1],min(k∗2+e[i−1],d[i−1])]);

其中好多是小细节部分:

for(int j=w[i];j>=0;--j)
for(int k=0;k<=j;k++)   // 这下面的f[i][j-k] (用"*"标注的)1其实并没有变动,因为是反着来求的,所以还是表示 单单是 2^i 的物品,大家可以写一下式子,就明白了!
f[i][j] = getmax(f[i][j] , **f[i][j-k]** + f[i-1][min(w[i-1],(k<<1)|(m>>(i-1)&1))]);

贴上代码 , 有明显的小错误的哦!

#pragma GCC optimize(3)
#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <cstring>
#include <queue>
#include <vector>
#define N 200005
#define INF 0x3f3f3f3f
using namespace std;
int n,m,len,x,y,f[35][1150],w[35];
inline int read()
{
int x=0;
char c=getchar();
while(c<'0'||c>'9'){c=getchar();}
while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}
return x;
}
vector<int> G[35],V[35];
inline int getmax(int x,int y){return (x>y)?x:y; }
inline void init()
{
len = 0;
memset(G,0,sizeof G);
memset(V,0,sizeof V);
memset(w,0,sizeof w);
memset(f,0,sizeof f);
}
int main(int argc,char const* argv[])
{
while(~scanf("%d%d",&n,&m) && n!=-1)
{
init();
for(int i=1;i<=n;++i)
{
x = read(); int j = 0;
while(!(x&1)) { x>>=1; j++; }
G[j].push_back(x); w[j] += x;
len = getmax( len , j);
scanf("%d",&y); V[j].push_back(y);
}
for(int i=0;i<=len;i++)
for(int j=0;j<(int)G[i].size();++j)
for(int k=w[i];k>=G[i][j];k--)
f[i][k] = getmax(f[i][k],f[i][k-G[i][j]]+V[i][j]);
while(m>>len) len++; len--;
for(int i=1;i<=len;i++)
{
w[i] += (w[i-1]+1)/2;
for(int j=w[i];j>=0;--j)
for(int k=0;k<=j;k++)
f[i][j] = getmax(f[i][j] , f[i][j-k] + f[i-1][min(w[i-1],(k<<1)|(m>>(i-1)&1))]);
}
printf("%d\n",f[len][1]);
}
return main();
}