MIP启发式求解:局部搜索 (local search)
*本文主要记录和分享学习到的知识,算不上原创。
*参考文献见链接。
本文讲述的是求解MIP问题的启发式算法。
启发式算法的目的在于短时间内获得较优解。
个人认为局部搜索(local search)几乎包括所有的求解MIP的启发式算法的核心框架,从简单的爬山算法(Hill-climbing)到复杂的禁忌搜索(Tabu search),从一个初始解出发的爬山算法(Hill-climbing)到一群初始解出发的遗传算法(Genetic algorithm),其核心框架都是local search。
所以说,学习求解MIP的启发式算法,一定要掌握到local search的思想和算法过程。
目录
引言
local search的过程
local search的要素
引言
求解一个MIP问题的困难之一在于,当问题的解空间很大时,我们想从中获取最优解是很耗时间的。但是,我们可以做到的是,从一个小的解空间短时间获得最优解。
就像一个标准型的线性规划问题,我们并不想逐一判断解空间中的所有解的大小,从而确定出最优解,而是希望有一种更加高效的方法。由于我们知道一个标准型的线性规划问题,其解空间可以表示成一个凸多边形(polyhedron),而最优解只会发生在凸多边形的顶点(vertex)处,所以我们可以仅仅判断顶点上的解的大小,从而获得最优解。
那么如何找到“最有可能成为最优解”的解呢?
我们知道,一个全局最优解一定是局部最优解,而一个解如果都不是局部最优解,那也不可能是全局最优解。所以我们认为那些局部最优解是最有可能成为全局最优解的解。
那么能不能找到所有的局部最优解呢?
只有判断所有的局部最优解的大小,才能确定全局最优解。但是由于解空间很大,想找到所有的局部最优解也是很耗时间的。从时间的角度上看,我们只能获得一部分的局部最优解,找到一个问题的较优解。
那么如何找到局部最优解,实现当前解的迭代呢?
local search寻找局部最优解的思想和过程和“梯度下降法”类似。
local search 的过程
(1)生成初始解:算法从一个初始解或若干个初始解开始;
从一个初始解出发,如:Tabu search,
从若干个初始解出发,如:genetic algorithm
(2)定义邻域,获得候选解:定义解的邻域,并产生若干个候选解;
不同的算法会考虑不同的邻域结构。邻域结构的定义对算法的性能影响比较大。
(3)确定新界:从候选解中确定新解;
最简单的选择策略是将候选解中的最优解作为新解。不过不同的算法会考虑不同的选择策略。
(4)迭代:重复上述搜索过程,直至满足终止条件。期间可能伴随着参数的调整。
终止条件可能是时间、迭代次数等。
local search 的要素
(1)目标函数:用于判断解的优劣;
(2)初始解的产生方法;
初始解的好坏可能会对结果造成很大的影响,与具体算法有关。
(3)邻域的定义;
(4)新解的确定规则;
(5)终止条件。
*local search,包括以local search为核心框架的启发式算法,一般在设计算法时都必须考虑算法的两个方面:
(1)Intensification:深度搜索。尽可能找到局部最优解。
(2)Diversification:广度搜索。尽可能搜索到更多的解空间,以及避免陷入局部最优解。
参考文献
http://www.cnblogs.com/JiePro/p/Metaheuristics_0.html
MIP启发式求解:局部搜索 (local search)的更多相关文章
- 【智能算法】迭代局部搜索(Iterated Local Search, ILS)详解
迭代局部搜索(Iterated Local Search, ILS) 源代码下载请关注微信公众号[程序猿声],在后台回复:[ILS],不包括[]即可下载. 00 目录 局部搜索算法 简单局部搜索 迭代 ...
- 代码 | 自适应大邻域搜索系列之(7) - 局部搜索LocalSearch的代码解析
前言 好了小伙伴们我们又见面了,咳咳没错还是我.不知道你萌接连被这么多篇代码文章刷屏是什么感受,不过,酸爽归酸爽.今天咱们依然讲代码哈~不过今天讲的依然很简单,关于局部搜索LocalSearch的代码 ...
- 对《禁忌搜索(Tabu Search)算法及python实现》的修改
这个算法是在听北大人工智能mooc的时候,老师讲的一种局部搜索算法,可是举得例子不太明白.搜索网页后,发现<禁忌搜索(Tabu Search)算法及python实现>(https://bl ...
- 线搜索(line search)方法
在机器学习中, 通常需要求某个函数的最值(比如最大似然中需要求的似然的最大值). 线搜索(line search)是求得一个函数\(f(x)\)的最值的两种常用迭代方法之一(另外一个是trust re ...
- 集束搜索beam search和贪心搜索greedy search
贪心搜索(greedy search) 贪心搜索最为简单,直接选择每个输出的最大概率,直到出现终结符或最大句子长度. 集束搜索(beam search) 集束搜索可以认为是维特比算法的贪心形式,在维特 ...
- 选择性搜索(Selective Search)
1 概述 本文牵涉的概念是候选区域(Region Proposal ),用于物体检测算法的输入.无论是机器学习算法还是深度学习算法,候选区域都有用武之地. 2 物体检测和物体识别 物体识别是要分辨出图 ...
- java中什么是局部内部类Local inner class?
5.局部内部类Local inner class 马克-to-win:什么叫局部内部类?内部类声明位置:1.它的外部类的范围之内.2.在几个程序块的范围之内.例如,由方法定义的块中或甚至在for循环体 ...
- 采用梯度下降优化器(Gradient Descent optimizer)结合禁忌搜索(Tabu Search)求解矩阵的全部特征值和特征向量
[前言] 对于矩阵(Matrix)的特征值(Eigens)求解,采用数值分析(Number Analysis)的方法有一些,我熟知的是针对实对称矩阵(Real Symmetric Matrix)的特征 ...
- [SAP ABAP开发技术总结]搜索帮助Search Help (F4)
声明:原创作品,转载时请注明文章来自SAP师太技术博客( 博/客/园www.cnblogs.com):www.cnblogs.com/jiangzhengjun,并以超链接形式标明文章原始出处,否则将 ...
随机推荐
- Python相对导入导致SystemError的解决方案(译)
原文出处: http://stackoverflow.com/ 译文出处:yibohu1899 这个问题是如何解决在相对导入的时候,如果出现’System Error’的时候的解决方案.顺带一提, ...
- JavaMailSender怎么发送163和qq邮件
https://blog.csdn.net/Tracycater/article/details/73441010 引入Maven依赖包 <dependency> <groupId& ...
- matlab各向异性扩散滤波
主要是用来平滑图像的,克服了高斯模糊的缺陷,各向异性扩散在平滑图像时是保留图像边缘的(和双边滤波很像). 通常我们有将图像看作矩阵的,看作图的,看作随机过程的,记得过去还有看作力场的. 这次新鲜,将图 ...
- kafka java api生产者
import java.util.HashMap; import java.util.List;import java.util.Map;import java.util.Properties; im ...
- CF1079C Playing Piano
思路: dp. 实现: #include <bits/stdc++.h> using namespace std; ], dp[][]; int main() { int n; while ...
- html5.0学习记录(一)——可拖动视频播放器
最近自己在重新学习html5新特性,了解到有视频标签和拖动标签,于是自己用这两个特性写了一个小demo,主要功能就是可以通过拖动视频来直接播放.效果图如下: 页面使用了<video>标签和 ...
- selenium的定位
id定位 find_element_by_id()方法通过id来定位元素 例如: find_element_by_id("kw") find_element_by_id(&quo ...
- hihoCoder #1162 : 骨牌覆盖问题·三 (矩阵快速幂,DP)
题意:有一个k*n的棋盘,要求用1*2的骨牌来铺满,有多少种方案?(k<8,n<100000001) 思路: 由于k是比较小,但是又不那么小,可以专门构造这样的一个矩阵M,使得只要我们有一 ...
- -[UPAInitViewController startAPPay] in libUPAPayPlugin.a(UPAInitViewController.o)
问题 Undefined symbols for architecture arm64: "_PKPaymentNetworkChinaUnionPay", referenced ...
- Java面试题全集(下)
这部分主要是开源Java EE框架方面的内容,包括hibernate.MyBatis.spring.Spring MVC等,由于Struts 2已经是明日黄花,在这里就不讨论Struts 2的面试题, ...