运筹(第四章整数规划)
运筹学-整数规划与分配问题PPT
但 z=13 不是最优。实际问题的
最优解为(4 , 1)这时 z*= 14。
逻辑(0-1)变量在建立数学模型中的作用
1. m 个约束条件中只有 k 个起作用
设 m 个约束条件可以表示为:
n
aijxj bi (i1, ,m)
j1
定义逻辑变量
1,假定第 i 个约束条件不起作用 yi 0,假定第 i 个约束条件起作用
第四章 整数规划与分配问题
整数规划的特点及作用 分配问题与匈牙利法 分枝定界法 割平面法 应用举例
1 整数规划的特点及应用
在实际问题中,全部或部分变量取值必须是整数。比如人 或机器是不可分割的,选择地点可以设置逻辑变量等。
在一个线性规划问题中要求全部变量取整数值的,称纯整
数线性规划或简称纯整数规划;只要求一部分变量取整 数值的,称为混合整数规划。
如果完成任务的效率表现为资源消耗,考虑的是如何分配 任务使得目标极小化;如果完成任务的效率表现为生产效 率的高低,则考虑的是如何分配使得目标函数极大化。
在分配问题中,利用不同资源完成不同计划活动的效率常
用表格形式表示为效率表,表格中数字组成效率矩阵。
例2. 有一份说明书,要分别翻译成英、日、德、俄 四种文字,交甲、乙、丙、丁四个人去完成。因各人专长 不同,使这四个人分别完成四项任务总的时间为最小。效 率表如下:
又设 M 为任意大的正数,则约束条件可以改写为:
n
aijxj
bi Myi
j1
y1 y2 ym mk
2. 约束条件的右端项可能是 r 个值中的某一个
n
即
aijxj b1或b2或或br
j1
定义逻辑变量:
yi 10, ,假 其定 它约束右端项b为 i
运筹第四章整数规划与分配问题
i=1,2
则问题可以表示为
4 用以表示含固定费用的函数 总费用
K j + c j x j ( x j > 0) Cj(xj ) = ( x j = 0) 0
则上述条件可以表示成
r n ∑ aij x j ≤ ∑ b; y + ... + y = 1 m 2 1
3、 两组条件中满足其中的一组 、
若 x1 ≤ 4, 则 x2 ≥ 1
若 x1 > 4, 则 x2 ≤ 3
定义
1 第i组条件不起作用 yi = 0 第i 组 条件 起作 用
0 0 X = 1 0 0 0 1 1 0 0 0 0 0 0 1 0
用矩阵形式表示为: 用矩阵形式表示为: 解矩阵
一般分配问题 设有n项任务 需有n个人去完成 项任务, 个人去完成, 设有 项任务,需有 个人去完成,每个人只能完成一 项任务,每项任务只能由一个人去完成,设第i人完成 项任务,每项任务只能由一个人去完成,设第 人完成 项任务需要的时间是a 第j 项任务需要的时间是 ij , 问如何分配才能使完成任 务的总时间最少? 务的总时间最少? 设
2. 整数规划问题的特征与性质
特征—变 特征 变量整数性要求 来源 问题本身的要求 引入的逻辑变量的需要 性质—可 性质—可行域是离散集合
3. 整数规划的分类
纯整数规划 要求全部决策变量的取值都为整数, 要求全部决策变量的取值都为整数 则称为纯整数规划 (All IP); ; 混合整数规划 仅要求部分决策变量的取值为整数,则称为混合整数规 仅要求部分决策变量的取值为整数, 划(Mixed IP); ; 0-1整数规划 整数规划 要求决策变量只能取0或 值 则称为0-1规划 规划(0-1 要求决策变量只能取 或1值,则称为 规划 Programming)。 。
运筹学基础及应用第4章-整数规划与分配问题
整数规划的特点及应用
解:对每个投资项目都有被选择和不被选择两种可能,因此 分别用0和1表示,令xj表示第j个项目的决策选择,记为:
j投 资 1 对 项 目 xj ( j 1,2,..., n) j不 投 资 0 对 项 目
投资问题可以表示为:
max z
c
j 1
n
j
xj
n a j x j B j 1 x2 x1 s .t x 3 x4 1 x5 x6 x7 2 ) x j 0或者1 (j 1, 2, L n
B1 B2 B3 B4 年生产能力
A1
A2 A3 A4 年需求量
2
8 7 4 350
9
3 6 5 400
3
5 1 2 300
4
7 2 5 150
400
600 200 200
工厂A3或A4开工后,每年的生产费用估计分别为1200万或1500万元。 现要决定应该建设工厂A3还是A4,才能使今后每年的总费用最少。
0-1型整数线性规划:决策变量只能取值0或1的整数线性 规划。
整数规划的特点及应用
整数规划的典型例子
例4.1 工厂A1和A2生产某种物资。由于该种物资供不应求,故需要 再建一家工厂。相应的建厂方案有A3和A4两个。这种物资的需求地 有B1,B2,B3,B4四个。各工厂年生产能力、各地年需求量、各厂至各 需求地的单位物资运费cij,见下表:
例4.3 设整数规划问题如下
max Z x1 x 2 14x1 9 x 2 51 6 x1 3 x 2 1 x , x 0且 为 整 数 1 2
首先不考虑整数约束,得到线性规划问题(一般称为松弛问 题)。
运筹学基础及应用_(第四章_整数规划与分配问题)
(d) 8
(e)1号、
4号、6号、9号开采时不能超过两个,试表示上
述约束条件。
Next
基础教研室
(a)当x8=1 当x8=0 ∴ x8 x6
x6=1,x6≠0 x6=1,x6=0
(b)当x5 =1 当x5 =0 ∴ x5 + x3 1
x3=0, x3 ≠1 x3=0, x3 =1
基础教研室
【例1】求下述整数规划的最优解
Max z= 3x1 + 2x2 st . 2x1 + 3x2 14 x1 + 0.5x2 4.5 x10,x20,且为整数
基础教研室
x2 x1+0.5x2=4.5
4
(3.25, 2.5) 2 2x1+3x2=14
2
4
6
x1
3x1+2x2=6
二、整数规划的求解方法
1 -选择电网供应 设 y1 0 -不选择电网供应
10 d j x j f (1 y1 ) M j 1 10 0.3d j x j p (1 y2 ) M j 1 y1 y2 1 y1 , y2 0或1
基础教研室
【例3】投资决策问题 某公司准备1000万元资金在10个地点中选择若干个建立 工厂(工厂名称用地点名来命名),有关数据如下:
由于各个工厂之间有配套和协作关系,因此必须满足条件: 1、 建工厂1就必须同时建工厂2; 2、 若建工厂2就不允许建工厂3; 3、 工厂4和工厂5至少建一个; 4、 工厂6,7,8恰好建2个; 5、 工厂8,9,10最多建2个; 6、 建工厂4或者建工厂6,就不能建工厂8,反过来也一样; 7、 条件2,3,5最多满足2个。 问选择哪几个地点建厂最有利? Next
运筹学 第4章 整数规划与分配问题
匈牙利法思路:若能在 [Cij] 中找出 n 个位于
不同行不同列的0元素(称为独立0元素),则
令解矩阵[xij]中对应这n个独立0元素的元素
取值为 1 ,其他元素取值为 0 ,则它对应目
标函数zb=0是最小的。这就是以[Cij]为系数
矩阵分配问题的最优解,也得原问题的最
优解。
定理1 若从分配问题效率矩阵[cij]的每一行元素中分别减去 (或加上)一个常数ui(称为该行的位势),从每一列分别减去 (或加上)一个常数vj(称为该列的位势),得到一个新效率矩阵 [bij],若其中bij=cij-ui-vj,则[bij]的最优解等价于[cij]的最优解
第1步:找出效率矩阵每行的最小元素,并分别从每行
中减去。
第2步:再找出矩阵每列的最小元素,并分别从各列中 减去。
2 10 9 7 2 15 4 14 8 4 13 14 16 11 11 4 15 13 9 4
0 8 7 5 11 0 10 4 0 3 5 0 0 11 9 5
表明m个约束条件中有(m-k)个的右端项为( bi+M ),不起约 束作用,因而,只有k个约束条件起作用。 ② 约束条件的右端项可能是r个值b1 , b2 ,, br 中的某一个 即: 定义:
n
aij x j b1 或b2或或br
j 1
1 假定约束右端项为 bi yi 否则 0
现用下例来说明: max z=40x1+90x2 9x1+7x2≤56 7x1+20x2≤70 x1,x2≥0 x1,x2整数 ① ② ③ ④ ⑤
解:先不考虑条件⑤,即解相应的线性规划B,①~④(见图5-2), 得最优解x1=4.81,x2=1.82,z0=356
运筹学 整数规划( Integer Programming )
检查所有分枝的解及目标函数值,若某分枝的解是整数并且目标函数 值大于(max)等于其它分枝的目标值,则将其它分枝剪去不再计算,若 还存在非整数解并且目标值大于(max)整数解的目标值,需要继续分枝, 再检查,直到得到最优解。
割平面法的内涵:
Page 18
通过找适当的割平面,使得切割后最终得到这样的可行域( 不一定一次性得到), 它的一个有整数坐标的顶点恰好是 问题的最优解.
-Gomory割平面法
例: 求解
max z x1 x2 s.t. x1 x2 1
3x1 x2 4 x1 , x2 0, 整 数
1 x1 3/4 1 0 -1/4 1/4 0
1 x2 7/4 0 1 3/4 1/4 0
0 x5 -3 0 0 -3 -1 1
0 0 -1/2 -1/2 0
由对偶单纯形法, x5为换出变量, x3为换入变量, 得Page 29
cj CB XB b 1 x1 1 1 x2 1 0 x3 1
1 100 0 x1 x2 x3 x4 x5 1 0 0 1/3 1/12 0 1 0 0 1/4 0 0 1 -1 -1/3 0 0 0 -1/2 -1/6
收敛性很慢. 但若下其它方法(如分枝定界法)配合使用,
也是有效的.
分支定界法
Page 33
分支定界法的解题步骤:
1)求整数规划的松弛问题最优解; 若松弛问题的最优解满足整数要求,得到整数规划的最优解,否则转下
一步; 2)分支与定界:
任意选一个非整数解的变量xi,在松弛问题中加上约束: xi≤[xi] 和 xi≥[xi]+1
运筹学习题解答(chap4 整数规划与分配问题)
第四章 整数规划与分配问题一、建立下列问题的数学模型1、P143, 4.1 利用0-1变量对下列各题分别表示成一般线性约束条件 (a) 221≤+x x 或53221≥+x x ; (b) x 取值0,3,5,7中的一个; (c) 变量x 或等于0,或50≥; (d) 若21≤x ,则12≥x ,否则42≤x ; (e) 以下四个约束条件中至少满足两个:6225433121≥+≥≤≤+x x x x x x ,,,。
解:(a) 设⎩⎨⎧=否则。
,个条件起作用;第1i ,0y i (i=1,2),M 为任意大正数。
则有 ⎪⎩⎪⎨⎧=+≥++≤+1y y My -5x 3x 2My 2x x 21221121(b) 设⎩⎨⎧=≠=ix i x y i ,1,0,7,5,3,0=i ,则原条件可表示为⎩⎨⎧=++++++=1753075307530y y y y y y y y x(c) 设⎩⎨⎧≥==50,10,0x x y ,则原条件可表示为⎪⎩⎪⎨⎧≥--≥≤0)1(50x M y x yM x(d)⎩⎨⎧=否则。
,组条件起作用;第1i ,0y i (i=1,2),M 为任意大正数。
则有⎪⎪⎪⎩⎪⎪⎪⎨⎧=++≤->-≥+≤.1,4,2,1,22122211211y y My x My x My x My x (e)设⎩⎨⎧=个条件不成立第个条件成立第i ,1i ,0y i ,4,3,2,1i =,则原条件可表示为:⎪⎪⎪⎩⎪⎪⎪⎨⎧≤+++-≥+-≥+≤+≤+2y y y y My 6x x My 2x M y 2x M y 5x x 43214433321121 2、P143, 4.2 某钻井队要从以下10个可供选择的井位确定5个钻井探油,目的是使得总的钻探费用最小。
若10个井位代号为101S ,...,S ,相应的钻探费用为101C ,...,C ,并且井位的选择要满足下列条件:(1)或选择1S 和7S ,或选择8S ;(2)选择了3S 或4S 就不能选择5S ,反过来也一样; (3)在10962S ,S ,S ,S 中最多只能选两个。
第四章_整数规划整数规划数学模型运筹学基础及其应用胡运权第五版
表4-3
人员 工作 译成英文
译成日文
甲 2
15
乙 10
4
丙 9
14
丁 7
8
译成德文
译成俄文
13
4
14
15
16
13
11
9
【解】此工作分配问题可以采用枚举法求解,即将所有分配方案 求出,总分最大的方案就是最优解。本例的方案有 4!=4×3×2×1=24种,当人数和工作数较多时,方案数是人数 的阶乘,计算量非常大。用0-1规划模型求解此类分配问题显得非 常简单。
【例4.1 】某人有一背包可以装10公斤重、0.025m3的物品。他准备 用来装甲、乙两种物品,每件物品的重量、体积和价值如表5-1所 示。问两种物品各装多少件,所装物品的总价值最大? 表4—1 重量 体积 价值 物品 (公斤/每件) (m3/每件) (元/每件)
甲 乙
1.2 0.8
0.002 0.0025
max Z 4 x1 3 x 2 (a) 1.2 x1 0.8 x 2 10+My2 1.8 x 0.6 x 12 My (b) 1 2 1 (c ) 2 x1 2.5 x 2 25 My2 (d ) 1.5 x1 2 x 2 20 My1 y y 1 2 1 x1 , x 2 0, 且均取整数, y 0或1
要求每人做一项工作,约束条件为
x11 x12 x13 x14 1 x x x x 1 21 22 23 24 x31 x32 x33 x34 1 x 41 x 42 x 43 x 44 1
每项工作只能安排一人,约束条件为
【例4.3 】企业计划生产2000件某种产品,该种产品可利用A、B、 C设备中的任意一种加工。已知每种设备的生产准备结束费用、生 产该产品时的单件成本以及每种设备限定的最大加工数量(件)如 下表所示,试建立总成本最小的数学模型。 设备 生产准备结束 费(元) 生产成本 (元/件) 限定最大加工 数(件)
运筹整数规划素材
OPERATIONS RESEARCH
2024/10/20
1
第四章 整数规划与分配问题
整数规划的有关概念及特点 整数规划的应用 指派问题及匈牙利解法 整数规划的求解方法:分枝定界法、割平面法
2024/10/20
2
§1 整数规划的有关概念及特点
§1.1 概念
整数规划: 要求决策变量取整数值的规划问题。 (线性整数规划、非线性整数规划等)
解:设 xij 表示学生i在周j的值班时间。
0, 学生i在周j不值班 yij 1, 学生i在周j值班
a表ij 示学生i在周j的最多可值班时间。
65
则目标函数:
min z
ci x ij
i1 j1
2024/10/20
12
约 (1) 束 条 (2) 件
(3)
(4)
(5)
(6)
6
xij 14,
15
x6 400
x5 x6 850
x4 x5 x6 1750
x3 x4 x5 x6 2450 x2 x3 x4 x5 x6 3000
x1 x2 x3 x4 x5 x6 3500 x j 0, j 1,...6
2024/10/20
16
例3(固定成本问题) 高压容器公司制造小、中、大三种尺寸的金属容器, 所用资源为金属板、劳动力和机器设备,制造一个 容器所需的各种资源的数量如表所示。每种容器售 出一只所得的利润分别为 4万元、5万元、6万元, 可使用的金属板有500吨,劳动力有300人/月,机 器有100台/月,此外不管每种容器制造的数量是多 少,都要支付一笔固定的费用:小号是l00万元, 中号为 150 万元,大号为200万元。现在要制定一 个生产计划,使获得的利润为最大。
运筹学课件第四节0-1型整数规划
目录
CONTENTS
• 0-1型整数规划概述 • 0-1型整数规划的数学模型 • 0-1型整数规划的求解算法 • 0-1型整数规划的案例分析 • 0-1型整数规划的软件实现
01 0-1型整数规划概述
CHAPTER
定义与特点
定义
0-1型整数规划是一种特殊的整数规 划,其中决策变量只能取0或1。
解决方案通常采用动态规划或混合整数线性规 划方法,通过迭代和优化算法来找到最优解。
05 0-1型整数规划的软件实现
CHAPTER
Excel求解工具
适用范围
适用于简单的0-1型整数规划问题。
优点
操作简单,易学易用,适合初学者。
使用方法
利用Excel的Solver插件,设置目标函数、 约束条件和决策变量,进行求解。
其他约束
除了资源和需求约束外,还可能 存在其他类型的约束,如数量约 束、时间约束等,这些约束条件 都对决策变量的取值范围进行了 限制。
决策变量
离散变量 0-1型整数规划中的决策变量通常 是离散的,只能取0或1两个值。 这些决策变量代表了不同的策略 或选择。
最优解 最优解是指在所有可行解中使目 标函数达到最优值的决策变量的 取值组合。
缺点
对于大规模问题求解能力有限,可能存在精 度问题。
Python求解库
适用范围
适用于各种规模的0-1型整数规 划问题。
使用方法
利用Python的优化库,如PuLP 或CVXPY,编写目标函数和约束 条件,进行求解。
优点
功能强大,可处理大规模问题 ,精度高。
缺点
需要一定的编程基础,学习成 本较高。
MATLAB求解工具
运筹学chap.5整数规划
利用已知的数据,年运输成本为: TC1=5x11+2x12+3x13+4x21+3x22+4x23+9x31+7x32 +5x33+10x41+4x42+2x43+8x51+4x52+3x53
14
建新工厂的年固定成本为:
TC2=175y1+300y2+375y3+500y4; 总成本为:TC=TC1+TC2; 生产能力的约束条件为: 从新工厂B1运到A1,A2,A3三个城市销售中心的总量应小于 等于B1的生产能力,所以约束条件为: x11+x12+x13≤10y1 B1的生产能力; 同理可得: x21+x22+x23≤20y2 B2的生产能力; x31+x32+x33≤30y3 B3的生产能力; x41+x42+x43≤40y4 B4的生产能力; x51+x52+x53≤30 B5的生产能力(已建);
17
例3:投资决策问题
某部门在今后五年中可用于投资的资金总额为
B 万 元 , 有 n(n 2)个 可 以 考 虑 的 投 资 项 目 , 假
定 每 个 项 目 最 多 投 资 一 次 ,第 j ( j 1 , , n ) 个 项 目 所 需 的 资 金 为 b j 万 元 , 将 会 获 得 的 利 润 为c
90 x1 40 x 2 10 x 3 37 x 4
同理得到后三年的约束 条件。
26
0-1整数规划的标准型和解法
标准型
min z
必须大于 等于0
n
c
j 1
运筹学-4-整数规划ppt课件
.
8
第四章 整数规划 0-1规划
解:设xi
1 0
带第 i件物品
不带第 i件物品 数学模型:
Z表示所带物品的总价值
m
Z ci 带第i件
ci xi
i 1
m
携带物品的总重量 bi x i
i 1
m
max Z ci xi
m i1
s.t
i1
bi xi
b
xi 0,1,
i 1, 2, m
i1
1, 2,..., m
i1
s.t. xij bj j 1, 2 , n
i1
xij
0
,
yi 0,1
混合型整数规划
.
11
第四章 整数规划
例 工厂A1和A2生产某种物资。由于该种物资供不应求,故需要再 建一家工厂。相应的建厂方案有A3和A4两个。这种物资的需求地有 B1,B2,B3,B4四个。各工厂年生产能力、各地年需求量、各厂至各需 求地的单位物资运费cij,见下表:
.
10
第四章 整数规划
解:设 xij表示A 工 i运厂 往B 商 j的店 运量
m
n
则总运费为
c ij x ij
i1 j 1
数学模型:
mn
m
设yi
1 0
则总建厂费为
在第 i个地点建m厂in Z
不在第 i个地点建厂 n
m
fi yi
j1 m
xij
i1
j
ai
1
yi
cij xij
i
fi yi
1 若 建 工 厂 yi 0 若 不 建 工 厂(i3,4)
再设xij为由Ai运往Bj的物资数量,单位为千吨;z表示总费用, 单位万元。
运筹学 第四章 整数规划与分配问题
第四章 整数规划与分配问题
冯大光制作
(4)
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
第二节 分配问题与匈牙利法
在实际中经常会遇到这样的问题,有n 项不同 的任务,需要n 个人分别完成其中的一项,但由 于任务的性质和各人的专长不同,因此各人去 完成不同的任务的效率(或花费的时间或费用) 也就不同。于是产生了一个问题,应指派哪个 人去完成哪项任务,使完成 n 项任务的总效率 最高(或所需时间最少),这类问题称为指派 问题或分配问题。
种下料方式可以得到各种零件的毛坯数以及每种
零件的需要量,如表所示。问怎样安排下料方式, 使得即满足需要,所用的原材料又最少?
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
设:xj 表示用Bj (j=1.2…n) 种方式下料根数模型:
x1 … xn
零件 方 个数 式 零件
A1 b1 Am am1 amn bm
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
逻辑变量的应用
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
(3)两组条件满足其中一组
若 x1 4,则 x2 1 ;否则(即 x1 4 时) 2 3 x
列的零元素,则只要令这些零元素位置的 xij 1 ,其 n n 余的 xij 0 ,则 z aij xij 就是问题的最优解.
i 1 j 1
沈阳农业大学
第四章 整数规划与分配问题
冯大光制作
如效率 矩阵为
管理运筹学第四章整数规划与指派问题
货物运输路线选择案例
案例描述
某物流公司需要为其客户提供从起点到终点的货物运 输服务。在运输过程中,有多种可能的路线可以选择 ,每条路线都有不同的运输成本和时间。此外,客户 对货物的运输时间和成本也有一定的要求。
整数规划应用
该案例可以通过整数规划来解决。首先,将每条路线的 选择定义为整数决策变量,1表示选择该路线,0表示 不选择。然后,根据每条路线的运输成本和时间,构建 目标函数,即最小化总运输成本和时间。接下来,根据 客户的要求和路线的特点,构建约束条件,如运输时间 限制、成本限制和路线连通性等。最后,使用整数规划 求解算法,找到满足所有约束条件的最优路线组合,即 最小化总运输成本和时间的路线选择方案。
展望
未来,整数规划与指派问题将在更多领域得到应用和推广 ,为实际问题的解决提供更加有效的方法和工具。同时, 随着相关技术的不断发展,整数规划与指派问题的求解方 法将更加高效和精确,为相关领域的发展提供更加有力的 支持。
THANKS
感谢观看
要点一
Xpress
Xpress是一款功能强大的数学优化求 解器,适用于线性规划、整数规划等 多种问题。它提供了丰富的算法和工 具,支持大规模问题的求解和分析。
要点二
LINGO
LINGO是一款易于使用的数学优化建 模工具,具有直观的语法和丰富的函 数库。它可以帮助用户快速构建和求 解线性规划、整数规划等问题,并提 供详细的解决方案和报告。
原理
通过添加割平面约束条件,逐 步缩小问题的可行域,从而找 到整数最优解。
添加割平面
根据松弛问题的最优解,构造 一个割平面约束条件,添加到 原问题中。
迭代
重复添加割平面和求解新问题 的步骤,直到找到整数最优解 或确定无整数最优解为止。
管理运筹学4 整数规划
甲
乙
丙 丁
39
34 24
38
27 42
26
28 36
20
40 23
33
32 45
x ij 0或1 ,i、j 1,2,3,4
整数规划的特点及应用
整数规划问题的求解方法: 分支定界法和割平面法
Page 19
匈牙利法(指派问题)
分配问题与匈牙利法
指派问题的数学模型的标准形式:
Page 20
设n 个人被分配去做n 件工作,规定每个人只做一件工作, 每件工作只有一个人去做。已知第i个人去做第j 件工作的效率 ( 时间或费用)为Cij(i=1.2…n;j=1.2…n)并假设Cij ≥0。问应 如何分配才能使总效率( 时间或费用)最高? 设决策变量
每项工作只能安排一人,约束条件为:
x11 x 21 x 31 x 41 x12 x 22 x 32 x 42 x13 x 23 x 33 x 43 x14 x 24 x 34 x 44 1 1 1 1
Page 18
变量约束:
0-1型整数线性规划:决策变量只能取值0或1的整数线性 规划。
整数规划的特点及应用
如
Page 5
1. 变量是人数、机器设备台数或产品件数等都要求是整数 2. 对某一个项目要不要投资的决策问题,可选用一个逻辑变 量 x,当x=1表示投资,x=0表示不投资; 3. 人员的合理安排问题,当变量xij=1表示安排第i人去做j工作,
整数规划的特点及应用
min z c ij x ij [1200y1 1500y 2 ]
i 1 j 1 4 4
运筹学整数规划
2 0 5 0 3
0 0 7 0 6
2 0 2 4 5
3。重复。依行,不考虑划去的0,只有一个0的 对0打圈,划去列 2
1
5 2 0 9 0
0 3 10 8 6
2 0 5 0 3
0 0 7 0 6
2 0 2 4 5
3
4。重复。依列,不考虑划去的0,只有一个0的 对0打圈,划去行 2
7 8 11
任务 人员
分派情况
甲 乙 丙
E J G R
丁
4
15
所需时间
13
9
甲 1 1 乙 1 丙 1 丁
对应每个指派问题, 都有类似的表格,我们称之 为效率矩阵或系数矩阵,某元素 cij ( i , j = 1,2, · · · · · · , n ) 表示指派第 i 个人去完成 第 j 项任务时的效率(或
整数规划问题的求解要比一般的线性规划困难
本章将着重讨论 1。一类特殊的整数规划——指派问题,它的数 学模型和求解。 2。求解整数规划方法——分枝定界法。
一、指派问题的数学模型
1。数学模型
某单位需要指派 n 个人去完成 n 项任务,每个人 只做一项工作,同时,每项工作只由一个人完成。由 于各人的专长不同,每个人完成各项任务的效率也不 同。于是产生了应指派哪一个人去完成哪一项任务, 使完成 n 项任务的总效率最高(如所用的时间为最 少)。我们把这类问题称之为指派问题或分派问题 (Assignment Problem)。
二、匈牙利法
指派问题的效率矩阵的每一个元素aij≥0
解矩阵是每行或每列只能有一个元素为1,其余 均为 0 的 n 阶方阵。如:
0 0 ( xij ) 1 0 1 0 0 0 0 1 0 0 0 0 0 1
运筹学-第4章--整数规划习题
第四章 整数规划4.1 某工厂生产甲、乙两种设备,已知生产这两种设备需要消耗材料A 、材料B,有关数据如下,问这两种设备各生产多少使工厂利润最大?(只建模不求解)解:设生产甲、乙这两种设备的数量分别为x 1、x 2,由于是设备台数,则其变量都要求为整数,建立模型如下:2123max x x z +=⎪⎪⎩⎪⎪⎨⎧≥≤+≤+为整数21212121,0,5.45.01432x x x x x x x x4.2 2197max x x z +=⎪⎩⎪⎨⎧≥≤+≤+-且为整数0,35763.212121x x x x x x t s 割平面法求解。
(下表为最优表)线性规划的最优解为:63max ,0,2/7,2/94321=====z x x x x由最终表中得:27221227432=++x x x ﻩ④ 将系数和常数项分解成整数和非负真分式之和,上式化为;2132********+=++x x x移项后得:①②③④①②③即:21221227212212274343-≤--→≥+x x x x 只要把增加的约束条件加到B 问题的最优单纯形表中。
表4-4由x1行得:7327171541=-+x x x 将系数和常数项分解成整数和非负真分数之和:74476715541+=+-+x x x x得到新的约束条件: 74767154-≤--x x747671654-=+--x x x 在的最优单纯形表中加上此约束,用对偶单纯形法求解:则最优解为3,421==x x ,最优目标函数值为z *=55。
4.3 m ax z =4x1+3x 2+2x 3⎪⎪⎩⎪⎪⎨⎧=≥+≥++≤+-10,,13344352.32132321321或x x x x x x x x x x x t s 隐枚举法解:(1)先用试探的方法找出一个初始可行解,如x 1=x2=0,x 3=1。
满足约束条件,选其作为初始可行解,目标函数z 0=2。
运筹四 整数规划
– 若松弛问题最优解为整数解,则其也是整数规划的 最优解
2、分枝过程
– 若松弛问题最优解中某个 xk=bk 不是整数,令 bk 为 bk 的整数部分 – 构造两个新的约束条件 xk bk 和 xk bk +1,分 别加于原松弛问题,形成两个新的整数规划
0 3 1 ( 0) 0 1 2 0*
逐 列 检 查
2 3 ( 0) 1
( 0) 2 2 2
0 * 1 0 * 2
3、重复1、2后,可能出现三种情况; a. 每行都有一个 (0),显然已找到最优解,令对应(0)位置的 xij=1; b. 仍有零元素未标记,此时,一定存在某些行和列同时有多个零, 称为僵局状态,因为无法采用 1. 2 中的方法继续标记。 4、打破僵局。令未标记零对应的同行同列上其它未标记零的个数为 该零的指数,选指数最小的先标记 ( );采用这种方法直至所有零都 被标记,或出现 情况 a,或 情况 c 。 10
表4.2.1 分枝问题解可能出现的情况
序号 问题 1 问题 2 无可行解 无可行解 1 无可行解 整数解 2 无可行解 非整数解 3 整数解 整数解 4 非整数解 5 整数解,目标函 数优于问题 2 整数解 非整数解,目标 6 函数优于问题 1 说 明 整数规划无可行解 此整数解即最优解 对问题 2 继续分枝 较优的一个为最优解 问题 1 的解即最优解 问题 1 停止分枝(剪 枝), 其整数解 为 界, 对问题 2 继续分枝
9
清华算法的步骤:例4.6.1
2、逐列检查,若该列只有一个未标记的零,对其加( )标记,将( )标 记元素同行同列上其它的零打上*标记。若该列有二个以上未标记的 零,暂不标记,转下一列检查,直到所有列检查完;
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2020/5/30
28
(3)设法使每一行都有一个打括号的“0”元素。 按定理1继续对矩阵进行变换:
ⅰ)从矩阵未被直线覆盖的元素中找出最小者k,
ⅱ)对矩阵中无直线覆盖的行,令 ui k ,有直 线覆盖的列,令 v j k 。其余为0。
ⅲ)对矩阵的每个元素计算 aij ui v j ,得到 一个新矩阵,转第三步重复进行,直至每一行都有 一打括号的0元素。
(2)从第一列开始,若该列只有一个“0”元素, 则对该“0”元素打括号( ),并划去该“0”元 素所在的行;若该列无“0”元素或有两个以上的 “0”元素(不含划去的0),则转下一列;Βιβλιοθήκη 2020/5/3026
(0) 8
2
5
11 (0) 5
4
2
3 (0) 0
0
11
4
5
完成上述步骤后可能出现下列情况:
ⅰ)效率矩阵的每一行都有一个打括号的0元素, 则按照打括号的0元素位置指派任务,即是最优解;
2020/5/30
14
解:设x j 表示代号为j的包装箱的订做数量。
0, 不订j包装箱 y j 1, 订j包装箱
目标函数
6
min z 1200 y j 5x1 8x2 10x3 j 1
12.1x4 16.3x5 18.2x6
约束条件
x j My j , j 1,...6
2020/5/30
2020/5/30
27
ⅱ)打括号的0元素个数小于m,但未被划去的0元
素之间存在闭回路,则沿此闭回路,每隔一个0元
打一括号,然后对打括号的0元素所在行或所在列
画直线;
0 0
(0) 0
0
0 0
(0)
0 0
(0) 0
ⅲ)矩阵中所有0元素或被打括号,或被划去,但打
括号的0元素个数 m ,则进入下一步;
0 第i组条件起作用,i=1,2 yi 1 第i组条件不起作用
可将上述条件表示为 x1 4 y1M
x2 x1
1 4
y1M y2M
x2
3
y2M
其中:M是任意大正数 y1 y2 1
2020/5/30
9
4、表示含有固定费用的函数
例如: x j 表示产品 j 的生产数量,其生产费用函数
则 [aij ] 的最优解等价于 [bij ] 的最优解。
定理2:若矩阵A的元素可分为“0”元和“非0”元, 则覆盖“0”元的最少直线数等于位于不同行、不 同列的“0”元的最大个数。
2020/5/30
23
例:有一份说明书,要分别译成英、日、德、俄四种语言, 交给甲、乙、丙、丁四人去完成,各人的效率不同,如何 分配任务,可使总效率最高。
表中数据为完成任务所需时间(单位:小时)。
人甲
乙
丙
丁
任务
英文 2
10
9
7
日文 15
4
14
8
德文 13
14
16
11
俄文 4
15
13
9
2020/5/30
24
匈牙利解法步骤: 1、在效率矩阵每行减去该行最小元素; 2、在效率矩阵每列减去该列最小元素;
2 10 9 7 2 0 8 7 5
15
4
14
8
18
max Z 4x1 5x2 6x3 100 y1 150 y2 200 y3
2x1 4x2 8x3 500 2x1 3x2 4x3 300 x1 2x2 3x3 100 xj Myj x j 0, y j是0 -1变 量 ,j 1,2,3
2020/5/30
19
§3 指派问题及匈牙利解法
y
j
0或1
2020/5/30
10
§2.2 应用举例
例1 东方大学计算机实验室聘用4名大学生(代号
1,2,3,4)和2名研究生(代号5,6)值班。已知各学生从 周一至周五每天可安排的值班时间及每人每小时报酬见下 表所示。
学生 代号
1 2 3 4 5 6
酬金 (元/h) 10.0 10.0
9.9 9.8 10.8 11.3
周一 6 0 4 5 3 0
每天可安排的值班时间(h) 周二 周三 周四
0
6
0
6
0
6
8
3
0
5
6
4
0
4
6
6
2
4
周五 7 3 5 0 0 4
2020/5/30
11
实验室每天开放时间为8:00AM—10:00PM,共14小时。开 放时间内需要有一名学生值班。规定大学生每周值班时间 是8—15小时,研究生是7—12小时,每次值班不小于2小 时。又每名学生每周值班次数不得多于三次,每天值班人 员中至少有一名研究生,每天值班人数不超过3人。试为 该实验室安排一张人员值班表,使得总酬金支出为最少。
运筹学
OPERATIONS RESEARCH
2020/5/30
1
第四章 整数规划与分配问题
整数规划的有关概念及特点 整数规划的应用 指派问题及匈牙利解法 整数规划的求解方法:分枝定界法、割平面法
2020/5/30
2
§1 整数规划的有关概念及特点
§1.1 概念
整数规划: 要求决策变量取整数值的规划问题。 (线性整数规划、非线性整数规划等)
2020/5/30
3
§1.2 整数规划的求解特点
求整数解的线性规划问题,不是用四舍五入法或
去尾法对性规划的非整数解加以处理就能解决的,
用枚举法又往往会计算量太大,所以要用整数规
划的特定方法加以解决。
例: 求解下列整数规划:
max z 3x1 2x2
2x1 3x2 14
s.t
x1
0.5x2
每天不超过3人 每天有一研究生
(7) 2 yij xij aij yij i 1,...6, j 1,..., 5
值班不超过每人可安排的时间
2020/5/30
13
例2 红星日用化工厂为发运产品,下一年度需要6 种不同容积的包装箱,每种包装箱的需求量及生产 一个的可变费用如下表所示。
包装箱代号
4.5
x1, x2 0, 并取整数
2020/5/30
4
x2
分析: 若当作一般线性规划求解, 图解法的结果如下。
x1 0.5x2 4.5
1、非整数规划最优解 (3.25, 2.5) 显然不是整数规划的可行解。
2、四舍五入后的结果 (3, 3) 也不是整数规划的可行解。
(3.25, 2.5)
3、可行解是阴影区 域交叉点,可比较这 些点对应的函数值, 找出最优。(4, 1)
i 1
5
8 xij 15 ,i 1,...4
j 1
5
7 xij 12 , i 5,6
j 1 5
yij 3,
i 1,..., 6
j 1
6
yij 3,
i 1
j 1,...,5
y5 j y6 j 1 j 1,..., 5
每天开放14小时 大学生值班8-15小时 研究生值班7-12小时 每周不超过3次
2020/5/30
6
设M是任意大正数,则原约束中只有k个真正 起作用的情况可表示为:
n
aij xij bi Myi , i 1,2,...,m
j1 y1 y2 ... ym m k
2020/5/30
7
2、约束条件右端项是r个可能值中的一个
n
aij xij b1, 或b2 ,... 或br
容积(m3) 需求量(个)
可变费用(元/个)
1
0.08 500 5.0
2
0.10 550 8.0
3
0.12 700 10.0
4
0.15 900 12.1
5
6
0.20 0.25 450 400 16.3 18.2
由于生产不同容积包装箱时需进行专门的准备、下 料等,生产每一种包装箱的固定费用都是1200元。 又若某容积的包装箱数量不够时,可用比它大的代 替。试问该厂应订做哪几种代号的包装箱各多少个, 可使得费用最省?
xij 0 否则
于是建立模型如下:
mm
min z
aij xij
i1 j 1
m
xij 1,
j 1
i 1,...m
m
xij 1,
i 1
xij 0或1,
j 1,...m i, j 1,...m
2020/5/30
21
§3.1 指派问题的匈牙利解法
该指派问题可当作运输问题解决,但匈牙利解法更 有效。
§3.1 指派问题与模型
m项任务分配给m个人去完成,每人只能完成其中 一项,每项任务只能分给一人完成,应如何分配 使得效率最高?
aij是第j个人完成第i项任务的效率(如 时间)。
人
任务
1
2
…
m
1
a11
a12
…
a1m
2
a21
a22
…
a2m
…
…
…
…
…
m
am1
am2
amm
2020/5/30
20
设
1 第j个人完成第i项任务
2020/5/30
17
资源 金属板(吨) 劳动力(人月) 机器设备(台月)
小号容器 2 2 1
中号容器 4 3 2
大号容器 8 4 3
解:设 x1, x2, x3 分别为小号容器、中号容器和大号容 器的生产数量。
0, 不生产j型号容器 y j 1, 生产j型号容器
建立如下的数学模型: