变邻域搜索(VNS)原理梳理和应用细节-附求解VRPTW问题C++代码

1、简介-introduction

\qquad变邻域搜索(Variable Neighborhood Search)是一种求解组合优化和全局优化问题的元启发式算法,它的基本思想是:在局部搜索(LS)算法的框架中系统地变化邻域结构的选择。
\qquad下面本文将介绍以下几部分内容,第二部分本文将介绍VNS的基本机制;第三部分将介绍VNS的相关拓展;第四部分将介绍VNS和其他启发式算法的混杂应用;第五部分将介绍VNS在求解VRPTW中的具体思路,最后将给出使用VNS求解VRPTW问题的代码。

2、VNS的基本机制

\qquad首先定义N k , ( k = 1 , . . . , k m a x ) N_k,(k= 1,...,k_{max})Nk,(k=1,...,kmax)表示事先定义好的邻域结构集合,定义N k ( x ) N_k(x)Nk(x)表示第k kk个邻域结构中构造的解的集合。VNS基于三个简单的事实:
①某一个邻域结构的局部最优解不一定是另外一个邻域结构的局部最优解;
②对于所有邻域结构来说,全局最优解均是其局部最优解;
③对于许多问题情境,一个或者多个邻域的局部最优解应该彼此相近。
\qquad最后一个事实,虽然有点经验主义,但是暗示一个局部最优解通常会为全局最优解提供有效信息。

2.1 变邻域下降方法-Variable Neighborhood Descent

\qquad在VND中,事先选定的邻域结构按照一个确定的顺序进行执行,其流程如下图所示:
在这里插入图片描述
\qquad首先进行初始化,选择邻域结构N k , ( k = 1 , . . . , k m a x ) N_k,(k= 1,...,k_{max})Nk,(k=1,...,kmax),构建一个初始解,选定终止准则;之后循环执行VND,直到达到终止准则。其中VND分为两步,第一步为重置计数变量k = 1 k=1k=1,第二步为顺序执行邻域结构,直到k = k m a x k=k_{max}k=kmax。在每一次选定邻域结构进行执行的时候又可以分为两步(a)在第k kk个邻域结构生成的解方案中随机选择一个解方案x ′ x'x-shaking;(b)若选择的解方案x ′ x'x优于当前解,则使用x ′ x'x替换当前解,同时重新从第一个邻域结构开始执行优化;否则选择下一个邻域结构进行解的优化。
\qquad除了上述顺序执行的方法,VND还可以进行一定的改进,使用嵌套策略来执行邻域结构。假设有3个邻域结构,k m a x = 3 k_{max}=3kmax=3,则嵌套结构执行方法为在外层只使用第3个邻域结构,而对于第3个邻域结构生成的每一个解方案x ′ x'x,顺序执行邻域结构1和2组成的VND。
\qquad缩减VNS(Reduced Variable Neighborhood Search)方法同上述VND方法类似,但在shaking之后,不是对x ′ x'x执行后续的邻域操作,而是对当前最优解进行后续的邻域操作。

2.2基本变邻域搜索方法-Basic Variable Neighborhood Search

\qquadVNS综合了邻域选择的确定性和随机性因素,其流程如下图所示:
在这里插入图片描述
\qquad基本VNS与VND的区别主要在于(2.b)中在执行完邻域操作随机选定一个解x ′ x'x之后,VNS不是直接判断解x ′ x'x的优劣,而是对其进行进一步的局部搜索(LS)探索,将LS得到的解x ′ ′ x''x与当前解进行判断,保留较优的解。

2.3 终止准则

\qquad终止准则的选择可以有但不限于以下三种方式:
①规定一个最大的CPU运算时间;
②规定一个总的最大的迭代次数;
③规定一个最大迭代次数m a x I t e maxItemaxIte,若在m a x I t e maxItemaxIte之内解没有发生优化,则终止算法。

3、VNS的拓展和混杂算法

3.1 VNS的拓展

\qquad基本的VNS是一种基于descent,first improvement method和randomization的算法。在基本VNS流程的(2.c)中,可以加入一定的概率机制,当解x ′ ′ x''x次于当前解时,以一定的几率接受退化的解,来实现一种descent-ascent 机制;在选择邻域结构时,也可以加入轮盘赌的机制,优先选择优化效果较好的邻域结构进行执行;在(2.a)中可以选出最优的解x ′ x'x来进行后续的优化;还可以规定一个k m i n k_{min}kmink s t e p k_{step}kstep,替换k ← 1 k←1k1k ← k m i n k←k_{min}kkmink ← k + 1 k←k+1kk+1k ← k s t e p k←k_{step}kkstep

