vrp 节约算法 c++_优化 | 浅谈旅行商问题（TSP）的启发式算法

最新推荐文章于 2024-06-24 14:08:22 发布

weixin_39586915

最新推荐文章于 2024-06-24 14:08:22 发布

阅读量721

点赞数

文章标签： vrp 节约算法 c++

↑↑↑↑↑点击上方蓝色字关注我们！

『运筹OR帷幄』原创

作者：高源

编者按：作为经典NP-hard问题，TSP和VRP问题一直是业界和学界的关注重点。本文以业界为主要视角，通过VRP和TSP的联系做切入点来引入TSP问题，分类讨论了启发式算法在对称TSP与非对称TSP中的运用。一、TSP简介

TSP全称为Travelling Salesman Problem(旅行商问题)，通俗而言，它是指对于给定的一系列城市和每对城市之间的距离，找到访问每一座城市仅一次并回到起始城市的最短回路。

TSP问题在运筹学发展史上有重要的意义，1952年，Danzig, Fulkerson和Johnson成功的解决了美国本土分数不同州的48个城市和哥伦比亚特区共49个城市的TSP实例，使更多的人初次了解了组合优化研究的意义，也感受到了离散问题求解的准度。

因为TSP是NP-C问题，因此其没有多项式时间内可以求最优解的算法。(如果对于NP-C的定义不了解，可以简单理解为对于某些实例，随着规模的增加，求解最优解所需的时间是爆炸式的指数增长)。因此研究TSP的启发式算法就显得很重要了。

二、由难入简 - VRP与TSP的联系

在深入TSP前首先简单介绍车辆路径问题(VRP)。事实上，在实际生活中，大多数我们所遇见的问题大都是VRP问题。1959年，《The truck dispatching problem》的作者Dantzig和Ramser在书中写到，TSP可被看成是VRP的某一类子问题。因此，理解TSP将是研究更深层路径优化算法的基石。

VRP的数学描述即为对于图G = (V, A, E) 找到成本最小的闭迹(无重复边的闭合回路)组合遍历已知的点集合(V)和弧(单向路)(A)或边(双向路)(E)集合。看似抽象，但在实际生活中，这样的问题其实无处不在。

比如在供应链领域，最经典的VRP即卖家对客户的补货策略：车队从配送中心(0)出发，遍历所有客户(U)，满足客户货物需求后回到配送中心{0}，在一定约束下，达到路程最短，成本最低等目的。此时A⋃E=∅，闭迹仅需遍历已知点集合(V=U⋃{0})。

而进一步当没有容量和时间等约束且仅有一辆货车进行一次性补货时，该问题就是TSP问题，即找到成本最小的一条闭迹遍历已知的点集合(V)，且在最优解中所有点仅会出现一次。

三、TSP模型

那么首先还是先抛出TSP的数学模型，对于有向图，其中为弧ij的成本， TSP可以被分解为以下四个约束和其目标函数。首先最后一个约束指变量为二进制变量，其等于1当仅当弧ij属于解回路。

约束(1)和(2)用于定义闭合回路，即任意顶点j/ i，解回路中必有且仅有一条弧以其作为起点的同时有且仅有一条弧以其作为终点。

但是仅仅有约束(1)和(2)是不足以定义TSP，因为虽然保证了解为闭合回路却没有保证解仅含一条闭合回路，因此我们需要约束(3)以防止小圈的出现。

那么约束(3)可以理解为对于任意点集合S(非V且含2个顶点及以上)，必定存在至少有一个S中的顶点与S外的顶点相连接。综上，我们可以得到以下数学模型。

其中约束(3)尽管可被简化但仍使得该问题成为了NP问题，因此基于上述讨论，本文将简单介绍可用于解决TSP的启发式算法。

四、ATSP的启发式算法(基于指派问题)

TSP问题依据边的成本正反向是否相等可以分为对称 (STSP) 或非对称 (ATSP)。而STSP与ATSP的性质有许多不一样，因此以下将分开进行讨论。

注释：

以下讨论条件均为满足三角不等式下的度量TSP，即

该文中对于ATSP的启发式算法的讨论相对简单，因为可以发现当我们如果忽略约束 (3) ，该问题便成了简单的指派问题 (Assignment Problem (AP))。

而因为AP具有最优解整数性质(全幺模矩阵)，因此约束 (4) 可被改写为

综上，上述模型改写成了可快速求解的线性规划问题。

那么，当得到AP问题的最优解(p个遍历所有点的有向闭迹 )后，如果发现此时p=1，即仅有一个闭合回路，那么则为ATSP的最优解，如若不然则对解进行优化改写。

注释：

对于ATSP而言，如果成本矩阵非对称性非常强，那么当解 满足

时，我们可以认为当前解为相对最优；

而对于STSP，该值将在30%附近或更高，因为STSP如用AP作为下界，那么解中将会有很多成对出现的子迹。

基于初始解进行启发式算法的流程如下：

记p个有向子闭迹为, 若p=1，停止，该为ATSP的最优解，否则进入步骤2；
找到两个拥有最多覆盖点的子迹；
合并，使得合并增加的成本最小，另p = p - 1，如果p = 1，此时解为相对最优解，否则返回步骤2；

我们可以通过以下具体例子[2]进一步理解，假设有6个点，每个点之间的距离成本矩阵如下

通过AP模型，可以得到为

此时p=2，进行启发式算法合并两个子迹，得到

在上述简单例子中，由于AP初始解仅有2个闭合回路，因此我们仅需要通过一次循环就可以找到启发式算法的最优解，更一般而言，当初始解回路为M时，我们需要重复M-1次以得到最优解。

在供应链领域，ATSP问题通常是针对城内规划，路径多样且不对称，而对于STSP，多为城际规划，道路单一且大致对称。因此接下来我们将着重描述适用于更大型系统的STSP的启发式算法。

五、STSP的启发式算法一(Nearest Neighbor)

因为对称TSP性质的特殊性，模型将建立在无向图上，因此我们可以改写其数学模型为

约束(6)即等效于ATSP的约束(1)和(2)，而约束(7)即为防止子迹的出现。

现在我们先介绍简单的最近邻点算法(nearest neighbor heuristic)的具体步骤：

任选点作为起点，令
找到下一邻点使得 , 将k加入C中
当|C|=|V|则算法停止，否则记h=k重复步骤2

下面用一个最近邻点算法的具体实例来进一步阐释：

选择最左端点为起点(1)
(1)其可与最右端点(9)，第(2)和(3)端点相连，选择有最小成本1的(3)与(1)相连。

重复上述过程后，可知最近邻点发的最优解为红色路径，成本为10，而通过肉眼观察容易发现最优解即为直接依次从最左端至最右端，其最优长度为9+5ϵ。

基于上述理解，结合严格的数学证明可知，最近邻点算法的所得长度与最优长度相比可能随点集数量的增加而趋于无穷，但因为其在TSPLIB(a library of sample instances for the TSP (and related problems) from various sources and of various types)实例中所得到的结果对比平均值，即(Nearest Neighbor的最优解/ 目前已知最优解)，约在1.26附近且算法简单容易实现，因此也是广为实用的一个算法。