Проблема транспортировки целлюлозы по маршрутам, а не по количеству

Ваш код не очень ясен и неполон; мы понятия не имеем, как выглядят затраты. Пожалуйста, улучшите его, например, объясняя или используя более ясные имена переменных.

Чтобы добавить параметр, обозначающий, где вы используете маршрут, в терминах LP должно работать следующее условие:

Let c_r = cost of using route r
    r   = whether route r is being used
    d_r = total quantity shipped over route r
    M   = a very big number (at least the sum of all quantities or the capacity of a truck)

min  sum(c_r * r)
s.t. Mr >= d_r
     d_r >= 0
     r in {0, 1}

Здесь, если по маршруту r ничего не отправлено, тогда r будет равно нулю, чтобы минимизировать целевую функцию, а если d_r > 0, то r будет 1, Mr = M, что будет работать, если d_r <= M. Таким образом, все зависит от значения, которое вы выбираете для M.

В терминах Python:

prob = LpProblem("LP Problem", LpMinimize)

route_vars = LpVariable.dicts("Route",(VendorStores, deficitStores), 0, None, LpInteger)  # d_r in the example
route_used = LpVariable.dicts("Route",(VendorStores, deficitStores), 0, 1, LpInteger)  # d_r in the example
a_very_large_number = 10000 # replace or calculate

# Objective function
prob += lpSum([(route_used[w][b],1)*costs[w][b] for (w,b) in Routes])

for w in VendorStores:
    prob += lpSum([route_vars[w][b] for b in deficitStores]) <= excess[w], "Sum of Cases out of VendorStore {0}".format(str(w))
for b in deficitStores:
    prob += lpSum([route_vars[w][b] for w in VendorStores]) >= deficit[b]

for (w, b) in routes:
    prob += a_very_large_number * route_used[w][b] >= route_vars[w][b], "Route used"