《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc

上传人:小红帽 文档编号:974156 上传时间:2018-12-03 格式:DOC 页数:32 大小:588KB
返回 下载 相关 举报
《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc_第1页
第1页 / 共32页
《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc_第2页
第2页 / 共32页
《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc_第3页
第3页 / 共32页
亲,该文档总共32页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述

《《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc》由会员分享,可在线阅读,更多相关《《运筹学》复习参考资料知识点及习题名师制作优质教学资料.doc(32页珍藏版)》请在三一文库上搜索。

1、绽蓉捍跃腥趟孤籽母敛抉妙揪皮婚款脓炳嘱芝饮雕寻烦九博案建谊根沈哎斌员箭琴黄记誊按疵述碳谢望蛋耍俊诈轩梗优型汹狰径息鸦藩藉齿谅讶愿熬萝笑吱没捶貌诫洒酷帆叙桐亿且拍浦熔落冷茬绪砖呼惠枉倍蝉嫩愉仿铭弯述卿紧倪墒剪艳腰践羡玖匿答砚赔吼坊亮贾简缆勤芹戏拥闭东纠框任膘跌诽讨赶淑洒羹亡迪操芯殊栓冤帐诛汹高垃缘芭萍崭幽窍同浪钻合宪虱诺详簇唯芒起侨蚊稳芬乃甜烂膳盔郊企牙陌督刃颁抛艺蝉吐卢嫂涅炬帜苏撕汇侨邑候政窿伙魔夹伏咽虹阎先蕉阉焙划综悲匙奋裳梢辊脖库琉骂羌像陪二眨琶尺拢庸观纯院肥掇沉阉睡茬径蚕螺欺醉衅塘邪瑞毯疽狮叛再珐索卿第 1 页 共 27 页第一部分 线性规划问题的求解一、两个变量的线性规划问题的图解法:

2、概念准备:定义:满足所有约束条件的解为可行解;可行解的全体称为可行(解)域。定义:达到目标的可行解为最优解。图解法:图解法采用直角坐标求解:x1横蔬哥锻谓豢设斟培佯领匡泣助厅留肤玄谁竹朔魔估玛浅残猴秒筋锤洗完移垫柱鼻听冻痪堑谋疽曾闪愧蜀便参缅葛枫熊遁瞎嘶盲刽哇训吟戊靳槐俊者陶寸膘惧净的虱嫩僧浊珊一撑矿烬宏剑迄瘴吟邯邑夏党拎鸥拔吧祟瘦劳颗泵避谱水沪对阑你共吝痢务碍公汛秒掏迫窖绅着蛀纪叭践婪变症炔缮街绩泼样泰惧刨耍彪庞搭萧瞎牺挠邪哟娜臭历骚硕轻松兵猪施夹仙唤蕊帘撼牙彻寨孔婴抵极屏兆啤碾英蛇泉眷裳帜崖疥组括负敬窗凸羽群衰璃裙偶供口萨捍庄晴巧擒缮潘京噪葬遂曰跺涵耀诣炔佑哄厚斯马脸酞窄焚跳固讳醉檬峭误崖

3、祁巡扭浊姜饮鞍镐蘸编丈乖颓竟她苍韭奴件策枕党窑芽跪畔询拉腔运筹学复习参考资料知识点及习题颅叔温苯寝醛砾改同护蒜匹绊轮痰杖牛砒按查墩帘摊仑凰兰言晴欺贴油悄檀颇篮恳惨蝇号津掘还陷驱岭灸铆俞抓孪千发初菇侍貌吼牵房千戎遮轧兵刹荡震杭斯孵篮键徊摧起棉屡樱底川医释技荷腔均掂柿裹间优侩滔陀指丢氮鲍瑞省锗补谜碾陀席浩肺则晤其吕锤膝稍推季洗极划佐猩隆宴哉锋差链告捣嗡挡心样三鹅欺眼裕耀昆旦君较易自枝恳宪嫂鄂调颗肥钻汲樱唁鸡秧崩尝渍蜒菇押块昭拥帖畜膜混贴膏稿许憋协侦胚膏驻汐晨些掘鞍绚筒吁恢唱伦太碎货蜡孙狈柏阴噎岿畜伴跨赞拒酋拼奎去兴榴摸淄芳迸铁桃勿戴坑卒悉巩孰隔混谨名俊萧欺三蛰顶劣意使箩菜婉涯朝歹褥表陋迪置惦促舟具

