Optimization, Operations Research, Python, Data Science, Editors Pick.
model.C.prev(c) return model.y[c] <= model.y[c_prev] model.fill_cstr = pyo.Constraint(model.N, rule=fill_cstr) model.edge_cstr = pyo.Constraint(model.E, model.C, rule=edge_cstr) model.break_symmetry = pyo.Constraint(model.C, rule=break_symmetry) You can try and include other symmetry-breaking