1、承諾書我們仔細閱讀了全國大學生數學建模競賽章程和全國大學生數學建模競賽參賽規則（以下簡稱為競賽章程和參賽規則”，可從全國大學生數學建模競賽網站下載）。我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網上咨詢等）與隊外的任何人（包括指導教師）研究、討論與賽題有關的問題。我們知道，抄襲別人的成果是違反競賽章程和參賽規則的，如果引用別人的成果或其他公開的資料（包括網上查到的資料），必須按照規定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽章程和參賽規則，以保證競賽的公正、公平性。如有違反競賽章程和參賽規則的行為，我們將受到嚴肅處理。我們授權全國

2、大學生數學建模競賽組委會，可將我們的論文以任何形式進行公開展示（包括進行網上公示，在書籍、期刊和其他媒體進行正式或非正式發表等）。我們參賽選擇的題號是（從A/B/C/D中選擇一項填寫）：B我們的參賽報名號為（如果賽區設置報名號的話）：所屬學校（請填寫完整的全名）：許昌學院參賽隊員（打印并簽名）：1.徐晨曦2. 陳永牛3. 劉志寬指導教師或指導教師組負責人（打印并簽名）：（論文紙質版與電子版中的以上信息必須一致，只是電子版中無需簽名。以上內容請仔細核對，提交后將不再允許做任何修改。如填寫錯誤，論文可能被取消評獎資格。）日期：2016年8月27日賽區評閱編號（由賽區組委會評閱前進行編號）：編號專用

3、頁賽區評閱編號（由賽區組委會評閱前進行編號）：賽區評閱記錄（可供賽區評閱時使用）：評閱人評分備注全國統一編號（由賽區組委會送交全國前編號）：全國評閱編號（由全國組委會評閱前進行編號）：貨運列車編組運輸問題摘要對于這次我們需要求的貨車編組運輸，通過不同的情況制定最佳運送方案對于問題一，我們首先確定的是以運輸貨物最多，運輸總量最小為目標函數的雙目標優化問題，這里我們首先是將復雜的B類貨物單獨的分開來，看成是兩種類型的貨物，我們為了簡化運算我們先針對單個目標數量最多對其進行優化求解，用lingo軟件得出數量最多為24,分別有幾組數據，然后在以數量為最多的條件下為約束，求取另一個目標總重量最小，用li

4、ngo分析得出其中最小的總重量為179噸，然后再將兩者的求得結果相互結合得出，數量最多為24的情況下，總重量最小為179噸。對于問題二：問題二是下料問題，因此需要先確定可行的下料方式，即兩種車廂可行的貨物裝載方式。以每種裝載方式的使用次數為決策變量，總使用次數最少為目標函數，建立整數線性規劃模型求解。用MATLAB解得：要將貨物運輸完畢，B,C,E分別為68、50、41件時使用的最少車廂數量為25,B,C,E分別為48,42,52件時使用的最少車廂數量為21對于問題三給出了最近100天上午和下午需要運的集裝箱數目，根據所給的數據我們做出了散點圖根據散點圖并用MATLAB擬合我們發現最近100天

5、需要運的集裝箱數目符合正態分布。然后我們算出上午和下午的日利潤，再把他們相加R=R1+R2,得到每天的利潤之和。其中上午的利潤我們把它分為集裝箱可以全部運完和集裝箱運不完兩種情況分別計算，下午的同上午的，但是若上午的集裝箱沒有運完要加到下午需要運的集裝箱數目上。關鍵詞：lingo線性規劃雙目標優化Matlab正態分布、問題重述列車編組問題貨運列車編組調度的科學性和合理性直接影響著貨物運輸的效率。請根據問題設定和相關數據依次研究解決下列問題：1、假設從甲地到乙地每天有5種類型的貨物需要運輸，每種類型貨物包裝箱的相關參數見附錄一。每天有一列貨運列車從甲地發往乙地，該列車由1節I型車廂和2節II型車

