00 前言

各位小伙伴大家好，相信大家已经看过前面column generation求解vehicle routing problems的过程详解。该问题中，子问题主要是找到一条reduced cost最小的合法路径，然后加入到Master Problem中。其实，子问题也是一个著名的NP-Hard问题，今天我们就来介绍一下。

01 ESPPRC

考虑图1.1中描述的网络。 除了每条边的成本c_ij之外，还存在经过边（i，j）的所消耗的资源t_ij，比如时间。 我们的目标是找到从开始节点到结束节点的最短路径，每个节点只能访问一次，同时使得资源消耗满足可用的资源约束，比如全程不能超过多少时间。[1]

当然上面描述问题只是ESPPRC中的一个例子，实际的资源约束可能有很多种，比如在VRPTW的子问题中：[2]

起始节点和结束节点一样，每个节点有固定的时间窗和固定的需求。车辆不能超过容量约束的要求等等。

ESPPRC vs SPPRC

SPPRC和ESPPRC一样，只不过SPPRC去掉了elementary的约束，允许最短路中一个节点被访问多次。

02 应用

我们知道，ESPPRC是可以应用在column generation中的算法框架中的。那么具体是怎么应用的呢？我们知道，在column generation中，subproblem每次迭代就是找一条reduced cost最小的路径，然后加入到Master Problem中。但是对于ESPPRC来说，每次的cost一样的，那不每次都求出同一条路径吗？？？

不！在column generation中，其子问题ESPPRC中边的cost是会随着Master Problem的求得到对偶变量改变而改变的。还记得reduced cost是怎么计算的吗？

其中，式子(22)可以表达成式子(23)，(23)是什么意思呢？b_ijk意思是边ij是否在路径k中。那么，在每一次迭代中，我们就可以利用Master Problem的对偶变量，来更新ESPPRC中每条边的cost，最终求得的路径cost就是column generation中的reduced cost。

03 常见的算法

ESPPRC的建模如下：[4]


求解SPPRC和ESPPRC常见的算法主要有以下几种：[3]

  • Dynamic programming and labeling algorithms
  • Lagrangean relaxation
  • Constraint programming(建模)
  • Heuristics

04 Pulse Algorithm

这一节介绍一个ESPPRC的精确算法Pulse Algorithm[5]，算法的伪代码如下：

其中：

  • r：表示到达当前节点时的cost
  • q：表示到达当前节点时的装载量
  • t：表示到达当前节点所需的总时间，早到需要等待，不能晚到

大体的思想是通过bound算法确定到达每个节点的最低cost，然后pulse进行路径搜索，而之前bound求出来的最低cost就可以在pulse搜索的过程中起到定界的作用，去掉一些不好的路径。

bound的算法如下：

每个节点的最低cost（相当于一个lower bound）是怎么算出来的呢？简单来说是通过限定每个节点最晚到达时间。在每个节点都给定最晚到达时间t，然后计算在这个最晚到达时间下，每个到达每个节点的最低cost。然后让t递减，直到t减少到临界值。最终得到的是每个节点关于各个最晚到达时间的最低cost矩阵。

最后来看看pulse算法：

pulse是找路的过程，在该过程中：

  • isFeasible检查到达节点时路径是否满足各种资源约束。
  • checkBounds检查在当前到达时间下，到达节点的cost是否小于等于之前bound算法求出来的那个最晚时间下的最低cost，如果不是就丢弃该支路径。
  • rollback进行如下检查：

在每一个节点（开始节点除外），比如上图中节点j，如果实线路径的cost < 虚线路径的cost。那么砍掉实线的路径（已经有人的cost比你更低，你可以滚了），这就相当于一个回滚的操作。

以上>>>>>>>>>>>>>>>>>>>>>就是整个pulse算法。这里只是起到一个抛砖引玉的过程。可能讲的不是很详细，详细的过程请大家去阅读文献。

05 算法代码

关于整个pulse算法伪代码和讲解已经够详细了，这里给出一个C++的实现代码，是我远房的一个学长的学姐的男朋友写的（真话）。关于编译运行上面也说明得够详细了。

代码下载请关注我们的公众号哦！在后台回复【espprc】不包括【】即可下载。

reference

  • [1] A PRIMER IN COLUMN GENERATION, Jacques Desrosiers Marco E. L ubbecke
  • [2] A tutorial on column generation and branch-and-price for vehicle routing problems, Dominique Feillet
  • [3] SHORTEST PATH PROBLEMS WITH RESOURCE CONSTRAINTS, Stefan Irnich Guy Desaulniers
  • [4] An Exact Algorithm for the Elementary Shortest Path Problem with Resource Constraints: Application to Some Vehicle Routing Problems, Dominique Feillet, Pierre Dejax, Michel Gendreau, Cyrille Gueguen
  • [5] An Exact Algorithm for the Elementary Shortest Path Problem with Resource Constraints, Leonardo Lozano, Daniel Duque, Andrés L. Medaglia

