算法模型知識點(diǎn)第02章整數(shù)

第二章整數(shù)規(guī) 于整數(shù)線性規(guī)劃。目前還沒有法能有效地求解一切整數(shù)規(guī)劃。1o變量全限制為整數(shù)時(shí)，稱純（完全）1原線性規(guī)劃為 zx12x14x2 x1 x2 x10x24minz42原線性規(guī)劃為 zx12x14x2 x1 x2 x10,x22minz2。 對有約束條件的最優(yōu)化問題（其可行解為有限數(shù)）Landg和Dakin等人。由于這方法靈活且便于用計(jì)算機(jī)求解，所以現(xiàn)在它已是zAz*z。分枝定界法就Bzzz*?，F(xiàn)用下例來3 z40x19x17x27x20x x14.8092,x21.8168,zzAz*z即0z*356x1x2當(dāng)前均為非整數(shù)，故不滿足整數(shù)要求，任選一個(gè)進(jìn)行分枝。設(shè)選x1[4.8092]4，x1[4.8092]1問題B1： z40x19x17x27x20x 問題B2： z40x19x17x27x20x

再定界：0z*349。對問B1再進(jìn)行分枝得問題B11B12B11B12

x14,x22,z11x11.43,x23.00,z12對問B2再進(jìn)行分枝得問題B21B22B21 x15.44,x21.00,z22將B21,B22剪枝。x14,x22,z*A，將與它相應(yīng)的線性規(guī)劃問題稱為問題B。Axj0,j1,Ln，試探，求得其目標(biāo)函數(shù)值，并記作zz*表示問題A的最優(yōu)目標(biāo)函數(shù)值；這時(shí)有zz*以[bj]表示小于bj的最大整數(shù)。構(gòu)造兩個(gè)約束條件xjbj]和xjbj將這兩個(gè)約束條件，分別加入問題B，求兩個(gè)后繼規(guī)劃B1B2。不考慮整數(shù)條件值為最大者作為新的下界z，若無作用z不變。jz，且不符合整數(shù)條件，則重復(fù)第一步驟。一直到最后得到z*z為止。得最優(yōu)整數(shù)解x*,j1,Ln。j§301為01變量，或稱二進(jìn)制變量。xj01這個(gè)條件可由下述約束條件：0xj1先介紹引入01變量的實(shí)際問題，再研究解法。引入01變量的實(shí)際例 某公司擬在市東、西、南三區(qū)建立門市部。擬議中有7個(gè)位置（點(diǎn)Ai(i1,2,L,7)在東區(qū)。A1A2A3三個(gè)點(diǎn)中至多選兩個(gè)；在西區(qū)。由A4A5兩個(gè)點(diǎn)中至少選一個(gè)；在南區(qū)，由A6A7兩個(gè)點(diǎn)中至少選一個(gè)。i點(diǎn)，設(shè)備投資估計(jì)為bi元，每年可獲利潤估計(jì)為ci元，但投資總額不能B元。問應(yīng)選擇哪幾個(gè)點(diǎn)可使年利潤為最大？解題時(shí)先引入01xi(i令xxi

i1,2,L,7

zcibixix1x2x3xx x6x7

xi0或5x14x224或7x13x2455x14x2247x13x245(1

x10或500x1500yx1800y0或ai1x1Lainxn i1,2,L,為了保證這m個(gè)約束條件只有一個(gè)起作用，我們引入m個(gè)01yi(i1,2,Lai1x1LainxnbiyiMy1Lymm1

i1,2,L, 代入（1），就只有ii*的約束條件起作用，而別的式子都是多余的。關(guān)于固定費(fèi)用的問題（FixedCost5某工廠為了生產(chǎn)某種產(chǎn)品，有幾種不同的生產(chǎn)方式可供選擇，如選定的生產(chǎn)Pkjcjxj 當(dāng)xj j 當(dāng)xj在構(gòu)成目標(biāo)函數(shù)時(shí)，為了統(tǒng)一在一個(gè)問題中討論，現(xiàn)引入01yjyj

當(dāng)采用第j種生產(chǎn)方式，即xj0時(shí)，當(dāng)不采用第j種生產(chǎn)方式，即xj0時(shí)

minz(k1y1c1x1)(k2y2c2x2)(k3y3c3x3yjxjyjM j 01型整數(shù)規(guī)劃解法之一（過濾隱枚舉法解01變量取值的2nn較大（n100方法稱為隱枚舉法（ImplicitEnumeration下面舉例說明一種解01型整數(shù)規(guī)劃的隱枚舉法。例 z3x12x2x12x2x3x4xx x1x24xx x1x2x30或個(gè)可行解，且z3?！?法（隨機(jī)取樣法7Maxzx2x23x24x22x 8x12x23x3x40xi (ixxxxx x2x2xx6x8002x x3x45x5用去隨機(jī)計(jì)算106個(gè)點(diǎn)，便可找到滿意解，那么這種方法的可信度究竟怎樣下面就分析隨機(jī)取樣106個(gè)點(diǎn)計(jì)算時(shí)，應(yīng)用概率理論來估計(jì)一下可信度假設(shè)目標(biāo)函數(shù)落在高值區(qū)的概率分別為0.01，0.00001，則當(dāng)計(jì)算 個(gè)點(diǎn)后，1 0.99L99(100多位10.999991000000 function[f,g]=mengte(x);fori=1:10^6ifp0<=fifsum(g<=0)==4ifa=111111221216000115;§5直接利用的函數(shù)，必須利用編程實(shí)現(xiàn)分枝定界解法和割平面解法。但對于指派問題等01整數(shù)規(guī)劃問題，可以直接利用的函數(shù)bintprog進(jìn)行求解。8874222978423697 555 c=[382103;87297;642784235;9106910];fori=1:5a(i,(i-x15x23x32x44x51121。求解的LINGO程序如下：c=3821087296427842391069§69某公司用兩種原油（AB）混合加工成兩種汽油（甲和乙。甲、乙兩種汽油含原油的最低比例分5060％，每噸售價(jià)分別為48005600元。該公司現(xiàn)有原油AB的庫存量分別為5001000還可以從市場上買到不超過1500噸的原油A。原油A的市場價(jià)為：量不超過500噸時(shí)的單價(jià)為10000元/噸；量超過500噸單不超過1000噸時(shí)，超過500噸的部分8000元/噸；量超過1000噸采購價(jià)利潤為銷售汽油的收入與原油A的支出之差這里的難點(diǎn)在于原油A的采購：c(x)

0x500x1000x

設(shè)原油A用于生產(chǎn)甲、乙兩種汽油的數(shù)量分別為x11x12，原油B用于生產(chǎn)甲、乙兩種汽油的數(shù)量分別為x21和x22，則總的收入為4.8(x11x215.6(x12x22（千 z4.8(x11x21)5.6(x12x22) 限制，以及兩種汽油含原油A的比例限制，它們表示為x11x12500xx21x221000x

x11,x12,x21,x22,x

)且，對于這樣用分段函數(shù)定義的c(x)，一般的非線性規(guī)劃軟件也難以輸入和求解。能一個(gè)自然的想法是將原油A的采x分解為三個(gè)x1x2x3分別表示以價(jià)10千元噸、8千元噸、6千元噸采購的原油A的噸數(shù)，總支出為10x18x26x3xx1x2maxz4.8(x11x21)5.6(x12x22)(10x18x26x3

應(yīng)該注意到，只有當(dāng)以10千元/噸的價(jià)格x1500（噸）時(shí)，才能以8千元/ x2(0)，這個(gè)條件可以表示為(x1500)x2 x3(0)(x2500)x3此外x1x2x3的取值范圍0x1,x2,x3

15（16型輸入LINGO軟件如下：c=108可以用菜單命令“LINGO|Options”在“GlobalSolver”選項(xiàng)卡上啟動(dòng)全局優(yōu)化選型，并運(yùn)行上述程序求得全局最有解：1000噸原油A，與庫存的500噸原油A和1000B一起，共生產(chǎn)2500噸汽油乙，利潤為5000（千元。z11z21z3110千元/噸、8千元/噸、6千元/噸的價(jià)格采購原油A，則約束（15）和（16）可以替換為500z2x1500z1 500z3x2500z2 x3500z3 z1,z2,z30或 （7）～（14@for(var1(i)|i#lt#3:500*z(i+1)<x(i);x(i)<500*z(i));c=108直接處理分段線性函數(shù)c(x。式（5）表示的函數(shù)c(x)1x軸上的分點(diǎn)為b10，b2500，b31000x屬于第1個(gè)小區(qū)間[b1b2時(shí)xw1b1w2b2w1w21w1w20，因?yàn)閏(x在[b1b2上是線性1c(x所以c(xw1c(b1w2c(b2。同樣，當(dāng)x屬于第2個(gè)小區(qū)間[b2b3時(shí)， [b3b4]時(shí)，xw3b3w4b4，w3w41，w3w40x在第k個(gè)小區(qū)間時(shí)，zk1，否則zk0。這樣w1w2w3w4z1z2z3應(yīng)足w1z1，w2z1z2，w3z2z3，w4w1w2w3w41，wk0（k1,2,3,4z1z2z31z1z2z30或x和c(xxw1b1w2b2w3b3w4b4500w2 c(x)w1c(b1)w2c(b2)w3c(b3)w4c(b4)5000w2 （6）～（13z=,,,0增加的虛擬變量z(4)=0;@for(var(i)|i#ne#1:w(i)<z(i-1)+z(i));設(shè)一個(gè)nf(x的分點(diǎn)為b1Lbnbn1wkxf(x為xkf

)

f(b

w1 kwk和0-1zk滿

z1，z1z2，…，wnzn1zn，wn1 zzLz1，z0 w1w2Lwn11，wk0（k1,2,L,n Maxzx1x9x 14 2xx x1x20,x1x2試將下述非線性的01規(guī)劃問題轉(zhuǎn)換成線性的01規(guī)劃問題 zx1x1x2x32x13x2x3x xj105個(gè)鉆井探油，使總的鉆探費(fèi)用為最小。若10個(gè)井位的代號為s1s2L,s10，相應(yīng)的鉆探費(fèi)用為c1c2L,c10，并且或選s1s7s9選擇s3s4就不能s5s5s6s7s8成十五個(gè)小區(qū)，每個(gè)小區(qū)內(nèi)標(biāo)有三個(gè)數(shù)字，分別表示該區(qū)域的海拔高度h(m)、面積S(km2和被完全淹沒時(shí)土地、房屋和財(cái)產(chǎn)等損失總數(shù)k（百萬元各小區(qū)間有相對高度為1.2m的小堤相互，例如第一區(qū)和第二區(qū)之間事實(shí)上有海拔5.2m小堤。當(dāng)洪水淹沒一個(gè)小區(qū)且洪水高于該小區(qū)高度pmf(k,kpf(k,p)

0pp求整個(gè)區(qū)域全部受損失的最小洪水量Q當(dāng)洪水量為Q6Q3