6、廂編組。I型車廂為單層平板車，II型車廂為雙層箱式貨車，這兩種車廂的規格見附錄二。貨物在車廂中必須按占用車廂長度最小方式放置（比如：A類貨物占用車廂長度只能是米，不能是3米；再比如：一節車廂中B類貨物裝載量為2件時，必須并排放置占用長度米，裝載量為3件時，占用長度米），不允許貨物重疊放置；R型箱式車廂下層裝載貨物后剩余長度小于等于米，才能在上層放置貨物。試設計運輸貨物數量最多的條件下，運輸總重量最小的裝運方案。2、如果現有B,C,E三種類型的貨物各69、50、51件，試設計一個使用車廂數量最少的編組方案將貨物運輸完畢。由于整個鐵路系統I型車廂較多，要求在編組中I型車廂的數量多于II型車廂數量，

7、II型箱式車廂下層裝載貨物后剩余長度小于等于5米，才能在上層放置貨物，貨物裝車其它規則同問題1。若B,C,E三種類型的貨物各有58,42,62件，請重新編組。3、從甲地到乙地每天上午和下午各發送一列由I型車廂編組的貨運列車，每列火車開行的固定成本為30000元，每加掛一節車廂的可變成本為1500元。為了裝卸的方便，鐵路部門擬將貨物放置到長、寬、高分別為4米，3米及米的集裝箱中運輸，每個集裝箱的總重量不超過18噸，集裝箱的運費為1000元/個。每天需要運輸的集裝箱數量是隨機的，附錄三給出了過去最近100天上午和下午分別需要運輸的集裝箱的數量。上午的需求如果不能由上午開行列車運輸，鐵路部門要支付5

8、0元/個的庫存費用；下午列車開行后如果還有剩余集裝箱，鐵路部門將支付200元/個的賠償，轉而利用其它運輸方式運輸。試制定兩列火車的最佳編組方案。、問題分析問題一分析對于問題一，我們首先確定的是以運輸貨物最多，運輸總量最小為目標函數的雙目標優化問題，這里我們首先是將復雜的B類貨物單獨的分開來，看成是兩種類型的貨物，我們為了簡化運算我們先針對單個目標數量最多對其進行優化求解，用lingo軟件得出數量最多為24,分別有幾組數據，然后在以數量為最多的條件下為約束，求取另一個目標總重量最小，用lingo分析得出其中最小的總重量為179噸，然后再將兩者的求得結果相互結合得出，數量最多為24的情況下，總重量

9、最小為179噸。問題二分析問題二為求解在所有貨物都能運走的條件下使用車廂最少的情況。可以看出此題為最優化問題，也就是整數規劃問題。針對此問題可以建立模型使用matlab和lingo取得最優值。貨物類型為B,C,E,根據貨物要以占用車廂長度盡可能小的要求可知，擺放貨物C和E只有只有一種方式。由于貨物C,E寬為3m恰好等于車廂寬度，所以根據要求只能使CE的寬的方向和車寬度的方向平行，這樣才能使貨物占用長度最小。針對貨物B,已知尺寸為2.22m1.5m,寬度為1.5m,所以要使占用長度最小就要分情況而定了。當貨物B的數量為偶數時可以兩兩配對豎放，為奇數時取其中一個橫放，這樣占用長度最小。由于貨物不能

10、重疊放置，我們可以將貨物車廂中的裝載問題抽象為二維矩形件的排樣問題，只是增加了貨物總重量的上限約束。如果將一節I車廂和兩節R車廂一起進行分析，情況較為復雜，為減少計算負荷，我們先對兩種車廂各自的可行裝載方式進行分析，再將其進行組合。也就是在滿足車廂空間和重量要求的前提下，列出I車廂和R車廂所有裝載的可能情況問題三分析題目中給出了最近100天上午和下午需要運的集裝箱數目，根據所給的數據我們做出了散點圖根據散點圖并用MATLAB擬合我們發現最近100天需要運的集裝箱數目符合正態分布。然后我們算出上午和下午的日利潤，再把他們相加R=R1+R2,得到每天的利潤之和。其中上午的利潤我們把它分為集裝箱可以