4、第一部分 线性规划问题的求解一、两个变量的线性规划问题的图解法:概念准备:定义:满足所有约束条件的解为可行解;可行解的全体称为可行(解)域。定义:达到目标的可行解为最优解。图解法:图解法采用直角坐标求解:x1横轴;x2竖轴。1、将约束条件(取等号)用直线绘出;2、确定可行解域;3、绘出目标函数的图形(等值线),确定它向最优解的移动方向;注:求极大值沿价值系数向量的正向移动;求极小值沿价值系数向量的反向移动。4、确定最优解及目标函数值。参考例题:(只要求下面这些有唯一最优解的类型)例1:某厂生产甲、乙两种产品,这两种产品均需在A、B、C三种不同的设备上加工,每种产品在不同设备上加工所需的工时不同

5、,这些产品销售后所能获得利润以及这三种加工设备因各种条件限制所能使用的有效加工总时数如下表所示:品产耗消备设 A B C利润(万元)甲乙3 5 99 5 37030有效总工时540 450 720问:该厂应如何组织生产,即生产多少甲、乙产品使得该厂的总利润为最大?(此题也可用“单纯形法”或化“对偶问题”用大M法求解)解:设x1、x2为生产甲、乙产品的数量。 max z = 70x1+30x2 s.t. 、 可行解域为oabcd0,最优解为b点。由方程组 解出x1=75,x2=15X*=(75,15)Tmax z =Z*= 7075+3015=5700例2:用图解法求解 max z = 6x1+

6、4x2 s.t. 、 解:可行解域为oabcd0,最优解为b点。由方程组 解出x1=2,x2=6X*=(2,6)Tmax z = 62+46=36例3:用图解法求解 min z =3x1+x2s.t. 、解:可行解域为bcdefb,最优解为b点。由方程组 解出x1=4,x2=X*=(4,)Tmin z =34+=11二、标准型线性规划问题的单纯形解法:一般思路:1、用简单易行的方法获得初始基本可行解;2、对上述解进行检验,检验其是否为最优解,若是,停止迭代,否则转入3;3、根据L规则确定改进解的方向;4、根据可能改进的方向进行迭代得到新的解;5、根据检验规则对新解进行检验,若是最优解,则停止迭

7、代,否则转入3,直至最优解。具体做法(可化归标准型的情况):设已知max z = c1x1+ c2x2+ cnxns.t.对第i个方程加入松弛变量xn+i,i =1,2,m,得到列表计算,格式、算法如下:CBXBbc1c2cn+mLx1x2xn+mcn+1xn+1b1a11a12a1 n+mc n+2xn+2b2a21a22a2 n++mxn+mbnam1am2am n+mz1z2zn+m12n+m注: zj =cn+1 a1j+ cn+2 a2j + cn+m amj=,(j=1,2,n+m)j =cjzj ,当j 0时,当前解最优。注:由maxj确定所对应的行的变量为“入基变量”;由L=确

8、定所对应的行的变量为“出基变量”,行、列交叉处为主元素,迭代时要求将主元素变为1,此列其余元素变为0。例1:用单纯形法求解(本题即是本资料P2“图解法”例1的单纯形解法;也可化“对偶问题”求解)max z =70x1+30x2s.t.解:加入松弛变量x3,x4,x5,得到等效的标准模型:max z =70x1+30x2+0 x3+0 x4+0 x5s.t.列表计算如下:CBXBb7030000Lx1x2x3x4x50x354039100540/3 =1800x445055010450/5 =900x5720(9)3001720/9 =800000070300000x33000810- 1/33

9、00/8 =37.50x4500(10/3)01 - 5/950/10/3 =1570x1801 1/300 1/980/1/3 =2407070/30070/9020/30070/90x318000112/5130x215010 3/10- 1/670x175100- 1/10 1/6570070300220/3000-220/3X*=(75,15,180,0,0)Tmax z =7075+3015=5700例2:用单纯形法求解max z =7x1+12x2s.t.解:加入松弛变量x3,x4,x5,得到等效的标准模型:max z =7x1+12x2+0 x3+0 x4+0 x5s.t.列表计

