旅行商问题，车辆路径问题的数学规划模型

根据书籍 Snyder, Lawrence V., and Zuo-Jun Max Shen. Fundamentals of supply chain theory. John Wiley & Sons, 2019 中的相关章节整理。

1. 旅行商问题

旅行商问题（Travelling Salesman Problem): 找到一条将所有地点走完，除了起点外，每个地点只经过一次，再回到起点的最短路线。

下图是一个美国几十个城市的 TSP 路线。

参数

• N = { 1 , 2 , … , n } N=\{1,2,\dots,n\}\quad N={1,2,…,n} 地点的集合
• $c_{ij}$ 从地点 $i$ 到地点 $j$ 的距离
• $x_{ij}$ 0-1决策变量，表示路线是否经过从地点 $i$ 到地点 $j$

由于从 i 到 j 到距离与从 j 到 i 的距离一般是一样的（$c_{ij}=c_{ji}$,对称的），因此 $x_{ij}=x_{ji}$.

数学模型

构建数学规划模型如下：

min ⁡ ∑ i ∈ N ∑ j ∈ N c i j x i j s.t. ∑ i ∈ N x i h + ∑ j ∈ N x h j = 2 ∀ h ∈ N ∑ i , j ∈ S x i j ≤ ∣ S ∣ − 1 ∀ S ⊂ N , 2 ≤ ∣ S ∣ ≤ n − 1 x i j ∈ { 0 , 1 } ∀ i , ∀ j \begin{align*} \min \quad &\sum_{i\in N}\sum_{j\in N} c_{ij}x_{ij}\\ \text{s.t.} \quad & \sum_{i\in N}x_{ih}+\sum_{j\in N}x_{hj}=2\qquad \forall h \in N\\ & \\ &\sum_{i,j\in S}x_{ij}\leq |S|-1\qquad \forall S\subset N, 2\leq |S|\leq n-1\tag{1} \\ &x_{ij}\in\{0,1\}\quad\forall i,\forall j \end{align*} mins.t.​i∈N∑​j∈N∑​cij​xij​i∈N∑​xih​+j∈N∑​xhj​=2∀h∈Ni,j∈S∑​xij​≤∣S∣−1∀S⊂N,2≤∣S∣≤n−1xij​∈{0,1}∀i,∀j​(1)​

第一个约束条件表示任何一个地点都通过（到达这个地点，从这个地点出发），第二个约束条件为消除子环条件（sub-tour elimination constraint) ，其中 $S$ 为 $N$ 的子集，$|S|$ 表示 $S$ 中的元素个数。

消除子环条件也可表示为：

$$\sum _{i\in S,j\in \overline{S}\text{ or }i\in \overline{S},j\in S}{x}_{ij}\le |S|-1\forall S\subset N,2\le |S|\le n-1$$

子环表示所有地点中的某几个地点的路线形成了一个闭环。例如，下图中，本来有 7 个地点，但是却存在两个子环：

2. 车辆路径问题

车辆路径问题（VRP）与 TSP 问题的区别是，路线中不仅有一辆车，有好几辆车。所有车辆的路线将所有地点走完，例如下图：

参数

因为 VRP 中不仅一辆车，涉及到分配车辆的问题，所以决策变量也多了些：

• N = { 0 , 1 , 2 , … , n } N=\{0,1,2,\dots,n\}\quad N={0,1,2,…,n} 地点的集合, 0 表示起点
• N − = { 1 , 2 , … , n } N^-=\{1,2,\dots,n\}\quad N−={1,2,…,n} 地点的集合, 0 表示起点
• $K$ 一共 $K$ 辆车
• $C$ 每个车辆的容量为 $C$
• $d_i$ 地点 $i$ 的需求量
   
• $c_{ij}$ 从地点 $i$ 到地点 $j$ 的距离
• $x_{ijk}$ 0-1决策变量，表示分配车辆 $k$ 从地点 $i$ 到地点 $j$
• $y_{ik}$ 0-1决策变量，表示车辆 $k$ 是否通过地点 $i$

数学模型

构建数学规划模型如下：

min ⁡ ∑ k = 1 K ∑ i ∈ N ∑ j ∈ N c i j x i j k s.t. ∑ k = 1 K y i k = 1 ∀ i ∈ N −   ∑ k = 1 K y 0 k = 1 ∑ i ∈ N x i h k = ∑ j ∈ N x h j k = y h k ∀ h ∈ N , ∀ k ∑ i ∈ N d i y i k ≤ C ∑ i , j ∈ S x i j k ≤ ∣ S ∣ − 1 ∀ k , ∀ S ⊆ N − , ∣ S ∣ ≥ 2 x i j k ∈ { 0 , 1 } ∀ i , ∀ j , ∀ k y i k ∈ { 0 , 1 } ∀ i , ∀ k \begin{align*} \min \quad &\sum_{k=1}^{K}\sum_{i\in N}\sum_{j\in N} c_{ij}x_{ijk}\\ \text{s.t.} \quad & \sum_{k=1}^Ky_{ik}=1\qquad \forall i \in N^-\ \\ & \sum_{k=1}^Ky_{0k}=1\\ &\sum_{i\in N} x_{ihk}=\sum_{j\in N} x_{hjk}=y_{hk}\quad \forall h\in N,\forall k\\ &\sum_{i\in N} d_iy_{ik}\leq C\\ & \\ &\sum_{i,j\in S}x_{ijk}\leq |S|-1\qquad \forall k, \forall S\subseteq N^-, |S|\geq 2 \\ &x_{ijk}\in\{0,1\}\quad\forall i,\forall j,\forall k\\ &y_{ik}\in\{0,1\}\quad\forall i,\forall k \end{align*} mins.t.​k=1∑K​i∈N∑​j∈N∑​cij​xijk​k=1∑K​yik​=1∀i∈N− k=1∑K​y0k​=1i∈N∑​xihk​=j∈N∑​xhjk​=yhk​∀h∈N,∀ki∈N∑​di​yik​≤Ci,j∈S∑​xijk​≤∣S∣−1∀k,∀S⊆N−,∣S∣≥2xijk​∈{0,1}∀i,∀j,∀kyik​∈{0,1}∀i,∀k​

第一个约束条件表示除了起点，每个地点只分配给一辆车；第二个约束条件表示起点出发 K 辆车；
第三个约束条件表示 $x$ 与 $y$ 的关系；第四个约束条件表示车辆的容量约束；第五个约束条件是子环消除约束。

第三个约束条件也可以写成：

$$\sum _{i\in N}{x}_{ihk}+\sum _{j\in N}{x}_{hjk}=2y_{hk}$$