MIP с использованием пульпы не приближается к результату - PullRequest
Новая
       6

MIP с использованием пульпы не приближается к результату

0 голосов
/ 31 декабря 2018

Я пытаюсь решить проблему MIP.Я пытаюсь найти количество экзаменов, которые должны быть выполнены каждым специалистом на дату в течение недели, сводя к минимуму общее количество используемых технологий.У меня есть спрос, время, затрачиваемое каждым техником, список техников и т. Д. В отдельных фреймах данных.

Изначально я использовал функцию стоимости, чтобы минимизировать общее время, затрачиваемое на удовлетворение спроса, которое @kabdulla помогло мне решить, ссылка здесь !

Теперь, с новой функцией стоимости, скрипт застревает и, кажется, не сходится, и я не могу определить причину.

Нижемой код до сих пор: 

# Instantiate problem class
model = pulp.LpProblem("Time minimizing problem", pulp.LpMinimize)

capacity = pulp.LpVariable.dicts("capacity",
                                 ((examdate , techname, region) for examdate, techname, region in tech_data_new.index),
                                 lowBound=0,
                                 cat='Integer')

tech_used = pulp.LpVariable.dicts("techs",
                             ((examdate,techname) for examdate,techname,region in tech_data_new.index.unique()),
                             cat='Binary')

model += pulp.lpSum(tech_used[examdate, techname] for examdate,techname in date_techname_index.index.unique())

for date in demand_data.index.get_level_values('Exam Date').unique():
    for i in demand_data.loc[date].index.tolist():
        model += pulp.lpSum([capacity[examdate,techname,region] for examdate, techname, region in tech_data_new.index if (date == examdate and i == region)]) == demand_data.loc[(demand_data.index.get_level_values('Exam Date') == date) & (demand_data.index.get_level_values('Body Region') == i), shiftname].item()

for examdate, techname,region in tech_data_new.index:
    model += (capacity[examdate, techname, region]) <= tech_data_new.loc[(examdate,techname,region), 'Max Capacity']*tech_used[examdate, techname]

# Number of techs used in a day should be less than 8    
for examdate in tech_data_new.index.get_level_values('Exam Date').unique():
    model += pulp.lpSum(tech_used[examdate, techname] for techname in tech_data_new.index.get_level_values('Technologist Name').unique()) <=8


# Max time each tech should work in a day should be less than 8 hours(28800 secs)    
for date in tech_data_new.index.get_level_values('Exam Date').unique():
    for name in tech_data_new.loc[date].index.get_level_values('Technologist Name').unique():
        #print(name)
        model += pulp.lpSum(capacity[examdate,techname,region] * tech_data_new.loc[(examdate,techname,region), 'Time taken'] for examdate, techname, region in tech_data_new.index if (date == examdate and name == techname)) <= 28800

Последнее условие кажется проблемой, если я ее уберу, проблема сходится.Однако я не могу понять проблему.

Пожалуйста, дайте мне знать, чего мне не хватает в моем понимании.Спасибо.

...