Tengo un programa entero lineal que quiero resolver. Instalé solver glpk (gracias aesta respuesta) y pyomo. Escribí un código como este:

from pyomo.environ import *
from pyomo.opt import SolverFactory

a = 370
b = 420
c = 2

model             = ConcreteModel()
model.x           = Var([1,2], domain=NonNegativeIntegers)
model.Objective   = Objective(expr = a * model.x[1] + b * model.x[2], sense=minimize)
model.Constraint1 = Constraint(expr = model.x[1] + model.x[2] == c)
# ... more constraints

opt = SolverFactory('glpk')

results = opt.solve(model)

Esto produce una solución para presentarresults.yaml.

Tengo muchos problemas que quiero resolver usando el mismo modelo pero con diferentesa, byc valores. Quiero asignar diferentes valores aa, byc, resolver el modelo, obtener la solución demodel.x[1] ymodel.x[2]y tener una lista dea, b, c, model.x[1] ymodel.x[2]. Yo leodocumentación pero los ejemplos solo escriben soluciones a archivos comoresults.yaml.

¿Hay alguna forma de acceder a los valores de la solución desde el código?

Gracias,