Wolsey“强整数规划模型”经典案例之一单源固定费用网络流问题
Wolsey“强整数规划模型”经典案例之一单源固定费用网络流问题
阅读本文可以理解什么是“强”整数规划模型。
单源固定费用网络流问题见文献[1]第13.4.1节(p229-231),是"强整数规划建模“的极好案例。
本文是本博客原创,本博客不转贴他人作品。
单源固定费用网络流问题(The Signle Source Fixed Charge Network Flow Problem)
单源固定费用网络流问题:给定一个有向网络 (边数为m,节点数n),网络上只有一个流量流入节点(标记为1),但有若干个流量流出节点, 节点 i 的流出流量标记为 B[i], 每个边 e (e =1,...,m) 有流量限制(此例中不考虑流量限制,因此流量限制取为-B[1]) 。网络情况如下图:
网络的节点数n=4, 边数m=8。
对应节点1,2,3,4, 流(出)量B={ -6 2 3 4}。
对应边e=1,2,3,4,5,6,7,8, 边由节点对定义,数据是 E={(1 4) (1 2) (4 2) (1 3) (2 3) (3 2) (4 3) (3 4)}.
对应边e=1,2,3,4,5,6,7,8, 边上的一次性费用是C={5 2 7 3 2 4 9 12}.
非“强”整数规划模型
设非负变量 x[e] 是边e (e =1,...,8) 上的流量, 又设0-1变量 y[e] 表示边e上是否有流量。于是模型的目标是极小化一次性费用和,约束无外乎节点上的流量平衡约束和x[e]-y[e]之间的关联约束。
用+Leapms写出模型,其PDF摘录如下:
上述模型用+Leapms中的solve命令求松弛解,可以看到y变量非0-1:
+Leapms>solve
The LP is solved to optimal.
找到线性规划最优解.非零变量值和最优目标值如下:
.........
x1*=
x2*=
x4*=
y1*=0.166667
y2*=0.333333
y4*=0.5
.........
Objective*=
.........
+Leapms>
若要获得整解则必须使用mip命令对问题进行分支定界/割平面求解:
+Leapms>mip
relexed_solution=; number_of_nodes_branched=; memindex=(,)
The Problem is solved to optimal as an MIP.
找到整数规划的最优解.非零变量值和最优目标值如下:
.........
x1* =
x2* =
x5* =
y1* =
y2* =
y5* =
.........
Objective*=
.........
+Leapms>
结果在网络上表现:
“强”整数规划模型(Strong Integer Formulation)
“强”整数规划模型对上述模型采用了Multicommodity改写。方法是引入一个新非负变量 z[e][k], 其含义是流过e边最终贡献给k节点流出的流量(显然此处k=2,...,n, k$\neq 1$)。
新模型的目标不会改变,约束中的节点平衡条件逻辑上也不改变,即对任何节点 i , 流出和流入之差应该为0 或者 当i==k时等于k的流出量。
用+Leapms写出模型,其PDF摘录如下:
上述所谓“强”模型比之前的模型的好处在于在+Leapms中直接使用solve命令就可以求出整数解,即不需要分支定界/割平面过程!
不十分严格地:此即是“强”整数规划模型和含义,强模型更容易求解,或者说对求解器更友好。
+Leapms的求解过程:
+Leapms>solve
The LP is solved to optimal.
找到线性规划最优解.非零变量值和最优目标值如下:
.........
x1*=
x2*=
x5*=
y1*=
y2*=
y5*=
z1_4*=
z2_2*=
z2_3*=
z5_3*=
.........
Objective*=
.........
+Leapms>
结果在网络上表现:
强整数规划模型的详细解释 及 “强”建模原理
本案例的“强”建模原理来源于文[1]中的三个观察(Oberservation 13.2-13.4, page 230) 。
关于“强”模型的详细解释,见文[2]。
其他
文[1]并未在13.4.1中给出此案例的所有数据,其余数据(主要是边上费用数据C是从13.1, page 222)的XPRESS MP模型中读出的。贴在这里,供看官与+Leapms建模语言作对比:
结论
“强”整数规划建模概念是经典建模方法中的重要内容,如果说在本科《运筹学》教学中只要讲述0-1变量的应用即可,那么在研究生层次的《高级运筹学》教学中最好增加“强”整数规划建模内容。
使用本土完全自主知识产权的+Leapms建模语言和+Leapms求解器,可以有效辅助教学,比之传统的舶来建模语言和求解器有优势。
参考文献
[1] Wolsey L A. Integer Programming. New York: Jonh Wiley & Sons, 1998 / ISBN 978-0-471-28366-9
[2] Wolsey L . Strong formulations for mixed integer programming: A survey[J]. Mathematical Programming, 1989, 45(1-3):173-191.
Wolsey“强整数规划模型”经典案例之一单源固定费用网络流问题的更多相关文章
- Wolsey "强整数规划“ 建模的+Leapms实践——无产能批量问题
Wolsey "强整数规划“ 建模的+Leapms实践——无产能批量问题 <整数规划>[1]一书作者L. A. Wolsey对批量问题(Lot-sizing Problem)做了 ...
- Altera OpenCL用于计算机领域的13个经典案例(转)
英文出自:Streamcomputing 转自:http://www.csdn.net/article/2013-10-29/2817319-the-application-areas-opencl- ...
- (zhuan) 资源|TensorFlow初学者必须了解的55个经典案例
资源|TensorFlow初学者必须了解的55个经典案例 2017-05-27 全球人工智能 >>>>>>欢迎投稿:news@top25.cn<<< ...
- javascript的理解及经典案例
js的简介: JavaScript是一种能让你的网页更加生动活泼的程式语言,也是目前网页中设计中最容易学又最方便的语言. 你可以利用JavaScript轻易的做出亲切的欢迎讯息.漂亮的数字钟.有广告效 ...
- jQuery基础的工厂函数以及定时器的经典案例
1. jQuery的基本信息: 1.1 定义: jQuery是JavaScript的程序库之一,它是JavaScript对象和实用函数的封装, 1.2 作用: 许多使用JavaScript能实现的交 ...
- Linux运维之道(大量经典案例、问题分析,运维案头书,红帽推荐)
Linux运维之道(大量经典案例.问题分析,运维案头书,红帽推荐) 丁明一 编 ISBN 978-7-121-21877-4 2014年1月出版 定价:69.00元 448页 16开 编辑推荐 1 ...
- 经典案例:那些让人赞不绝口的创新 HTML5 网站
在过去的10年里,网页设计师使用 Flash.JavaScript 或其他复杂的软件和技术来创建网站.但现在你可以前所未有的快速.轻松地设计或创造互动的.有趣好看的网站.如何创建?答案是 HTML5 ...
- php中foreach()函数与Array数组经典案例讲解
//php中foreach()函数与Array数组经典案例讲解 function getVal($v) { return $v; //可以加任意检查代码,列入要求$v必须是数字,或过滤非法字符串等.} ...
- 阿里云资深DBA专家罗龙九:云数据库十大经典案例分析【转载】
阿里云资深DBA专家罗龙九:云数据库十大经典案例分析 2016-07-21 06:33 本文已获阿里云授权发布,转载具体要求见文末 摘要:本文根据阿里云资深DBA专家罗龙九在首届阿里巴巴在线峰会的&l ...
随机推荐
- 在Eclipse上Maven环境配置使用
1. 安装配置Maven: 1.1 从Apache网站 http://maven.apache.org/ 下载并且解压缩安装Apache Maven. Maven下载地址: http://maven. ...
- CountDownLatch和CyclicBarrier 区别
CountDownLatch : 一个线程(或者多个), 等待另外N个线程完成某个事情之后才能执行. CyclicBarrier : N个线程相互等待,任何一个线程完成之前,所有的线程都 ...
- Eureka的功能特性及相关配置
1.服务提供者1.1服务注册服务提供者启动时,会通过rest请求的方式将自己注册到Eureka Server上,同时带上了自身服务的一些元数据信息.Eureka Server接收到请求后,将元数据信息 ...
- Java工程师修炼之路(校招总结)
Java工程师修炼之路(校招总结) 前言 在下本是跨专业渣考研的985渣硕一枚,经历研究生两年的学习积累,有幸于2019秋季招聘中拿到几个公司的研发岗offer,包括百度,阿里,腾讯,今日头条,网易, ...
- 介绍几款 Python 类型检查工具
近日,微软在 Github 上开源了一个 Python 静态类型检查工具:pyright ,引起了社区内的多方关注. 微软在开源项目上的参与力度是越来越大了,不说收购 Github 这种大的战略野心, ...
- 教程:关于如何通过Maven仓库安装Spire系列的 Java产品
Spire系列库中已发布的Java产品目前有三个,即Spire.PDF for Java.Spire.Presentation for Java.Spire.Barcode for Java.使用该J ...
- JavaScript对象类型判断注意点
注意点 不要使用 new Number() . new Boolean() . new String() 创建包装对象:用 parseInt() 或 parseFloat() 来转换任意类型到numb ...
- Android项目实战(五十四):zxing 生成二维码图片去除白色内边距的解决方案
目录:zxing->encoding->EncodingHandler类 中修改 createQRCode方法 private static final int BLACK = 0xff0 ...
- Android Studio 受不了了
今天遇到一个问题,受不了了,来吐槽一下: Android Studio 里面创建两个工程:app, baselib.一个是App 一个是库.编译不过去. targetSdkVersion 26 imp ...
- Golang随机数生成
需求:10以内随机生成2个幸运数字 package main import ( "fmt" "math/rand" ) func main() { fmt.Pr ...