分枝定界法,割平面法_共51页
6、最大的骄傲于最大的自卑都表示心灵的最软弱无力。——斯宾诺莎 7、自知之明是最难得的知识。——西班牙 8、勇气通往天堂，怯懦通往地狱。——塞内加 9、有时候读书是一种巧妙地避开思考的方法。——赫尔普斯 10、阅读一切好书如同和过去最杰出的人谈话。——笛卡儿
Thank you
33、如果惧怕前面跌宕的山岩，生命 就永远 只能是 死水一 潭。 34、当你眼泪忍不住要流出来的时候 ，睁大 眼睛， 千万别 眨眼!你会看到 世界由 清晰变 模糊的 全过程 ，心会 在你泪 水落下 的那一 刻变得 清澈明 晰。盐 。注定 要融化 的，也 许是用 眼泪的 方式。
35、不要以为自己成功一次就可以了 ，也不 要以为 过去的 光荣可 以被永 远肯定 。
分枝定界法,割平面法_
31、别人笑我太疯癫，我笑他人看不 穿。(名 言网) 32、我不想听失意者的哭泣，抱怨者 的牢骚 ，这是 羊群中 的瘟疫 ，我不 能被它 传染。 我要尽 量避免 绝望， 每天的 成功， 避免以 失败收 常在别 人停滞 不前时 ，我继 续拼搏 。
整数规划算法
⑵
先求（LP3）,如图所示。 此时D 在点取得最优解。
x2
A 3 B C
⑴
(18/11,40/11)
D ⑶
即 x1＝12/5≈2.4, x2 =3,
Z(3)＝-87/5≈-17.4<Z≈-19.8
但x1＝12/5不是整数，可继 续分枝。即 3≤x1≤2。
（三）、整数规划与线性规划的关系
从数学模型上看整数规划似乎是线 性规划的一种特殊形式，求解只需在线 性规划的基础上，通过舍入取整，寻求 满足整数要求的解即可。但实际上两者 却有很大的不同，通过舍入得到的解 （整数）也不一定就是最优解，有时甚 至不能保证所得倒的解是整数可行解。 举例说明。
例：设整数规划问题如下
用 图 解法求出最优解 x1＝3/2, x2 = 10/3 且有Z = 29/6
x2
3
⑴
⑵
(3/2,10/3)
现求整数解（最优解）： 如用“舍入取整法”可得 到4个点即(1，3) (2， 3)(1，4)(2，4)。显然， 它们都不可能是整数规划 的最优解。
3
x1
按整数规划约束条件，其可行解肯定在线性规划问题 的可行域内且为整数点。故整数规划问题的可行解集 是一个有限集，如图所示。
各分枝的目标函数值中，若有小于Z 者，则剪掉此 枝，表明此子问题已经探清，不必再分枝了;否则继续 分枝。
如此反复进行，直到得到Z＝Z*＝Z 为止，即得最优解 X* 。
（二）、例题 例一：用分枝定界法求解整数规划问题（用图解法计算） min Z x1 5 x2 x1 x2 2 5 x 6 x 30 记为（IP） 1 2 x 4 1 x1 , x2 0且全为整数 解：首先去掉整数约束，变成一般线性规划问题 min Z x1 5 x2 x1 x2 2 5 x 6 x 30 记为（LP） 1 2 x 4 1 x1 , x2 0
第五章整数规划
第五章 整数规划主要内容：1、分枝定界法； 2、割平面法； 3、0－1型整数规划； 4、指派问题。
重点与难点：分枝定界法和割平面法的原理、求解方法，0－1型规划模型的建立及求解步骤，用匈牙利法求解指派问题的方法和技巧。
要 求：理解本章内容，熟练掌握求解整数规划的方法和步骤，能够运用这些方法解决实际问题。
§1 问题的提出要求变量取为整数的线性规划问题，称为整数规则问题（简称IP ）。
如果所有的变量都要求为（非负）整数，称之为纯整数规划或全整数规划；如果仅一部分变量要求为整数，称为混合整数规划。
例1 求解下列整数规划问题211020m ax x x z +=⎪⎪⎩⎪⎪⎨⎧≥≤+≤+为整数21212121,0,13522445x x x x x x x x 如果不考虑整数约束，就是一个线性规划问题（称这样的问题为原问题相应的线性规划问题），很容易求得最优解为：96m ax ,0,8.421===z x x 。
用图解法将结果表示于图中画“+”号的点都是可行的整数解，为满足要求，将等值线向原点方向移动，当第一次遇到“+”号点（1,421==x x ）时得最优解为1,421==x x ，最优值为z=90。
由上例可看出，用枚举法是容易想到的，但常常得到最优解比较困难，尤其是遇到变量的取值更多时，就更困难了。
下面介绍几种常用解法。
§2 分枝定界法分枝定界法可用于解纯整数或混合的整数规划问题。
基本思路：设有最大化的整数规划问题A ，与之相应的线性规划问题B ，从解B 开始，若其最优解不符合A 的整数条件，那么B 的最优值必是A 的最优值*z的上界，记为z ；而A 的任意可行解的目标函数值是*z的一个下界z ，采取将B 的可行域分枝的方法，逐步减少z 和增大z ，最终求得*z 。
现举例说明： 例2 求解A219040m ax x x z +=⎪⎪⎩⎪⎪⎨⎧≥≤+≤+为整数21212121,0,702075679x x x x x x x x 解：先不考虑条件⑤，即解相应的线性规划B （①--④），得最优解=1x 4.81, =2x 1.82,①② ③ ④ ⑤=0z 356（见下图）。
分支定界法和割平面法
分支定界法和割平面法在上学期课程中学习的线性规划问题中，有些最优解可能是分数或消失，但现实中某些具体的问题，常要求最优解必须是整数，这样就有了对于整数规划的研究。
整数规划有以下几种分类：（1）如果整数规划中所有的变量都限制为（非负）整数，就称为纯整数规划或全整数规划；（2）如果仅一部分变量限制为整数，则称为混合整数规划；（3）整数规划还有一种特殊情形是0-1规划，他的变量取值仅限于0或1。
本文就适用于纯整数线性规划和混合整数线性规划求解的分支定界法和割平面法，做相应的介绍。
一、分支定界法在求解整数规划是，如果可行域是有界的，首先容易想到的方法就是穷举变量的所有可行的整数组合，然后比较它们的目标函数值以定出最优解。
对于小型问题，变量数量很少，可行的整数组合数也是很小时，这个方法是可行的，也是有效的。
而对于大型的问题，可行的整数组合数很大时，这种方法就不可取了。
所以我们的方法一般是仅检查可行的整数组合的一部分，就能定出最有的整数解。
分支定界法就是其中一个。
分枝定界法可用于解纯整数或混合的整数规划问题。
在二十世纪六十年代初由Land Doig 和Dakin 等人提出。
由于这方法灵活且便于用计算机求解，所以现在它已是解整数规划的重要方法。
目前已成功地应用于求解生产进度问题、旅行推销员问题、工厂选址问题、背包问题及分配问题等。
设有最大化的整数规划问题A ，与它相应的线性规划为问题B ，从解问题B 开始，若其最优解不符合A 的整数条件，那么B 的最优目标函数必是A 的最优目标函数z *的上界，记作z ；而A 的任意可行解的目标函数值将是z *的一个下界z 。
分枝定界法就是将B 的可行域分成子区域再求其最大值的方法。
逐步减小z 和增大z ，最终求到z *。
现用下例来说明：例1 求解下述整数规划 219040Maxx x z +=⎪⎩⎪⎨⎧≥≥+≤+且为整数0,702075679212121x x x x x x解 （1）先不考虑整数限制，即解相应的线性规划B ，得最优解为：124.81, 1.82,356x x z ===可见它不符合整数条件。
求解整数规划常用的方法有分枝定界法和割平面法
寻找割平面方程
（1）由单纯形最终表得到决策变量非整数解方程，设 1 x1 ai k xk bi 为
k
其中bi是基变量的非整数解。 （2）将aik和bi分解为整数N和正真分数f 两部分之和
a ik N ik f ik , bi N ni f bi
2
将（2）代入(1)中,然后将整数置于方程左边，分 数置于方程右变，即
xi
N
k
ik x k
N bi f bi
f
k
ik xk
0
(3)得割平面方程
f bi
f
k
ik x k
0
3
整数线性规划模型的求解——分枝定界法
基本思想 通过分枝枚举来寻找最优解。首先不考虑对变 量的整数要求，求解相应的线性规划模型，如求得 最优解不符合整数要求，则把原模型分解为两部分， 每一部分都增加新的约束条件以减少相应线性规划 模型的可行域。通过不断分解，逐步逼近满足要求 的整数最优解，在这个过程中包括了“分枝”和 “定界”两个关键步骤。
1 1 1 0
利用这一性质，可以使原系数矩阵(cij)变换成含有
很多0元素的新系数矩阵
11 c ij ，而最优解保持不变。
匈牙利法是针对目标要求极小化问题提出的 基本原理：为了实现目标极小，在系数矩阵 元素cij≥0条件下，如果能使矩阵具有一组处于 不同行不同列的零元素cij’＝0，画上圈符号 “◎”，表示对应该元素的决策变量xij=1，未画 圈元素对应的决策变量xij＝0，那么目标的数值 z’＝0为最小，这样的组合解x就是最优解。所以 匈牙利法又称画圈法。 画圈法的关键是如何实现系数矩阵具有一组 处于不同行又不同列的0元素（独立零），并保 证所画的圈的个数等于矩阵的阶数。
整数规划 割平面法 分枝定界法
割平面法的关键在于如何确定切割方程，使之能对可行域进行 真正的切割，而且切去部分不含有整数解点。
下面讨论切割方程的求法。 设与整数规划相对应的线性规划最优解中基变量XBi=（B-1b)i不 是整数，将最优单纯形表中该基变量对应的行还原成约束方程，即 XBi +ΣaijXj=（B-1b)i ⑴ 将（B-1b)i，aij都分解成整数与非负真分数之和的形式，即 （B-1b)i=Ni+fi 其中0＜ fi ＜1 ⑵ aij=Nij+fij 其中0≤ fij ＜1 ⑶ 这里Ni、Nij是整数，将⑵、 ⑶代入⑴，得 XBi +Σ（Nij+fij）Xj=Ni+fi 即 XBi +ΣNijXj-Ni=fi-ΣfijXj ⑷ 当诸Xi都是整数时， ⑷式左端是整数，所以右端亦应是整数，但右 端是两个正数之差，且∵0＜ fi ＜1，∴ fi-ΣfijXj是小于1的整数，从
9x1+ 7x2=56 Z=40x1+90x2 D1
4
7x1+20x2=70
D2
6
10
x1
求解线性规划L1、L2 得最优解为： 问题L1： L+ x1≤4 x1=4.00 x2=2.10 Z1=349 问题L2： L+ x1≥5 x1=5.00 x2=1.57 Z2=341
因为没有得到整数解，所以继续对L1进行分解，增加约束： x2≤2，x2≥3将L1分解成问题L3与L4，并求得最优解如下： 问题L3： L1+ x2≤2 问题L4：L1+x2≥3
例2 求解下面整数规划
x2 8
maxZ=40x1+90x2 ⑴ 9x1+ 7x2≤56 ⑵ 7x1+20x2≤70 ⑶ 4 x1,x2≥0 ⑷ x1,x2 整数 ⑸ 解：先不考虑条件⑸，求解相 0 应的线性规划问题L，得最优解 x1=4.81，x2=1.82，Z0=356（见图） 该解不是整数解。选择其中一个 非整数变量，如x1=4.81，对问题 L分别增加约束条件： x1≤4，x1≥5 将问题L分解为两个子问题L1，L2 （分枝），也就是去掉问题L不含 整数解的一部分可行域，将原可 行域D变为D1、D2两部分（如图）。
分枝定界法
束——缩小可行域;将原整数规划问题分枝——分为两个子 规划,再解子规划的伴随规划……通过求解一系列子规划的 伴随规划及不断地定界 .最后得到原整数规划问题的整数最 优解 . 下面结合一个极大化例题来介绍分枝定界法的主要思路 .
例2 某公司计划建筑两种类型的宿舍.甲种每幢占地0.25 ×103m2, 乙种每幢地0.4×103m2.该公司拥有土地3×103m2. 计划甲种宿舍不超过 8 幢,乙种宿舍不超过4幢.甲种宿舍每 幢利润为10万元，乙种宿舍利润为每幢20万元.问该公司应
x2 3 x1, x2 0
问题 B4
max f 20 x1 10 x2
5x1 8x2 60
x1 8
s.t
x2 4 x1 6
x2 4 x1, x2 0
它们的可行域分别为 K3, K4 ( ). 见图3。
第21页/共34页
x2
因为 K4 ，问题 B4
4
无可行解，此问题已
3
作出问题 A1, 的A2伴随规划 B则1, 问B2题, 的可B1行, B2, 域为 K1, K见2图, 2(b). 以下我们将由同一问题分解出的两
个分枝问题称为"一对分枝".
第15页/共34页
x2
4
x2
3
2 1
O
246
8 x1
O
12 4
6
8
x1
(a)
(b)
图2 ( a )
4. 分别求解一对分枝
在一般情况下,对某个分枝问题(伴随规划)求解时,可能出现 以下几种可能:
x1, x2 0, 整数
（1）
第3页/共34页
若暂且不考虑 x1, x取2 整数这一条件.则(1)就变为下列 线性规划 :
4.3-分枝定界法和割平面法
剪枝 x 3 再分枝： 2
不是问题A解 而z (12 ) z
B1 : x1 4.00 x2 2.10 z (1) 349 z
z (12 ) 327
4
x2 2
B11 : x1 4.00 x2 2.00 z (11) 340
定界： z 340 z 341
z 340 z 341
分支定界的全过程： B2 : x1 5.00 B1 : x1 4.0 0
x 2 2 .1 0
x2 1.57
x2 2
z
(1)
B11 : x1 4.00 x2 2.00 z
(11)
xx 34 2 最优解： 1
3 49
z ( 2 ) 341
（3）求解
点
求解过程如表4－6所示。
过滤条件 约束 ④ × √ √ √ 4x1＋3x2＋2x3≥5 × √ × √ √ √ √ × √ √ 5 √ 2 ① ② ③ z值
4x1＋3x2＋2x3≥2
（0，0，0）T （0，0，1）T （0，1，0）T （0，1，1）T （1，0，0）T （1，0，1）T
§4 分枝定界法
第二步:定界
记A的目标函数最优值为z*,以z(0)作为z* 的上界,记为 z =z(0).再用观察法找的一个整数可 行解X′,并以其相应的目标函数值z′作为z*的下 界,记为z＝z′,也可以令z＝－∞,则有: *
zz z
§3 分枝定界法
第三步:分枝
在以上界 z 所对应的解 X (b1,, br ,, bm ,0,,0)T 中,任选一个不符合整数条件的变量,例如 br(不 为整数),以 [br ]表示不超过 br 的最大整数.构造 两个约束条件
分枝定界法
§3 分枝定界法 （2）
x2 10 1)不受整数限制，作为 9 8 普通线性规划求解，可 得出最优解为：x1=10/3， 7 6 x2=4/3，z=26/3（ 见 图 5 2-1）。 4 3 该解示如图2-4中的节 2 点①。 1
[解]
2x1+x2=8 最优解 x1+2x2=6 1 2 3 4 5 6 7 8 9 图2-1 x1
x1+2x2≥6
x1≤3 x2≥2 x1≥0，全取整。 其图解法见图2-5，最优解为x1=2，x2=2，z=10。
§3 分枝定界法 （9）
x1≤3 1 z =26/3 x1=10/3 x1≥4 x2=4/3 图2-4
2
z =9 x1=3 x2=3/2 3 不可行
x2
8 7 6 5 4 3 x2=2
目标函数 1 2 3 4 5 6 7 8 9 x1 图2-5
§3 分枝定界法 （10）
从对求解[例2-5]的过程，可归纳出求解整数规划的分枝 定界法有下述特点： (i) 既可求解全整数规划，亦可求解混合整数规划。 (ii)求解每个子集的最优整数解，都是首先放弃整数约束， 用线性规划解法求出无约束时的最优解，此时的目标函 数值即为该子集所有可行解的目标下界（对于求极小值 的规划而言）。
§3 分枝定界法 （3）
2）因为x1、x2当前均为非整数，故不满足整数要求，任 选1个进取分枝。设选x1进行分枝，把可行集分成2个子 集： x1≤[10/3]=3及x1≥[10/3]+1=4 3)x1≤3时 目标函数 min z=x1+4x2 约束条件 2x1+x2≤8
x1+2x2≥6
x1≤3 x1，x2≥0且为整数。
(iii)如果子集的非整数最优解的下界超过迄今已得到的最 好可行整数解目标值，或者子集无解，则这个子集将被 剪掉，又称剪枝。
算法讲稿5分枝定界法
Q.Delete(Ew); // 取下一扩展结 点
15
四、构造最优解
为了在算法结束后能方便 地构造出与最优值相应的 最优解，算法必须存储相 应子集树中从活结点到根 结点的路径。为此目的， 可在每个结点处设置指向 其父结点的指针，并设置 左、右儿子标志。
找到最优值后，可以根据 parent回溯到根节点，找到 最优解。
for (int i = 0; i <= n+1; i++) grid[i][0] = grid[i][m+1] = 1; // 左翼和右翼
for (int i = 0; i < NumOfNbrs; i++) {
}
}
11
6.3 装载问题
一、问题描述 二、队列式分支限界法 三、算法的改进 四、构造最优解 五、优先队列式分支限界法
12
一、问题描述
有 的一轮批船共，个其集 中装 集箱 装要 箱装i的上重2量艘为载w重i，量且分∑别wi为≤CC1+1和C2 C2 装载问题要求确定是否有一个合理的装载方案可将
优先队列式分支限界法：
[A] B, C => B(45), C(0) [B, C] D, E => E(45) [E, C] J, K => K(45) [1, 0, 0] [C] F, G => F(25), G(0) [F, G] L, M => L(50), [0, 1, 1] M(25) [G] N, O => N(25), O(0)
在优先队列式分支限界法中，一旦有一个叶结 点成为当前扩展结点，则可以断言该叶结点所 相应的解即为最优解。此时可终止算法。
17
6.4 布线问题
运筹学_分支定界法课件
m ax Z x1 5 x2
m ax Z x1 5 x2
x1 x2 2
( IP 1)
5
x1 x1
6
x2
30 4
x1
1
x 1 , x 2 0 且 为 整 数
x1 x2 2
( IP
2)
5
x1 x1
6
x2
30 4
x1
2
x 1 , x 2 0 且 为 整 数
现在只要求出（LP1）和（LP2）的最优解即可。
x1
x1 5x2 Z
运筹学_分支定界法
x1＝18/11, x2 =40/11
Z(0) =218/11≈(19.8)
x2
即Z 也是（IP）最大值的上限。
m ax Z x1 5 x2
3
x1 x2 2
5
x1 x1
6
x2
30 4
x 1 , x 2 0
⑵ 5x16x2 30 ⑴ x1 x2 2
＃
运筹学_分支定界法
例一：用分枝定界法求解整数规划问题（用图解法计算）
m ax Z x1 5 x2
x1 x2 2
5
x1
6
x2
30
x1
4
x 1 , x 2 0 且 全 为 整 数
记为（IP）
解：首先去掉整数约束，变成一般线性规划问题
m ax Z x1 5 x2
x1 x2 2
5
x1 x1
6
x2
30 4
3
x1
1
x 1 , x 2 0 且 为 整 数
1
⑵ 5x16x2 30 ⑴ x1 x2 2
A
⑶ x1 4
1