如何在解决方案生成期间更改 Cplex 模型？

Question

我有一个包含多个约束和解决方案池的 Cplex 模型。我的限制之一是：

 R_alt=[i for i in R if i not in SetAlt]
        model.add_constraints((model.sum(x[i, j] for j in R2 ) == 2 for i in R_alt),"6C" )
        model.add_constraints((x[i, n1-4] ==x[i, n1-2]  for i in R_alt ),"7C" )

SetAlt 是一组 2 个值，将在进行约束之前从 R_alt 中删除。我需要 cplex 为每个解决方案随机选择这两个值。换句话说，我需要 cplex 在解决方案池生成期间更改此约束的模型。

例如，如果我在一个解决方案中的 R_alt=[0,1,2,3,6,7,8] 上有 6C，我会得到 R_alt=[0,2,3,4,5,6,8]在另一个解决方案中。

之前，我使用 python random 来选择这个 SetAlt，但问题是我在所有解决方案中都有相同的 SetAlt。

Answer 1

在https://github.com/AlexFleischerParis/zoodocplex/blob/master/zoomontecarlooptimization.py

import random
import math

random.seed(1)

from docplex.mp.model import Model

# original model

nbKids=300
mdl = Model(name='buses')
nbbus40 = mdl.integer_var(name='nbBus40')
nbbus30 = mdl.integer_var(name='nbBus30')
costBus40=500.0;
costBus30=400.0;
mdl.add_constraint(nbbus40*40 + nbbus30*30 >= nbKids, 'kids')
mdl.minimize(nbbus40*costBus40 + nbbus30*costBus30)

nbSamples=20
nbMaxKidsAbsent=30;

nbKidsLess=[random.randint(0,nbMaxKidsAbsent) for i in range(0,nbSamples)]
nbKidsOptions=[nbKids-nbKidsLess[i] for i in range(0,nbSamples)]

#Monte Carlo optimization

totalCost=0.0;
for i in range(0,nbSamples):
    
   mdl.get_constraint_by_name("kids").rhs=nbKidsOptions[i]
   mdl.solve()
   cost=mdl.solution.get_objective_value()
   totalCost+=cost
   print("if we need to bring ",nbKidsOptions[i]," kids  to the zoo");
   print("cost = ",cost)

print()   
averageCost=1/nbSamples*totalCost


print("------------------------------");
print("average cost = ",math.ceil(averageCost));

您可能会找到一个更改约束并再次求解的示例：

mdl.get_constraint_by_name("kids").rhs=nbKidsOptions[i]