• bzoj1190 [HNOI2007]梦幻岛宝珠


    传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=1190

    【题解】

    首先,我们把所有物品都分解成$a imes 2^b$的形式,然后把物品按照$b$分组;

    我们按$b$从高到低考虑。$f(i,j)$表示考虑到$2^i$,当前还剩余$j imes 2^i$的空间,所能取到的最大值。

    每层先从上一层传递$f$数组,然后再更新。每次就是一个背包转移了。

    考虑这个$j$可能随着$b$减小越来越大,我们需要优化。

    对于每一层的物品,最多只能表示成$na imes 2^b$的形式,即最多$1000 imes 2^b$。

    所以对于上一层传下来的值,如果它乘了2之后大于1000,实际上这些多于1000的部分是没有用的,我们可以直接舍弃!

    这样的话就能保证复杂度了!

    # include <vector>
    # include <stdio.h>
    # include <string.h>
    # include <iostream>
    # include <algorithm>
    
    using namespace std;
    
    typedef long long ll;
    typedef unsigned long long ull;
    typedef long double ld;
    
    const int M = 5e5 + 10;
    const int mod = 1e9 + 7;
    
    int n, W; 
    
    vector<int> w[32], v[32];
    ll f[32][1010];
    int to[32];
    
    bool Main() {
      cin >> n >> W;
      if(n == -1 && W == -1) return 0;
      for (int j=0; j<=31; ++j) {
        w[j].clear(), v[j].clear(), to[j] = 0;
        for (int i=0; i<=1000; ++i) f[j][i] = 0;
      } 
      for (int j=1; j<=n; ++j) {
        int ww, vv;
        cin >> ww >> vv;
        for (int i=30; ~i; --i) {
          if(ww % (1 << i) == 0) {
            w[i].push_back(ww / (1 << i));
            v[i].push_back(vv);
            break;
          }
        }
      }
      
      
      for (int j=30, curW=0; ~j; --j) {
        bool flag = 0;
        to[j] = to[j+1] * 2;
        if(W & (1<<j)) flag = 1, to[j] ++;
        
        if(to[j] <= 1000) {
          for (int i=0; i<=to[j]; ++i) 
            if(flag == 0) {
              if(i%2 == 0) f[j][i] = f[j+1][i/2];
            } else if(i%2 == 1) f[j][i] = f[j+1][i/2];
        } else {
          for (int i=0; i<=1000; ++i) 
            if(flag == 0) {
              if(i%2 == 0) f[j][i] = f[j+1][i/2];
            } else if(i%2 == 1) f[j][i] = f[j+1][i/2];
          for (int i=500; i<=to[j+1]; ++i)
            f[j][1000] = max(f[j][1000], f[j+1][i]); 
          to[j] = 1000;
        }
        
        
        for (int k=0, kto = w[j].size(); k<kto; ++k) {
          for (int i=0; i<=to[j]; ++i) {
            int tw = w[j][k], tv = v[j][k];
            if(i - tw >= 0) f[j][i-tw] = max(f[j][i-tw], f[j][i]+tv);
          }
        }
        
      }
      
      ll ans = 0;
      for (int i=to[0]; ~i; --i) ans = max(ans, f[0][i]);
      cout << ans << endl; 
      return 1;
    }
    
    int main() {
      while(Main());
      return 0;
    }
    View Code

     

  • 相关阅读:
    转】Apache解决高并发和高可用
    Kafka学习(一)配置及简单命令使用
    unity3d教程动态创建简单平面地形
    LeetCode: Unique Binary Search Trees [095]
    德惠也有星巴克
    一个css和js结合的下拉菜单,支持主流浏览器
    【图像处理】人类视觉成像原理
    windows使用nginx+memcached实现负载均衡和session或者缓存共享
    OpenCV基础篇之画图及RNG随机数对象
    在阿里云上布置git server
  • 原文地址:https://www.cnblogs.com/galaxies/p/bzoj1190.html
Copyright © 2020-2023  润新知