GAMS ord(i) 函数的 Python (gurobipy) 等价物

Question

我正在使用 gurobipy 在 Python 中制作作业车间调度模型，因为免费版本的 GAMS 无法处理这个特定问题的维度，但我完全是新手。

我的部分目标函数如下：

(sum((i,ii,c)$(ord(ii)>ord(i)),w(i,ii,c)*u(i,ii))*R)

ord() 给出了集合中元素的索引，因此w * u 的求和只针对那些索引大于i 的ii 元素。

我想知道是否有一种优雅的方式可以在 python 中编写这段代码，到目前为止我已经看到我可以使用 enumerate 函数来循环一维参数并保留索引，例如：

bs = [7, 6, 5, 4, 1]
cs = [2, 8, 3, 9, 0]

for i, b in enumerate(bs):
    for j, c in enumerate(cs):
        if i > j:
            print(b, " * ", c, " = ", b * c)

但我不知道如何使用多维参数和变量来做到这一点？

Answer 1

最优雅的方法是在一组元组上定义决策变量 w 和 u，其中这些元组仅包括 ord(i)<ord(ii) 的有效组合。然后你只需要遍历 u 或 w 的索引器。

使用 Gurobi，我会执行以下操作：

u = m.addVars(((i,ii) for m,i in enumerate(bs)
                      for n,ii in enumerate(cs) if m<n),
              name='u')