11、全部運完和集裝箱運不完兩種情況分別計算，下午的同上午的，但是若上午的集裝箱沒有運完要加到下午需要運的集裝箱數目上。算出每天利潤之和，再根據我們對最近100天上午和下午需要運的集裝箱數目分析利用它符合正態分布算出需要運輸的集裝箱數量是r1的概率為f(r1),然后把它3sl們相乘，得到上午的利潤之和為B0(1000r11500s130000)P(r1)dr13s.(1650s150r130000)P(r1)dr1同理可得下午的利潤之和，然后求出利潤之和最大時所上午需要運輸的集裝箱數和下午需要運輸的集裝箱數。三、模型假設1 .貨物不能重疊放置，且不能直立放置符號名稱xijligiSiujVjRRQS

12、S2ri2 .上午運不完的集裝箱，歸到下午需要運的集裝箱的范疇3 .出于利潤最大化的考慮，發出的列車車廂數達到最大編組量且每個車廂中裝滿三個集裝箱4 .超過需求量的集裝箱，鐵路部門收不到相應的運費四、符號說明第i種貨物放入第j號車廂的數量第i種貨物占用車廂總長度第i種貨物重量第i種貨物的總數I型車廂第j種裝載方式的使用次數口型車廂第j種裝載方式的使用次數上午、下午的利潤總和上午的利潤下午的利潤上午發出的車廂數下午發出的車廂數上午需要運輸的集裝箱數下午需要運輸的集裝箱數五、模型的建立與求解問題一基于上述分析，對問題一進行模型的建立和求解。首先確定的是在運輸數量最多的條件下，我們求的是運輸的重量最

13、小，這樣我們建立的目標函數就是雙目標類型了，這里我們為了簡化模型，分別先確定數量最多的情況,然后再求解重量最小。1 .由于貨物的不能重疊放置我們這里將1節I型和2節II型分別計算各自的可以裝載的運行方案，在進行組合。這里對于B類型貨物相較于其他的復雜所以我們這里采用的方法是將其看成兩種不同的貨物具體如下兩種：然后我們分析各個車廂內的分類情況如下圖所示：圖1各個車廂內的分類情況Xij如上圖中i表示的是貨物的類型A,B1,B2,C,D,Ej表示的是車廂數量IIIIIIIII2 .考慮單個車廂的情況時，如下條件：1）貨物占用車廂的高度"廂高度考慮實際情況以及題中所給的例子，我們假設貨物不能

14、豎直放置。此時只需考慮貨物實際高度與車廂高度的關系，得到n型車廂的第二層不能放置a類和B類貨物的結論。2）貨物按占用車廂長度最小方式放置對于A,C,D,E類的貨物，他們占用車廂的最小長度就是他們的實際長度。對于B類貨物，需要進行分類討論：B1的最小長度是，B2的最小長度為3）貨物占用車廂的寬度"廂寬度貨物按占用車廂長度最小的方式放置時，恰使得A,C,D,E類貨物占用車廂的寬度等于車廂寬度，而對B類貨物進行分類討論時，已經考慮到了車廂寬度的限制，因此這一條件可以不單獨列出。4）R型車廂下層裝載貨物后剩余長度小于等于米時，才能在上層放置貨物5）貨物占用車廂總長度"廂長度6）貨物

15、總重量可廂載重量57）由于每種貨物數量有限所斗闈iSGn為第i種貨物能放的車廂的車廂數。3.對兩種車廂可行的貨物裝載方式進行組合得到目標函數：在lingo軟件中編程（源程序見附錄四）得到各種情況下的裝載數量最多方式。數量最多條件下，求總重量最小得到目標函數：由以上的目標函數在lingo中得到數量最多能裝載是24,在數量最多的境況下即24時，由lingo編程可以得出總重量最小的裝載重量，最小為179噸。（具體源程序可見附錄表五）具體裝載方案如表5-1表5-1.最優裝載方案A232NANAB010NANAC40100D00023E02310問題二模型的建立問題二為求解在所有貨物都能運走的條件下使用

