Gurobi Python == 和变量

Question

== 在 Gurobi Python 中如何工作？

在我创建了所需的变量 (x[s,d,r]) 之后，我正在尝试创建一个类似于 x[s,d,r] == 1 的 if 语句。问题是，这个陈述总是正确的，所以不管我把它设置为x[s,d,r]==0 还是x[s,d,r]==2，它总是一样的。

x[s,d,r] 是由x[s,d,r] = m.addVar(0.0,1.0,1.0,GRB.BINARY,"x_"+ s+"_"+ d +"_" + r) 定义的二进制变量。

z = 0.0
for s in students:
    for d in dates:
        if (s,d) in preferences:        
            if preferences[s,d]!=0:
                for r in rooms:
                    if (d,r) in tutorials:
                        if x[s,d,r]>0.001:
                            print('%s:%s:%s:%s '%(s,d,r,preferences[s,d]))
                            z = z + preferences[s,d]
                        else:
                            print('no')                                 

 m.setObjective((z), GRB.MAXIMIZE)

所以如果 x[s,d,r]>0.001: 是始终评估为“真”的部分。

Answer 1

grb.Expr 的== 运算符返回一个 grb.Constraint 对象。正如answer 中提到的，约束对象总是真实的，所以你总是会得到'if'部分。

假设您已经优化，您想要优化解决方案中 x 值的实际值。这是由 grb.Variable 上的 X 属性给出的。 所以你应该用

替换你的代码

                        if x[s,d,r].X>0.001:

它应该可以按您的预期工作。

Answer 2

要对此进行调查，您应该在 Python shell 中单独评估 x[s,d,r] == 1。相信你会发现结果绝不是布尔值，而是另一个 Gurobi 对象。 （我希望像LinExpr 这样的东西。）

Gurobi（明智地）重载变量的所有算术运算，以便更容易构造约束。这会产生非常人类可读的模型构建代码。 （恕我直言，几乎与 OPL 或 AMPL 一样好，具有 Python 的所有其他优点）。

您是否要检查该变量的解值是否为 1？为此，我认为您想要x[s,d,r].x == 1，它将评估为布尔值。 （但前提是模型有解决方案）。