会议筹备问题.doc

.
郭继萍
郭 雪
姬娜娜
会议筹备问题
一、摘要
如何为与会代表预定宾馆客房，租借会议室，并租用客车的问题是本文重点讨论和分析的。经过仔细分析与讨论，可将本文分文三个问题进行解决并运用MATLAB、LINGO软件对模型进行辅助求解，建立优化模型对本文进行进一步求解：
关于问题一预测本届与会代表数及预订客房数，这首先预测本届实际与会人数及预订客房数，通过建立随机模型计算求解可以算出实际与会人数m=668，而且与会代表对住房的要求分为六类从而来确定需要预定客房的代表的数量即宾馆客房数为：合住一为104间；合住二为69间；合住三为23间；独住一为148间；独住二为86间；独住三为54间；筹备组在预定客房时需要根据往届所统计与会人员的数据信息这可以使本次会议订房时损失尽量降低到最小，并且可以降低与会代表的不满意度。
关于问题二制定宾馆客房的合理预定方案，根据问题的要求应该从管理、方便、费用合理、住房条件优越等角度来考虑问题 ，由于预订宾馆会造成一些损失；损失一：空房费(预定的客房数多于参加与会代表人数),损失二：房间不足(预定的客房数少于参加与会代表人数)带来损失由筹办方负责，由于客房房费是由与会代表自付，所以在确定宾馆的客房预订间数时要使宾馆总数为最小，建立最优目标函数来求解，得到预定的宾馆和各类宾馆中客房类型和数量见附录表。
关于问题三对租用会议室和客车的研究，要求建立整数规划模型计算租用会议室和客车所用费用达最小，在会议室的选定上，因为我们事先无法知道参加与会的代表是否在下榻的宾馆的会议室参加会议，所以假设每个与会代表参加任何一个会议的概率是相等的即[m/6] 所以会议室的规模至少也应该是 [m/6] ，再有问题（2）得到的四家宾馆中去除不合规格的会议室；由于租用会议室和客车时还有不在本宾馆开会需要移动的人数不确定，分析可建立优化模型使得租用客车和会议室费用达最小，用LINGO软件求解得到的最优解为实际应租会议室为4个，租用11辆45座的客车.
在解决本模型时，从数学角度、实际问题出发，建立优化模型及整数规划模型对问题进行分析，并运用软件进行求解。
关键词： 优化模型 线性规划 整数规划模型 随机模型 二项分布 MABLTB软件 LINGO软件
一、 问题重述
某市的一家会议服务公司负责承办某专业领域的一届全国性会议，会议筹备组要为与会代表预订宾馆客房，租借会议室，并租用客车接送代表。由于预计会议规模庞大，而适于接待这次会议的几家宾馆的客房和会议室数量均有限，所以只能让与会代表分散到若干家宾馆住宿。为了便于管理，除了尽量满足代表在价位等方面的需求之外，所选择的宾馆数量应该尽可能少，并且距离上比较靠近。
筹备组经过实地考察，筛选出10家宾馆作为备选，它们的名称用代号①至⑩表示，相对位置见附图，有关客房及会议室的规格、间数、价格等数据见附表1。
根据这届会议代表回执整理出来的有关住房的信息见附表2。从以往几届会议情况看，有一些发来回执的代表不来开会，同时也有一些与会的代表事先不提交回执，相关数据见附表3。附表2，3都可以作为预订宾馆客房的参考。
需要说明的是，虽然客房房费由与会代表自付，但是如果预订客房的数量大于实际用房数量，筹备组需要支付一天的空房费，而若出现预订客房数量不足，则将造成非常被动的局面，引起代表的不满。
会议期间有一天的上下午各安排6个分组会议，筹备组需要在代表下榻的某几个宾馆租借会议室。由于事先无法知道哪些代表准备参加哪个分组会，筹备组还要向汽车租赁公司租用客车接送代表。现有45座、36座和33座三种类型的客车，租金分别是半天800元、700元和600元。
请你们通过数学建模方法，从经济、方便、代表满意等方面，为会议筹备组制定一个预订宾馆客房、租借会议室、租用客车的合理方案。
二、 问题分析
本题主要是研究会议筹备组为与会代表预定宾馆客房，租借会议室，并租用客车接送代表的问题。由于规模庞大且宾馆的客房数和会议室数量均有限，与会人员只能分散到几家宾馆，综合：管理、代表满意、经济、方便等要求，需要解决以下三个问题：
（1） 首先要预测本届与会代表的数量及预订客房数；
（2） 其次就是制定各个宾馆客房数量的方案；
（3） 最后制定出租用客车和会议室的最优方案。
2.1对问题（1）的分析
要想预测本届与会代表的数量，就要根据以往每届与会代表的数量与以往每届发来回执的代表数量的比值进行预测。（以往每届与会代表数=发来回执的代表数量—发来回执但未与会的代表数量+未发回执但与会的代表数量；比值=与会代表数/发来回执的代表数）用他们的最大值来预测本届与会代表数量的最大值M。
预测本届与会代表m的客房数，这与客房费用有关，我们应参考往届参加与会代表人对客房住宿价位的居中值，来策划本届居住情况；这样可以使得筹备组在预订客房时以费用最少为目标函数进行建模求解，同时还应考虑如果预订客房的数量大于实际用房数量，筹备组需要支付一天的空房费，经计算估计为g=193.01，和预订客房数量不足时，引起代表不满的费用你 n=0.9560（由假设知每位代表是否参加会议是相互独立的，因此可以由以往每届与会代表的数量与与会代表数量的预测值的比例，来计算出每位代表与会的概率P）用MATLAB求解（见附录）。
2.2对问题（2）的分析
指定宾馆客房数量应该考虑管理、方便、代表满意等方面。同时在满足代表合住、独住、价位、以及独住客房不够被安排在双人房，及宾馆间的距离等要求，应以订宾馆客房数最少为目标函数建立优化模型，用软件进行求解。
2.3对问题（3）的分析
会议室的选定，因为我们事先无法知道参加与会的代表是否在下榻的宾馆的会议室参加会议，所以假设每个与会代表参加任何一个会议的概率是相等的即[m/6] 所一会议室的规模至少也应该是 [m/6] ，再有问题（2）得到的宾馆中去除不合规格的会议室，以会议室的费用最少为目标函数建立优化模型。
在确定租用客车时，应该考虑需要移动的代表的数量，这与与会代表参加任何一个会议的概率是相等的即[m/6]，以租车时的费用最少为目标函式建立优化模型。
三、 符号约束
本届的最大预测值；
以往每届与会代表的最大预测值；
:以往发来回执人数；
:以往发来回执但未与会的人；
:以往未发来回执但来与会的人；
:以往实际与会人数；
:本届发来回执的人；
实际与会代表数量与发来回执的代表数量比例；
以往每届与会代表的最大预测值与每人实际与会的概率；
：代表各宾馆之间距离的系数值；
：为选择宾馆的变量，(不选该宾馆，选择该家宾馆）
：表示第家宾馆的第类客房的合住预订间数；
：表示第家宾馆的第类客房的独住预订间数；
：表示第家宾馆的第类客房独住的客房数；
：表示第家宾馆的第类客房合住的间数；
：表示与会代表发来回执要求独住所需要预定的客房数（其中包含与会代表个人住双人间）；
：表示与会代表发来回执要求合住所需预定的客房数；
：租用类会议室的数量。
四、 模型假设
1、假设本届与前几届与会代表人数比例没有发生大的变化；
2、要求单住的与会代表被安排在基本设施相同的双人间不影响代表的满意程度， 而且代表入住时不更改回执中的住房要求；
3、与会代表合住男女不能同住；
4、假定每位来宾参加六个会议几率相同；
5、每个会议在分组中讨论的内容都是随机的；
6、由于宾馆之间的距离不远，不因接送代表而耽误开会的时间，租车在用的过程中一切正常；
7、参加会议的人数不因天气原因而影响出勤人数；
8、未发来回执前来参与会代表与发来回执参家与会代表的住房要求相同；
9、客车在运行过程中不出现任何意外。
五、 模型的建立与求解
5.1对问题（1）的模型建立与求解
（1）预测本届实际与会代表的最大值，应该根据以往每届与会代表的实际人数，与发来回执的人数比，来计算本届实际与会代表的最大值：
以往每届与会代表数=发来回执的代表数量—发来回执但未与会的代表数量+未发回执但与会的代表数量即：

比值f=与会代表数/发来回执的代表数即( i=1,2,3,4)

本届与会代表的最大值记做M=本届与会代表发来回执的人数e/以往每届与会代表与发来回执的代表数的比值即：

由题中给出的数据及用MATLAB求解得最大预测值为M=678.
（2）要想确定预订客房的数量我们应该对预测的到达与会代表的人数进行 析，由于预订人数不确定我们会面临的损失有：客房数大于参加与会代表人数而多掏客房费,设本届实际与会代表人数为,‚客房数小于与会代表人数引起不满造成的损失费；建立损失函数为;用表示需要预定的客房数量、g表示一间客房空一天的客房费：
(空房损失费)=*g
‚（缺房损失费）=*g
以往每次与会代表的最大预测值与每人实际与会的概率符合二项分布

损失函数的期望：
=+
设至少有r位代表团预定客房不足而不能入住的概率为即
=
由于不能入住的概率越大对与会代表声誉影响越大于是我们设定<=0.035
所以筹备组预订客房量m满足约束条件<=0.035的情况下用求得最小值用软件进行求解：求得最优解
（3） 确定需要预定各类客房的数量：
由题意知：应首先根据表10-2的数据统计出本届与会代表对合住及单住占总住房的比例，同时由假设知男女不能同住，根据题中住房要求的信息知即：男住房（i）总数比上男女住房总数，女住房（i）总数比上男女住房总数；（i=1,2,3,4,5,6，）
合住的房间数=男同志合住比例*m/2+女同志合住比例*/2
独住的房间数=男同志独住比例*m+女同志独住比例*m
六类客房的预定数量如下表： 表（1）

合住1
合住2
合住3
独住1
独住2
独住3
数量（单位：间）
104
69
23
148
86
54
2、对问题（2）的模型建立与求解根据表‚即每天每间120~160元、161~200元、201~300元三种不同价格的房间分类，以及对问题（1）的求解可以将原表转化为如下：
价位段
类型
价格（元/天）
宾馆代号
间数
120~160
单人间
150
③
27
160
⑥
40
160
⑦
40
双人间
140
②
50
④
50
⑤
35
150
③
50
⑦
50
160
②
35
⑤
35
⑧
40
161~200
单人间
180
①
30
180
⑥
30
180
⑧
45
161~200
双人间
180
①
50
180
②
30
200
②
35
180
③
24
200
④
45
200
⑤
40
170
⑥
40
180
⑧
40
201~300
单人间
220
①
20
300
⑦
30
260
⑨
40
280
⑨
40
双人间
220
①
30
220
⑥
30
260
⑨
30
280
⑨
30
260
⑩
55
280
⑩
45
为了使在住房、管理、经济、方便，距离等问题达到与会代表人员的满意，由此我们建立了如下目标函数，并规定代表各宾馆之间距离的系数值（可以使得计算的结果更完善），为选择宾馆的变量，(不选该宾馆，选择该家宾馆）表示第家宾馆的第类客房的合住预订间数；表示第家宾馆的第类客房的独住预订间数；表示第家宾馆的第类客房独住的客房数；表示第家宾馆的第类客房合住的间数；表示与会代表发来回执要求独住所需要预定的客房数（其中包含与会代表个人住双人间）；表示与会代表发来回执要求合住所需预定的客房数；通过计算使得所建立目标函数取得最小值
；

或1，i=1,2，…10，
，，0，且为整数
用软件求的结果为最少应预定四家宾馆，这四家宾馆分别为，‚，ƒ，⑦运算过程见附录：
各个宾馆的客房安排为：
宾馆
合住
独住
合住1
合住2
合住3
独住1
独住2
独住3
①

50
23
单间

单间
30
单间
20
双间

双间

双间
4
②
85
19

单间

单间

单间

双间
0
双间
46
双间

③
19

单间
27
单间

单间

双间
31
双间
24
双间

⑦
0

单间
40
单间

单间
30
双间
50
双间

双间

合住
104
69
23
148
86
54
（3） 对问题（3）的建立与求解
在租用会议室时，由于无法确定与会代表参加哪项会议，所以假设参加任何一个会议的概率相等即1/6，所以在租用客车时，应该考虑需要移动的人员（需要移动的人员=与会的实际人员—所在各个宾馆会议室开会的人员）
需要租车的人数为：
假设每位与会代表是否在自己的宾馆参加会议的概率是相等的即1/6的可能，从筹办方的角度考虑建立最优模型，是会议室费用和租车费用达到最少而且能按计划完成任务等方面考虑，我们规定用表示满足规模要求的第类会议室的单价，表示的类会议室所在宾馆下榻的代表数量，表示第类会议室的间数，这些数据有问题（2）以得出，表示第类客车的单价，为800元，700元，600元。表示第类客车可以容纳的人数，分别为45人,36人,33人，以会议室数和客车数为约束条件建立优化模型为：
会议室费和客车费为 ： ；
租用会议室的费用如表所示：
会议室编号
所在宾馆
规模
间数
价格（半天）
2
①
150人
1
1200元
3
②
130人
2
1000元
7
⑦
140人
2
800元
8
⑦
200人
1
1000元
六、 结果分析
在对问题（1）进行求解时用的是根据以往四届的数据进行预测，其数据较少、预测不精确，这些放面都有待与提高，同时在预订宾馆客房数时是以总费用最少为目标函数，这一点从实际出发比较经济适用，且满意度高。
解决问题（2）时，首先考虑的是各个宾馆之间的距离，这样可以在解决第三个问题时避免了在租车费用上的大量计算，同时以宾馆数最少为目标函数，实际可行。在租车方面考虑到如果顺道的话可以搭顺风车的情况，是租车费用最少以及根据问题（1）、(2)关于人数和宾馆的求解结果可以是本问题更简单易求得出的数据更实际。
七、 模型评价与推广
优点：在解决模型时运用了MATLAB、LINGO软件，简单明了且结果有一定的可靠性， 利用软件求解不仅方便快捷还可以达到事半功倍的效果！并且对于宾馆人数及客房数都以表格的形式直观明了，在对客房的预订时是以费用最少计算的，这样不仅是代表满意，且经济，也少支付空房费，比较实际，在确定宾馆客房数时考虑到各个宾馆之间的距离远近问题，在计算时是结果一举两得，是本文的一个亮点。
缺点：本模型建立时有些预测的数据不够精确，且在建模时忽略了某些实际问题对模型的影响，使得模型的结果不够精确，这在以后的建模中应尽量避免。同时个人认为从经济、实用、资源合理利用等方面考虑，由于客房费用是与会代表自付，建议最好不要个人独住双人房。同时考虑到与会人员是否为夫妻，所以筹备组不应过于强求男女不能同住。
推广：本模型是为了解决在开大型会议时的代表住房、会议室的安排、以及接送代表的问题。在现实中比较常见所以以后如果有类似的问题都可以使用本模型的建模。
八、 参考文献
【1】姜启源等主编《数学模型》，高等教育出版社，2011
【2】石生明等主编《高等代数》，高等教育出版社，2003
【3】魏振军等主编《概率论与数理统计》，中国铁道出版社，2008
【4】MATLAB软件使用教程
【5】LINGO软件使用教程
十、附录
附录1
附表1 10家备选宾馆的有关数据
宾馆代号
客房
会议室
规格
间数
价格（天）
规模
间数
价格(半天)
①
普通双标间
50
180元
200人
1
1500元
商务双标间
30
220元
150人
2
1200元
普通单人间
30
180元
60人
2
600元
商务单人间
20
220元

②
普通双标间
50
140元
130人
2
1000元
商务双标间
35
160元
180人
1
1500元
豪华双标间A
30
180元
45人
3
300元
豪华双标间B
35
200元
30人
3
300元
③
普通双标间
50
150元
200人
1
1200元
商务双标间
24
180元
100人
2
800元
普通单人间
27
150元
150人
1
1000元

60人
3
320元
④
普通双标间
50
140元
150人
2
900元
商务双标间
45
200元
50人
3
300元
⑤
普通双标间A
35
140元
150人
2
1000元
普通双标间B
35
160元
180人
1
1500元
豪华双标间
40
200元
50人
3
500元
⑥
普通单人间
40
160元
160人
1
1000元
普通双标间
40
170元
180人
1
1200元
商务单人间
30
180元

精品双人间
30
220元

⑦
普通双标间
50
150元
140人
2
800元
商务单人间
40
160元
60人
3
300元
商务套房（1床）
30
300元
200人
1
1000元
⑧
普通双标间A
40
180元
160人
1
1000元
普通双标间B
40
160元
130人
2
800元
高级单人间
45
180元

⑨
普通双人间
30
260元
160人
1
1300元
普通单人间
30
260元
120人
2
800元
豪华双人间
30
280元
200人
1
1200元
豪华单人间
30
280元

⑩
经济标准房（2床）
55
260元
180人
1
1500元
标准房（2床）
45
280元
140人
2
1000元
附表2 本届会议的代表回执中有关住房要求的信息（单位：人）

合住1
合住2
合住3
独住1
独住2
独住3
男
154
104
32
107
68
41
女
78
48
17
59
28
19
说明：表头第一行中的数字1、2、3分别指每天每间120~160元、161~200元、201~300元三种不同价格的房间。合住是指要求两人合住一间。独住是指可安排单人间，或一人单独住一个双人间。
附表3 以往几届会议代表回执和与会情况

第一届
第二届
第三届
第四届
发来回执的代表数量
315
356
408
711
发来回执但未与会的代表数量
89
115
121
213
未发回执而与会的代表数量
57
69
75
104
（其中500等数字是两家宾馆的间距，单位为米。）
⑤（与⑦间距300米） ⑦
⑦（与⑦间距300米） ⑦
⑧（与⑦间距300米） ⑦
①（与⑦间距300米） ⑦
②（与⑦间距300米） ⑦
④（与⑦间距300米） ⑦
③（与⑦间距300米） ⑦
1000
500
300
150
200
300
⑨（与⑦间距300米） ⑦
⑥（与⑦间距300米） ⑦
⑩（与⑦间距300米） ⑦
300
700
150

附录2
附表4
回执情况
第一届
第二届
第三届
第四届
本届
发来回执的代表数量
315
356
408
711
755
发来回执但未与会的代表数量
89
115
121
213

未发回执而与会的代表数量
57
69
75
104

实际来人
283
310
362
602
p
实际来人与发来回执
的人数比例
0.8984
0.8708
0.8873
0.8467

最大预测值
283
320
367
639
678
预测本届与会代表数量的最大值MATLAB如下：表5
>> %预测本届与会代表数量的最大值
>> a=[315,356,408,711];
>> b=[89,115,121,213];
>> c=[57,69,75,104];
>> d=a-b+c
d =
283 310 362 602
>> aman=[154,104,32,107,68,41];
>> awoman=[78,48,17,59,28,19];
>> e=aman+awoman
e =
232 152 49 166 96 60
>> f=d./a
f =
0.8984 0.8708 0.8873 0.8467
>> f1=max(f)
f1 =
0.8984
for i=1:6
w1=e(1)+e(2)+e(3)+e(4)+e(5)+e(6)
w1 =
755
End
>> disp('本届与会代表数量的最大预测值')
本届与会代表数量的最大预测值
>> M=floor(f1*w1)
M =678

确定需要预订客房的代表数量MATLAB求解如下：表6
>> %确定需要预订客房的代表数量
>> p=0.9761;g=193.005;
>> M=678;+
>> for m=200:678
i=0;es=0;
while i<=m
es=es+(m-i)*g*binopdf(i,M,p)/2;
i=i+1;
end
j=i
while j<=M
es=es+(j-m)*g*binopdf(i,M,p)/2;
j=j+1;
end
m
es
k=m+20;
pj=0;
while k<=M
pj=pj+binopdf(k,M,p);
k=k+1;
end
pj
end
六类客房的预定数量：表7
六类客房的预定数量MATLAB求解如下：
>> manzhufang=[154,104,32,107,68,41];
womanzhufang=[78,48,17,59,28,19];
t=sum(manzhufang)+sum(womanzhufang);
r1=manzhufang./t;
r2=womanzhufang./t;
C=668;
for i=1:3
zhufang(i)=ceil(r1(i)*C/2)+ceil(r2(i)*C/2);
end
for i=4:6
zhufang(i)=ceil(r1(i)*C)+ceil(r2(i)*C);
end
zhufang %六类客房的预定数量
zhufang =
104 69 23 148 86 54
每位代表实际与会的概率MATLAB计算如下：表8
>> %每位代表实际与会的概率
>> a=[315,356,408,711];
>> f1=0.8984;
>> for i=1:4
N=a(i)*f1
end
N =
282.9960
N =
319.8304
N =
366.5472
N =
638.7624
>> N=[283,320,367,693];
>> d=[283,310,362,602];
>> for i=1:4
p=d(i)/N(i)
end
p =
1
p =
0.9688
p =
0.9864
p =
0.8687
>> p1=1.0000,0.9688,0.9864,0.8687
>> for i=1;4
q=p1(1)+p1(2)+p1(3)+p1(4)
end
ans =
4
q =
3.8239
>> q/4
ans =
0.9560
平均房价：表9
>>%平均房价计算
a=[50*180+30*220+30*180+20*220+50*140+35*160+30*180+35*200+50*150+24*180+27*150+50*140+45*200+35*140+35*160+40*200+40*160+40*170+30*180+30*220+50*150+40*160+30*300+40*180+40*160+45*180+30*260+30*260+30*280+30*280+55*260+45*280]
a =
229870
>>b=[50+30+30+20+50+35+30+35+50+24+27+50+45+35+35+40+40+40+30+30+50+40+30+40+40+45+30+30+30+30+55+45]
b =
1191
>> c=a/b
c =
193.0059
各个宾馆之间距离的系数值：表10
各个宾馆之间距离的系数值： A=[0,150,900,650,600,600,300,500,650,1300;150,0,750,500,750,750,450,650,800,1450;900,750,0,250,1500,1500,1200,1000,1150,2200;650,500,250,0,1250,1250,950,1150,1300,1950;600,750,1500,1250,0,600,300,500,650,1300;600,750,1500,1250,600,0,300,500,350,700;300,450,1200,950,300,300,0,200,350,1000;500,650,1000,1150,500,500,200,0,150,1200;650,800,1150,1300,650,350,350,150,0,1050;1300,1450,1200,1950,1300,700,1000,1200,1050,0];
>> [V,D]=eig(A)%%求矩阵的最大特征值和特征向量
V =
0.2481 0.0239 0.2360 0.2437 -0.3029 -0.2242 -0.0056 0.7266 -0.0171 -0.0123
0.2673 0.1494 0.1580 0.2945 -0.3237 -0.2027 -0.3175 -0.6163 0.0104 -0.0302
0.4103 0.5681 -0.1935 -0.3790 0.0950 -0.0823 -0.1124 0.0891 0.4778 -0.4095
0.3655 0.5215 -0.1784 0.1772 0.0845 0.3349 0.2358 -0.0037 -0.3903 0.3590
0.3222 -0.2586 0.3235 0.1133 0.8068 -0.0428 -0.2962 0.0615 0.0844 -0.1032
0.2818 -0.3327 0.0924 0.1522 -0.3003 0.7779 -0.2461 0.0876 0.2976 -0.6357
0.2286 -0.2315 0.3235 -0.0004 0.0067 0.0422 0.6351 -0.1874 0.4143 0.1330
0.2590 -0.2060 0.3334 -0.5276 -0.0111 -0.3518 0.3499 -0.1589 -0.5274 -0.1412
0.2803 -0.2691 0.2349 -0.5285 -0.2113 0.0245 -0.3665 0.0896 0.0959 0.4278
0.4288 -0.1906 -0.6826 0.2919 -0.0381 -0.2380 0.1498 -0.0657 -0.2536 0.2588
D =
1.0e+003 *
7.8476 0 0 0 0 0 0 0 0 0
0 -3.2315 0 0 0 0 0 0 0 0
0 0 -2.0673 0 0 0 0 0 0 0
0 0 0 -1.0271 0 0 0 0 0 0
0 0 0 0 -0.6779 0 0 0 0 0
0 0 0 0 0 -0.4349 0 0 0 0
0 0 0 0 0 0 -0.1993 0 0 0
0 0 0 0 0 0 0 -0.1211 0 0
0 0 0 0 0 0 0 0 -0.0211 0
0 0 0 0 0 0 0 0 0 -0.0673
>> a=V(:,1)
a =
0.2481
0.2673
0.4103
0.3655
0.3222
0.2818
0.2286
0.2590
0.2803
0.4288
>> a./sum(a)
ans =
0.0802
0.0865
0.1327
0.1182
0.1042
0.0911
0.0739
0.0838
0.0906
0.1387

定会议室与客车数量：表11
model:
sets:
num_i/1..8/:L,Q1,h,W;
num_j/1..3/:t,Q2,n;
endsets
data:
Q1=1500,1200,1000,1500,1200,1000,800,1000;
h=1,2,2,1,1,1,2,1;
W=200,200,226,226,122,122,132,132;
Q2=800,700,600;
n=45,36,33;
m=668;
enddata
[OBJ]min=@sum(num_i(i):Q1(i)*L(i))+@sum(num_j(j):t(j)*Q2(j));
@sum(num_i(i):L(i))=6;
@for(num_i(i):L(i)<=h(i););
@sum(num_j(j):t(j)*n(j))>=m-@sum(num_i(i):L(i)*W(i)/6);
@for(num_i(i):L(i)>=0;@gin(L(i)););
@for(num_j(j):t(j)>=0;@gin(t(j)););
end

Global optimal solution found.
Objective value: 14600.00
Extended solver steps: 1
Total solver iterations: 54
Variable Value Reduced Cost
M 668.0000 0.000000