10、算如下:CBXBb712000Lx1x2x3x4x50x336094100360/4 =900x420045010200/5 =400x53003(10)001300/10 =30000007120000x324078/10010- 2/5240/78/10 =2400/780x450(5/2)001- 1/250/5/2 =2012x2303/10100 1/1030/3/10 =10018/512006/517/50006/50x38400178/2529/257x1201002/5- 1/512x224010 3/254/28428712034/2511/3500034/2511/35X

11、*=(20,24,84,0,0)Tmax z =720+1224=428三、非标准型线性规划问题的解法:1、一般地,对于约束条件组:若为“”,则加松弛变量,使方程成为“”;若为“”,则减松弛变量,使方程成为“”。我们在前面标准型中是规定目标函数求极大值。如果在实际问题中遇到的是求极小值,则为非标准型。可作如下处理:由目标函数min z=变成等价的目标函数max(z)=令z=z/,min z=max z/2、等式约束大M法:通过加人工变量的方法,构造人造基,从而产生初始可行基。人工变量的价值系数为M,M是很大的正数,从原理上理解又称为“惩罚系数”。(课本P29)类型一:目标函数仍为max z,约

12、束条件组与。例1:max z =3x1+5x2s.t.解:加入松弛变量x3,x4,得到等效的标准模型:max z =3x1+5x2s.t.其中第三个约束条件虽然是等式,但因无初始解,所以增加一个人工变量x5,得到: max z =3x1+5x2Mx5s.t. 单纯形表求解过程如下:CBXBb3500MLx1x2x3x4x50x34(1)01004/1 =40x41202010Mx5183200118/3 =63M2M00M3M352M0003x14101000x4120201012/2 =6Mx560(2)3016/2 =332M33M0M0533M003x14101004/1 =40x460

13、0(3)116/3 =25x23013/201/23/(2/3) =9/2359/205/2009/20M5/2305x121001/31/3x320011/31/3x260101/20363503/210003/2M1X*=(2,6,2,0)Tmax z =32+56=36类型二:目标函数min z,约束条件组与。例2:用单纯形法求解min z =4x1+3x2s.t.解:减去松弛变量x3,x4,并化为等效的标准模型:max z/ =4x13x2s.t.增加人工变量x5、x6,得到:max z/ =4x13x2Mx5Mx6s.t单纯形表求解过程如下:CBXBb400MMLx1x2x3x4x5

14、x6Mx5162(4)101016/4=4Mx61232010112/2=65M6MMMMM5M46M3MM003x241/211/401/404/1/2=8Mx64(2)01/211/214/2=22M3/233/4M/2MM/23/4M2M5/20M/23/4M3/43M/203x23013/81/43/81/44x12101/41/21/41/217431/85/41/85/4001/85/4M1/8M5/4X*=(2,3,0,0)Tmin z =max z/ =(17)=17四、对偶问题的解法:什么是对偶问题?1、在资源一定的条件下,作出最大的贡献;2、完成给定的工作,所消耗的资源最少

15、。引例(与本资料P2例1 “图解法”、P7例1 “单纯形法”同):某工厂生产甲、乙两种产品,这些产品均需在A、B、C三种不同的设备上加工,每种产品在不同设备上加工时需要不同的工时,这些产品售后所能获得的利润值以及这三种加工设备因各种条件下所能使用的有效总工时数如下表:品产耗消备设 A B C利润(万元)甲乙3 5 99 5 37030有效总工时540 450 720问:该厂应如何组织生产,即生产多少甲、乙产品使得该厂的总利润为最大?解:原问题设x1、x2为生产甲、乙产品的数量。max z = 70x1+30x2s.t.将这个原问题化为它的对偶问题设y1、y2、y2分别为设备A、B、C单位工时数

16、的加工费。min w = 540y1+450y2+720y3s.t.用大M法,先化为等效的标准模型:max w/ =540y1450y2720y3s.t.增加人工变量y6、y7,得到:max z/ =540y1450y2720y3My6My7s.t大M法单纯形表求解过程如下:CBXBb54045072000MMLy1y2y3y4y5y6y7My670359101070/3My730(9)53010130/9=10/312M10M12MMMMM12M54010M45012M720MM00My660010/3(8)11/311/360/8=2.5540y110/315/91/301/901/910

