子行程消除约束pyomo

Question

给定一个图G（V，A），其中节点1是仓库，我正努力为Pyomo中类似TSP的问题制定以下子行程消除约束：

其中x_ij和y_h是我之前定义为二进制变量的二进制约束。

首先，我创建了一个包含所有可能子集S的字典，以便始终包含节点1：subsets_s。

然后，我一直在尝试这样的事情，但是我遇到了错误：

model1.con3=ConstraintList()
for h in model1.V:
    if h is not 1:    
         for i in model1.S:
             if h not in subsets_s[i]['nodes_subset']:
                 S=subsets_s[i]['nodes_subset']
                 for v in S:
                     print(v)
                     model1.con3.add(sum(sum(model1.x[v,j]) for j in  
                     model1.V if j not in S)>=model1.y[h])                                     

您有什么建议吗？

谢谢