背包问题学习笔记-多重背包问题

news/2024/10/23 5:31:04/

题意描述:

有 N 种物品和一个容量是 V 的背包。第 i 种物品最多有 si 件,每件体积是 vi,价值是 wi。求解将哪些物品装入背包,可使物品体积总和不超过背包容量,且价值总和最大。
输出最大价值。输入格式
第一行两个整数,N,V (0<N≤1000, 0<V≤20000),用空格隔开,分别表示物品种数和背包容积。接下来有 N 行,每行三个整数 vi,wi,si,用空格隔开,分别表示第 i 种物品的体积、价值和数量。输出格式
输出一个整数,表示最大价值。

多重背包的问题根据数据范围的大小划分为了三个难度层次,分别对应三种解法。对应的数据范围分别是:

暴力求解
0<N,V100
0<vi,wi,si≤100
二进制优化
0<N1000
0<V2000
0<vi,wi,si≤2000
单调队列优化
0<N1000
0<V20000
0<vi,wi,si≤20000

示例:

4 5
1 2 3
2 4 1
3 4 3
4 5 2

10


解题思路:
Alice: 这是一道题拆成了三道题 ?这题这么难
Bob: 丰俭由人,😁
Alice: 多重背包和 01 背包的区别就是每类物品有了数量的限制,01 背包是只能选或者不选,完全背包是由无数种可以选,01 背包是有 k 种可以选。
Bob: 是的,那直接改一下状态转移公式得了,直接按照完全背包的方式换一下最内层的状态转移好了,dp[j] = max(dp[j], dp[j-vi * k] + wi * k) k 取值 0,1,2,3 …
Alice: 为啥要按照完全背包,从最大体积开始求解 ?
Bob: 看状态转移方程,我们还用一维度的 dp 数组的话,对于第 i 个物品,还需要用到第 i-1 个物品的 j-vi * k 的状态,从最大体积往小了计算才是对的。
Alice: 这个理解起来还不太难。
Bob: 是的,这里可以关注一下计算量,暴力去算实际上是 O^3 的,数据范围都是 100,最大也就是 10^6,这样不会超时的。
Alice: 二进制优化讲的是啥 ?
Bob: 二进制优化的前提是,把多重背包转换成 01 背包问题,但是转换的方式有很多种,二进制优化是其中一种的优化方法。举个例子,第 0 种物品有 7 个,我们该如何拆分呢,拆成 7 个1 ?每个体积和价值都是 v0,w0。这样的计算量大概,O^2 也就是双重循环,1000 * 2000 * 2000 大概 2* 10^9, 1s 的时间限制大概能完成10^6~10^7,这样会超时的。
Alice: 是的,直接拆,拆出来的物品数量太多了
Bob: 这样其实就转换成了拆分的问题,把 7 拆成若干个数字,且这些数字能组成 0-7 的所有数字。
Alice: 哦哦,就是二进制啊,二进制是能够表达所有整数的,0到7也就是,2^1 , 2^1, 2^2 这些。然后其实就是二进制的表示一样,7 就是 111,对应到三个物品,就是全选 ?
Bob: 对,但是还有一个问题,如果你的数字是 8 呢,拆成 1,2, 4,8 ?那样的话,就有可能选出 1+2 + 4 + 8 == 15 的选择,但是原来最多也就是 8个。
Alice: 可以这样子,8 的数字就用 7 的那套,剩下的差值额外补齐,0-7 的范围再加个 1,所能表达的范围就是 0-8。
Bob: 应该就是这样,然后具体拆分的时候还可以直接从 1,2,4 累乘,然后不断地减下来,最后剩下的额外加一套。这样应该比较方便而且很快。
Alice:nice,拆外之后直接用 01 背包就行了。这样算的计算量应该是 log 2000 * 1000 * 2000 == 2.2 * 10^7 ,应该没问题啦。
Bob:单调队列优化呢 ?
Alice: 这个看起来很麻烦的样子,是啊,这个需要先知道单调队列是啥。可以先看一下 这篇文章
Bob:单调队列优化的思路大概是这样的,还得从完全背包的递推公式讲起,考虑第 i 种物品 dp[j] = max(dp[j], dp[j-vi * k] + wi * k) k 取值 0,1,2,3 …。这里实际的计算过程是什么样的呢,j-vi*k 减到最后,无论 j 的取值是啥,一定剩下的是 vi 的各种余数,1,2,3 … vi-1,然后从 1+v,1+2v,1+3v … 的状态转移计算过程会相互影响,而 2+v,2+2v,2+3v 会相互影响,两个余数之间的计算过程互相不影响,这样我们就能把对 i 物品的计算划分为互相独立的 vi-1 个,然后单独计算。
Alice: 然后的,拆分成 多个计算过程就能变快 ?并行处理吗 ?
Bob: 然后就是难以理解的地方了,我先举个例子吧,假设背包体积是 100, 第 i 种物品的体积是 5,价值是 4,数量是 3,考虑第 i 个物品时候的状态的计算,dp[j] = max(dp[j- 5*0] + 4*0, dp[j- 5*1] + 4*1, dp[j-5*2] + 4*2), dp[j - 5*3] + 4*3,具体的计算过程,从底向上就是 求 dp[0 + 5*0], dp[0 + 5*1] + 4*1, dp[0 + 5*2] + 4*2, dp[0 + 5*3] + 4*3 之间的最大值,然后再根据放 1 个,2个,3个求出 dp[0+20]dp[1+20] 看的是 dp[1 + 5*0] + 4*0, dp[1 + 5*1] + 4*1, dp[1 + 5*2] + 4*2, dp[1 + 5*3] + 4*3 的最大值。明白了吗 ? 0,1,2 在这里就是不同的计算序列,每个计算序列都要根据 4 的滑动窗口求前面的最大值,然后再计算当前的最大值。
Alice: 滑动窗口就在这呢,原来是对第 i 种物品的体积余数的每个计算序列里面滑动,滑动窗口的大小就是物品的最大数量。
Bob: 其实滑动窗口的大小不一定是物品的最大数量,k 的实际取值范围是 math.min (si + 1, maxVolum / vi),不过这个可以在代码层直接实现掉,可以暂时认为是最大数量,不影响理解。
Alice: 然后单调队列里面维护的时候什么呢 ?队首和队尾都是怎么维护的呢 ?
Bob: 单调队列里面维护的是当前窗口里面最大价值所对应的体积,这样我们应该能够比较轻松的写出 dp[j] 的更新,dp[j] = lastRowDp[queue[0]] + (j - queue[0]) / vi * wi,想一下,queue[0] 里面是当前窗口的最大价值对应的体积,我们在这个体积之上更新 dp[j]。队首的维护其实还是窗口的大小,只不过这里窗口的大小是通过体积的计算来校验的。
Alice: 这些都还好理解,那队尾的维护呢 ?我记得滑动窗口最大值里面是直接计算窗口里面的数字之和,这里应该不是吧。
Bob: 确实不是,这里还有点不太好理解。这里还是按照最大价值来计算的,只是比较的是第 i-1 个物品对应的体积和 第 i 个物品对应的体积所能给 dp[j] 带来的价值收益。要知道我们在 queue[0] 队首的位置维护的是在窗口中能给 dp[j] 带来最大价值的体积,而单调队列的维护正式通过队首和队尾维护的,所以队尾的维护逻辑实际和 dp[j] 的更新逻辑是一致的。
Alice: 更新逻辑是一致的 ?!我好像有点明白了。还有一些细节问题,滑动窗口的大小不定是怎么实现的 ?
Bob: 这个好说,你从余数 r 开始,r 就是 r + 0 * vi 然后每次给 r += vi,让 r 不要超过最大体积就可以了。
Alice: 这题真难。
Bob:这题真难,我看别人的题解看了半天才明白滑动窗口在哪滑呢,看别人代码看了半天,单调队列维护的代码都快背下来了,也没看明白怎么维护的,还是得实际举个例子。
Alice: 还有一个小问题,这里为啥不能把状态压缩成一个一位数组,为啥还要一个 lastRowDp 呢 ?
Bob:简单,第 i 个物品的 dp[j] 的更新需要依赖于体积 j 的前 si 个状态,如果直接用 dp[j - vi],那用的就是更新过的值了,就不是 i-1 个物品的状态了。