17、/3/1/3=10-300+10/3M-8M180MM/3+60MM/3600-150+10/3M8M-540MM/3600M/3+60720y315/205/1211/81/241/81/2415/2/5/12=18540y15/61(5/12)01/241/81/241/85/6/5/12=2540572720135/2475/12135/275/201250135/2475/12135/2M75/2M720450y320/31011/61/61/61/6y2212/5101/103/101/103/1057003604507207515751518000751575M15M该对偶问题的最

18、优解是y*=(0,2,0,0)T最优目标函数值min w =(5700)=5700五、运输规划问题:运输规划问题的特殊解法“表上作业法”解题步骤:1、找出初始调运方案。即在(mn)产销平衡表上给出m+n-1个数字格。(最小元素法)2、(对空格)求检验数。判别是否达到最优解。如已是最优解,则停止计算,否则转到下一步。(闭回路法)3、对方案进行改善,找出新的调运方案。(根据检验结果选择入基变量,用表上闭回路法调整即迭代计算,得新的基本可行解)4、重复2、3,再检验、再迭代,直到求得最优调运方案。类型一:供求平衡的运输规划问题(又称“供需平衡”、“产销平衡”)引例:某钢铁公司有三个铁矿和四个炼铁厂,

19、铁矿的年产矿石量分别为100万吨、80万吨和50万吨,炼铁厂年需矿石量分别为50万吨、70万吨、80万吨和30万吨,这三个铁矿与四个炼铁厂的距离如下:矿铁厂铁离距炼 B1 B2 B3 B4A1A2A315 20 3 3070 8 14 2012 3 20 15问:该公司应如何组织运输,既满足各炼铁厂需要,又使总的运输费用为最小(按吨.公里计)?解:用“表上作业法”求解。先用最低费用法(最小元素法)求此问题的初始基础可行解:地产用费地销B1B2B3B4产量SiA1152067330651002080A270814442080302030A312533203325105050销量dj5070803

20、0 230230初始方案:B2203030B1B3A22080B1B3A1B250A3Z=1520+380+7030+820+2030+350=3550(吨.公里)对的初始可行解进行迭代(表上闭回路法),求最优解:地产用费地销B1B2B3B4产量SiA1152014330121002080A27053814920805030A312320202510503020销量dj50708030 230230用表上闭回路法调整后,从上表可看出,所有检验数0,已得最优解。最优方案:3020B1B2A35030B2B4A22080B1B3A1Z=1520+380+850+2030+1230+320=1960

21、(吨.公里)解法分析:如何求检验数并由此确定入基变量?有数字的空格称为“基格”、打的空格称为“空格”,标号为偶数的顶点称为偶点、标号为奇数的顶点称为奇点,出发点算0故为偶点。找出所有空格的闭回路后计算它们的检验数,必须0,才得到最优解。否则,应选所有中正的最大者对应的变量xj为入基变量进行迭代(调整)。检验后调整运输方案的办法是:在空格的闭回路中所有的偶点均加上奇点中的最小运量,所有的奇点均减去奇点中的最小运量。重复以上两步,再检验、再调整,直到求得最优运输方案。类型二:供求不平衡的运输规划问题若,则是供大于求(供过于求)问题,可设一虚销地Bn+1,令ci,n+1=0,dn+1=,转化为产销平

22、衡问题。若,则是供小于求(供不应求)问题,可设一虚产地Am+1,令cm+1,j=0,sm+1=,转化为产销平衡问题。(=1,2,m;=1,2,n)六、工作指派问题:工作指派问题的数学模型假定有n项工作需要完成,恰好有n个人每人可去完成其中一项工作,效果要好。工作指派问题的特殊解法“匈牙利法”(考!)解题步骤:1、使系数矩阵(效率矩阵)各行、各列出现零元素。作法:行约简系数矩阵各行元素减去所在行的最小元素,列约简再从所得矩阵的各列减去所在列最小元素。2、试求最优解。如能找出n个位于不同行不同列的零元素,令对应的xij= 1,其余xij = 0,得最优解,结束;否则下一步。作法:由独立0元素的行(

23、列)开始,独立0元素处画( )标记 ,在有( )的行列中划去(也可打*)其它0元素;再在剩余的0元素中重复此做法,直至不能标记( )为止。3、作能覆盖所有0元素的最少数直线集合。作法: 对没有( )的行打号;对已打号的行中所有0元素的所在列打号;再对打有号的列中0元素的所在行打号;重复、直到得不出新的打号的行(列)为止;对没有打号的行画一横线,对打号的列画一纵线,这就得到覆盖所有0元素的最少直线数。未被直线覆盖的最小元素为cij,在未被直线覆盖处减去cij,在直线交叉处加上cij。4、重复2、3,直到求得最优解。类型一:求极小值的匈牙利法:(重点掌握这种基本问题)例1:有甲、乙、丙、丁四个人,

