ortool 整数规划之车间调度提前停止搜索

最新推荐文章于 2023-12-06 11:34:12 发布

石头的小软

最新推荐文章于 2023-12-06 11:34:12 发布

阅读量596

点赞数 1

分类专栏： python库文章标签： jsp

本文链接：https://blog.csdn.net/weixin_42158546/article/details/117962006

版权

python库专栏收录该内容

2 篇文章

订阅专栏

这篇博客介绍了如何使用ortools解决大规模作业车间调度问题，并在短时间内找到解决方案。通过设置求解个数或时间限制，可以在满足条件时提前停止搜索。示例代码展示了如何创建解决方案回调函数以限制求解数量或时间，从而优化搜索过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

ortool解决作业车间调度问题

当问题规模比较大时短时间无法求得解，可以设置条件提前停止求解

1，以设置的求解个数来作为终止条件

"""Code sample that solves a model and displays a small number of solutions."""

from ortools.sat.python import cp_model


class VarArraySolutionPrinterWithLimit(cp_model.CpSolverSolutionCallback):
  """Print intermediate solutions."""

  def __init__(self, variables, limit):
    cp_model.CpSolverSolutionCallback.__init__(self)
    self.__variables = variables
    self.__solution_count = 0
    self.__solution_limit = limit

  def on_solution_callback(self):
    self.__solution_count += 1
    for v in self.__variables:
      print('%s=%i' % (v, self.Value(v)), end=' ')
    #   print('  objective value = %i' % self.ObjectiveValue())
    print()
    if self.__solution_count >= self.__solution_limit:
      print('Stop search after %i solutions' % self.__solution_limit)
      
      self.StopSearch()

  def solution_count(self):
    return self.__solution_count


def StopAfterNSolutionsSampleSat():
  """Showcases calling the solver to search for small number of solutions."""
  ###设置停止条件
  stop_early_count=2
  # Creates the model.
  model = cp_model.CpModel()
  # Creates the variables.
  num_vals = 3
  x = model.NewIntVar(0, num_vals - 1, 'x')
  y = model.NewIntVar(0, num_vals - 1, 'y')
  z = model.NewIntVar(0, num_vals - 1, 'z')
  model.Add(x != y)
  model.Maximize(x + 2 * y + 3 * z)
  # Create a solver and solve.
  solver = cp_model.CpSolver()
  solution_printer = VarArraySolutionPrinterWithLimit([x, y, z], stop_early_count)
  # Enumerate all solutions.
  solver.parameters.enumerate_all_solutions = True
  # Solve.
  status = solver.SolveWithSolutionCallback(model, solution_printer)
  print('Status = %s' % solver.StatusName(status))
  print('Number of solutions found: %i' % solution_printer.solution_count())
  assert solution_printer.solution_count() == stop_early_count
  print(solver.ObjectiveValue())


StopAfterNSolutionsSampleSat()

2.设置时间为停止条件

from ortools.sat.python import cp_model
from time import time

class VarArraySolutionPrinterWithLimit(cp_model.CpSolverSolutionCallback):
  """Print intermediate solutions."""

  def __init__(self,  limit):
    cp_model.CpSolverSolutionCallback.__init__(self)
    self.__solution_count = time()
    self.__solution_limit = limit

  def on_solution_callback(self):
    t1=time()
    if (t1-self.__solution_count) >= self.__solution_limit:
      self.StopSearch()



def StopAfterNSolutionsSampleSat():
  """Showcases calling the solver to search for small number of solutions."""
  ###设置停止条件
  stop_early_time=10
  # Creates the model.
  model = cp_model.CpModel()
  # Creates the variables.
  num_vals = 3
  x = model.NewIntVar(0, num_vals - 1, 'x')
  y = model.NewIntVar(0, num_vals - 1, 'y')
  z = model.NewIntVar(0, num_vals - 1, 'z')
  model.Add(x != y)
  model.Maximize(x + 2 * y + 3 * z)
  # Create a solver and solve.
  solver = cp_model.CpSolver()
  solution_printer = VarArraySolutionPrinterWithLimit( stop_early_time)
  # Enumerate all solutions.
  solver.parameters.enumerate_all_solutions = True
  # Solve.
  status = solver.SolveWithSolutionCallback(model, solution_printer)
  print('Status = %s' % solver.StatusName(status))
  print(solver.ObjectiveValue())


StopAfterNSolutionsSampleSat()