1、线性规划和非线性规划实验目的 1)了解最优化问题的基本结构和基本建模方法;2)线性规划的求解方法;3）非线性规划的求解方法.一,优化问题的普遍性以及引例1，无处不在的优化每一个人，高致总统首相，总裁经理，平民百姓，无不在做决策：该做什么，该怎么做，才能有最好的效果?甚至自然中的动植物，也时刻面临这样的问题.类似的问题，还广泛的存在于无机世界中.一,优化问题的普遍性以及引例看看下面的例子分别属于哪一类？a)证券的投资组合；b)国家经济发展战略；c)产品规格、性能设计；d)球形的水滴；e)狼群的集体捕食；f)好的购物方案；g)物质分子结构；h)生物的身体构造；i)乘务组排班表;j)光传播路径:直线
2、,反射,折射课堂作业：和你的同桌讨论还有什么方面需要优化的。一,优化问题的普遍性以及引例2,一些成功的优化例子：“最优人员安排”为美国航空每年节约两千万美元.“改进的出货流程”每年为Yellow Freight 公司节约一千七百多万美元.“改进的卡车分派”为 Reynolds 公司每年节约七百万美元.最优全局供应链为数字设备行业节约超过三亿美元.重建的 North America Operations,Proctor and Gamble减少 20%的工厂,每年节约两亿美元.大阪的Hanshin高速的最优安排每年节约一千七百万人小时.为说明最优化的价值,建立了专门的网站,列举了哪些公司的什么问
3、题,运用最优化方法节约和增加了多少金额.有可选的行业,考察的方面,受益的方式,希望同学们各选择其中的一个,提一份报告,以说明最优化的价值.一,优化问题的普遍性以及引例Google上相关搜索的结果：Search phrase number of hits（英文）短语 点击数（中文）“optimize the supply chain”1,160,000 优化供应链414,000“optimize(the)return”2,490,000 优化回报 453,000“optimal experience”32,400,000 最优经历“optimal investment”8,320,000 优化投
4、资 8,250,000“optimal system”84,200,000 优化系统 13,800,000“optimal decision”28,800,000 最优决策 2,890,000“optimize your PC”3,300,000 优化你的PC“optimal choice”25,800,000 最优选择10,900,000“optimal design”77,300,000 优化设计 1,270,000“optimal health”31,900,000 优化健康还有如:优化产业结构 2,830,000 优化人员结构 3,110,000同学们有没有发现,英文和中文短语间有很大
5、的不同,原因可能是什么?一,优化问题的普遍性以及引例3,相关的几句格言：Waste neither time nor money,but make the best use of both.-Benjamin Franklin Obviously,the highest type of efficiency is that which can utilize existing material to the best advantage.-Jawaharlal Nehru It is more probable that the average man could,with no injury
6、 to his health,increase his efficiency fifty percent.-Walter Scott请同学翻译上面的句子，你喜欢那一句？你有什么好的表述？引例1,动物饲料配置问题 美国一家公司以专门饲养并出售一种实验用的动物而闻名。这种动物的生长对饲料中的三种营养成分特别敏感，即蛋白质、矿物质和维生素。需要的营养量蛋白质：70克矿物质：3克维生素：9.1毫克 现有五种饲料，公司希望找出满足动物营养需要使成本达到最低的混合饲料配置。饲料蛋白质(克)矿物质(克)维生素(毫克)1(x1)2(x2)3(x3)4(x4)5(x5)需要量0.302.001.000.601.
7、80700.100.050.020.200.0530.050.100.020.200.089.1每一种饲料每磅所含的营养成分每种饲料每磅的成本饲料12345成本(美元)0.020.070.040.030.05引例2:供应与选址a1.258.750.55.7537.25b1.250.754.7556.57.75d3547611二,优化问题建模的基本步骤介绍 在我们的生活中,始终有这样的问题:为了一定的目的做一些事情,我们可能要考虑有哪些重要的因素,这些因素和要完成的目标之间有什么样的关系.也就是说,我们在做一个决定时,会注意下面的三个要点:目的是什么?有哪些重要的因素?这些因素和你的目标之间有什
8、么样的关系?二，优化问题的表述 目标函数 对应决策者而言,对其有利的程度必须定量的测度,在商业应用中，有效性的测度经常是利润或者成本,但对于政府，更经常的使用投入产出率来测度.表示有效性测度的经常称为目标函数.目标函数要表出测度的有效性,必须说明测度和导致测度改变的变量之间的关系.系统变量分为决策变量和参数.决策变量是指能由决策者直接控制的变量.而参数是指不能由决策者决定的量.实际上，数学模型很少有能表达变量和有效性测度之间的精确关系的.实际上，运筹学分析者的任务就是找出对测度有最重要影响的变量 然后找出这些变量和测度之间的数学关系.这个数学关系也就是目标函数.二，优化问题的表述决策变量和参数
9、 我们称对应决策者可控的量称为决策变量,决策变量的取值确定了系统的最终性能,也是决策者采用决策的依据.在系统中还有一些量,它不能由决策者所控制,而是由系统所处的环境所决定,我们称之为参数.二，优化问题的表述约束条件 约束条件就是决策变量和参数之间的关系.约束集界定决策变量可以取某些值而不能取其他的值.比如对应生产问题,任何活动中,时间和物品不能为负数.当然,也有一些优化问题不带约束条件,我们称之为无约束优化问题.而在实际问题中,决策变量带有约束是普遍的.三,优化问题的分类优化问题的分类可以从几个方面进行:1,从变量取值的连续和离散可以分成:连续优化,离散优化和混合优化2,从问题的线性非线性可以
10、分为:线性规划和非线性规划3,从变量是确定性和随机性可以分为:随机规划和确定性问题.John Von Neumann 约翰冯诺依曼(19031957),美藉匈牙利人.20世纪最杰出的数学家之一,被誉为”计算机之父”,”博弈论之父”.被认为是数学规划的三大创始人之一.以下的三个人物和线性规划的出现有重要的关系.George B.DantzigGeorge B.Dantzig George B.Dantzig(1914-2005),美国人,线性规划单纯形法的创始人,被誉为”线性规划之父”.美国科学院三院院士,美国军方数学顾问,教授.并以其名字设立Dantzig奖.数学规划的三大创始人之一.发现算法
11、时非常年轻,以至到日本时,人们以为”线性规划之父”是个老人,而对他无人问津.Leonid Vitalyevich Kantorovich Kantorovich(1912-1986)苏联人,著名数学家和经济学家,教授,年仅18岁获博士学位.因在经济学上提出稀缺资源的最优配置获诺贝尔奖.线性规划对偶理论的提出者,数学规划的三大创始人之一.非线性规划问题在实践中也是及其常见的.标志着这一学科的产生的奠基性工作由美国的数学家Tucker和Kuhn在1952年的一篇文章.该文章给出了非线性规划问题的必要条件和充分条件,后来成为Kuhn-Tucker条件.这为非线性规划问题的求解算法的提出提供了理论基础
12、和算法的基本思路.相关的规划问题,比如多目标规划,决策论等等.美国一家公司以专门饲养并出售一种实验用的动物而闻名。这种动物的生长对饲料中的三种营养成分特别敏感，即蛋白质、矿物质和维生素。需要的营养量蛋白质：70克矿物质：3克维生素：9.1毫克 现有五种饲料，公司希望找出满足动物营养需要使成本达到最低的混合饲料配置。四,线形规划问题的解法及举例饲料蛋白质(克)矿物质(克)维生素(毫克)1(x1)2(x2)3(x3)4(x4)5(x5)需要量0.302.001.000.601.80700.100.050.020.200.0530.050.100.020.200.089.1每一种饲料每磅所含的营养成
13、分每种饲料每磅的成本饲料12345成本(美元)0.020.070.040.030.05建立数学模型 决策变量：在混合饲料中，每天所需第j种饲料的磅数xj，j=1，2，3，4，5；约束条件：蛋白质：0.30 x1+2x2+x3+0.6x4+1.8x570 矿物质：0.10 x1+0.05x2+0.02x3+0.2x4+0.05x53 维生素：0.05x1+0.1x2+0.02x3+0.2x4+0.08x510自然约束条件：xi0 确定目标：混合饲料的成本最低 0.02x1+0.07x2+0.04x3+0.03x4+0.05x5min完整的线性规划模型：min 0.02x1+0.07x2+0.04
14、x3+0.03x4+0.05x5s.t.0.30 x1+2x2+x3+0.6x4+1.8x570 0.10 x1+0.05x2+0.02x3+0.2x4+0.05x53 0.05x1+0.1x2+0.02x3+0.2x4+0.08x510 xj0 j=1,2,3,4,5;min c cTx xs.t.A Ax xb x x010370,08.02.002.01.005.005.02.002.005.01.08.16.0123.005.0,03.0,04.0,07.0,02.0bAcT归纳：返 回linprog min c cTx x s.t.A Ax xb Aeqx x beq lb x x
15、ubSolve a linear programming problem where c,x,b,beq,lb,and ub are vectors and A and Aeq are matrices.调用格式:x=linprog(f,A,b,Aeq,beq)x=linprog(f,A,b,Aeq,beq,lb,ub)x=linprog(f,A,b,Aeq,beq,lb,ub,x0)x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options)x,fval=linprog(.)x,fval,exitflag=linprog(.)x,fval,exitflag,output
16、=linprog(.)x,fval,exitflag,output,lambda=linprog(.)原油生产计划原油生产计划原油类别原油类别 买入价买入价(元元/桶桶)买入量买入量(桶桶/天天)辛烷值辛烷值(%)硫含量硫含量(%)A455000120.5B35500062.0C25500083.0汽油类别汽油类别卖出价卖出价(元元/桶桶)需求量需求量(桶桶/天天)辛烷值辛烷值(%)硫含量硫含量(%)甲甲703000101.0乙乙60200082.0丙丙50100061.01:1 加工费:4元/桶 能力:=14000桶/天I:安排生产计划，在满足需求的条件下使利润最大 决策变量：目标：甲甲(3
17、000)乙乙(2000)丙丙(1000)A/45X1X2X3B/35X4X5X6C/25X7X8X9fxxxxxxxxxzmin35600025354560004100050200060300070max987654321约束：987654321253545minxxxxxxxxxf总利润最大 需求限制；原料限制；含量限制；非负限制 含量限制非负限制 原料限制需求限制约束1000325.020002325.03000325.01000686122000886123000108612500050005000100020003000963852741963852741987654321963852
18、741xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx0 x总盈利：总盈利：126000元 c=454545353535252525;a1=100100100;010010010;001001001;a2=111000000;000111000;000000111;-1200-600-800;0-1200-600-80;00-1200-600-8;0.500200300;00.50020030;0 00.5002003;b1=3000 20001000;b2=5000 50005000-30000-16000-6000300040001000;v1=zeros(1,9)
19、;x f=linprog(c,a2,b2,a1,b1,v1)z=356000-f甲甲(3000)乙乙(2000)丙丙(1000)A/452400800800B/35000C/256001200200II:通过广告增加销售(1元广告费：增加10桶销售)决策变量：目标：甲甲(3000+)乙乙(2000+)丙丙(1000+)A/45X1X2X3B/35X4X5X6C/25广告广告销售销售X7X103000+10X10X8X112000+10X11X9X121000+10X12fxxxxxxxxxxxxxxziiiimin3800004253545)101000(50)102000(60)103000
20、(70max121091987654321121110约束：121110987654321499599699293949minxxxxxxxxxxxxf总利润最大需求限制；原料限制；产量限制；含量限制；非负限制 含量限制非负限制 产量限制原料限制需求限制约束129631185210741)129631185210741987654321129631185210741101000325.0)102000(2325.0103000325.0101000(68612)102000(88612)103000(108612500050005000101000102000103000 xxxxxxxxxx
21、xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx0 x8001210iix总盈利：总盈利：287750元 c=49 49 49 39 39 39 29 2929-699-599-499;a1=1 00 100100 -100 0;0 10 010010 0 -10 0;0 01 001001 00 -10;a2=1110000 0 0 0 0 0;0001110 0 0 0 0 0;0000001 1 1 0 0 0;-1200-600-8 0 0 1000 0;0-1200-600 -8 0 080 0;00-1200-60 0-8 00 60;0.5002003 0
22、0-100 0;00.500200 3 0 0-20 0;000.50020 0 3 00 -10 000000 0 001 1 1;b1=300020001000;b2=500050005000-30000-16000-6000 3000 40001000800;v1=zeros(1,12);x f=linprog(c,a2,b2,a1,b1,v1)z=380000-f甲甲(3000)乙乙(2000)丙丙(1000)A/452121.82185.3692.9B/35695.54036.8267.6C/25广告广告182.703277.9750 39.40某公司有某公司有6 6个建筑工地，位置
23、坐标为个建筑工地，位置坐标为(a(ai i,b,bi i)()(单位：公单位：公里里),),水泥日用量水泥日用量d di i(单位：吨）单位：吨）ia1.258.750.55.7537.25b1.250.754.7556.57.75d35476111)现有 2 料场，位于 A(5,1),B(2,7),记(xj,yj),j=1,2,日储量 ej各有 20 吨。假设：料场假设：料场和工地之间和工地之间有直线道路有直线道路目标：制定每天的供应计划，即从 A,B 两料场分别向各工地运送多少吨水泥，使总的吨公里数最小。五,非线形规划问题的解法及举例用例中数据计算，最优解为i1234561ic（料料场场
24、A）3507012ic（料料场场 B）00406102,1,6,.,1,.)()(min612121612/122jecidctsbyaxcjijiiijjjiijijij线性规划模型线性规划模型决策变量：决策变量：c ci ji j (料场料场j j到到工地工地i i的运量）的运量）1212维维Shili084lin.mShili084lin.m选址问题：选址问题：NLPNLP2 2）改建两个新料场，需要确定新料场位置）改建两个新料场，需要确定新料场位置(x xj j,y,yj j)和运量和运量c cijij ，在其它条件不变下使总吨公里数最小。，在其它条件不变下使总吨公里数最小。2,1,6
25、,.,1,.)()(min612121612/122jecidctsbyaxcjijiiijjjiijijij决策变量：决策变量：c ci ji j，(x xj j,y,yj j)16)16维维非线性规划模型非线性规划模型结果：总吨公里数结果：总吨公里数为为85.385.3，但局部最，但局部最优解问题严重优解问题严重Shili084.m:Shili084.m:shili084fun.mshili084fun.m2611612/12222612/12121)(,6,1,.)()()()(minecdecidctsbyaxcdbyaxciiiiiiiiiiiiiii决策变量：决策变量：c ci i
26、，(x xj j,y,yj j)10)10维维 计算方法的改善计算方法的改善局部最优解问题有所改进i123456新新料料场场位位置置（),jjyx1ic304760（3.2552 5.6528）2ic0500011（7.2497 7.7499）+为工地,数字为用量;*为新料场,数字为供应量。012345678901234567835476112016约束非线性规划情形约束非线性规划情形;,0)(,0)(.)(min2121GGGULGGtsfnRxxxxx标准模型标准模型fminconfmincon语句的具体用法语句的具体用法ULbeqAeqbAGGtsfnRxxxxxx,x,0)(,0)(.
27、)(min21x,fval,exitflag,output=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,x,fval,exitflag,output=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,con)con)建立建立m m文件函数文件函数 function f,G=fun(x)function f,G=fun(x)f=f(x);f=f(x);G=G G=G1 1(x),G(x),G2 2(x);(x);选项选项optionsoptions 输出参数输出参数options(8)=options(8)=目标函数最优值；目标函数最优值；输入参数输入参数o
28、ptions(13)=options(13)=等式约束的个数等式约束的个数(m)(m)；L L，U U是决策变量的下界和上界；是决策变量的下界和上界；2 2、fminconfmincon语句的具体用法语句的具体用法maxmax f f(x x)=x=x1 12 2+x x2 22 2-x x1 1x x2 2-2-2x x1 1-5-5x x2 2 s.t.-(s.t.-(x x1 1 1)1)2 2+x x2 2 0 0 2 2 x x1 1-3-3x x2 2+60,+60,x x0 0=0,=0,11例2fminconfmincon语句的具体用法语句的具体用法转化成标准形转化成标准形mi
29、nmin f f(x x)=-x=-x1 12 2-x x2 22 2+x x1 1x x2 2+2+2x x1 1+5+5x x2 2 s.t.(s.t.(x x1 1 1)1)2 2-x x2 2 0 0 -2 -2 x x1 1+3+3x x2 2 60,60,x x0 0=0,1=0,1 function f=fun22(x)function f=fun22(x)f=-x(1)2-x(2)2+x(1)f=-x(1)2-x(2)2+x(1)*x(2)+2x(2)+2*x(1)+5x(1)+5*x(2);x(2);function G,Geq=cont2(x)function G,Geq=
30、cont2(x)G=(x(1)-1)2-x(2);G=(x(1)-1)2-x(2);Geq=;Geq=;x0=0 1;x0=0 1;A=-2,3;b=6;A=-2,3;b=6;Aeq=;beq=;Aeq=;beq=;lb=;ub=;fun22 lb=;ub=;fun22x,fval=fmincon(fun22,x0,A,b,Aeq,beq,lb,ub,cont2)x,fval=fmincon(fun22,x0,A,b,Aeq,beq,lb,ub,cont2)x=1.0e+008 x=1.0e+008*-0.0006 -2.7649-0.0006 -2.7649 fval=-7.6432e+01
31、6 fval=-7.6432e+016fun22=fun22=-x(1)2-x(2)2+x(1)-x(1)2-x(2)2+x(1)*x(2)+2x(2)+2*x(1)+5x(1)+5*x(2)x(2)fminconfmincon语句的具体用法语句的具体用法例179.22,662.26708.101.24986.14.721.24914.5054.34164.17089.610,.,1,10,.,1,101212222.)/ln()(min61098737654106321101101，是常数，其值分别为：其中jcjxxxxxxxxxxxxxxxtsxxcxxfjjiijjjjfminconfm
32、incon语句的具体用法语句的具体用法function fm=ex02(x)function fm=ex02(x)n=10;n=10;c=-6.089,-17.164,-34.054,-5.914,-24.721,c=-6.089,-17.164,-34.054,-5.914,-24.721,-14.986,-24.1,-10.708,-26.662,-22.179;-14.986,-24.1,-10.708,-26.662,-22.179;sx=0;sx=0;for i=1:nfor i=1:n sx=sx+x(i);%sx=sx+x(i);%endendfm=0;fm=0;for i=1:
33、nfor i=1:n fm=fm+x(i)fm=fm+x(i)*(c(i)+log(x(i)/sx);(c(i)+log(x(i)/sx);endendniix1计算1)1)建立目标函数建立目标函数 ex02.mex02.m供应与选址a1.258.750.55.7537.25b1.250.754.7556.57.75d3547611供应与选址范例：供应与选址(a1,b1),d1=3(a2,b2),d2=5(a6,b6),d6=11A(x1,y1)e1=20B(x2,y2)e2=20 x11x12x16x21x26x22目标：使得吨公里min建立规划模型02,1,6,1,.)()(min6121
34、222161ijjiijijijijijjiijxjexidxtsbyaxxf20吨吨其中其中a a，b b，d d，e e 都是常数向量（数组）。都是常数向量（数组）。当料场位置为当料场位置为(x(xj j,y,yj j)为已知时为已知时,上述模型为线上述模型为线性规划模型。从料场性规划模型。从料场j j向工地运送量为向工地运送量为x xijij（决策变（决策变量）。量）。02,1,6,1,.)()(min6121222161ijjjijiiijijijijijxjexidxtsbyaxxf思考：矩阵形式 当料场位置为当料场位置为(x(xi i,y,yi i)为未知时为未知时,上述上述模型为
35、非线性规划模型。对决策变量模型为非线性规划模型。对决策变量(x(xi i,y yi i)而言。而言。线性规划求解结果：线性规划求解结果：最优目标值最优目标值f=136.23f=136.23（吨千米）（吨千米）j123456x1j(料场A)350701x2j(料场B)0040610程序程序liaocx.mliaocx.mfunction f,g=liaoch(x)function f,g=liaoch(x)a=1.25 8.75 0.5 5.75 3 7.25;a=1.25 8.75 0.5 5.75 3 7.25;b=1.25 0.75 4.75 5 6.5 7.75;b=1.25 0.75
36、4.75 5 6.5 7.75;d=3 5 4 7 6 11;d=3 5 4 7 6 11;e=20 20;e=20 20;f1=0;f1=0;for j=1:6for j=1:6s(j)=sqrt(x(13)-a(j)2+(x(14)-b(j)2);s(j)=sqrt(x(13)-a(j)2+(x(14)-b(j)2);f1=f1+s(j)f1=f1+s(j)*x(j);%A(xx(j);%A(x1313,x,x1414)到各工地的吨公里数endendf2=0;f2=0;for i=7:12for i=7:12s(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2);s
37、(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2);f2=f2+s(i)f2=f2+s(i)*x(i);%B(xx(i);%B(x1515,x,x1616)到各工地的吨公里数endendf=f1+f2;%f=f1+f2;%总的目标函数for i=1:6for i=1:6 g(i)=x(i)+x(i+6)-d(i);%g(i)=x(i)+x(i+6)-d(i);%第一组约束条件endendg(7)=sum(x(1:6)-e(1);g(7)=sum(x(1:6)-e(1);g(8)=sum(x(7:12)-e(2);%第二组约束条件6,1,21idxiiij2,1,61j
38、exjjijx0=ones(1,12);x0=ones(1,12);L=zeros(1,12);L=zeros(1,12);opt(13)=6;opt(13)=6;x,opt=constr(liaoch,x0,opt,L)x,opt=constr(liaoch,x0,opt,L)f=opt(8),n=opt(10)f=opt(8),n=opt(10)MATLABMATLAB运行程序运行程序(liaochj.m)(liaochj.m)i123456ABCi1(料场A)0507085.906x3xCi2(料场B)3040635.073y6.5y计算结果：最优目标值计算结果：最优目标值f=93.14
39、f=93.14（吨千米）（吨千米）新料场位置的改变，目标值比原来减少了新料场位置的改变，目标值比原来减少了43.143.1吨千米。吨千米。原最优目标值原最优目标值f=f=136.23136.23六,实验布置实验目的 1)了解最优化问题的基本结构和基本建模方法;2)线性规划的求解方法;3）非线性规划的求解方法.1,(糖果生产计划)某糖果公司用三种原料:糖/果仁/巧克力生产两种糖果A和B.所有重量单位为公斤,所有价格单位为元.假定原料混合过程中重量既不增加也不减少.三种原料分别只有100/20/30公斤,并要求糖果A中的果仁和巧克力含量不少于10%,糖果B中的果仁含量不少于20%.糖果A和B的单价
40、分别为25元和20元.问如何安排生产才能得到最高的收益.试建立相应的数学模型并进行求解.2,(运输问题)从Toronto和Detroit两市分别有两批货物途径Chicago和Buffalo最终到达New York、Phila.和St.louis市.之间的路线表述如下图:Toronto1Detroit2Chicago3Buffalo4NewYork 5Phila.6St.louis7SourceTransshipmentPointDestinationToFromChicagoBuffaloSupplyTorontoDetroit$4$5$7$7600500ToFromNew YorkPhila
41、.St.louisChicagoBuffaloDemand$3$1450$2$3350$2$4300 其中Toronto和 Detroit 分别有600和500的货物需要运出,New York、Phila.和St.louis的货物需求分别是450、350和300.每一段上的运输单价如下面两表：问：如何进行运输安排使整个的运输费用最少？试建立问题的数学模型并求出最有解。3 3、求解非线性规划、求解非线性规划5,4,3,2.33.22,1,3.23.20105;10.)exp(min3231543251254321ixixxxxxxxxtsxxxxxfiiii4 4、桃李花园服务中心选址、桃李花园
42、服务中心选址初始点x0=20，20设设(a ai i,b,bi i)(i=120)(i=120)为第为第i i 栋住宅楼的坐标；栋住宅楼的坐标；a=29.74 4.9 69.32 65.0 98.3 a=29.74 4.9 69.32 65.0 98.3 55.27 40.0 19.8 55.27 40.0 19.8 62.5 73.3 37.58 0.98 41.98 62.5 73.3 37.58 0.98 41.98 75.37 79.38 92.0 75.37 79.38 92.0 84.47 36.77 62.08 73.13;84.47 36.77 62.08 73.13;b=19.39 90.48 56.92 63.18 23.44 b=19.39 90.48 56.92 63.18 23.44 54.88 93.16 54.88 93.16 33.5 65.5 39.19 62.73 69.9 39.72 33.5 65.5 39.19 62.73 69.9 39.72 41.37 65.5241.37 65.52 83.75 37.16 42.52 59.46 56.58;83.75 37.16 42.52 59.46 56.58;感谢下感谢下载载
