数学建模社区-数学中国
标题:
python代码解决0-1背包问题
[打印本页]
作者:
2744557306
时间:
2023-11-7 11:17
标题:
python代码解决0-1背包问题
一、首先介绍一下0-1背包问题:
. h& L5 X; z; J: `
0-1背包问题(0-1 Knapsack Problem)是一个经典的组合优化问题,通常在计算机科学和运筹学中讨论。这个问题涉及到一个背包和一组物品,每个物品都有一个特定的重量和价值。问题的目标是在给定背包的最大容量下,选择一组物品放入背包,以使得所选物品的总重量不超过背包的容量,同时最大化这些物品的总价值。
9 u- I6 g9 X% A6 m5 k
0-1背包问题的名称中的“0-1”表示每个物品要么完全放入背包(选择)要么完全不放入背包(不选择),不能部分放入。这是问题的一个关键特征,与分数背包问题不同,分数背包问题允许部分放入物品。
+ u9 E* F) y' `" X
问题的形式描述如下:
7 ~1 X5 I% O9 `
1.给定一个固定容量的背包,通常表示为一个正整数W(背包的最大承载重量)。
: w- G$ h5 U3 u/ x. N# G8 `
2.给定一组物品,每个物品都有两个属性:重量(weight)和价值(value)。
. d" d/ N! _# @/ Z1 r' @7 f: J
3.对每个物品,你可以选择将其放入背包(选择)或不放入背包(不选择)。
$ J4 Y1 N" a: x* s6 Z! X1 D0 Z$ `
4.每个物品只能选择一次,即要么放入背包,要么不放入。
3 U$ t1 L9 \; T7 g1 _ \
5.目标是选择一个物品组合,使得它们的总重量不超过背包容量W,同时使它们的总价值最大化。
( ?! [' r2 M/ |/ P% M6 v! `
! x: w$ z9 K: B* E7 O* f, b
解决0-1背包问题的一种常见方法是使用动态规划(Dynamic Programming)算法。这个问题有广泛的应用,包括资源分配、排程问题、投资组合优化等领域。它还是计算复杂性理论中的一个经典问题,通常被用来说明NP难问题的概念。
/ z6 I; Y/ D9 E9 ` z4 B
9 |3 \" D- `+ o# Z$ C0 J7 V
, X/ e3 ^, n, K4 p, N) Q2 M
二、 介绍代码
; d7 v" t' i$ n5 W- P
这段代码是一个Python实现的0-1背包问题的解决方法,使用了动态规划算法来找到最优解。以下是对代码的详细解释:
def knapsack(v, w, n, capacity):
' v4 s$ q5 G0 u0 o
i = 0
. C4 m, u- r# y x0 [9 b
capacity = capacity + 1 # 初始化背包容量最大值
- l( ~5 {4 }2 o* {! Q
m = np.zeros((n, capacity)) # 初始化
; ^% C/ j. y1 S; f* A \
x = np.zeros(n)
复制代码
1.v 是一个长度为 n 的列表,表示每个物品的价值。
. y. I9 N$ v8 r7 P" V
2.w 是一个长度为 n 的列表,表示每个物品的重量。
+ B0 M# ~( a7 Y7 ]- V# U* I& F; Q. D
3.n 表示物品的数量。
n. L: |* ^- Q. a: W7 T
4.capacity 表示背包的容量。
6 e X" A+ v/ q- D: E; i
0 f% a# S* ?2 d- c1 h) M& ^
代码首先初始化了一个二维数组 m 作为动态规划表,其中 m
[j] 表示在考虑前 i 个物品时,背包容量为 j 时可以获得的最大总价值。数组 x 用来存储最终的解,x
表示是否选择第 i 个物品。
for i in range(n):
0 u/ P5 ]2 l3 k) z& K, X
for j in range(capacity):
) R R+ u% F' i6 R1 [! b
if (j >= w[i]):
6 H& }4 R6 _ N) `" ~; {2 @8 s& A' `
m[i][j] = max(m[i - 1][j], m[i - 1][j - w[i]] + v[i])
\+ q& u. p e2 m* ?. ~$ u
else:
. r3 `6 M3 n% x
m[i][j] = m[i - 1][j]
复制代码
在这个部分,代码使用了一个嵌套的循环,遍历了所有物品和不同的背包容量。对于每个物品 i 和容量 j,代码计算了两种情况下的最大总价值:
# W0 o x/ z2 {/ I( |$ V
5.如果当前物品的重量 w
小于等于当前容量 j,那么可以选择将第 i 个物品放入背包,此时总价值为 m[i-1][j-w
] + v
,或者选择不放入,此时总价值为 m[i-1][j]。代码选择其中较大的值作为 m
[j]。
. @; g. c" e% i: ^4 q4 C, N8 B
6.如果当前物品的重量 w
大于当前容量 j,则无法放入物品,所以总价值等于上一行的值 m[i-1][j]。
" h+ n6 c& w0 @6 L% ~' V
6 J0 S5 P+ p, A; @) h1 h, C: `( \
这个循环填充了动态规划表 m,最终 m[n-1][capacity] 包含了问题的最优解,即在给定容量下可以获得的最大总价值。
capacity = capacity - 1
' u1 @9 b, L' o7 K& Z$ ]9 f) |% }6 P
for i in range(n - 1, 0, -1):
' Z7 ?! b9 W6 ^! I ?
if (m[i][capacity] == m[i - 1][capacity]):
X5 a$ l2 n9 ~. `
x[i] = 0
2 g8 U s5 [. h7 W5 N
else:
& R/ ~ G- f& J1 `
x[i] = 1
' w9 q& E7 g! Q: j G# `0 }, i
capacity -= w[i]
. d5 L+ H# A0 q/ A
x[0] = 1 if (m[1][capacity] > 0) else 0
复制代码
在这一部分,代码反向遍历动态规划表,从最后一行向前找到解的路径。如果 m
[capacity] 等于 m[i-1][capacity],表示第 i 个物品没有放入背包,否则放入背包,并更新剩余容量 capacity。
weight = 0
! K9 @2 c6 u6 f- y$ j
value = 0
/ ?: M( g; w2 c2 B. [# j
print('装载的物品编号为:')
9 z- w$ n% K$ \" ]2 @7 [' Q6 I0 I
for i in range(len(x)):
B4 g* r) O0 L( L Y- Z7 l
if (x[i] == 1):
# b+ B0 U; b/ |/ m& V
weight = weight + w[i]
$ A4 A3 {. j! l* ?6 n# C
value = value + v[i]
7 _5 g* d, o( r5 b0 L3 e
print(' ', i + 1)
4 R- V, y# U' K! p
print('装载的物品重量为:')
* j3 {0 B5 D0 g' G' v
print(weight)
( ^2 W4 S8 R' R$ a) W% E
print('装入的物品价值为:')
. K; r" |7 e9 ]" N9 }/ p
print(value)
: ~4 |* H% t* x$ o0 T
return m
复制代码
最后,代码计算了被选择的物品的总重量和总价值,并将它们打印出来。函数返回动态规划表 m。
( a1 x% ]- u+ z# \5 w5 |
这段代码实现了0-1背包问题的解决方法,它通过动态规划算法找到最优解,即在给定背包容量下可以获得的最大总价值,以及选择哪些物品放入背包。
* w: n, {: Y0 P) e# Y6 j) U
2 R2 ~3 u2 |! A' G) ^: [$ J
最后在函数的输入文件中,如下例,第一行物品数量为:5 背包载重量为:10物品的重量列表为:[2, 2, 6, 5, 4] 物品的价值列表为: [6, 3, 5, 4, 6]
' f7 z1 I4 r" F! ?' |) X9 o- ]$ f- j; O
2023-11-7 11:13 上传
下载附件
(8.27 KB)
1 W* P+ a+ M8 X i: t" y# l. A3 ?
5 a9 s1 ?- r: f) Q
接下来展示我们的输出结果:
$ @: y o! ?# x$ B
2023-11-7 11:16 上传
下载附件
(126.79 KB)
$ R4 B, i5 q* X, c0 p2 d3 B/ [
/ S' k) t0 I9 Z: |* {. z) X
1 H- _( z# s" Y4 e; t$ p
具体代码如下:
% y; D: w9 s; k* v$ S
( {9 K2 ?( i; _" j8 g- U' W( r3 M0 e$ u
( D; F: P+ D) C6 b' n9 N
3 F$ c U$ c% L6 W1 u+ l
: K$ b( j8 A# t+ S W! f; K
5 L+ W9 C% z" ^" d4 }( Q
基础0-1背包问题(动态规划).rar
2023-11-7 11:17 上传
点击文件名下载附件
下载积分: 体力 -2 点
3.58 KB, 下载次数: 0, 下载积分: 体力 -2 点
售价:
2 点体力
[
记录
] [
购买
]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5