0

我有一个模型，它有一个二进制变量 x [i] [j] [k]。我需要添加一个满足这个条件的约束：


if x[i][j][k1]==1  and  x[j][i][k2]==1 

Then:

 k2-k1>8

我有这段代码，但我认为它不正确：

  mdl.add((y[(i,j,k)]+y[(j,i,k1)]==2),(k1-k>8) )

我也，把这个：

 mdl.add(mdl.if_then(y[(i,j,k1)]+y[(j,i,k2)]==2, k2-k1>8))

但我收到了这个错误：

    raise DOcplexException(resolved_message)

DOcplexException: Expecting linear constraint, got: False

如何使用 cplex python API 对此进行建模？

4

2 回答 2

0

让我分享一下 if then 的例子

https://www.linkedin.com/pulse/making-optimization-simple-python-alex-fleischer/

from docplex.mp.model import Model

mdl = Model(name='buses')
nbbus40 = mdl.integer_var(name='nbBus40')
nbbus30 = mdl.integer_var(name='nbBus30')
mdl.add_constraint(nbbus40*40 + nbbus30*30 >= 300, 'kids')
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()

for v in mdl.iter_integer_vars():
   print(v," = ",v.solution_value)

print()
print("with if nb buses 40 more than 3  then nbBuses30 more than 7")

#if then constraint
mdl.add(mdl.if_then(nbbus40>=3,nbbus30>=7))
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()



for v in mdl.iter_integer_vars():
    print(v," = ",v.solution_value)

如果你想看到并且在 if

from docplex.mp.model import Model

mdl = Model(name='buses')
nbbus40 = mdl.integer_var(name='nbBus40')
nbbus30 = mdl.integer_var(name='nbBus30')
mdl.add_constraint(nbbus40*40 + nbbus30*30 >= 300, 'kids')
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()

for v in mdl.iter_integer_vars():
   print(v," = ",v.solution_value)

print()
print("with if nb buses 40 more than 3 and less than 7 then nbBuses30 more than 7")

#if then constraint
mdl.add(mdl.if_then((nbbus40>=3) + (nbbus40<=7)>=2,nbbus30>=7))
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()



for v in mdl.iter_integer_vars():
    print(v," = ",v.solution_value)

你也可以重写

from docplex.mp.model import Model

mdl = Model(name='buses')
nbbus40 = mdl.integer_var(name='nbBus40')
nbbus30 = mdl.integer_var(name='nbBus30')
mdl.add_constraint(nbbus40*40 + nbbus30*30 >= 300, 'kids')
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()

for v in mdl.iter_integer_vars():
   print(v," = ",v.solution_value)

print()
print("with if nb buses 40 more than 3 and less than 7 then nbBuses30 more than 7")

#if then constraint
mdl.add((((nbbus40>=3) + (nbbus40<=7)>=2))<=(nbbus30>=7))
mdl.minimize(nbbus40*500 + nbbus30*400)

mdl.solve()



for v in mdl.iter_integer_vars():
    print(v," = ",v.solution_value)
于 2020-06-15T16:01:48.493 回答
0

Model.if_then从左到右链接两个线性约束。如果第一个满意，第二个也会满意。

在您的情况下，我了解其中k2-k1>8不涉及决策变量。所以这纯粹是依赖于数据的，在模型构建时就知道了。

在这种情况下，因果关系反过来：如果 k2-k1>8，则 x[i,j,k1] 和 x[i,j,k2] 都必须等于 1。

那么最简单的代码是：

if k2-k1>8:
   m.add(x[i,j,k1] == 1)
   m.add(x[i,j,k2] == 1)
于 2020-06-16T14:31:31.823 回答