【bzoj 1190】梦幻岛宝珠（DP）

这题是在01背包问题的基础上，扩充了重量，需要用时间换空间。

思路：

1.仔细看题，注意到重量wi为a*2^b（a<=10，b<=30），很容易想到要按 b 分开做背包的DP。接下来的重点就是怎么使DP从b-1继承到b。

2.再仔细看题，发现只有一次询问，那么就可以在这个W上做文章——依W的大小进行所有的DP。由于是按b的大小进行DP，所以把数都看成二进制位来继承和转移，每对一组 b 相同的数做完后，就略去这一二进制位。

实现：

1.“去位”：若W的二进制数上第 i 位为1，那么这时背包里已经放的物品的重量的第 i 位为0或1都无所谓，装得了；若W的二进制数上第 i 位为0，那么背包已有的重量的第 i 位为1就要从二进制的前一位借一个“1”拿来用。

2.对每组 b 进行DP时，直接用最大的背包容量N*wi（mW），之后转移时再按W进行修改就行了。

其中，以上实现可行的原因是：

对于每个第 i 位，它处于一个很“尴尬”的境地，它的前一位一给就是给“1”，对于它就是“2”了，多了也无需退回去。因为对于每个第 j 位，它的下一位根本就帮不了它，贡献“1”对于它也只是“0.5”，是没有用的。因此要略去某一位时，只需干脆地按W是否能满足它这一位的需求来对自己的前一位进行调整。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long LL;
const int N=110,M=30,mW=1000;
int n,W;
struct node{int w,v,a,b;}s[N];
LL f[M+10][mW+10];

bool cmp(node x,node y)
{   return x.b<y.b;   }
LL mmax(LL x,LL y)
{   return x>y?x:y;   }

int main()
{
    while (1)
    {
      scanf("%d%d",&n,&W);
      if (n==-1&&W==-1) break;
      for (int i=1;i<=n;i++)
      {
        scanf("%d%d",&s[i].w,&s[i].v);
        s[i].a=s[i].w,s[i].b=0;
        while (s[i].a%2==0) s[i].a/=2,s[i].b++;
      }
      sort(s+1,s+1+n,cmp);
      int l,r;
      l=1;
      memset(f,0,sizeof(f));
      for (int i=0;i<=M;i++)
      {
        r=l-1;//
        while (r<n && s[r+1].b==i) r++;//<n
        for (int j=l;j<=r;j++)
         for (int k=mW;k>=s[j].a;k--)
           f[i][k]=mmax(f[i][k],f[i][k-s[j].a]+s[j].v);
        for (int j=0;j<=mW;j++)//or (int j=mW;j>=0;j--)
          if (W&(1<<i)) f[i+1][j>>1]=mmax(f[i+1][j>>1],f[i][j]);
          else f[i+1][(j+1)>>1]=mmax(f[i+1][(j+1)>>1],f[i][j]);
        l=r+1;
      }
      printf("%lld
",f[31][0]);
    }
    return 0;
}