什么时候使用约束求解而不是机器学习

做者|Antoine Champion
编译|VK
来源|Towards Data Science

机器学习和深度学习一直是业界的热门话题。品牌领先于功能，致使深度学习在许多人工智能应用中被过分使用。

这篇文章将提供对约束求解的快速理解，这是一个强大但未被充分利用的方法，能够解决人工智能和其余计算机科学领域的大量问题，例如物流和调度时间推理和图形问题。

解决现实问题

让咱们来考虑一个事实性的和高度话题性的问题。

病人人数正在上升。医院必须迅速组织起来治疗病人。

世界上须要一种算法，在疾病严重程度、患者年龄和位置、医院容量和设备等多个标准下，将感染者和医院匹配起来。

许多人会说，神经网络将是最适合它的：它能够有不一样的配置，普遍的参数范围，能够根据须要减小到一个独特的解决方案。

然而，也有一些不利因素会破坏这个方案：

• 模型须要训练，所以须要之前案例的历史数据，

• 清理和整合数据集会浪费不少时间，

• 各类体系结构都须要经过冗长的训练而且要进行测试。

另外一方面，若是用一个布尔可知足性问题来描述，在不肯定多项式时间(NP彻底问题)中仍然给出次优解，而且不须要任何历史数据的状况下，不会有上述任何缺点。

这篇文章帮助你快速一览CSPs。理论和问题的表述将被忽略。有关更严格的方法，请参考论文，论文在文章的末尾

抽象问题

这篇文章将介绍约束编程，旨在解决这个案例。上面那张图说明了咱们算法的输出，应该该算法将感染者与医院匹配。现有几个框架用于约束求解。Google Optimization Tools(又称Tools)是一个用于解决组合优化问题的开源软件套件。咱们的问题将使用Python中的这个框架进行建模。

from ortools.sat.python import cp_model

colab:https://colab.research.google.com/drive/1vFkt5yIQtyelqvCh2TsJ9UDeM5miXqui

参数

如今，让咱们将问题简化为4个参数(1)：

• 感染者所在地

• 感染者的严重程度

• 医院位置

• 每家医院的床位数

让咱们用python定义这些参数：

# 医院数量
n_hospitals = 3
# 感染者人数
n_patients = 200
# 每家医院的床位数
n_beds_in_hospitals = [30,50,20]
# 病人位置，tuple (x,y)
patients_loc = [(randint(0, 100), randint(0, 100)) for _ in range(n_patients)]
# 医院位置，tuple (x,y)
hospitals_loc = [(randint(0, 100), randint(0, 100)) for _ in range(n_hospitals)]  
# 病人严重等级 1~5
patients_severity = [randint(1, 5) for _ in range(n_patients)]

变量

约束知足问题由一组变量组成，这些变量必须以知足一组约束。

• 令I为医院的集合

• 令\(J_i\)为医院i的床位集合

• 令\(K\)为病人集合

定义变量的索引族：

若是在医院i中，床j由病人k取走，则\(x_{ijk} = 1\)。为了将医院的每一张床与一个病人联系起来，咱们的目标是找到一组知足全部约束条件的变量。

咱们能够将这些变量添加到模型中：

model = cp_model.CpModel()
x = {}
for i in range(n_hospitals):
  for j in range(n_beds_in_hospitals[i]):
    for k in range(n_patients):
      x[(i,j,k)] = model.NewBoolVar("x(%d,%d,%d)" % (i,j,k))

硬约束

硬约束定义了模型的目标。它们是必不可少的，若是它们得不到解决，就没法解决问题：

• 每张床上最多只能有一我的
• 每一个人最多只能有一张床

让咱们关注第一个硬约束。对于每家医院的每一张床，我：

• 要么有一个惟一的病人k，

• 要么床是空的。

所以，能够用如下方式表示：

咱们的求解器是一个组合优化求解器，它只能处理整数约束。所以，必须转化为一个整数方程：

这个不等式能够加到咱们的模型中。

# 每张床最多只能住一我的
for i in range(n_hospitals):
  for j in range(n_beds_in_hospitals[i]):
    model.Add(sum(x[(i,j,k)] for k in range(n_patients)) <= 1)

接下来，第二个硬约束：对于每一个患者k：

• 要么他在一个惟一的病床上j在一个惟一的医院i
• 要么他在家。

