0-1整数规划与隐枚举法-感受剪枝的魅力
0-1整数规划与隐枚举法-感受剪枝的魅力
整数规划是线性规划的特殊情况,即当约束条件是变量为整数时,线性规划就变成了整数规划。若要求所有变量都为整数,即为纯整数规划;若允许存在一部分变量不一定为整数,则称为混合整数规划。而本文要讨论的0-1整数规划则是纯整数规划的特殊情况,即所有变量要么等于0,要么等于1,故这种变量又成为逻辑变量。
0-1整数规划在生活中还是很常见的,通常可以总结为“是”“否”问题。例如,有n个产品销地x1,...,xn可供选择,为使得利润最大,那么每一个销地都面临是否选择的问题,通常还会有一些限制条件,由于销地xi与销地xj距离较近,所以规定若选择xi就不能选择xj等。那么如何求解0-1规划问题?最朴素的方法是枚举,即将所有销地是否被选择的情况都考虑,那么就是从{0, ... ,0}枚举到{1, ... ,1},需要2的n次方的枚举次数。显然,当n较大时,这种方式的效率就非常低。本文要介绍的隐枚举法就可以提高求解出最优解的效率。
所谓隐枚举法,从字面上理解,就是隐去一些不需要枚举的情况,下面从一个例子出发,来给出隐枚举法的步骤。
【例】求解下列规划问题
max z = 8*x1 + 2*x2 - 4*x3 - 7*x4 - 5*x5; s.t.{ 3*x1 + 3*x2 + x3 + 2*x4 + 3*x5 <= 4 5*x1 + 3*x2 - 2*x3 - x4 + x5 <= 4 xi = 0或1,i = 1, ..., 5 }
1. 预处理
首先需要对原问题进行预处理,至于为什么后文将会解释。预处理的步骤如下:
1) 将目标函数统一为求最小值,即"min", 同时将约束条件都化为">="。
- 若原约束条件为"<=",则不等式左右同乘-1;
- 若原约束条件为"Ai * X = bi",则化为"Ai * X >= bi" 和 "-Ai * X >= -bi",其中Ai为系数行向量,X为变量列向量。
min z' = -8*x1 - 2*x2 + 4*x3 + 7*x4 + 5*x5; s.t.{ -3*x1 - 3*x2 - x3 - 2*x4 - 3*x5 >= 4 -5*x1 - 3*x2 + 2*x3 + x4 - x5 >= 4 xi = 0或1,i = 1, ..., 5 }
2) 将目标函数中系数为负的变量xi化为系数为正的变量xi',其中 xi = 1 - xi’ (若xi = 0 则 xi' = 1; 若xi = 1则xi' = 0)。
故针对本问题,在目标函数中x1和x2前的系数为负,故令x1 = 1 -x1', x2 = 1 - x2',代入1)中化简得
min z' = 8*x1' + 2*x2' + 4*x3 + 7*x4 + 5*x5 - 10; s.t.{ 3*x1' + 3*x2' - x3 - 2*x4 - 3*x5 >= 2 5*x1' + 3*x2' + 2*x3 + x4 - x5 >= 4 xi或xi' = 0或1,i = 1, ..., 5 }
3) 重新排列变量在目标函数和约束条件的先后顺序,使其在目标函数中的系数递增。
min z' = 2*x2' + 4*x3 + 5*x5 + 7*x4 + 8*x1' - 10; s.t.{ 3*x2' - x3 - 3*x5 - 2*x4 + 3*x1' >= 2 (a) 3*x2' + 2*x3 - x5 + x4 + 5*x1' >= 4 (b) xi或xi' = 0或1,i = 1, ..., 5 }
2. 隐枚举
隐枚举的思想是首先枚举找到一个可行解,并得到目标函数值z0,之后的枚举若目标函数值没有z0优,那么就一定不是最优解。
现在说明预处理的作用:
预处理使得目标函数是求最小值,变量的系数都为正且由小到大排列,所以有如下规律:
- 从xi = 0开始枚举是使目标函数最优的,此时得到的函数值也就是最优解的下界;
- 只要按照目标函数中变量的顺序枚举也就是二进制数位从小到大(0...0到1...1)就能尽量较早的枚举出使得目标函数取最小值(最优值)的可行解z0。
- 若解形如0..0xj...x1的目标函数z1取值大于已得到的可行解,那么只要是以xj...x1结尾和只将xj...x1中某些位由0变为1的解的目标函数取值一定大于z1,当然也就大于z0,故一定不会是最优解,可以直接剪枝,即不考虑上述两种情况,直接按照二进制数码顺序枚举其他形式。
隐枚举步骤如下:
(1) 先忽略除" xi或xi' = 0或1 "以外的约束条件,从xi = 0也就是0...0开始枚举。
(2) 计算出枚举出的目标函数值。
- 若小于已有可行解的函数值,或者还无可行解,则执行(3);
- 若大于已有可行解的函数值,剪枝,再进行枚举。
(3) 检查枚举的解是否满足除去的约束条件。(只要检查出一个约束条件不满足就无需再检查)
- 若不满足,则此时的枚举值不是可行解,继续枚举;
- 若满足,则更新可行解和目标函数值z0。可行解0..0xj...x1(前面'0'的个数可能为0),那么只要是以xj...x1结尾和只将xj...x1中某些位由0变为1的解的目标函数取值一定大于z0,剪枝,再进行枚举。
对于本问题,从xi = 0 (i = 1到5)开始枚举,得到z' = -10,所以-10便是最优解的下界(所以10便是原问题的上界)。枚举过程列表如下('-'代表没有判断):
x1' x4 x5 x3 x2' |
z' |
是否(Y/N)满足约束条件 (a) (b) |
是否(Y/N)为可行解 |
---|---|---|---|
0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 1 1 |
-10 -8 -6 -4 |
N - Y N N - Y Y |
N N N Y 剪枝 |
0 0 1 0 0 0 0 1 0 1 0 0 1 1 0 0 1 0 0 0 0 1 1 0 0 |
-5 -3 -1 -3 2 |
N - 函数值-3大于已知可行解的函数值-4,一定不会是最优解 -1 > -4 -3 > -4 2 > -4 |
N 无需判断约束条件,该分支也无需再枚举,即剪枝 同上 同上 同上 |
由表可以看出,我们在第4次枚举得到了一个较优的可行解,其目标函数值z0 = -4,之后的枚举要么是不满足约束条件,要么是函数值大于-4,剪枝。最后我们只枚举了9次就完成了整个过程(比直接枚举的2^5 = 32次快了很多),得到最优解为{0,0,0,1,1},min z' = -4,将其还原成原问题,最优解为{1, 0, 1, 0, 0},max z = 4.
总结:
在解决很多问题的时候,枚举(搜索)似乎是一种直接了当的方式。但是,当解空间较大时,枚举的效率可能就很低,无法达到目的。此时,不妨想想是否在枚举过程中有一些解可以在枚举之前就判断它一定不满足要求,直接不考虑它们(剪枝),这样就可以缩小解空间,提高效率。
- Vuejs和其他前端框架的对比
- 详述 IntelliJ IDEA 中恢复代码的方法
- C# 通过IEnumberable接口和IEnumerator接口实现自定义集合类型foreach功能
- 微信小程序之picker组件
- 详述 IntelliJ IDEA 中恢复代码的方法「进阶篇」
- mac环境下mongodb的安装和使用
- C# 终极基类Object介绍
- EF基础知识小记五(一对多、多对多处理)
- 字符串的方法汇总
- Kotlin和anko融合进行Android开发
- EF基础知识小记六(使用Code First建模自引用关系,常用于系统菜单、文件目录等有层级之分的实体)
- Vue.js简介
- Kotlin之Elvis 操作符
- C# 文件操作系列一
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 使用 Django + Vue.js 开发个人博客网站(完整版附源码)—— Python-课程设计-期末项目
- LVM添加硬盘并扩容至已有分区
- 新增项目时如何修改Jenkins
- 使用封锁协议解决事务的并发问题
- 命令小技巧总结
- Java 中 HashMap 数据结构分析(语言无关)
- Java8中日期的特性,别再用Date来处理日期啦~
- 数据库系统概论考试心得+复习指南
- Docker的基础知识和简单应用
- Spring 注解开发之 @ComponentScan
- 事务并发调度的可串行性及两段锁协议
- 面试不易,控场情况下,感觉少要了1万!
- Java 操作 ElasticSearch
- Docker部署常见应用
- Spring 中实现自动注入的几种方式