Python 線形計画問題 その２

# -*- coding: utf-8 -*-

from pulp import LpProblem, LpVariable, LpMaximize, LpStatus

problem = LpProblem('sample', LpMaximize)

x1 = LpVariable('x1', lowBound=0,upBound=None,cat='Continuous')
x2 = LpVariable('x2', lowBound=0,upBound=None,cat='Continuous')

problem += (15*x1 + 20*x2, 'Objective function')

problem += (3*x1 + x2 <= 156, 'Constraint1')
problem += (x1 + 3*x2 <= 136, 'Constraint2')
problem += (8*x1 + 7*x2 <= 486, 'Constraint3')

print(problem)

result = problem.solve()

print(LpStatus[result])
print('objective function value', problem.objective.value())
print('x1 = ',x1.value())
print('x2 = ',x2.value())