How to use a variable as a divisor in PuLP

问题

I was trying to resolve a LP problem with a constraint that is calculated by dividing variable A by variable B.

The simple version of the problem is as below:

The product is made by two materials (A and B)

% of A should be greater than 50%

% of B should be less than 40%

Total amount of A and B are 100

Objective: What's the minimum amount of A?

The code is like:

from pulp import *

prob = LpProblem('Simple problem', LpMinimize)
x = LpVariable('x', 0, None, 'Integer')
y = LpVariable('y', 0, None, 'Integer')
prob += x
prob += x / (x + y) > 0.5  # <== Where the error happens
prob += y / (x + y) < 0.4
prob += x + y == 100
prob.solve()

print 'Result: %s' % LpStatus[prob.status]
print 'Amount of A: %s' % value(prob.objective)

However I'm getting an error message saying:

TypeError: Expressions cannot be divided by a non-constant expression

It looks like PuLP does not support variable as divisor. https://github.com/coin-or/pulp/blob/master/src/pulp/pulp.py#L800

Any idea? If PuLP is not the right library to use, I'm happy to switch to any library that fits in.

UPDATE 27 Nov 2015

For some reason, the sample above does not make sense (not working as expected). I am very newbie to this library. Maybe it's not the right one to solve my problem at all. So if any one has suggestions of other libraries, it'd be appreciated.

BTW, Koen Peters's advice below is great. The error is gone after taking his advice. Thank you.

回答1:

Linear Programming doesn't understand divisions, hence the error :) You have to reformulate it so that the division is formulated linearly. In this case:

prob += x / (x + y) > 0.5  
prob += y / (x + y) < 0.4

is equivalent to:

prob += x > 0.5 * (x + y)
prob += y < 0.4 * (x + y)

Which are linear constraints. Good luck!

来源：https://stackoverflow.com/questions/33929353/how-to-use-a-variable-as-a-divisor-in-pulp

标签

python

pulp