我认为您在 f 中创建二进制“开-关”变量是正确的。
f = LpVariable('f',0,1,cat='Integer')
如果我理解正确的话,只要 A > 20,就需要 B > 5,对吧?
然后你需要调整你的代码,只要 A 大于 20,f 就设置为 1。
由于 f 是二进制的,因此您可以这样做:
prob+= f>= (A-20)/80
prob+= B>= 6*f
在第一行中,对于从 0 到 19 的 A 值,(A-20)/80 将为负值,并且当 A 为 20 时为零。这将确保对于这些 A 值,f 为零。
当 f 为零时,第二个约束仅意味着 B 必须至少为零,无论如何这是它的下界。
然而,如果 A 为 21 及以上,即 A>20,则 (A-20)/80 变为正数,但绝不会大于 1(稍后会详细介绍)。当 A 为 21 及以上时,这会强制 f 至少为 1。由于 f 只能为 1 或 0,因此 f 被设置为 1。
这将导致第二个约束在 B 为 1(即 A 大于 20)时强制 B 至少为 6。简而言之,只要 f 为 1(即 A 大于 20),B 就大于 5。
希望这可以帮助!如果不起作用,请告诉我。我自己一直在研究 puLP 问题,并使用这种方法编写了一些我的约束。
笔记:
我们除以 81 以确保 (A-20)/80 永远不会大于 1。例如,如果 A 为 21,则 (A-20)/80 的计算结果为 1/80。由于 A 只能大到 100,因此 (A-20)/80 只能大到 (100-20)/80,即 1。如果我们将其更改为 (A-20)/X,其中这里的 X 是低于 (A-20) 最大值的任何其他值,则分数 (A-20)/X 可以大于 1。并且由于 f 是二进制(1 或 0),因此约束 f>= (A-20)/X 意味着我们实际上会迫使 A 比其他情况更小。