16、車廂最少的情況。可以看出此題為最優化問題，也就是整數規劃問題。針對此問題可以建立模型使用matlab和lingo取得最優值。問題二中的貨物類型為B,C,E,根據貨物要以占用車廂長度盡可能小的要求可知，擺放貨物C和E只有只有一種方式。由于貨物C,E寬為3m恰好等于車廂寬度，所以根據要求只能使CE的寬的方向和車寬度的方向平行，這樣才能使貨物占用長度最小。針對貨物B,已知尺寸為2.22m1Bm,寬度為1.5m,所以要使占用長度最小就要分情況而定了。當貨物B的數量為偶數時可以兩兩配對豎放，為奇數時取其中一個橫放，這樣占用長度最小。具體放置方式如圖1.圖2.貨物B兩種放置方式由于貨物不能重疊放置，我們可

17、以將貨物車廂中的裝載問題抽象為二維矩形件的排樣問題，只是增加了貨物總重量的上限約束。如果將一節I車廂和兩節n車廂一起進行分析，情況較為復雜，為減少計算負荷，我們先對兩種車廂各自的可行裝載方式進行分析，再將其進行組合。也就是在滿足車廂空間和重量要求的前提下，列出I車廂和R車廂所有裝載的可能情況。約束條件由于重量和空間的約束條件大致和問題一相同，只改變了II型箱式車廂下層裝載貨在問題一中我們將貨物物后剩余長度小于等于5米，才能在上層放置貨物在此不再進行討論B分類為兩種類型，即b/DB2,這兩種貨物分別為B貨物單個放置和成對放置。所以B貨物放置的長度條件重新建立。求解結果在MATLAB中編程求解（源

18、程序見附錄六）,通過Excel對結果進行數據分析和整合，排除明顯劣解，得到只考慮B,C,E時I型車廂可行的貨物裝載方式22種、R車廂可行的貨物裝載方式125種（附錄七）。問題二模型的建立用uj表示只考慮B,C,E時I型車廂第j種裝載方式的使用次數，用Vj表示只考慮B,C,E時R型車廂第j種裝載方式的使用次數，則Uj（j1,2,.,22）、Vj（j1,2,.,125）是模型二的決策變量，均為非負整數。確定約束條件題目規定將一定數量的貨物B、C、E運輸完畢，既是運走的各類貨物數量大于等于現有數量，由此可得到：其中，UjVj分別是I型車廂R型車廂第j種裝載方式的使用次數，§為第i種貨物各方

19、案中貨物的運載量，Qi為第i種貨物的現有數量。I型車廂的數量多于II型車廂確定目標函數問題二要求使用車廂數量最少，即各裝載方式使用次數之和最少，所以目標函數為：根據目標函數利用Lingo求解模型（附錄八），得到B,C,E分別為69、50、51件時使用的最少車廂數量為26,B,C,E分別為58,42,62件時使用的最少車廂數量為24。具體編組方案如表5-2,表5-3。表,C,E分別為69、50、51件時使用車廂數量最少的編組方案方案次數根據表格IX1514222NANA計算此種方案裝載Y29203230貨物IIY100322400B,C,E的Y123750110數量分別為69,53,51,比現有

20、貨物數量多出3件C表,C,E分別為58、42、62件時使用車廂數量最少的編組方案格計算此種方力殺次數貨物IX15132數量分Y7521IIY1076358,43,62,Y12134貨物數根據表件C.22NANAB,C,E的3211另U為0500比現有3000量多出1案裝載問題三模型建立數據處理根據過去最近100天上午和下午需要運的集裝箱數量的數據，做出散點圖。系列一為上午的，系列二為下午的。圖一，最近100天上午和下午需要運的集裝箱數量散點圖圖3最近100天需要運的集裝箱數量沒有發現數據的明顯規律，用MATLAB進行數據分布擬合，發現兩組數據均服從正態分布，接受概率分別為、。模型的建立確定目標

