[HNOI2007]梦幻岛宝珠

https://www.zybuluo.com/ysner/note/1332315

题面

给你(n)颗宝石，每颗宝石都有重量和价值。要你从这些宝石中选取一些宝石，保证总重量不超过(m)，且总价值最大。

• (nleq100,mleq2^{30},val_i=a*2^b,aleq10,bleq30)

解析

既然都告诉了重量是(2)的幂次方带系数，那我们就可以用系数来代表一种重量。
就可以设状态为(f[i][j])表示只取重量形如(a2^i)的物品，系数和小于等于(j)的最大价值
（因为系数转移起点可以大于(0)）。
然后可以在各个(i)内跑背包。

接下来考虑合并(i)不同的状态。
注意到一个问题，(i)从小往大转移时，(j)是会缩水的，因为(j/2)会强制向下取整。

假定(f[i-1])系列对(f[i])系列的系数贡献为(k)。
此时，如果二进制中(m)的第(i-1)位为(1)，说明对应的(f[i-1])系列中的(j)可以是(k*2+1)。否则就只能是(k*2)。
这两个值(/2)后都是(k)，但是(f[i-1][k*2+1]geq f[i-1][k*2])。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;++i)
#define fq(i,a,b) for(re int i=a;i>=b;--i)
using namespace std;
const int N=5e5+100;
int n,m,f[35][1010],V[35][1010],mx,t[35],W[35][1010];
il ll gi()
{
  re ll x=0,t=1;
  re char ch=getchar();
  while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
  if(ch=='-') t=-1,ch=getchar();
  while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  return x*t;
}
il int max(re int x,re int y){return x>y?x:y;}
il void init()
{
  mx=0;
  memset(f,0,sizeof(f));memset(V,0,sizeof(V));
  memset(W,0,sizeof(W));memset(t,0,sizeof(t));
}
int main()
{
  while(233)
    {
      n=gi();m=gi();
      if(n==-1&&m==-1) break;
      init();
      fp(i,1,n)
    {
      re int x=gi(),y=gi(),j=0;
      while(!(x&1)) {x>>=1;++j;}
      mx=max(mx,j);
      t[j]+=x;
      V[j][++V[j][0]]=x;//volume
      W[j][++W[j][0]]=y;//val
    }
      fp(i,0,mx)
    fp(j,1,V[i][0])
    fq(k,t[i],V[i][j])
    f[i][k]=max(f[i][k],f[i][k-V[i][j]]+W[i][j]);
      while(m>>mx) ++mx;--mx;
      fp(i,1,mx)
      {
    t[i]+=(t[i-1]+1)/2;
    fq(j,t[i],0)
    fp(k,0,j)
      f[i][j]=max(f[i][j],f[i][j-k]+f[i-1][min(t[i-1],(k<<1)|((m>>i-1)&1))]);
      }
      printf("%d
",f[mx][1]);
    }
  return 0;
}