3.2 变邻域分解搜索-Variable Neighborhood Decomposition Search

\qquadVNDS方法将基本VNS拓展称为一种两级VNS,VNDS的步骤如下所示:
在这里插入图片描述
\qquadVNDS相对于基本VNS的主要区别在于(2.b)中,VNDS不是直接对于解x ′ x'x进行局部搜索LS操作,而是在(2.a)时选出一个解集合y = x ′ /   x y=x'/\ xy=x/ x,在新的解空间y yy(全部解空间的子集)中进行局部寻优,记寻找到的最优解为y ′ y'y,整个解空间中的最优解记为x ′ ′ x''x,之后再执行(2.c)更新最优解或者执行下一次邻域结构优化。

3.3 偏变邻域搜索-Skewed VNS

\qquadSVNS用来探索距离当前解较远的另外的求解空间,以求得到不同的解结构,寻找到更优的解,SVNS的流程如下图所示:
在这里插入图片描述
\qquad从上述流程可以看出,SVNS相对于基本VNS,主要将(2.c)分解成为了两部分,新的(2.c)用来判断是否进行更新最优解,新的(2.d)用来判断是否进行重新开始邻域搜索还是继续进行邻域搜索,其中利用了函数ρ ( x , x ′ ′ ) \rho(x,x'')ρ(x,x)来判断当前解x xx和局部最优解x ′ ′ x''x之间的距离,(本人理解:两个解之间的差异程度,如在VRPTW中,可以指两个路径方案用的车辆数,访问客户点的顺序差异程度等等),当ρ ( x , x ′ ′ ) \rho(x,x'')ρ(x,x)的值较小的时候,参数α \alphaα取值应该较大,来使求解空间转移到差距更大的地方进行寻优。α ρ ( x , x ′ ′ ) \alpha\rho(x,x'')αρ(x,x)的选取还可以采取学习机制或者自适应的机制进行更新。

3.4 并行变邻域搜索-Parallel VNS

\qquad在(2.b)时采用多线程的编程方式,对于某一个邻域产生的多个解方案来说,给每一个解方案开一个进程进行局部搜索LS,最终保留所有进程中求解结果最好的结果。

4、VNS和其他启发式算法的混合使用-hybrids

4.1 VNS 和 TS

\qquad有两种方式来进行VNS和TS的混合使用,第一种是在TS框架之中使用VNS,第二种是在VNS框架中使用TS。
\qquad在VNS框架中使用TS时,只需要将禁忌表机制加入到VNS搜索的过程即可,为所有的邻域搜索算子建立一个禁忌表,规定解禁步长和宽恕准则,可以有效避免backward搜索。
\qquad在TS框架中使用VNS时,只需要将TS中解的优化步骤中的LS优化变为VNS优化即可。

4.2 VNS 和 GRASP

\qquadGRASP是一种两阶段元启发式算法,在第一阶段使用贪婪随机的思想进行初始解的构建,在第二阶段使用LS或者枚举的思想进行解的优化。所以,一种顺其自然的方式来结合VNS和GRASP是在GRASP的第二阶段使用VNS来进行优化。

5、利用VNS框架求解VRPTW

\qquad首先选定邻域操作算子,本文在实现VNS时使用3种邻域操作算子,分别是2-opt inter-route算子,2-opt intra-route 算子和1-opt inter-route算子;在实现时使用基本VNS框架进行实现;终止准则本文使用了规定最大的迭代次数的终止准则。在基本VNS框架上,本文进行了些许改进,在内循环依次执行邻域操作算子时,本文只使用2-opt inter-route算子和2-opt intra-route 算子两种算子进行当前最优解的优化,在内循环结束之后,使用1-opt inter-route算子进行LS局部搜索优化,寻找更优的解。为了避免陷入局部最优,本文使用了“整条路径”破坏的思想,来跳出局部最优解,以求找到全局最优解。

6、求解代码已经上传值GitHub

https://github.com/Dragon-wang-fei-long/Huristic-code/tree/Dragon-wang-fei-long-Heuristic-VNS

NOT THE END

各位大佬,求点赞,求分享,求关注,嘻嘻嘻~~


版权声明:本文为weixin_43160744原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。