同理，能够转化为一个整数不等式:

最后，能够将此约束添加到模型中。

# 每一个人最多只能睡一张床
for k in range(n_patients):
  inner_sum = []
  for i in range(n_hospitals):
    inner_sum.append(sum(x[(i,j,k)] for j in range(n_beds_in_hospitals[i]))) 
  model.Add(sum(inner_sum) <= 1)

软约束

接下来是软约束。这些都是很是须要的：咱们的解决方案必须尽量知足它们，但它们不是找到解决方案的必要条件：

• 每一个病人都应该躺在床上

• 每一个人都应该由最近的医院处理

• 病床不足时，应先处理病情严重的病人

当硬约束被建模为等式或不等式时，软约束是最小化或最大化的表达式。

设Ω为知足硬约束的全部解的集合。

每个病人都应该被安排在一张床上，这意味着最大限度地增长被占用的床的数量。

每一个人都应该由最近的医院处理，以尽可能减小每一个病人与其指定医院之间的距离。

若是没有足够的床位，应首先处理病情严重的病人，以最大限度地提升全部处理病人的总严重程度。经过表示sev(k)患者k的严重程度：

而后咱们能够将全部软约束简化为一个目标：

须要注意的是：这些软约束没有相同的定义域。

• 患者最大化约束范围从0到n，其中n是患者数，

• 病情严重性限制范围从0到5n

• 距离约束范围从0到全部i和k的最大欧几里得距离。

考虑到全部这些约束具备相同的优先级，咱们必须定义惩罚因子来平衡不一样的约束。

下面是相应的代码：

# 整数的距离函数
idist = lambda xy1, xy2: int(((xy1[0]-xy2[0])**2 + (xy1[1]-xy2[1])**2)**0.5)

gain_max_patients = 140
gain_severity = int(140/5)
gain_distance = -1
#最大化的目标
soft_csts = []
for i in range(n_hospitals):
  for j in range(n_beds_in_hospitals[i]):
    for k in range(n_patients):
      factor = \
        gain_max_patients \
        + gain_distance * idist(hospitals_loc[i], patients_loc[k]) \
        + gain_severity * patients_severity[k]
      soft_csts.append(factor * x[(i,j,k)])
model.Maximize(sum(soft_csts))

求解

如今咱们能够启动求解器了。它将试图在指定的时间限制内找到最优解。若是没法找到最优解，则返回最近的次优解。

solver = cp_model.CpSolver()
solver.parameters.max_time_in_seconds = 60.0
status = solver.Solve(model)

在咱们的例子中，求解器在2.5秒内返回一个最优解。

结论

要建立这个解决方案，只须要1小时的研究和30分钟的编程。

若是使用深度学习，要进行几天的数据清理，至少一天测试不一样的架构，另外一天进行训练。

此外，若是模型良好，CP-SAT模型是很是健壮的。下面是不一样模拟参数的结果。结果在许多不一样的状况下仍然是一致的，随着模拟参数的增长（3000名患者，1000张病床），解决方案推断只需不到3分钟。

固然，csp几乎不适用于计算机视觉和NLP等主题，在这些主题中，深度学习有时是最好的方法。然而，在物流、调度和计划方面，这每每是能够实现的方法。

深度学习的炒做激发了一些人尝试一些疯狂的举动来得到承认。有时，最好仍是经过阅读几篇关于你正在研究的问题的调查报告再想一想你应该如何解决。

引用

[1] Jingchao Chen, Solving Rubik’s Cube Using SAT Solvers, arXiv:1105.1436, 2011.

[2] Biere, A., Heule, M., and van Maaren, H. Handbook of satisfiability, volume 185. IOS press, 2009a

[3] Knuth, D. E., The art of computer programming, Volume 4, Fascicle 6: Satisfiability. Addison-Wesley Professional, 2015

[4] Vipin Kumar, Algorithms for constraint-satisfaction problems: a survey, AI Magazine Volume 13, Issue 1, 1992.

原文连接：https://towardsdatascience.com/where-you-should-drop-deep-learning-in-favor-of-constraint-solvers-eaab9f11ef45

欢迎关注磐创AI博客站：
http://panchuang.net/

sklearn机器学习中文官方文档：
http://sklearn123.com/

欢迎关注磐创博客资源汇总站：
http://docs.panchuang.net/