Networkx 节点可达性 - 如果条件

如何解决Networkx 节点可达性 - 如果条件

我正在尝试在 Python 中实现集合覆盖问题的基本版本。我的二元决策变量是否在节点处定位设施是 x[j]。如果两个节点之间的边的权重小于 4000，我正在努力解决一个约束，即一个节点 i 只能覆盖另一个节点 j。

我实施的约束是

#add constraint
for i in G.nodes():
    model.addConstr(sum(a[i,j]*x[j] >= 1 for j in G.nodes))

但我不知道如何定义 a(i,j) == 1 仅当节点之间的距离是“可覆盖”时。

在定义模型并将其用作参数之前，我曾尝试定义类似的东西，但它不起作用...

for i,j in G.nodes():
    if(nx.shortest_path_length(G,i,j) <= 4000):
        return a[i,j] == 1