我是Google OR-Tools的新手(通常是约束编程),我试图为Jobshop示例添加截止日期,但实际上并没有用。
我在此求职的例子可以在这里找到:https://developers.google.com/optimization/scheduling/job_shop#entire-program
我改变了几件事:
我向Jobs_data和名为tuple的task_type添加了截止日期参数。
然后,我向模型添加了一个截止日期变量。
最后,我向模型添加了一个约束,其中task.end必须小于或等于截止日期。 model.Add(all_tasks[job_id, task_id].end <= all_tasks[job_id, task_id].deadline)
但这是行不通的。它会做一些事情,因为时间表输出更改,而当我输入不可能的截止日期时,它将返回0个结果。但这没有适当考虑截止日期。我究竟做错了什么?
这是修改后的版本的输出:
Optimal Schedule Length: 11
Machine 0: job_0_0 job_1_0
[0,3] [3,5]
Machine 1: job_2_0 job_0_1 job_1_2
[0,4] [4,6] [6,10]
Machine 2: job_1_1 job_0_2 job_2_1
[5,6] [6,8] [8,11]
如您所见,作业0的截止日期为7,但在时间表中,截止日期为8。
这是我完整的修改示例:
from __future__ import print_function
import collections
# Import Python wrapper for or-tools CP-SAT solver.
from ortools.sat.python import cp_model
def MinimalJobshopSat():
"""Minimal jobshop problem."""
# Create the model.
model = cp_model.CpModel()
jobs_data = [ # task = (machine_id, processing_time).
[(0, 3, 7), (1, 2, 7), (2, 2, 7)], # Job0
[(0, 2, 12), (2, 1, 12), (1, 4, 12)], # Job1
[(1, 4, 12), (2, 3, 12)] # Job2
]
machines_count = 1 + max(task[0] for job in jobs_data for task in job)
all_machines = range(machines_count)
# Computes horizon dynamically as the sum of all durations.
horizon = sum(task[1] for job in jobs_data for task in job)
# Named tuple to store information about created variables.
task_type = collections.namedtuple('task_type', 'start end deadline interval')
# Named tuple to manipulate solution information.
assigned_task_type = collections.namedtuple('assigned_task_type',
'start job index duration')
# Creates job intervals and add to the corresponding machine lists.
all_tasks = {}
machine_to_intervals = collections.defaultdict(list)
for job_id, job in enumerate(jobs_data):
for task_id, task in enumerate(job):
machine = task[0]
duration = task[1]
deadline = task[2]
suffix = '_%i_%i' % (job_id, task_id)
start_var = model.NewIntVar(0, horizon, 'start' + suffix)
end_var = model.NewIntVar(0, horizon, 'end' + suffix)
interval_var = model.NewIntervalVar(start_var, duration, end_var,
'interval' + suffix)
deadline_var = model.NewIntVar(deadline, deadline,
'deadline' + suffix)
all_tasks[job_id, task_id] = task_type(
start=start_var, end=end_var, deadline=deadline_var, interval=interval_var)
machine_to_intervals[machine].append(interval_var)
# Create and add disjunctive constraints.
for machine in all_machines:
model.AddNoOverlap(machine_to_intervals[machine])
# Precedences inside a job.
for job_id, job in enumerate(jobs_data):
for task_id in range(len(job) - 1):
model.Add(all_tasks[job_id, task_id].end <= all_tasks[job_id, task_id].deadline)
model.Add(all_tasks[job_id, task_id +
1].start >= all_tasks[job_id, task_id].end)
# Makespan objective.
obj_var = model.NewIntVar(0, horizon, 'makespan')
model.AddMaxEquality(obj_var, [
all_tasks[job_id, len(job) - 1].end
for job_id, job in enumerate(jobs_data)
])
model.Minimize(obj_var)
# Solve model.
solver = cp_model.CpSolver()
status = solver.Solve(model)
if status == cp_model.OPTIMAL:
# Create one list of assigned tasks per machine.
assigned_jobs = collections.defaultdict(list)
for job_id, job in enumerate(jobs_data):
for task_id, task in enumerate(job):
machine = task[0]
assigned_jobs[machine].append(
assigned_task_type(
start=solver.Value(all_tasks[job_id, task_id].start),
job=job_id,
index=task_id,
duration=task[1]))
# Create per machine output lines.
output = ''
for machine in all_machines:
# Sort by starting time.
assigned_jobs[machine].sort()
sol_line_tasks = 'Machine ' + str(machine) + ': '
sol_line = ' '
for assigned_task in assigned_jobs[machine]:
name = 'job_%i_%i' % (assigned_task.job, assigned_task.index)
# Add spaces to output to align columns.
sol_line_tasks += '%-10s' % name
start = assigned_task.start
duration = assigned_task.duration
sol_tmp = '[%i,%i]' % (start, start + duration)
# Add spaces to output to align columns.
sol_line += '%-10s' % sol_tmp
sol_line += '\n'
sol_line_tasks += '\n'
output += sol_line_tasks
output += sol_line
# Finally print the solution found.
print('Optimal Schedule Length: %i' % solver.ObjectiveValue())
print(output)
MinimalJobshopSat()
最佳答案
这是行不通的,因为您在for
循环中为优先级引入了约束。
创建一个新循环并从-1
中删除for task_id in range(len(job) - 1)
。
您还可以通过限制end_var的上限来设置最后期限。
此外,此github问题还有一些您可以使用的想法:https://github.com/google/or-tools/issues/960
关于python - 如何在Google OR-Tools作业示例中添加截止日期?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/58267491/