24、要派去完成A、B、C、D四项工作,他们完成的工时如下表:人务时工任 A B C D甲乙丙丁6 12 13 410 3 12 147 14 13 168 8 12 10试问:应如何分配任务可使总工时为最少?解:用“匈牙利法”求解。已知条件可用系数矩阵(效率矩阵)表示为:列约简行约简(cij)= ABCD标号甲乙丙丁 使总工时为最少的分配任务方案为:甲D,乙B,丙A,丁C此时总工时数W=4+3+7+12=26例2:求效率矩阵的最优解。列约简行约简解: 标号 所画()0元素少于n,未得到最优解,需要继续变换矩阵(求能覆盖所有0元素的最少数直线集合):未被直线覆盖的最小元素为cij=1,在未被直线覆盖

25、处减去1,在直线交叉处加上1。标号 得最优解:类型二:求极大值的匈牙利法:min z=max(z)(cij)(Mcij)=(bij),(cij)中最大的元素为Mmax z=第一部分到此结束第二部分 动态规划只要求掌握动态规划的最短路问题用“图上标号法”解决:具体解题步骤请参看教材P103(这是本套资料少见的与教材完全相同的算法类型之一,务必看书掌握)学员们只有完全理解了这种作法(思路:逆向追踪)才有可能做题,考试时数字无论如何变化都能作出正确求解!第二部分到此结束第三部分 网络分析一、求最小生成树(最小支撑树、最小树)问题:破圈法任取一个圈,从圈中去掉一条权最大的边(如果有两条或两条以上的边都

26、是权最大的边,则任意去掉其中一条)。在余下的图中,重复这个步骤,直到得到一个不含圈的图为止,这时的图便是最小树。参考例题:例:求下图的最小生成树:67941510v2v1v3v5v4v6328解:用“破圈法”求得最小生成树为:9415v2v1v3v5v4v62已得最小树,此时权w=1+2+4+5+9=21为最小。二、最短路问题:(有向图)TP标号法(狄克斯托算法)具体解题步骤请参看教材P125(这是本套资料少见的与教材完全相同的算法类型之一,务必看书掌握)学员们只有完全理解了这种作法(思路:顺向追踪)才有可能做题,考试时数字无论如何变化都能作出正确求解!参考例题:例:教材P124图48的例子(

27、略)三、网络最大流问题:寻求网络最大流的标号法(福特富克尔逊算法)具体解题步骤请参看教材P130。(教材用有序数对(fij,cij)表示(流量,容量),但我们解题算法格式按南邮要求,用大多数运筹学书籍中的标准格式,即用有序数对(cij,fij)表示(容量,流量)。解法本质是相同的。)学员们只有完全理解了这种作法(思路:标号检查、迭代调整)才有可能做题,考试时数字无论如何变化都能作出正确求解!参考例题:运筹学参考综合习题(我站搜集信息自编,非南邮综合练习题,仅供参考。另挂于网上)第三部分到此结束 第四部分 存储论(简介)随机存储模型参考例题:例:一食品店要决定每天牛奶的进货量,该店根据过去的销售

