blob: 72448623b95771f43026a2f87ba4a4985b603b03 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
from pyomo.environ import *
A = ['hammer', 'wrench', 'screwdriver', 'towel']
b = {'hammer': 8, 'wrench': 3, 'screwdriver': 6, 'towel': 11}
w = {'hammer': 5, 'wrench': 7, 'screwdriver': 4, 'towel': 3}
W_max = 14
model = ConcreteModel()
model.x = Var(A, within=Binary)
model.value = Objective(
expr=sum(b[i] * model.x[i] for i in A), sense=maximize)
model.weight = Constraint(
expr=sum(w[i] * model.x[i] * model.x[i] for i in A) <= W_max)
opt = SolverFactory('ipopt')
result_obj = opt.solve(model, tee=True)
model.pprint()
|
