Python 使用pybnb包进行分支和绑定

我想用分枝定界法解决以下类型的整数线性规划问题：
最大（Ax）
s、 t.Cx

import pyomo.environ as pmo
from pyomo.opt import SolverFactory

M = pmo.ConcreteModel()
N = [1, 2]
M.x = pmo.Var(N, bounds=(0, None), domain=pmo.PositiveIntegers)
M.o = pmo.Objective(expr=100*M.x[1] + 150*M.x[2], sense=pmo.maximize)
M.c1 = pmo.Constraint(expr=8000*M.x[1] + 4000*M.x[2] <= 40000)
M.c2 = pmo.Constraint(expr=15*M.x[1] + 30*M.x[2] <= 200)
model = M
solver = SolverFactory('glpk')
solver.solve(model, report_timing=True, tee=True)