代码:

暴力

const solve = (count, maxVolum, volumAndWeight) => {const dp = new Array(maxVolum + 1).fill(0);for(let i=0; i<count; ++i){// 第 i 个物品的体积和价值const [ivolum, iweight, itotal] = volumAndWeight[i];for(let j=maxVolum; j>=ivolum; --j) {const candiantes = [];for(let k=0; k<=itotal; ++k) {j - ivolum * k >= 0 && candiantes.push(dp[j - ivolum * k] + k * iweight);}dp[j] = Math.max(...candiantes);}}console.log(dp[maxVolum]);
}

二进制优化

const solve = (count, maxVolum, volumAndWeightAndCount) => {// 二进制拆分为 01 背包const volumnAndWeight = [];volumAndWeightAndCount.forEach(item => {let [v, w, s] = item;for (let k=1; k <= s; k*=2) {s -= k;volumnAndWeight.push([k*v, k*w]);}if(s > 0) {volumnAndWeight.push([s*v, s*w]);}});// 01 背包解法const dp = new Array(maxVolum + 1).fill(0);for(let i=0; i<volumnAndWeight.length; ++i){// 第 i 个物品的体积和价值const [ivolum, iweight] = volumnAndWeight[i];for(let j=maxVolum; j>=ivolum; --j) {dp[j] = Math.max(dp[j], dp[j - ivolum] + iweight);}}console.log(dp[maxVolum]);
}

单调队列优化

const fs = require('fs');
let buffer = '';process.stdin.on('readable', () => {const chunk = process.stdin.read();if (chunk) {buffer += chunk.toString()}
});// 输入的字符串转换为数字
const convert = (inputString) => {const list = [];inputString.split('\n').forEach((line) => {const tokens = line.split(' ');list.push(tokens.map(num => parseInt(num, 10)));});return list;
}// 批量调用
const batchCall = (list, solve) => {// 划分数据const data = [];let countAndVolumIndex = 0;while(countAndVolumIndex < list.length) {const [count, volum] = list[countAndVolumIndex];data.push({volum: volum,count: count,volumAndWeightAndCount: list.slice(countAndVolumIndex + 1, countAndVolumIndex + 1 + count)});countAndVolumIndex += count + 1;}data.forEach(item => {if(solve && item && item.count && item.volum) {solve(item.count, item.volum, item.volumAndWeightAndCount);}});
}const solve = (count, maxVolum, volumAndWeightAndCount) => {// 单调队列优化方法const dp = new Array(maxVolum + 10).fill(0);// 对于每种物品 for (let i=0; i<count; ++i) {// 状态压缩const lastRowDp = [...dp];// 取出第 i 种物品的体积,价值,数量const [vi, wi, si] = volumAndWeightAndCount[i];// 对于每种可能剩余的体积,0,1,2, ... vi-1 for (let r=0; r<vi; ++r) {// 单调队列求解每种可能的最大值,滑动窗口大小是,math.min (si, maxVolum / vi) 下取整// 0 + 0v, 0 + 1v, 0 + 2v ... 0 + kv 的数组中滑动,每次一步// 最大价值对应的体积的单调队列,双端队列const queue = [];for(let j=r; j<=maxVolum; j+=vi) {// 维护队首// i 物品的体积超了,注意这里是大于而不是大于等于,要把 r+0*vi也包括进来while(queue.length && j-queue[0] > vi*si) {queue.shift();}// 维护队尾while(queue.length && lastRowDp[queue[queue.length-1]] + (j - queue[queue.length-1]) /vi * wi <= lastRowDp[j]) {queue.pop();}// 入队queue.push(j);// 更新 dpdp[j] = lastRowDp[queue[0]] + (j-queue[0]) / vi * wi;}}}console.log(dp[maxVolum]);
}process.stdin.on('end', function() {batchCall(convert(buffer), solve)
});

参考:

  • 题目链接-多重背包1
  • 题目链接-多重背包2
  • 题目链接-多重背包3
  • 参考题解

http://www.ppmy.cn/news/1138568.html

相关文章

自媒体工作内容管理助手

内容助手 访问地址&#xff1a;editor.yunwow.cn 背景介绍 最近在学习流量运营&#xff0c; 流量运营的第一站是内容创作&#xff0c; 我试过不少原创内容&#xff0c;都是跟生活相关的例如&#xff1a;录一段联琴的视频、录一段秋天的风景、写一段生活感悟、发一段小宠物的生…

数据结构与算法-(7)---栈的应用-(4)后缀表达式求值

&#x1f308;write in front&#x1f308; &#x1f9f8;大家好&#xff0c;我是Aileen&#x1f9f8;.希望你看完之后&#xff0c;能对你有所帮助&#xff0c;不足请指正&#xff01;共同学习交流. &#x1f194;本文由Aileen_0v0&#x1f9f8; 原创 CSDN首发&#x1f412; 如…

tkinter中如何隐藏/去掉最大化窗口

在Tkinter中&#xff0c;去掉最大化窗口的按钮是无法实现的&#xff0c;因为这是Tkinter内置的特性。 #我的Python教程 #官方微信公众号&#xff1a;wdPython然而&#xff0c;你可以通过禁用窗口的最大化和最小化按钮来隐藏它们。 import tkinter as tk from tkinter impor…

Python3操作Redis最新版|CRUD基本操作(保姆级)

Python3中类的高级语法及实战 Python3(基础|高级)语法实战(|多线程|多进程|线程池|进程池技术)|多线程安全问题解决方案 Python3数据科学包系列(一):数据分析实战 Python3数据科学包系列(二):数据分析实战 Python3数据科学包系列(三):数据分析实战 Win11查看安装的Python路…

acwing算法基础之基础算法--高精度加法算法

目录 1 知识点2 模板 1 知识点 大整数 大整数&#xff0c;它们的长度都为 1 0 6 10^6 106。大整数是指长度为 1 0 6 10^6 106的整数。 大整数 - 大整数 大整数 * 小整数 大整数 / 小整数 把大整数存储到向量中&#xff0c;需要考虑高位在前还是低位在前&#xff0c;低位在前…

简单两步实现离线部署ChatGPT,ChatGPT平替版,无需GPU离线搭建ChatGPT

文末附主程序安装包和大模型参数文件~ 演示效果如下图所示&#xff1a; 一、使用方法 软件主要分为两个部分&#xff1a;GPT4ALL软件主体&#xff08;主程序&#xff09;模型参数&#xff08;离线模型&#xff09;&#xff0c;如果使用API Key的话则不需要下载模型参数。 可以…

高效使用Python的subprocess模块:基本用法与进阶技巧

文章目录 Python的subprocess模块参数和方法的用法subprocess.run、subprocess.check_call 和 subprocess.Popen的区别如何使用`subprocess`实现一些功能执行系统命令并获取输出通过管道执行多个命令在后台运行进程一些额外的`subprocess`模块技巧超时功能错误处理输入/输出重定…

Go语言中的指针介绍

Go语言中的指针 文章目录 Go语言中的指针一、Go语言中的指针介绍1.1 指针介绍1.2 基本语法1.3 声明和初始化1.4 Go 指针的3个重要概念1.4.1 指针地址&#xff08;Pointer Address&#xff09;1.4.2 指针类型&#xff08;Pointer Type&#xff09;1.4.3 指针取值&#xff08;Poi…