2009-7-26 1
运 筹 学 Operations Research
§ 4.3 割平面法

njx
bAxts
xcz
IP
j
T
,,2,1,,0
..
m a x
:)(
整数给定整数规划对应的松弛线性规划问题为
0
..
m a x
:)(
x
bAxts
xcz
LP
T
2009-7-26 2
运 筹 学 Operations Research
割平面法的基本思想:
.
)(
)(
.)(
)(
)(
)(
.)(
解为止的最优到继续上述步骤,直到得
,,设得最优解为继续求解新的的任一可行解可行解,而不割去的一部分在内的平面,割去含割增加一个约束条件否则,给的最优解;即为为整数解,则若
,设得最优解利用单纯形法求解
IP
xLP
IP
LPx
LP
IPxx
xLP



2009-7-26 3
运 筹 学 Operations Research
割平面的选取:
相应的单纯形表为
,最终得最优基设利用单纯形法求解 ),,,()(...,21 mPPPBLPGOLW




0
1
0
1
,,2,1,
zxrz
mibxbx
n
mj
jj
i
n
mj
jiji?




n
mj
jj
n
mj
jijii
xrzz
mixbbx
1
0
1
0,,2,1,?
2009-7-26 4
运 筹 学 Operations Research
单纯形表为
.)0,,0,0,,,,()( 002010 zbbbxLP Tm,最优值为的最优解为
2009-7-26 5
运 筹 学 Operations Research
.)0,,0,0,,,,()( 020100 Tmi bbbxIPZb 的最优解为,则若的方程为行对应,则单纯形表的第否则,不妨设 kmkZb k )1(0



n
mj
jkjkkk
n
mj
jkjk xbbxbxbx
1
00
1
0
[ ] [ ]
0,[ 0,1 ),
k j k j k j k j k j k j k j
k k j
b b f b b f b
ff


令,其 中 为 的 整 数 部 分,为 的 小 数 部 分,

.][][
)]([)][(
1
0
1
0
1
00
1
0






n
mj
jkjk
n
mj
jkjk
n
mj
jkjkjkk
n
mj
jkjkk
xffxbb
xfbfbxbbx
2009-7-26 6
运 筹 学 Operations Research
柯莫利割平面( Gomory cutting plane),
01
1
1
1
0
1
0
0
0
kn
n
mj
jkj
n
n
mj
jkjk
n
mj
jkjk
fxxf
xxff
xff






为加快计算的速度,常令
}0{m a x 010 imik ff
2009-7-26 7
运 筹 学 Operations Research
割平面的两个性质:
.)(1h 的非整数最优解割平面割去 LPT
.)0,,0,0,,,,( 02010 不满足割平面即可证:仅需证 Tmbbbx
00
1
0

k
n
mj
jkjk fxff?
.)0,,0,0,,,,( 02010 不满足割平面Tmbbbx ▌
2009-7-26 8
运 筹 学 Operations Research
.)(h2 的任一可行解割平面不割去 IPT
.)( 平面即可的任一可行解都满足割证:仅需证 IP
).(?)()?,,?,?(? 21 LPKxIPKxxxx Tn,则?
满足典式又 x?
.][][
)]([)][(
1
0
1
0
1
00
1
0






n
mj
jkjk
n
mj
jkjk
n
mj
jkjkjkk
n
mj
jkjkk
xffxbb
xfbfbxbbx
,?][][?
1
0 Zxbbx
n
mj
jkjkk

,?,?
1
0 Zxff
n
mj
jkjk

2009-7-26 9
运 筹 学 Operations Research
,又 1? 0
0
0?10

k
f
x
n
mj
jkjk fxff
kj
j
0?
1
0

n
mj
jkjk xff ▌
2009-7-26 10
运 筹 学 Operations Research
增添割平面后得新的线性规划



.1,,,2,1,0
..
m a x
01
1
nnjx
fxxf
bAxts
xcz
j
kn
n
mj
jkj
T
增添割平面后的新单纯形表:
对松弛线性规划问题
0
..
m a x
:)(
x
bAxts
xcz
LP
T
2009-7-26 11
运 筹 学 Operations Research
新单纯形表为新表仅是在原表中增加松弛变量的行、列而已!
显然,可继续利用对偶单纯形法求解,
2009-7-26 12
运 筹 学 Operations Research
根据以上讨论,设计算法如下:
2009-7-26 13
运 筹 学 Operations Research
例 1 利用割平面法求解整数规划



整数,0,
43
1..
m a x
:)(
21
21
21
21
xx
xx
xxts
xxz
IP
的松弛线性规划问题为解,)( IP



0,
43
1..
m a x
:)(
21
21
21
21
xx
xx
xxts
xxz
LP
2009-7-26 14
运 筹 学 Operations Research
243 ),( IPPB取初始可行基




4,3,2,1,0
43
1..
m a x
421
321
21
jx
xxx
xxxts
xxz
j
其标准形为
2009-7-26 15
运 筹 学 Operations Research
).,()( 21 PPBLP 的最优基为
100 4
3}
4
3,
4
3m a x { ff
k
0)( 41431310 xfxff作割平面
0)4143(43 43 xx即 434143 543 xxx
2009-7-26 16
运 筹 学 Operations Research
作新单纯形表,并利用对偶单纯形法继续求解:
2009-7-26 17
运 筹 学 Operations Research
.2)1,1()(
.2)0,0,1,1,1()(
,最优值为的最优解为故
,最优值为的最优解为
T
T
IP
LP? ▌
2009-7-26 18
运 筹 学 Operations Research
注:图解法割平面:
1
3)34(
)1(3
33
4
3
4
1
4
3
2
21
21
43
43





x
xx
xx
xx
xx
2009-7-26 19
运 筹 学 Operations Research
Ex.利用割平面法求解整数规划


,整数0,
52
1023..
m a x
:)(
21
2
21
21
xx
x
xxts
xxz
IP
.4)2,2(,最优值为最优解为 T
2009-7-26 20
运 筹 学 Operations Research
§ 4.3 over