干货 | 列生成VRPTW子问题ESPPRC（ Elementary shortest path problem with resource constraints）介绍附C++代码的更多相关文章

  1. 干货|技术小白如何在45分钟内发行通证(TOKEN)并上线交易（附流程代码

    https://blog.csdn.net/HiBlock/article/details/80071478

  2. 干货 | 10分钟带你彻底了解column generation（列生成）算法的原理附java代码

    OUTLINE 前言 预备知识预警 什么是column generation 相关概念科普 Cutting Stock Problem CG求解Cutting Stock Problem 列生成代码 ...

  3. 列生成算法（求解Cutting Stock问题）

    列生成是用于求解大规模线性优化问题的一种算法,其实就是单纯形法的一种形式.单纯性可以通过不断迭代,通过换基变量的操作,最终找到问题的最优解.但是当问题的规模很大之后,变量的个数就会增大到在有限时间内无 ...

  4. 根据oracle的主键列生成SQLserver的主键

    根据oracle的主键列生成MsSQLServer的主键列 select 'alter table  ' || cu.table_name ||'  add constraint  '||' PK_' ...

  5. pandas 选择列或者添加列生成新的DataFrame

    选择某些列 import pandas as pd # 从Excel中读取数据,生成DataFrame数据 # 导入Excel路径和sheet name df = pd.read_excel(exce ...

  6. Bootstrap Blazor Table 组件（四）自定义列生成

    原文链接:https://www.cnblogs.com/ysmc/p/16223154.html Bootstrap Blazor 官方链接:https://www.blazor.zone/tabl ...

  7. excel的列生成算法

    echo '<pre>'; $i = 1; while($i < 703){ $char1 = floor($i / 26); $char2 = $i % 26; if($i % 2 ...

  8. pandas 由其中几列生成新的列

    data是一个dataframe #data["x1"]=data[["a","b"]].apply(lambda x:x["a& ...

  9. 模板列传值到子窗体中，子窗体中多选gridview中checkbox保存数据多项到数据库中

    <html xmlns="http://www.w3.org/1999/xhtml"> <head runat="server"> &l ...

随机推荐

  1. CentOS7安装MySQL（完整版）

    在CentOS中默认安装有MariaDB,这个是MySQL的分支,但为了需要,还是要在系统中安装MySQL,而且安装完成之后可以直接覆盖掉MariaDB. 1 下载并安装MySQL官方的 Yum Re ...

  2. oracle基础知识语法大全

    ORACLE支持五种类型的完整性约束NOT NULL (非空)--防止NULL值进入指定的列,在单列基础上定义,默认情况下,ORACLE允许在任何列中有NULL值.CHECK (检查)--检查在约束中 ...

  3. (四) Docker 使用Let's Encrypt 部署 HTTPS

    参考并感谢 周花卷 https://www.jianshu.com/p/5afc6bbeb28c 下载letsencrypt镜像(不带tag标签则表示下载latest版本) docker pull q ...

  4. 单点logi,n

    using Newtonsoft.Json; using Newtonsoft.Json.Linq; using System; using System.Collections.Generic; u ...

  5. 神奇的外部嵌套（使用ROW_NUMBER()查询带条件的时候提示列名无效）

    declare @pageIndex int -- 第几页 declare @pageSize int -- 每页包含的记录数 --这里注意一下,不能直接把变量放在这里,要用select select ...

  6. Android 下拉列表Spinner 使用

    1.布局文件 <LinearLayout xmlns:android="http://schemas.android.com/apk/res/android" xmlns:t ...

  7. CentOS 7 使用 firewalld 打开关闭防火墙与端口

    1.firewalld的基本使用启动: systemctl start firewalld关闭: systemctl stop firewalld查看状态: systemctl status fire ...

  8. GNS3、Wireshark、SecureCRT 环境部署

    本次GNS3环境部署教程基于官方推荐的稳定版1.5.4.初次接触此软件,详细的使用方法不是很清楚,所以以此作为学习记录,仅供参考,后期补充. 软件介绍 GNS3 GNS3是一款具有图形化界面可以运行在 ...

  9. XAMPP + PhpStorm + Xdebug本地实验环境搭建

    参考: 下载合适的XDebug 点击这里,选择合适xdebug XAMPP配置 php_xdebug-xxxx.dll 拷贝dll至 D:\XAMPP\php\ext php.ini 文末追加 [XD ...

  10. Linux命令——pidof

    参考:Linux pidof Command Examples To Find PID of A Program/Command Linux pidof Command Tutorial for Be ...