[BZOJ 1190][HNOI2007]梦幻岛宝珠

1190: [HNOI2007]梦幻岛宝珠

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1057  Solved: 611
[Submit][Status][Discuss]

Description

给你N颗宝石，每颗宝石都有重量和价值。要你从这些宝石中选取一些宝石，保证总重量不超过W，且总价值最大为，并输出最大的总价值。数据范围：N<=100;W<=2^30,并且保证每颗宝石的重量符合a*2^b（a<=10;b<=30）

Input

输入文件中包含多组数据。每组数据的格式如下：第一行是两个正整数n和W，1≤n≤100,1≤W≤2^30，分别表示宝石的数目和最多能带走的宝石重量。接下来的n行，每行有两个正整数weighti和valuei，1≤weighti≤2^30, 0≤valuei≤2^30，分别表示第i颗宝石的重量和价值，且保证weighti能写成a*2^b(1≤a≤10,0≤b≤30)的形式。同一行的两个正整数之间用空格隔开。最后一组数据的后面有两个-1，表示文件的结束。这两个-1并不代表一组数据，你不需对这组数据输出结果。并且输入文件中数据的组数不超过20。

Output

对于输入的每组数据，输出一个整数C，表示小P最多能带走的宝石的总价值。每个结果整数C单独占一行，且保证C不会超过2^30。

Sample Input

4 10
8 9
5 8
4 6
2 5
4 13
8 9
5 8
4 6
2 5
16 75594681
393216 5533
2 77
32768 467
29360128 407840
112 68
24576 372
768 60
33554432 466099
16384 318
33554432 466090
2048 111
24576 350
9216 216
12582912 174768
16384 295
1024 76
-1 -1

Sample Output

14
19
1050650

HINT

 

Source

题解

这波0/1背包的容量可有点大...直接0/1背包打上去肯定是 $MLE+TLE$ ...

注意到数据保证物品质量都是2的整次幂的较小倍数, 我们可以把所有物品按照最接近的2的整次幂分组, 分别做0/1背包之后, 由于各组物品质量数量级相差较大, 最后直接把同一组的结果当成一个然后乱搞就好了...

然而交到BZOJ上直接 $AC$ , HZOJ上 $TLE$ $64\%$ ...这可真蠢...

参考代码

GitHub

 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <iostream>
 #include <algorithm>

 const int MAXN=1e3+;

 int n,w;
 long long ans=;
 long long dp[][MAXN];

 void FastRead(int&);

 int main(){
     FastRead(n);
     FastRead(w);
     while(n>=&&w>=){
         memset(dp,,sizeof(dp));
         ans=;
         for(int i=;i<=n;i++){
             int val,cnt=,wgh;
             FastRead(wgh);
             FastRead(val);
             while((wgh&)==){
                 wgh>>=;
                 cnt++;
             }
             for(int j=;j>=wgh;j--){
                 dp[cnt][j]=std::max(dp[cnt][j],dp[cnt][j-wgh]+val);
             }
         }
         for(int i=;i<=;i++){
             for(int j=;j<=;j++){
                 dp[i][j]=std::max(dp[i][j],dp[i][j-]);
             }
         }
         for(int i=;i<=std::min(,w);i++){
             ans=std::max(ans,dp[][i]);
         }
         for(int i=;i<=&&(<<i)<=w;i++){
             for(int j=std::min(,w>>i);j>=;j--){
                 for(int k=;k<=j;k++)
                     dp[i][j]=std::max(dp[i][j],dp[i][j-k]+dp[i-][std::min(k*+((w>>(i-)&)),)]);
                 ans=std::max(ans,dp[i][j]);
             }
         }
         printf("%lld\n",ans);
         FastRead(n);
         FastRead(w);
     }
     return ;
 }

 void FastRead(int& target){
     target=;
     bool inv=false;
     register char ch=getchar();
     while(!isdigit(ch)){
         if(ch=='-')
             inv=true;
         ch=getchar();
     }
     while(isdigit(ch)){
         target=target*+ch-'';
         ch=getchar();
     }
     if(inv)
         target=-target;
 }

Backup