问在python中编写非负流变量

EN

我正在用Python语言和Gurobi解决设施问题，其中一个决策变量是一个非负流变量g_ij > 0，用于表示设施之间的连接。也就是说，当且仅当y_i = 1和y_j = 1时为g_ij > 0，如果设施位于站点i，则为y_i = 1，否则为0。

我的约束之一是g_ij*(d_ij - k) >= 0

其中：

g_ij是非负流变量

d_ij是设施点i和j之间的距离(以公里为单位

k是Xkm

也就是说，设施点i和j之间的最大允许距离为Xkm

我将变量g_ij作为一个具有下界0和上界1的连续变量

m = Model("flowTest")
n = 10
g = {}

for i in range(n):
    for j in range(n):
        g[(i,j)] = m.addVar(lb=0, ub=1, vtype=GRB.CONTINUOUS, name="g%d%d" % (i,j))

以及约束条件：

for i in range(n):
    for j in range(n):
        m.addConstr(g[(i,j)] * (d[(i,j)] - k) >= 0)

但我得到了一个不可行的解决方案。我非常确定这是因为g_ij flow变量。

有谁能发现错误或建议我如何对g_ij变量进行编程？