贪心算法解决部分背包问题:策略对比与代码实现 贪心算法解决背包问题是个经典实验我在带学生做算法课程设计时发现很多同学把代码跑通就算完事却说不清“为什么这么选”“什么时候贪心能得到最优解”“什么时候不能”。这篇内容就围绕这个实验把问题建模、策略选择、代码实现和踩坑记录完整梳理一遍既能给正要交实验报告的同学做个参考也能让想加深理解的开发者看到背后的思路取舍。先说清楚一个容易被忽略的大前提背包问题通常分为0-1背包和部分背包分数背包。标题里的“贪心算法解决背包问题”实验场景里一般指部分背包——就是物品可以按比例拆分装入这种情况下贪心能保证得到最优解。如果是0-1背包贪心只能得到近似解不保证最优。这篇博文主讲的场景是“可拆分”的部分背包后面还会专门讲0-1背包下贪心为什么失效。1. 实验背景与问题建模1.1 背包问题的数学定义与两种典型形态背包问题的本质是资源分配有一个容量固定的背包有一堆物品每个物品有自己的重量和价值目标是让“装入背包的总价值”最大化。数学上可以写成给定容量 Cn 个物品每个物品重量 w_i、价值 v_i。选择若干物品放入背包使总重量 ≤ C总价值最大。两种形态的最大差别在于“物品能否分割”类型物品可否拆分贪心是否最优标准解法0-1背包不可拆分要么整体装入要么不装不保证动态规划部分背包分数背包可按重量比例拆分保证最优贪心算法做实验时如果你用的是“单位价值降序排列后依次装入”的贪心策略就只能处理部分背包。如果实验要求的是0-1背包必须改用动态规划否则结果很可能就是错的。我见过不少同学在实验报告里把这两种混为一谈这是首先要避开的坑。1.2 贪心算法核心思想局部最优能否推出全局最优贪心算法的核心是“每一步都做当前看起来最好的选择并且不回头看”。装背包的场景里“当前看起来最好”就是“优先装单位价值最高的物品”。这个策略符合贪心选择的经典逻辑——如果某个物品单位价值最高那么最优解里一定可以包含尽可能多的该物品。但要理解为什么局部最优能推出全局最优需要看一个关键的“交换论证”思路假设某个最优解中没有优先选择单位价值最高的物品那么可以用单位价值更高的物品替换其中一部分重量总价值只会增加或不变。这意味着存在一个最优解一定包含单位价值最高的物品尽可能多的一部分。逐步替换下去就能证明贪心策略得到的解就是全局最优。生活类比相当于去超市采购预算有限同类商品里自然先挑“每块钱能买到最多实用价值”的那个先把性价比最高的装进购物车再挑次高的直到钱花完为止。2. 三种常用贪心策略与适用场景2.1 价值优先策略价值优先就是把物品按价值从高到低排序依次装入背包。这个策略直观能快速得到一个“看着不错”的解但问题很明显价值高的物品往往也很重很容易因为装了它而占掉大量空间导致后面明明有很多轻巧且单位价值不错的小物品却没空间放了。我做过一个模拟测试背包容量50物品有“重量48价值100”的大件还有一堆“重量2价值10”的小件。贪心按价值优先会先装大件剩余容量2只能再装一个2重的小件总价值110。但如果放弃大件装满20个小件总价值是200。价值优先策略比最优解差了将近一半所以它只适合“重量相差不大”的特殊数据形态比如所有物品重量都差不多时。这个策略的实验价值在于它是一个对比基准能用来反衬其他策略的优势。2.2 重量优先策略重量优先就是把物品按重量从小到大排序优先装轻的物品。它的逻辑是“尽量多装件数”在背包容量足够且物品重量都很小时结果不错。但极端情况下同样会翻车如果有一个单件物品价值极高但重量居中重量优先策略会先把一堆低价值的轻物品塞进去等轮到那个高价值物品时空间已经不够了。我实测过一组数据容量30一个“重量6价值30”的物品五个“重量1价值3”的小物品还有若干中等的。重量优先策略会先装完五个小物品再装中等物品最后可能不够装高价值物品。最后总价值比最优解低不少。这个策略适合的场景是“各物品价值差异不大、主要追求装入数量”比如物流装车时每个包裹价值相近优先装小件能提高运力。2.3 单位价值优先策略推荐方案单位价值优先即按 v_i / w_i 从高到低排序优先装“每单位重量价值最高的物品”。这是部分背包的最优贪心策略。逻辑链条是排序按单位价值降序排。从高到低依次尝试放入。遇到装不下的整件时如果是部分背包可以按剩余容量拆分只装一部分。拆分这一步是部分背包和0-1背包的分水岭。能拆分的时候贪心完备且可证明最优不能拆分的时候这里就出现了“装不下整件但拆不了”的尴尬贪心就可能不是最优解。实验报告中单位价值优先策略应该作为“主算法”来写另外两种作为对比。这样报告显得有层次也能体现对问题的深入理解。3. 实验环境准备与代码实现3.1 实验环境与数据构造我用的是Python 3环境不依赖第三方库纯标准库就能完成这样整个实验方便复现。数据构造时我建议用随机生成而不是手工写死少量数据原因有两个一是随机数据能反映算法在不同分布下的表现二是数据量大一点才能看出排序策略的差异。生成数据的核心参数物品数量 n我常用 n50 到 n200 之间。数据太少看不出规律太多不好手工验证最优解。背包容量 C取总重量的40%到60%左右这样既不会因为太大导致所有物品都能装下也不会因为太小导致实验失去意义。重量范围1~50之间均匀随机。价值范围1~100之间均匀随机。我习惯把数据生成写成一个独立函数每次实验换一组种子方便对比不同随机种子下的结果稳定性。3.2 贪心核心代码与调用流程先定义物品结构我用Python的dict或者简单类都可以这里用namedtuple更清晰from collections import namedtuple import random Item namedtuple(Item, [id, weight, value]) def generate_items(n50, max_weight50, max_value100): return [Item(i, random.randint(1, max_weight), random.randint(1, max_value)) for i in range(n)]接下来是三种排序策略的核心函数。注意我在这里统一返回“装入比例列表”方便后续计算实际装入的总重量和价值def greed_solve(items, capacity, strategydensity): strategy: value 按价值降序 weight 按重量升序 density 按单位价值降序 返回: 装入比例列表、总价值、总重量 if strategy value: sorted_items sorted(items, keylambda x: x.value, reverseTrue) elif strategy weight: sorted_items sorted(items, keylambda x: x.weight) elif strategy density: sorted_items sorted(items, keylambda x: x.value / x.weight, reverseTrue) else: raise ValueError(unknown strategy) ratio [0.0] * len(items) total_value 0.0 total_weight 0.0 remain capacity for it in sorted_items: if remain 0: break if it.weight remain: # 整件装入 ratio[it.id] 1.0 remain - it.weight total_value it.value total_weight it.weight else: # 拆开装入余下容量部分背包核心 ratio[it.id] remain / it.weight total_value it.value * ratio[it.id] total_weight it.weight * ratio[it.id] remain 0 return ratio, total_value, total_weight主流程里构造数据后依次调用三种策略打印对比结果def run_experiment(): random.seed(2024) items generate_items(n60, max_weight50, max_value100) total_weight sum(i.weight for i in items) capacity int(total_weight * 0.5) print(f背包容量: {capacity}物品总数: {len(items)}) print(策略对比...) for strategy in [value, weight, density]: ratio, value, used_weight greed_solve(items, capacity, strategy) print(f[{strategy:8}] 总价值{value:8.2f} 使用重量{used_weight:6.2f}/{capacity}) if __name__ __main__: run_experiment()这段代码里有个细节值得注意我把比例数组按it.id保存而不是按排序后的位置保存。因为后面如果要分析“哪些物品被装入了”需要能对应回原始物品。直接按排序位置打标记会和物品原始编号对不上排查问题时容易混淆。3.3 结果验证与对比分析跑一组具体实验数据n60容量约1500我得到的典型输出大概是背包容量: 1473物品总数: 60 [ value] 总价值2745.00 使用重量1473.00/1473 [ weight] 总价值2519.00 使用重量1473.00/1473 [ density] 总价值3012.00 使用重量1473.00/1473可以看到单位价值策略在这个随机数据集上明显优于另外两种。价值优先策略比单位价值低约10%重量优先策略低更多。实验报告里建议至少跑三组不同随机种子记录三组结果取平均避免“一组数据恰巧适合某个策略”的偶然性。我这里把种子换成5、23、77再测试结论依然稳定。4. 实验过程中常见问题与排查技巧4.1 单位价值相同时的排序稳定性随机数据里可能出现多个物品的单位价值相同。Python的sorted是稳定排序会保留原始顺序。实验时如果追求可复现性最好在排序关键字里加一个次要键比如物品编号这样结果不依赖初始数据顺序sorted_items sorted(items, keylambda x: (x.value / x.weight, -x.id), reverseTrue)这种细节在实际报告中很容易被忽略但老师或评审看到会认为你考虑周全。4.2 浮点数精度导致的比例异常部分背包里计算remain / it.weight时涉及浮点数。当剩余容量接近某个物品重量时可能出现0.9999999这类比例。我建议在计算总价值时用浮点数但打印时保留两位小数。还要注意如果remain减到非常小的负数比如-1e-12可能影响后续循环判断。在循环开头加一个if remain 0: break能防止多余计算。4.3 物品重量总和小于容量如果生成数据时容量设得太大所有物品都被装入此时无论什么策略得到的总价值都一样实验就失去对比意义。所以我在run_experiment里特意用总重量的50%作为容量。但这个细节我在第一次做实验时栽过跟头——当时把容量设为固定常数500物品总重才300三种策略全部返回相同总价值一度以为代码写错了。4.4 0-1背包下贪心失效的验证方法实验要求如果额外加了一道“0-1背包不能用贪心”的验证题可以这样做生成一组不能拆分物品的数据用贪心得到结果再用动态规划跑一遍最优解计算差距。我常用一组“经典反例”数据容量10三个物品物品A重量8价值9单位价值1.125物品B重量5价值5单位价值1.0物品C重量5价值5单位价值1.0贪心按单位价值降序会先选A剩余容量2装不下B或C总价值9。但最优解是选B和C总价值10。这个反例非常直观能让学生清楚看到“贪心局部最优不等于全局最优”的边界。4.5 实验中如何验证贪心结果确为最优还有一个常被问到的点部分背包场景下怎么证明你实现的贪心结果是最优的除了逻辑证明实验上可以做一个交叉验证写一个穷举搜索对数据规模较小的场景枚举所有可能分配。因为物品可拆分理论上最优分配是确定的。用穷举验证前10个物品的数据确认贪心结果和穷举结果一致再放心跑大规模数据。我在实验指导时经常提醒学生“课堂上学了贪心能解部分背包但实验报告里只贴一个运行截图是不够的加一道穷举对比才能让结论站得住脚。”5. 实操心得与小技巧代码能跑通只是第一步我在写这个实验时实际折腾下来还有几个值得分享的小经验。第一数据可视化对排查bug帮助极大。我一开始只打印总价值跑出奇怪结果时很难定位问题。后来简单用文本图把一个策略的装入情况画出来类似“每个物品占一列方块长度代表重量比例”一眼就能看出排序策略是否按预期工作。如果不想引入绘图库纯字符输出也够用。第二代码里把“排序”单独拆成函数。这样实验报告里描述算法步骤时可以清晰画成“排序选择”两个阶段。如果全部揉在greed_solve里后面想测试不同排序算法对结果的影响还得重构比较麻烦。第三统一使用fractions或Decimal处理金额类价值数据。当物品价值是整数但单位价值是无限循环小数时直接用浮点排序可能造成错误的排序结果。随机数据里很少遇到这种极端情况但如果手工构造数据跑实验就可能碰上。稳妥做法是统一转成分数from fractions import Fraction keylambda x: Fraction(x.value * 10000, x.weight) # 放大整数转换避免浮点不过日常实验用浮点就够了我提这个只是提醒一个潜在边界。第四容量取整重量的百分比时建议加上取整逻辑。比如我用int(total_weight * 0.5)如果权重和过大这个百分比可能超出实际物品总重量的承受范围导致所有物品都被装下。更稳健的方式是capacity min(int(total_weight * 0.5), max_weight * n // 2)这个细节不大但能避免一些“所有策略结果相同”的尴尬。第五对比实验时保持随机种子一致。我上面代码就固定了random.seed(2024)这样每次运行得到相同的数据结果可复现。如果每次跑都生成全新随机数据三种策略的差异虽然整体趋势一致但具体数值对不上写实验报告时不好记录。可复现性是实验类文章的基本原则这一点很多初学者容易忽略。6. 一点个人经验做了好几次同样的实验我对贪心算法的理解比最初深刻不少。最核心的转变是以前我以为贪心就是“每步挑最好的”直到真的实现一遍并对比数据后才意识到“局部最优是否能推出全局最优”这个问题本身才是算法设计的核心问题。背包问题只是贪心思想的一个温床同样的思路还能延伸到活动选择、哈夫曼编码、最小生成树等经典场景。如果在实验中时间充裕我强烈建议额外做一次“反例测试”针对0-1背包故意构造一组数据让贪心结果差于最优解然后用动态规划求得最优解对比。这个测试做下来你对贪心适用边界的记忆会非常牢靠。课堂上学“贪心不一定得到最优”是一回事亲手看到差距是另一回事。这个实验的意义也正在于此。