问题描述
令 x
是 n 个变量的向量,定义为:x = M.addMVar(shape = n,vtype = GRB.BINARY,name = "x")
。令 A
为 n × n 矩阵。令 v
是一个由正整数组成的 n x 1 常数向量。我感兴趣的约束:
np.multiply(v,x) <= A @ x
但是,当我在 gurobi 中添加此约束时:
M.addConstr(np.multiply(v,x) <= A @ x,name = "c1")
它给出了错误:
文件“src/gurobipy/model.pxi”,第 3325 行,在 gurobipy.Model.addConstr 中
文件“src/gurobipy/model.pxi”,第 3586 行,在 gurobipy.Model.addMConstr 中
TypeError: 必须是实数,而不是 minexpr
知道为什么会这样吗?我已经研究了几个小时的解决方案。我目前的猜测是 Gurobi 对不等式两边出现的变量不满意。但是,我还没有想出解决方法。
解决方法
这行不通,因为 Gurobi Python 接口中的矩阵表示需要规范形式 Ax = b。所以你需要将 v 合并到 A 矩阵中,你的代码就变成了这样:
A2 = A-v*np.eye(n)
M.addConstr(A2 @ x >= 0)