[贪心+01背包] 能量石(01背包+贪心)

0. 前言

相关：

• [背包] 背包问题算法模板(模板)
• 01 背包

调整法证明贪心问题：

• [贪心] 耍杂技的牛(公式推导+贪心)
• 相关题目：114. 国王游戏

1. 贪心+01背包

734. 能量石

参考题解：

• mrk：重点参考题解

• 辰风：能量石 十分清晰！！！！

重点： 贪心证明、01 背包、思维、Google Kickstart2019 Round B Problem B

贪心证明：

• 暴力来做的话，就得枚举所有能量石的全排列，每次计算当前这种排列的最优值即可。但是全排列是 $n!$ 的
• 假设最优解能量石有 $k(1\le k\le n)$ 个，假设序列为 $a_1,a_2,...,a_k$，在此是将贡献为 0 的能量石剔除，保证每个 $a_i$ 均大于 0
• 那么交换相邻两个宝石一定不会使贡献变大，假设之前的总时间为 $t$，则有 $$E_i+E_{i+1}-t\times L_i-(t+S_i)\times L_{i+1}\ge E_i+E_{i+1}-t\times L_{i+1}-(t+S_{i+1})\times L_i$$ 化简：$$S_i\times L_{i+1}\le S_{i+1}\times L_i$$等价于：$${\frac{S_i}{L}}_i\le \frac{S_{i+1}}{L_{i+1}}$$
• 这个就说明，最优解一定满足这个性质，以该条件对宝石进行排序。所以我们就得到了一个包含最优解的序列，将全排列的 $O(n!)$ 缩小到了 $O(nlogn)$，史诗级优化！
• 接下来，就是在这个序列中，采用 01 背包，判断该石头是否拿取的问题了
  • 状态表示：
    • f[i][j] 从前 i 个物品中选择且总时间恰好为 j 的所有方案的最大价值。因为在此 j 的取值和最大价值 f[i][j] 精确对应。 具体原因参考文章末尾 mrk juju 的分析!!!
• 状态计算：
  • 状态转移方程为：f[i][j]=max(f[i-1][j], f[i-1][j-si]+ei-(j-si)*Li)。ei 是原能量，j 是当前吃完这个能量石的时刻，那么吃花费了 s 时刻，所以吃之前应该是处于 j-s 这个时刻，所以从起始时刻到 j-s 时刻，总共减少了 (j-s)*Li 能量

当然，空间优化为 1 维也是完全 ok 的。

引自 辰风：能量石 十分清晰！！！！

利用贪心来对所有的能量石排序，这样问题就变成了按照一个线性顺序取对每一个能量石判断是否拿取的问题了，就变成了经典的01背包问题

还有一个是总时间 m的取值问题，因为我们状态定义的是前i个能量石中取，经过了恰好 j 时刻的解，所以只要找出 j 最大可能是多少就可以了，就是所有能量石时间总和

令f[i][j]为从前i个能量石中取，经过了恰好 j 时刻的解
状态转移方程 f[i][j] = max(f[i][j], f[i - 1][j - s] + max(0, e - l * (j - s)));

解释一下里面的细节：首先吃第i个能量石的时刻已经过了 j-s 的时间了，因为能力石的能力不为负数，所以和0取一个max

由于我们不知道刚好进过哪一个时刻是解最大的时候 ，所以对每一个时间遍历，找出最大值


引自麻烦写注释大佬！！！首先，此题想到 01 背包应该是不难的，因为这题无非就是对每一块石头进行选与不选的抉择。但光光进行这一步，答案显然是错误的，因为选到后面，可能会导致他们现有的能量变成 0。前面所选的石头如果需要花费越长的时间吃，则后面石头所剩余的能量就越少，我们如果不进行贪心，则会造成前面的决策真的只是成了一个局部最优解，并且利用这个局部最优解，你无法构造出全局最优解。

为此我们需要利用贪心来缩小我们的决策范围。

这里要弄清楚 01 背包是解决了个什么问题：它是枚举了所有选法，并在其中找到了最大值。所以对于排序之后第一个能量石，我们会把所有吃它的方案，和不吃它的方案全部枚举到。

理解了hh

代码：

#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

const int N = 105, M = 25005;

int n;
int a[N];
int f[M];

int main() {
    int T;
    cin >> T;
    
    while (T --) {
        cin >> n;
        for (int i = 0; i < n; ++i) cin >> a[i];
        
        sort(a, a + n);
        
        memset(f, 0, sizeof f);
        
        f[0] = 1;
        
        int res = 0;
        for (int i = 0; i < n; ++i) {
            if (!f[a[i]]) res ++;                       // f数组记录了所有能转移到的数的方案数
            for (int j = a[i]; j <= a[n - 1]; ++j)      // 在此背包容量为最大的那个数
                f[j] += f[j - a[i]];
        }
        cout << res << endl;
    }
    return 0;
}

2020/11/24
针对体积，选恰好这块，还是不能太理解。可能一开始做背包问题都是最大不超过 j。自己也没举一反三的能力，在做买书、潜水员问题时初次接触到体积恰好、体积至少、这样的定义，也做了详细总结。但还是感到很别扭…在背包问题求方案数时，体积也为恰好，本题也是。在我的理解下，状态定义不同，初始化方式不同，最后返回结果需要遍历求最大值等等…目前能够理解本题，但还是抠不出来里面的更多细节。二刷回顾的时候，后续这类状态定义方式会很注意！！！加入 疑问。