21、函數因為每天上午、下午需要運輸的集裝箱數量都是隨機的，所以我們對上午、下午分別考慮，則鐵路部門的日利潤等于上午、下午的利潤之和，即目標函數為推導過程因為集裝箱和車廂的規格都固定，所以當上午發出的列車有si節車廂時，可運輸集裝箱的為3sio鐵路部門上午的利潤R1與上午需要運輸的集裝箱的數量ri有關，當ri3s1時，鐵路部門獲得最多的運費；當ri3si時，鐵路部門需要支付未被運走的集裝箱的庫存費用。即對于下午，需要運輸的集裝箱數量r2除了原來的需求，還可能包括上午剩余的集裝箱。則對上午的分析：假設上午需要運輸的集裝箱數量是ri的概率為f(ri),可以由過去的數據得到，用鐵路部門的利潤期望值來衡量利

22、潤，則3slR(S)(IOOOri1500sl30000)f(r1)(I65350r130000)f(r1)q0ri3si1即在f(ri)已知時，求si使彳3Ri最大。為了便于分析，將概率f(ri)轉化為概率密度函數P(ri),則對Ri(si)求導，并讓導數等于0,得到P(r)dr1因為0,所以將上式左右兩邊的分母都加上分子，得到由數據分析，已知ri服從正態分布，可以用正態分布的逆概率分布求解得到si對下午的分析：類似的，我們可以得到注意：下午需要運輸的集裝箱數量還包括上午未運輸完的集裝箱。模型的求解用MATLAB求解正態分布的逆概率分布(源程序見附錄四表1),解si=、s2=所以最佳編組方案

23、是上午發的列車帶41節I型車廂、下午發的列車帶38節I型車廂。六、模型評價優點1,模型把所述要求考慮的非常全面，能充分的利用數據。2,利用lingo和MATLAB把復雜的問題變得簡單化，減少了計算難度。3,通過數形結合和軟件擬合方法把看似毫無規律的數據總結出了規律。缺點1,模型涉及的問題太多，不易計算，難以理解。2,問題三中上午的數據雖然服從正態分布，但接受概率很低，對計算的結果帶來了誤差；七、參考文獻1姜啟源，謝金星，葉俊，數學模型（第四版），北京：高等教育出版社，20112張志勇，楊祖櫻，MATLAB教程R2010a，北京航空航天大學出版社，20133薛山，MATLAB2012簡明教程M,

24、北京清華大學出版社，2012;l-24BSEN581-3:2007Outdoorfurniture.Seatingandtablesforcamping,domesticandcontractuse.Mechanicalsafetyrequirementsandtestmethodsfortables5基于Matlab和Lingo的數學實驗（第一版）M,西北工業大學出版社，2009附錄1491001061329710297123124971031301461441081101061331449912898133101951001441111031061251051121501051449412

25、2148137103140121146148132120115117103931281271371001211491261301449311795911221251201359891附錄一：貨物包裝箱相關參數貨物類型長度（米）寬度（米）高度（米）重量（噸）數量A37B6C395D387E36附錄二：火車車廂相關參數車廂類型長度（米）寬度下層高度上層高度載重量（噸）I型355II型15370附錄三:近100天上午集裝箱數量:134107143143146115109139107971111411491121011111311401441309510813914211711512213612990

26、近1C)0大卜午集篇支箱數11:128137115106133569395113661551058910813110798122102102104109106971058786125124165738212182119618611362116738783136102751069312497121119103121688410811192881138578112908011675107889212511191991139811092807510185986961103851121281011029082111118128888547附錄四：Maxx11+x12+x13+x21+x22+x23+x

27、31+x32+x33+x41+x42+x43+x44+x45+x51+x52+x53+x54+x55+x61+x62+x63+x64+x65x11+x12+x13<=7x21+x22+x23+2x31+2x32+2x33<=6x41+x42+x43+x44+x45<=5x51+x52+x53+x54+x55<=7x61+x62+x63+x64+x65<=6+<=+>=+<=15+<=15+>=+<=15+<=15+21x31+9x41+8x51+<=55+21x32+9x42+8x52+9x44+8x54+<=

