【训练题8：(更新)大体积01背包变式】不开心的金明 | 洛谷 P3985

不开心的金明 | 洛谷 P3985

题外话

$2020/11/29，在ZJNU个人赛遇到了一样的题但是过不了。$
$于是发现了这题的蹊跷。$
这题所有能过的贪心策略的题解可能都有问题！

难度

$普及+/提高$

题意

物品 $N$ 种，背包容积 $V$
每样物品有体积 $v_i$ ， 重要度 $p_i$。
求01背包之后的重要度和最大为多少。

特殊条件：$\underset{i}{\min }(v_i)\le v_i\le \underset{i}{\min }(v_i)+3$
并且他装不下所有的物品。

数据范围

$1\le N\le 100$
$1\le W\le 10^9$
$1\le v_i\le 10^9$
$1\le p_i\le 10^7$

思路

• 物品种类虽然少，但是背包体积太大了，普通的 $O(NV)$ 背包肯定没法做。
• 考虑特殊条件，我们有两种做法：

---

• 【做法1】
• 因为最多就物品体积有四种，每种体积我们肯定优先选择重要度最大的。
• 因此我们枚举每种体积我们选几种物品。
• 时间复杂度 $O(N^4)$，经过优化可以 $O(N^3)$。
• 但是如果物品体积有 $x$ 种，时间复杂度 $O(N^x)$ 怎么优化都不是很优秀的样子。

---

• 【做法2】
• 因为物品种类很少，且物品体积差别也很小。我们考虑什么时候无论怎么选都只能选固定数量的物品
• 易得，当 $\lfloor \frac{V}{\min v_i}\rfloor =\lfloor \frac{V}{\max v_i}\rfloor$时成立。并且有隐含条件：
• $\max v_i=\min v_i+3$
• $V<N\times \max v_i$
• 化简可得：
  
• 综上：
• 当 $\min v_i<300$时，普通的01背包即可，时间复杂度 $O(NV)=O(N\times N\times \max v_i)=O(N^3)$
• 当 $\min v_i\ge 300$时，因为只能选择固定的 $\frac{V}{\min b_i}$个物品，每次贪心选择重要度最高的即可，时间复杂度 $O(N\log N)$

但是$上述策略并不正确！$

• 因为对于某一个 $\min =\frac{V}{k}$，易得 $\min =\max +1$，$\min \ne \max \text{min=max}\min \ne \max$
• 因此，$贪心策略得逞于题目数据氺了$

正解DP

• 因为每个物品的体积很大，我们缩小体积为 $v_i=v_i-\min v$
• 设 $dp[i][j]$ 表示 选了i个物品，新的背包体积为j时的最大收益
• 易得，如果我们选择了 $i$ 个物品，那么目前的 原来背包的体积为 $i\times \min v+k$
• 合法转移，我们需要现在的背包体积合法、原来的背包体积也合法。
• 那么，我们 修改后的背包体积就类似于新的背包的属性了，就变成了二维背包
• 二维01背包的每一个属性都需要 逆序枚举即可。

核心代码

时间复杂度为 $O(x{N}^3)$

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const int MAX = 3e3+50;
const int INF = 0x3f3f3f3f;

int w[MAX];
int v[MAX];
int dp[MAX][MAX];
int main()
{
    int N;
    ll M;
    cin >> N >> M;
    int mn = INF;
    for(int i = 1;i <= N;++i){
        cin >> w[i] >> v[i];
        mn = min(mn,w[i]);
    }
    for(int i = 1;i <= N;++i)
        w[i] -= mn;

    int ans = 0;
    for(int i = 1;i <= N;++i)
    for(int j = min(N,(int)M / mn);j >= 1;--j)
    for(int k = 3 * j;k >= w[i];--k){
        if((ll)j * mn + k <= M)dp[j][k] = max(dp[j][k],dp[j-1][k-w[i]]+v[i]);
        ans = max(ans,dp[j][k]);
    }
    cout << ans;
    return 0;
}