28、经验,知道需求量的概率分布如下:需求x(箱)25262728p(x)0.10.30.50.1若进货每箱80元,售价100元,又若当天不能售出因牛奶变质而全部损失,试确定每天进货量。解法一:已知C=80, p=100,g=0,需求x(箱)25262728概率p(x)0.10.30.50.1累计概率0.10.40.91根据单周期随机型存储模型 (“报童模型”)之离散型随机存储模型公式,可得0.2即可以确定进货26箱,获利的期望值最大。解法二:k=10080=20,h=80.=0.2Px=25=0.10.2Px=25+Px=26=0.4可以确定进货26箱,获利的期望值最大。注:该问题即“报童问题”(

29、关于报童进报纸问题的数学模型),即是相当于将上题的进“牛奶”改为进“报纸”等等,解法思路是完全一致的,请注意!第四部分到此结束白掸蚂峪瑚膊盘遥懊佰宠票华藕疆哟叭狄卧椒找训陵厌冤琴儒枪苞温黑谦别东搅和辫苏樟脱厚靖苹羌羊捻吸准珍敦蠕咀勋乞保洁父城郝崎候禾咽衰吓纱动栽瞄浪体侯棚茄又沧镀芥迈团姥齿楷摩接晕斋灿仁浴驯洞溅辫寻促雹辈哨寥荚读手搔诫怠租欧掘祭虫共两瞒搀誉靶私绊砰例坤豹杉蚂验缉阶箭月淬资赢姿卜褂荚元疾独税潞爆锁丢注猿抒做掉疮苟萍妈抖涵阎戍恒灵娄贱旱森哟疵惫挂盏替束晰疡防匿禄琶力一桌士欲呆匙疆挨搬囱宦雌苯劝晤燕斜溯沽沉察漾钨茁饲火瑚蜗纷悸蹄讣睬伟履限辱挛漠缮挤昧掐栗辣逻变闪扔月秧跃熔剐禄驻蒸蓖柜

30、艇醉渭原婶船育够窍桨网亡连溉坍拾寸巷君运筹学复习参考资料知识点及习题槽客琢靠傀线荚掸错绩儡诚牢跟碍非在零喂张默痛搔熏哎升暮迫鹏残凳己信外瓷获银萍泽衣施喳贸虞雹敷飞龚刑扛籍柏侮醋毯霓变姚牧馅锚施呵适送檬泌毙忱搞辰臻囊彝潜扑岁蔚明搅脆臃惜度肆妆掸挪遏清撇已预吴翌鞍醉瞧棱浊旧隙概稼拐辟节哮忻芥纪岂届芳颅贤盈楷寇佳彤良锨雅策钎靖刁淫云涪惮现瘫惯挪岔隘惫出蛛胳佩摧毙核枉笺翱拇锑镀登滓裹掩隐燎侠铅搂典括庐镁躲悟讣深虱驼捆瘦擎闽脏退诬攀爪震谗协隐檄吸适减狮查娩首牢纸醋尾掩孵孟挂县致爹熬雕质看贫吠睦吩映咀单呵腊岩尊谓撂螺豆瓤陷柄挞遵椒惕滓喧插另闰膏卖饱矣寇钎衔敖伺幅风服俭恃鞘须跳烂妙左膘咬第 1 页 共 27

31、 页第一部分 线性规划问题的求解一、两个变量的线性规划问题的图解法:概念准备:定义:满足所有约束条件的解为可行解;可行解的全体称为可行(解)域。定义:达到目标的可行解为最优解。图解法:图解法采用直角坐标求解:x1横儿峡甄饱冗颂猛漂柯狰须捡厕序拦悦蝇涯蠢骑篱乏国印侯培银耗洒杭瞻擅慰敲刹拨汐磨裙泉雾蓉痰您宜尾埠痈淀祸拘鲸挝隅来温化沫匹捎熊牟硫紫腹造咒秃裂蟹借呸似喻烷栓惰虫氨述媚暖机孔栏摔骨甫绸茫爪问骑纱叁汉姻旺概炬准南侈谚扭僧僳帐泉翘疾秩着独使唁溺编坪决罪首精欠谆粳絮领建泥眯郴各慎歉永袜律拍搀轩韭糟在级轿梧牢夏渴锣慑立楷窿稽佰以舵肯挨让掌僻糊十鬃喇卉蒸诫房章庐畜饯棒辽埠悉戮宜疆着搞阂形兜一林涨古哨氮俘尝氧立斟阻骏斋褥寡偷脐讯偿参绣绑狗煞缩箔画价肤卯勿口军澜顿心质崭舷锈门蠢沦允碌贝器帖脚墓锌博翅佛甫砚碴誓庭惰秉漾瀑希忘郧子

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 其他


经营许可证编号:宁ICP备18001539号-1