28、70+21x33+9x43+8x53+9x45+8x55+<=70ENDGIN24附錄五：min+21x31+21x32+21x33+9x41+9x42+9x43+9x44+9x45+8x51+8x52+8x53+8x54+8x55+x11+x12+x13<=7x21+x22+x23+2x31+2x32+2x33<=6x41+x42+x43+x44+x45<=5x51+x52+x53+x54+x55<=7x61+x62+x63+x64+x65<=6+<=+>=+<=15+<=15+>=+<=15+<=15+21x31

29、+9x41+8x51+<=55+21x32+9x42+8x52+9x44+8x54+<=70x11+x12+x13+x21+x22+x23+x31+x32+x33+x41+x42+x43+x44+x45+x51+x52+x53+x54+x55+x61+x62+x63+x64+x65=24ENDGIN24附錄六：對于I型車廂()p=;a11=0fora11=0:7;%aij中i表示貨物ABCDE,j表示層數fora21=0:6;fora31=0:5;fora41=0:7;fora51=0:6;ifa21=0;%對B進行分類l1=a11*+a31*+a41*+a51*;endifa21

30、=1,3,5;l1=a11*+(a21-1)*+a31*+a41*+a51*;endifa21=2,4,6;l1=a11*+a21*1/2*+a31*+a41*+a51*;endw1=a11*+a21*+a31*9+a41*8+a51*;%計算車廂載重ifl1<=&w1<=55&<=;p=p;a11a21a31a41a51a11+a21+a31+a41+a51w1;endendendendendendP對于II型車廂()p=;a11=0;fora11=0:5;%aij中i表示貨物ABCDE,j表示層數fora21=0:6;fora31=0:7;fora41=0

31、:5;fora51=0:5;ifa21=0;%對B進行分類l1=a11*+a31*+a41*+a51*;endifa21=1,3,5;l1=a11*+(a21-1)*+a31*+a41*+a51*;endifa21=2,4,6;l1=a11*+a21*1/2*+a31*+a41*+a51*;endfora32=0:7;fora42=0:5;fora52=0:5;I2=*a32+*a42+*a52;W2=*a11+*a21+9*(a31+a32)+8*(a41+a42)+*(a51+a52);if12<=15&w2<=70;p=p;a11a21a31a41a51a32a42a

32、52a11+a21+a31+a41+a51+a32+a42+a52w2;endendendendendendendendendendP附錄七：只考慮B,C,E時兩種車廂可行的貨物裝載方式I型車廂編PBCE10042013302340325041605170608103911310122111311214113203142121522216230173021831119320204012141022500n型車廂編pB下C下D下C上D上100405200414300422400431500440600504700513800521900530100140411014121201421130143

33、014015031501511160152017023031802312190232120023302102402220241123024202402501250251026032032703212280322129032303003302310331132033203303401340341035035003604202370421138042203904301400431041044004205102430511144052014505210460530047060024806010490610150062005107000521040353104125410421551043056105

34、025710511581052059113036011312611132162113306311402641141165114206611501671151068123026912311701232071124017212410731250074132027513211761330177133107813400791410280141108114201821430083150018415010851510086160008720402882041189204209020501912051092213029321311942140195214109621500972220298222109922

35、3011002240010123200102241001032500010430302105303101063040110730500108312011093121011031300111322001123310011334000114402011154021011640300117411011184111011941200120421001214300012250101123501101245020012551100附錄八:minx1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+x21+x22+Y1+

36、Y2+Y3+Y4+Y5+Y6+Y7+Y8+Y9+Y10+Y11+Y12+Y13+Y14+Y15+Y16+Y17+Y18+Y19+Y20+Y21+Y22+Y23+Y24+Y25+Y26+Y27+Y28+Y29+Y30+Y31+Y32+Y33+Y34+Y35+Y36+Y37+Y38+Y39+Y40+Y41+Y42+Y43+Y44+Y45+Y46+Y47+Y48+Y49+Y50+Y51+Y52+Y53+Y54+Y55+Y56+Y57+Y58+Y59+Y60+Y61+Y62+Y63+Y64+Y65+Y66+Y67+Y68+Y69+Y70+Y71+Y72+Y73+Y74+Y75+Y76+Y77+Y78+

