import networkx as nx
import pyomo.core as pyo

def mis(graph: nx.Graph) -> pyo.ConcreteModel:
model = pyo.ConcreteModel()
model.x = pyo.Var(graph.nodes, domain=pyo.Binary)

@model.Constraint(graph.edges)
def independent_rule(model, node1, node2):
return model.x[node1] + model.x[node2] <= 1

model.cost = pyo.Objective(expr=sum(model.x.values()), sense=pyo.maximize)

return model


This function generates a PYOMO model which formulates the Maximum Independent Set problem. The problem concerns finding the largest set of vertices in a graph such that no two vertices are adjacent. The model consists of:

• Index set declarations (model.Nodes, model.Arcs).
• Binary variable declaration for each node (model.x) indicating whether that node is chosen to be included in the set.
• Constraint rule - for each edge we require at least one of the corresponding node variables to be 0.
• Objective rule – the sum of the variables equals to the set size.