37、Y79+Y80+Y81+Y82+Y83+Y84+Y85+Y86+Y87+Y88+Y89+Y90+Y91+Y92+Y93+Y94+Y95+Y96+Y97+Y98+Y99+Y100+Y101+Y102+Y103+Y104+Y105+Y106+Y107+Y108+Y109+Y110+Y111+Y112+Y113+Y114+Y115+Y116+Y117+Y118+Y119+Y120+Y121+Y122+Y123+Y124+Y125SUBJECTTOb)+0丫36+0丫37+0丫38+0丫39+0丫40+0丫41+0丫42+0丫43+0丫44+0丫45+0丫46+0丫47+0丫48+0丫49+0丫50+

38、0丫51+1Y52+1Y53+1Y54+1Y55+1Y56+1Y57+1Y58+1Y59+1Y60+1Y61+1Y62+1Y63+1Y64+1Y65+1Y66+1Y67+1Y68+1Y69+1Y70+1Y71+1Y72+1Y73+1Y74+1Y75+1Y76+1Y77+1Y78+1Y79+1Y80+1Y81+1Y82+1Y83+1Y84+1Y85+1Y86+2Y87+2Y88+2Y89+2Y90+2Y91+2Y92+2Y93+2Y94+2Y95+2Y96+2Y97+2Y98+2Y99+2Y100+2Y101+2Y102+2Y103+3Y104+3Y105+3Y106+3Y107+3Y108+3

39、Y109+3Y110+3Y111+3Y112+3Y113+4Y114+4Y115+4Y116+4Y117+4Y118+4Y119+4Y120+4Y121+5Y122+5Y123+5Y124+5Y125+0x1+0x2+0x3+0x4+0x5+0x6+0x7+1x8+1x9+1x10+1x11+1x12+2x13+2x14+2x15+2x16+3x17+3x18+3x19+4x20+4x21+5x22>=69c)0Y1+1Y2+2Y3+3Y4+4Y5+0Y6+1Y7+2Y8+3Y9+1Y10+2Y11+3Y12+4Y13+1Y14+2Y15+3Y16+2Y17+3Y18+4Y19+5Y20

40、+2Y21+3Y22+4Y23+2Y24+3Y25+3Y26+4Y27+5Y28+6Y29+3Y30+4Y31+5Y32+3Y33+4Y34+3Y35+4Y36+5Y37+6Y38+4Y39+5Y40+4Y41+5Y42+6Y43+5Y44+6Y45+5Y46+6Y47+7Y48+6Y49+6Y50+7Y51+0Y52+1Y53+2Y54+3Y55+0Y56+1Y57+2Y58+1Y59+2Y60+3Y61+4Y62+1Y63+2Y64+3Y65+1Y66+2Y67+2Y68+3Y69+4Y70+2Y71+3Y72+2Y73+3Y74+4Y75+3Y76+4Y77+3Y78+4Y79+5Y80

41、+4Y81+4Y82+5Y83+6Y84+5Y85+6Y86+0Y87+1Y88+2Y89+0Y90+1Y91+1Y92+2Y93+1Y94+2Y95+1Y96+2Y97+3Y98+2Y99+2Y100+3Y101+4Y102+5Y103+0Y104+1Y105+0Y106+0Y107+1Y108+2Y109+1Y110+2Y111+3Y112+4Y113+0Y114+1Y115+0Y116+1Y117+2Y118+1Y119+2Y120+3Y121+0Y122+1Y123+0Y124+1Y125+0x1+1x2+2x3+3x4+4x5+5x6+6x7+0x8+1x9+2x10+3x11+4x12+0x13+1x14+2x15+3x16+0x17+1x18+2x19+0x20+1x21+0x