• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    兩級定位-路徑問題的變鄰域人工蜂群算法

    2014-07-07 01:49:00陳久梅龔英
    關(guān)鍵詞:中轉(zhuǎn)站鄰域蜂群

    陳久梅,龔英

    1.重慶工商大學(xué)電子商務(wù)及供應(yīng)鏈系統(tǒng)重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400067

    2.重慶工商大學(xué)商務(wù)策劃學(xué)院,重慶 400067

    兩級定位-路徑問題的變鄰域人工蜂群算法

    陳久梅1,2,龔英1,2

    1.重慶工商大學(xué)電子商務(wù)及供應(yīng)鏈系統(tǒng)重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400067

    2.重慶工商大學(xué)商務(wù)策劃學(xué)院,重慶 400067

    建立了兩級定位-路徑問題的數(shù)學(xué)模型,提出了一種求解該問題的人工蜂群算法。針對該算法容易出現(xiàn)早熟現(xiàn)象,將近年來國外出現(xiàn)的一種新穎的軌跡式啟發(fā)式算法——變鄰域搜索融入其中,由此提出三種變鄰域搜索策略?;诓煌冟徲蛩阉鞑呗缘娜斯し淙核惴ê腿斯~群算法的求解效果進(jìn)行對比仿真。實(shí)驗(yàn)結(jié)果表明,變鄰域人工蜂群算法能有效求解兩級定位-路徑問題。

    兩級定位-路徑問題;人工蜂群算法;變鄰域搜索;物流;配送

    1 引言

    近年來,物流進(jìn)城難、成本高、送錯(cuò)貨等問題屢見不鮮。為了解決這些問題,提高配送能力,越來越多的企業(yè)在城市周邊設(shè)立中轉(zhuǎn)站,用于進(jìn)城物品的分揀、拼裝和轉(zhuǎn)運(yùn)。由此形成由“配送中心-中轉(zhuǎn)站-顧客”為節(jié)點(diǎn)的兩級物流網(wǎng)絡(luò)。該網(wǎng)絡(luò)中涉及到的配送中心、中轉(zhuǎn)站的位置、數(shù)量,兩級網(wǎng)絡(luò)各節(jié)點(diǎn)間的任務(wù)分配,以及車輛的調(diào)度與配送路徑等問題的集成,即兩級定位-路徑問題(two-Echelon Location-Routing Problem,2E-LRP)[1]。該問題自Jacobsen等[2]提出以來,逐漸成為研究熱點(diǎn)。

    兩級定位-路徑問題是一個(gè)NP-Hard問題,目前的研究主要集中在啟發(fā)式求解算法方面。Ambrosino等[3]用大規(guī)模鄰域搜索算法和重新定位機(jī)制求解只有一個(gè)配送中心的帶異質(zhì)車隊(duì)的兩級定位-路徑問題。Sterle等[4]為設(shè)施點(diǎn)有容量限制的兩級定位-路徑問題建立整數(shù)規(guī)劃模型,使用禁忌搜索算法進(jìn)行求解。Jin等[5]使用遺傳算法和禁忌搜索算法的混合算法求解兩級定位-路徑問題。Boccia等[6]首先將兩級定位-路徑問題分解為兩個(gè)一級定位-路徑問題,每個(gè)一級定位-路徑問題進(jìn)一步分解為設(shè)施定位問題和車輛路徑問題,之后使用多階段迭代嵌入禁忌搜索算法進(jìn)行求解。王紹仁等[7]針對震后緊急響應(yīng)階段應(yīng)急物流系統(tǒng)優(yōu)化中的兩級定位-路徑問題,提出了一種基于兩階段分段思想的“三角”啟發(fā)式算法進(jìn)行求解。Nguyen等[8-9]采用貪婪隨機(jī)自適應(yīng)搜索、路徑重連及學(xué)習(xí)過程求解兩級定位-路徑問題。

    人工蜂群算法(Artificial Bee Colony,ABC)是由Karaboga[10]于2005年基于蜜蜂覓食原理提出的一種群智能優(yōu)化算法。它模擬了蜂群依各自分工不同協(xié)作采蜜,交換蜜源信息尋找最優(yōu)蜜源的群體行為。它已被廣泛應(yīng)用于求解旅行商問題、車輛路徑問題及設(shè)施定位等NP-Hard問題,均取得較好的求解效果。該算法具有良好的求取全局極值能力,并具有對初值、參數(shù)選擇不敏感、魯棒性強(qiáng)、簡單易實(shí)現(xiàn)等優(yōu)點(diǎn),但卻存在搜索早期過早收斂、搜索后期收斂較慢等問題[10]。近年來,國外出現(xiàn)了一種新穎的軌跡式啟發(fā)式算法——變鄰域搜索(Variable Neighborhood Search,VNS)。該算法的基本思想是在獲得局部最優(yōu)解的基礎(chǔ)上,通過系統(tǒng)地改變鄰域結(jié)構(gòu)集來拓展搜索范圍找到另一個(gè)局部最優(yōu)解。該算法思想簡單,容易實(shí)現(xiàn),算法結(jié)構(gòu)與問題無關(guān),適合于各類優(yōu)化問題,它可被方便地嵌入到其他啟發(fā)式算法中[11]。在此,針對人工蜂群算法存在的不足,設(shè)計(jì)嵌入變鄰域搜索的人工蜂群算法對兩級定位-路徑問題進(jìn)行求解。

    2 問題描述及數(shù)學(xué)建模

    2.1 問題描述

    兩級定位-路徑問題是傳統(tǒng)定位-路徑問題的拓展問題,是兩級物流網(wǎng)絡(luò)中的定位-路徑問題,其中涉及到配送中心和中轉(zhuǎn)站兩類物流設(shè)施的選址、設(shè)施與需求點(diǎn)之間的配給以及兩級物流網(wǎng)絡(luò)中車輛的行駛路徑。具體可描述如下:由一系列潛在配送中心、潛在中轉(zhuǎn)站及已知顧客所構(gòu)成的兩級物流網(wǎng)絡(luò)中,物品從配送中心經(jīng)中轉(zhuǎn)站到達(dá)顧客處,在滿足設(shè)施容量限制、車輛容量限制等約束條件下,確定哪些潛在設(shè)施點(diǎn)開放、開放的配送中心應(yīng)分別為哪些開放的中轉(zhuǎn)站提供物品、開放的中轉(zhuǎn)站應(yīng)分別為哪些顧客提供物品、兩級物流網(wǎng)絡(luò)中各需多少輛車及每輛車的行駛路線,以實(shí)現(xiàn)總成本最小,即兩級設(shè)施點(diǎn)的開放成本、車輛的使用成本及車輛的運(yùn)行成本之和最小。該問題的假設(shè)條件如下:

    (1)兩級物流網(wǎng)絡(luò)中所使用的車輛、兩級設(shè)施均有容量限制。

    (2)每個(gè)顧客只能接受來自一個(gè)已開放中轉(zhuǎn)站的物品,每個(gè)開放的中轉(zhuǎn)站只能接受來自一個(gè)已開放配送中心的物品。

    (3)二級物流網(wǎng)絡(luò)中每個(gè)顧客只能接受一輛車、一次服務(wù),一級物流網(wǎng)絡(luò)中每個(gè)開放的中轉(zhuǎn)站只能接受一輛車、一次服務(wù)。

    (4)物品從配送中心出發(fā),必須經(jīng)過中轉(zhuǎn)站中轉(zhuǎn)才能到達(dá)顧客處,不考慮配送中心直接送達(dá)顧客的情況,且同級設(shè)施之間無貨物流動。

    (5)物流網(wǎng)絡(luò)中每輛車僅完成一次旅行,且二級物流網(wǎng)絡(luò)中車輛從某中轉(zhuǎn)站出發(fā)完成行駛路徑上所有顧客需求服務(wù)后必須回到該中轉(zhuǎn)站,一級物流網(wǎng)絡(luò)中車輛從某配送中心出發(fā)完成行駛路徑上所有中轉(zhuǎn)站需求服務(wù)后必須回到該配送中心。

    (6)車輛行駛成本行駛距離成正比。

    2.2 數(shù)學(xué)建模

    2.2.1 符號說明

    集合:ND為潛在配送中心集合;NS為潛在中轉(zhuǎn)站集合;NC為顧客集合;E1為一級物流網(wǎng)絡(luò)邊的集合;E2為二級物流網(wǎng)絡(luò)邊的集合;P,R,T?ND∪NS∪NC,Z?E1∪E2。

    參數(shù):QDi為配送中心i的容量;FDi為配送中心i的開放成本;QSi為中轉(zhuǎn)站i的容量;FSi為中轉(zhuǎn)站i的開放成本;QVi為i級物流網(wǎng)絡(luò)中車輛的容量(i=1,2);FVi為i級物流網(wǎng)絡(luò)中車輛的固定使用成本(i=1,2);cij為邊(i,j)上車輛行駛成本;di為顧客i的需求量。

    決策變量:uij為一級物流網(wǎng)絡(luò)中邊(i,j)被使用一次為1,否則為0;vij為一級物流網(wǎng)絡(luò)中邊(i,j)被使用兩次為1,否則為0;xij為二級物流網(wǎng)絡(luò)中邊(i,j)被使用一次為1,否則為0;yij為二級物流網(wǎng)絡(luò)中邊(i,j)被使用兩次為1,否則為0;mij為中轉(zhuǎn)站i分配給顧客j為1,否則為0;nij為配送中心i分配給中轉(zhuǎn)站j為1,否則為0;oi為配送中心i開放為1,否則為0;pi為中轉(zhuǎn)站i開放為1,否則為0。

    2.2.2 模型建立

    上述模型中,目標(biāo)函數(shù)由八部分組成,其中,第一、二項(xiàng)分別為配送中心和中轉(zhuǎn)站的開放成本,第三、四項(xiàng)分別為一、二級網(wǎng)絡(luò)中車輛固定使用成本。第五、六項(xiàng)為一級網(wǎng)絡(luò)中車輛行駛成本,第七、八項(xiàng)為二級網(wǎng)絡(luò)中車輛行駛成本。

    式(2)~(27)為該問題的約束條件,式(2)表示二級物流網(wǎng)絡(luò)中對任意一個(gè)顧客,連接該顧客的邊有兩條;式(3)表示一級物流網(wǎng)絡(luò)中如果某中轉(zhuǎn)站開放,那么連接該中轉(zhuǎn)站的邊有兩條,否則沒有邊連接該中轉(zhuǎn)站;式(4)表示二級物流網(wǎng)絡(luò)中車輛運(yùn)輸量不能超過車輛的容量,當(dāng)K(NC′)=1時(shí),表示子回路消除;式(5)表示一級物流網(wǎng)絡(luò)中車輛運(yùn)輸量不能超過車輛的容量,當(dāng)K(NS′)=1時(shí),表示子回路消除;式(6)表示分配給某中轉(zhuǎn)站的顧客需求量之和不超過該中轉(zhuǎn)站的容量;式(7)表示分配給某配送中心的中轉(zhuǎn)站貨物量之和不超過該配送中心的容量;式(8)表示二級物流網(wǎng)絡(luò)中,如果某中轉(zhuǎn)站開放,那么連接該中轉(zhuǎn)站與某顧客的被使用一次和兩次的邊不能同時(shí)存在,否則沒有邊連接該中轉(zhuǎn)站;式(9)表示一級物流網(wǎng)絡(luò)中,如果某配送中心開放,那么連接該配送中心與某中轉(zhuǎn)站的用一次和兩次的邊不能同時(shí)存在,否則沒有邊連接該配送中心;式(10)表示顧客數(shù)不小于2的二級物流網(wǎng)絡(luò)中,從某中轉(zhuǎn)站出發(fā)的車輛,完成配送后必須回到該中轉(zhuǎn)站;式(11)表示中轉(zhuǎn)站數(shù)不小于2的一級物流網(wǎng)絡(luò)中,從某配送中心出發(fā)的車輛,完成配送后必須回到該配送中心;式(12)表示二級物流網(wǎng)絡(luò)中對任意顧客,連接該顧客與中轉(zhuǎn)站的被使用一次和兩次的邊的數(shù)量之和小于等于1;式(13)表示一級物流網(wǎng)絡(luò)中,如果中轉(zhuǎn)站開放,連接該中轉(zhuǎn)站與配送中心的被使用一次和兩次的邊的數(shù)量之和小于等于1,如果中轉(zhuǎn)站不開放,沒有邊連接;式(14)表示對任意一個(gè)顧客,有且僅有一個(gè)中轉(zhuǎn)站配給;式(15)表示對任意一個(gè)中轉(zhuǎn)站,若開放則有且僅有一個(gè)配送中心配給,否則無配送中心配給;式(16)表示如果某中轉(zhuǎn)站分配給某顧客,則該中轉(zhuǎn)站一定開放,如果某中轉(zhuǎn)站未開放,則一定不能分配給任何顧客;式(17)表示如果某配送中心分配給某中轉(zhuǎn)站,則該配送中心一定開放,如果某配送中心未開放,則一定不能分配給任何中轉(zhuǎn)站;式(18)表示若某顧客沒有分配給某中轉(zhuǎn)站,則兩者之間的邊不能被訪問,若某顧客分配給某中轉(zhuǎn)站,則連接兩者的被使用一次和使用兩次的邊不能同時(shí)存在;式(19)表示若某中轉(zhuǎn)站沒有分配給某配送中心,則該兩者之間的邊不能被訪問,若某中轉(zhuǎn)站分配給某配送中心,則連接兩者的被使用一次和使用兩次的邊不能同時(shí)存在;式(20)~(27)為決策變量的取值范圍。

    3 算法設(shè)計(jì)

    人工蜂群算法模擬蜂群的采蜜行為,通過不同角色蜜蜂間的信息交流、角色轉(zhuǎn)換和分工協(xié)作來找到問題的最優(yōu)解。人工蜂群由引領(lǐng)蜂、跟隨蜂和偵察蜂組成。引領(lǐng)蜂在蜜源處采蜜并提供所記憶的蜜源信息;跟隨蜂以一定的概率選擇某個(gè)引領(lǐng)蜂去蜜源處采蜜;偵察蜂負(fù)責(zé)尋找新蜜源。其中蜜源的位置代表優(yōu)化問題的可行解,蜜源的花蜜含量代表解的質(zhì)量,稱為適應(yīng)值。

    ABC算法工作原理如下:引領(lǐng)蜂尋找蜜源,跟隨蜂選擇并跟隨引領(lǐng)蜂在相應(yīng)的蜜源位置附近進(jìn)行鄰域搜索,若能搜索到適應(yīng)值更高的新蜜源,則更新當(dāng)前蜜源。如果某個(gè)蜜源的位置不能在預(yù)先設(shè)定的循環(huán)次數(shù)內(nèi)得到進(jìn)一步改進(jìn),則放棄該蜜源,將其對應(yīng)的引領(lǐng)蜂轉(zhuǎn)換為偵察蜂。被放棄的蜜源將由偵察蜂找到的新蜜源所取代,此時(shí)偵察蜂轉(zhuǎn)換為與新蜜源對應(yīng)的引領(lǐng)蜂。跟隨蜂再次選擇引領(lǐng)蜂跟隨。如此反復(fù)迭代,直到達(dá)到結(jié)束條件為止。

    3.1 算法總體框架

    針對兩級定位-路徑問題,本文在引領(lǐng)蜂及跟隨蜂所進(jìn)行的鄰域搜索過程中嵌入變鄰域搜索算法,即在鄰域搜索過程中系統(tǒng)地改變鄰域結(jié)構(gòu)集來拓展搜索范圍。具體步驟如下:

    步驟1初始化種群參數(shù)SN,limit,此時(shí)所有蜜蜂為偵察蜂。

    步驟2生成SN個(gè)初始解,按目標(biāo)函數(shù)值從小到大排序,前一半的解設(shè)為當(dāng)前解,此時(shí)一半偵察蜂轉(zhuǎn)換為引領(lǐng)蜂并各自攜帶一個(gè)當(dāng)前解,一半偵察蜂轉(zhuǎn)換為跟隨蜂。

    步驟3計(jì)算每個(gè)引領(lǐng)蜂攜帶解被選擇的概率Pi,引領(lǐng)蜂根據(jù)Pi招募跟隨蜂。

    步驟4每一只引領(lǐng)蜂及其跟隨蜂對其攜帶的當(dāng)前解進(jìn)行變鄰域搜索,若搜索到比當(dāng)前解更優(yōu)的解,則替換當(dāng)前解;否則,繼續(xù)進(jìn)行變鄰域搜索,直到達(dá)到一定次數(shù)limit,若此時(shí)仍然未搜索到更優(yōu)解,則該引領(lǐng)蜂轉(zhuǎn)換為偵察蜂。

    步驟5偵察蜂進(jìn)行全局搜索得到新解,并用新解替換當(dāng)前解,此時(shí)偵察蜂變?yōu)橐I(lǐng)蜂。

    步驟6當(dāng)所有引領(lǐng)蜂對其攜帶當(dāng)前解進(jìn)行變鄰域搜索后,將種群最優(yōu)解進(jìn)行保存;若滿足程序終止條件,則程序退出,輸出種群最優(yōu)解,否則轉(zhuǎn)步驟3。

    3.2 初始解的生成

    在此將兩級定位-路徑問題分為兩級,每級進(jìn)一步劃分為設(shè)施定位和車輛路徑兩個(gè)子問題。設(shè)施定位問題采用隨機(jī)最近鄰策略求解,車輛路徑問題采用經(jīng)典的C-W節(jié)約算法進(jìn)行求解。具體來說,在由中轉(zhuǎn)站與顧客構(gòu)成的二級網(wǎng)絡(luò)中,首先隨機(jī)選擇未開放的中轉(zhuǎn)站并將其開放;之后,從距離最近的r個(gè)未歸給任何中轉(zhuǎn)站的顧客中隨機(jī)選擇某個(gè)顧客,并將其歸給該中轉(zhuǎn)站,重復(fù)此步,直到所有顧客均已歸入某中轉(zhuǎn)站或超過該中轉(zhuǎn)站容量限制為止;若存在未歸入任意中轉(zhuǎn)站的顧客,則隨機(jī)選擇未開放的中轉(zhuǎn)站并將其開放;按上述方法將相應(yīng)顧客逐一歸入;重復(fù)此步,直到所有顧客均歸入某相應(yīng)中轉(zhuǎn)站為止,從而完成二級網(wǎng)絡(luò)中設(shè)施定位問題的求解。按類似的方法,完成由配送中心與中轉(zhuǎn)站構(gòu)成的一級網(wǎng)絡(luò)中設(shè)施點(diǎn)定位問題的求解。最后,分別針對開放的配送中心、中轉(zhuǎn)站,根據(jù)車輛容量限制,使用C-W算法求解兩級網(wǎng)絡(luò)中的車輛路徑問題。

    3.3 鄰域解的產(chǎn)生

    兩級定位-路徑問題,其產(chǎn)生鄰域解的考慮對象有兩級物流網(wǎng)絡(luò)中的顧客、中轉(zhuǎn)站及路徑三類。以顧客為對象的鄰域解采用單個(gè)顧客位置的移動(記為N1)及兩個(gè)顧客之間位置的交換(記為N2)兩種方式來產(chǎn)生;以路徑為對象的鄰域解采用常用的2-opt方式(記為N3)產(chǎn)生;以中轉(zhuǎn)站為對象的鄰域解采用單個(gè)中轉(zhuǎn)站開放與否的狀態(tài)改變(記為N4)及兩個(gè)狀態(tài)不同的中轉(zhuǎn)站之間進(jìn)行狀態(tài)交換(記為N5)兩種方式產(chǎn)生。

    據(jù)此提出三種變鄰域搜索策略:

    變鄰域搜索策略1(記為S1),順序執(zhí)行各鄰域搜索,若出現(xiàn)比當(dāng)前解更優(yōu)的鄰域解,則回到N1繼續(xù)執(zhí)行,一旦出現(xiàn)比該鄰域解更差的鄰域解,則繼續(xù)執(zhí)行下一步,否則,回到N1繼續(xù)執(zhí)行;否則,繼續(xù)執(zhí)行下一步,直到結(jié)束。

    變鄰域搜索策略2(記為S2),順序執(zhí)行各鄰域搜索,若出現(xiàn)比當(dāng)前解更優(yōu)的鄰域解,則繼續(xù)用該鄰域搜索模塊進(jìn)行搜索,直到比該鄰域解更差的解出現(xiàn),否則,繼續(xù)執(zhí)行下一步,直到結(jié)束。

    變鄰域搜索策略3(記為S3),順序執(zhí)行各鄰域搜索,若出現(xiàn)比當(dāng)前解更優(yōu)的鄰域解,則結(jié)束;否則,繼續(xù)執(zhí)行下一步,直到結(jié)束。

    3.4 選擇概率的確定

    在此采用人工蜂群算法常用選擇概率計(jì)算方法,即錦標(biāo)賽選擇策略。錦標(biāo)賽選擇策略是基于局部競爭機(jī)制的選擇過程,即在種群中隨機(jī)選擇q個(gè)個(gè)體(競賽規(guī)模)進(jìn)行比較,適應(yīng)度值大的被選中,并對其加權(quán)值1;如當(dāng)q=2時(shí),將種群中個(gè)體兩兩比較,并給適應(yīng)度值好的個(gè)體加權(quán)值1,對所有個(gè)體重復(fù)這一操作過程。根據(jù)個(gè)體權(quán)值,由此計(jì)算選擇概率:

    其中,ai為個(gè)體i的權(quán)值。

    該選擇策略中選擇概率的計(jì)算只涉及個(gè)體適應(yīng)度值之間的大小關(guān)系,不涉及適應(yīng)度值的具體取值,能較好避免陷入早熟。

    4 算例分析

    為分析變鄰域人工蜂群算法求解兩級定位-路徑問題的性能,需用算例進(jìn)行仿真測試。由于該問題目前缺乏公認(rèn)的國際算例,本文在文獻(xiàn)[12]一級定位-路徑問題算例25點(diǎn)、50點(diǎn)、75點(diǎn)基礎(chǔ)上,選擇前面部分設(shè)施點(diǎn)作為兩級定位-路徑問題的潛在配送中心、剩余設(shè)施點(diǎn)作為潛在中轉(zhuǎn)站、顧客點(diǎn)不變產(chǎn)生兩級定位-路徑問題的算例(算例中各類節(jié)點(diǎn)依次編號)。算例需補(bǔ)充和調(diào)整的參數(shù)如表1所示,其余參數(shù)詳見文獻(xiàn)[12]。

    表1 補(bǔ)充、調(diào)整的參數(shù)

    在此使用C sharp編程實(shí)現(xiàn)不同變鄰域搜索策略下的人工蜂群算法,在Windows XP系統(tǒng),Pentium?Dual-Core CPU E5700@3.00 GHz、2.00 GB內(nèi)存的計(jì)算環(huán)境下,對上述三組不同規(guī)模大小的算例進(jìn)行求解。種群規(guī)模SN=50,算法最大迭代次數(shù)Tmax=600,鄰域搜索最大迭代次數(shù)limit=5,初始解生成過程中所用參數(shù)r=4,車輛行駛成本與其行駛距離的比例系數(shù)取1。

    4.1 采用不同變鄰域搜索求解效果比較

    在此分別采用不同變鄰域搜索策略的人工蜂群算法求解兩級定位-路徑問題的上述算例。隨機(jī)運(yùn)行50次的求解結(jié)果見表2所示。

    由表2可知,采用鄰域搜索策略S1求解時(shí),取得的最好解和平均解較優(yōu);采用鄰域搜索策略S2求解時(shí),取得的最差解較優(yōu);采用鄰域搜索策略S3求解時(shí),平均計(jì)算時(shí)間較短。

    以下給出分別采用S1、S2、S3三種變鄰域搜索策略的人工蜂群算法求解算例“1-4-20”,隨機(jī)運(yùn)行50次種群最好解隨迭代次數(shù)變化趨勢圖,見圖1所示。

    從圖1可知,采用變鄰域搜索策略S1、S2、S3求解時(shí)均能較好避免早熟,隨著迭代次數(shù)的增加,目標(biāo)函數(shù)值不斷下降,算法具有較好的收斂性。

    4.2 不同算法的求解效果比較

    為了進(jìn)一步了解變鄰域人工蜂群算法求解兩級定位-路徑問題的求解效果,在此采用人工魚群算法進(jìn)行對比。借鑒文獻(xiàn)[13]的研究成果,將人工魚群算法的參數(shù)設(shè)置如下:最大試探次數(shù)try-number=5,擁擠度因子σ=0.4,可視范圍visual=(Ymax-Ymin)/4,其中Ymax、Ymin為當(dāng)代種群最好解、最差解對應(yīng)的目標(biāo)函數(shù)值。人工魚群算法的運(yùn)行環(huán)境與上述相同。隨機(jī)運(yùn)行50次,兩種算法求解結(jié)果見表3所示。

    表2 不同鄰域搜索策略求解效果比較

    圖1 種群最好解隨迭代次數(shù)的變化趨勢

    由表3可知,人工魚群算法在最差解及平均解方面較優(yōu),變鄰域人工蜂群算法在最好解及平均計(jì)算時(shí)間方面較優(yōu)。

    4.3 算例最好解網(wǎng)絡(luò)結(jié)構(gòu)

    在此,給出算例“3-7-65”最好解“3 220”對應(yīng)的網(wǎng)絡(luò)結(jié)構(gòu)圖,見圖2所示。

    從圖2可知,潛在的3個(gè)配送中心中開放了編號為“1”的1個(gè)配送中心,潛在的7個(gè)中轉(zhuǎn)站中開放了編號為“1”、“2”、“3”和“7”的4個(gè)中轉(zhuǎn)站,從配送中心出發(fā)到達(dá)各中轉(zhuǎn)站的路徑數(shù)為2,從各中轉(zhuǎn)站出發(fā)到達(dá)各顧客點(diǎn)的路徑數(shù)為12。根據(jù)該算例的相關(guān)數(shù)據(jù),可驗(yàn)證該最好解滿足假設(shè)條件(1)、(2)、(3)、(4)。由假設(shè)條件(5)可知,一級網(wǎng)絡(luò)車輛數(shù)即為2,二級網(wǎng)絡(luò)車輛數(shù)即為12。根據(jù)網(wǎng)絡(luò)中各節(jié)點(diǎn)的坐標(biāo),采用歐幾里德距離計(jì)算點(diǎn)與點(diǎn)之間的距離,求和取整得到各條路徑的長度,根據(jù)假設(shè)條件(6)可計(jì)算出各條路徑車輛行駛成本910(限于篇幅,詳細(xì)計(jì)算過程在此從略)。由此,可計(jì)算目標(biāo)函數(shù)=配送中心開放成本750×1+中轉(zhuǎn)站開放成本150×4+一級網(wǎng)絡(luò)車輛使用成本為240×2+二級網(wǎng)絡(luò)車輛使用成本40×12+兩級網(wǎng)絡(luò)中車輛行駛成本910=3 220。進(jìn)一步分析,由文獻(xiàn)[12]可計(jì)算出該算例的顧客需求量為1 168,結(jié)合表1所示參數(shù)可知,該最好解開放了最少數(shù)量的配送中心、最少的數(shù)量中轉(zhuǎn)站,使用了最小數(shù)量的車輛。因此,該解的質(zhì)量較好。

    表3 不同算法求解效果比較1)

    圖2 算例“3-7-65”最好解“3 220”的網(wǎng)絡(luò)結(jié)構(gòu)

    5 結(jié)束語

    針對城市物流配送中的兩級定位-路徑問題,建立了相應(yīng)的數(shù)學(xué)模型,并將求解NP-hard問題效果較好的人工蜂群算法應(yīng)用于包括多個(gè)NP-hard子問題的兩級定位-路徑問題的求解;針對人工蜂群算法容易出現(xiàn)早熟的問題,在算法的鄰域搜索階段,根據(jù)變鄰域搜索思想提出了三種變鄰域搜索策略。仿真實(shí)驗(yàn)結(jié)果表明,變鄰域人工蜂群算法能有效求解兩級定位-路徑問題;提出的三種變鄰域搜索策略各有優(yōu)勢,可根據(jù)實(shí)際應(yīng)用的需要及決策者的個(gè)人偏好,選擇相應(yīng)的搜索策略。因此,實(shí)際應(yīng)用時(shí),可根據(jù)決策者的風(fēng)險(xiǎn)偏好進(jìn)行選擇。若決策者為風(fēng)險(xiǎn)回避者,應(yīng)選擇基于變鄰域搜索策略S2的人工蜂群算法進(jìn)行求解,以期得到更優(yōu)的最差解;若決策者為風(fēng)險(xiǎn)中立者或風(fēng)險(xiǎn)追求者,應(yīng)選擇基于變鄰域搜索策略S1的人工蜂群算法進(jìn)行求解,以期得到更優(yōu)的平均解或最好解。此外,若在某些更關(guān)心決策效率的應(yīng)用情形,如應(yīng)急物流等,則應(yīng)選擇基于變鄰域搜索策略S3的人工蜂群算法進(jìn)行求解,以期在更短時(shí)間內(nèi)求得滿意解。

    [1]NagyG,Salhi S.Location-routing:issues,models and methods[J].European Journal of Operational Research,2007,177(2):649-672.

    [2]Jacobsen S K,Madsen O B G.A comparative study of heuristics for a two-level routing-location problem[J].European Journal of Operational Research,1980,5(6):378-387.

    [3]Ambrosino D,Sciomachena A,Scutellàb M G.A heuristic based on multi-exchange techniques for a regional fleet assignment location-routing problem[J].Computers&Operations Research,2009,36(2):442-460.

    [4]Sterle C.Location-routing models and methods for freight distribution and infomobility in city logistics[D].Universita Degli Studi di Napoli“Federico II”,Naples,Italy,2010.

    [5]Li Jin,Zhu Y,Shen H,et al.A hybrid genetic algorithm for two-layer location-routing problem[C]//4th International Conference on New Trends in Information Science and Service Science(NISS),2010:642-645.

    [6]Boccia M,Crainic T G,Sforza A,et al.A metaheuristic for a two Echelon Location-Routing Problem[C]//Lecture Notes in Computer Science.Berlin Heidelberg:Springer-Verlag,2010:288-301.

    [7]王紹仁,馬祖軍.震害緊急響應(yīng)階段應(yīng)急物流系統(tǒng)中的LRP[J].系統(tǒng)工程理論與實(shí)踐,2011,31(8):1497-1507.

    [8]Nguyen V P,Prins C,Prodhon C.Solving the two-echelon location routing problem by a GRASP reinforced by a learning process and path relinking[J].European Journal of Operational Research,2012,216(1):113-126.

    [9]Nguyen V P,Prins C,Prodhon C.A multi-start iterated local search with tabu list and path relinking for the two echelon location-routing problem[J].Engineering Applications of Artificial Intelligence,2012,25(1):56-71.

    [10]Karaboga D.An idea based on honey bee swarm for numerical optimization,TR06[R].Kaayseri:Erciyes University,2005.

    [11]暴勵.人工蜂群算法的混合策略研究[D].太原:太原科技大學(xué),2010.

    [12]胡大偉.設(shè)施定位和車輛路線問題模型及其啟發(fā)式算法研究[D].西安:長安大學(xué),2008:129-133.

    [13]李曉磊.一種新型的智能優(yōu)化方法——人工魚群算法[D].杭州:浙江大學(xué),2003.

    [14]董偉.變鄰域搜索算法研究及在組合優(yōu)化中的應(yīng)用[D].阜新:遼寧工程技術(shù)大學(xué),2011:1-2.

    CHEN Jiumei1,2,GONG Ying1,2

    1.Chongqing Key Laboratory of Electronic Commerce&Supply Chain System,Chongqing Technology and Business University,Chongqing 400067,China
    2.Strategical Planning Department,Chongqing Technology and Business University,Chongqing 400067,China

    A mathematical model of two-echelon location-routing problem is established.Artificial Bee Colony algorithm is put forward to solve this problem.Since this algorithm usually has premature convergence.Variable neighborhood search, a novel path heuristic algorithm appearing abroad in recent years,is blended in this algorithm.At the same time,three kinds of variable neighborhood search strategy are put forward.The comparison simulation between Artificial Bee Colony algorithm with different variable neighborhood search strategies and artificial fish swarm algorithm has been done.The experimental results show that Artificial Bee Colony algorithm with variable neighborhood search can effectively solve two-echelon location-routing problem.

    two-Echelon Location-Routing Problem(2E-LRP);Artificial Bee Colony(ABC)algorithm;Variable Neighborhood Search(VNS);logistics;distribution

    A

    F224.3

    10.3778/j.issn.1002-8331.1309-0335

    CHEN Jiumei,GONG Ying.Artificial Bee Colony algorithm with variable neighborhood search for two-Echelon Location-Routing Problem.Computer Engineering and Applications,2014,50(6):25-30.

    國家自然科學(xué)基金(No.71101159)。

    陳久梅(1976—),女,博士,副教授,研究領(lǐng)域?yàn)槲锪飨到y(tǒng)優(yōu)化、智能算法等;龔英(1968—),女,教授,研究領(lǐng)域?yàn)槲锪髋c供應(yīng)鏈。E-mail:chenjiumei@163.com

    2013-09-23

    2013-11-10

    1002-8331(2014)06-0025-06

    CNKI網(wǎng)絡(luò)優(yōu)先出版:2013-11-15,http://www.cnki.net/kcms/detail/11.2127.TP.20131115.1121.009.html

    猜你喜歡
    中轉(zhuǎn)站鄰域蜂群
    中亞是人類祖先關(guān)鍵“中轉(zhuǎn)站”?
    軍事文摘(2023年2期)2023-02-17 09:20:32
    高性能半柔性地坪在生活垃圾中轉(zhuǎn)站的應(yīng)用
    上海建材(2021年1期)2021-11-22 08:01:38
    “蜂群”席卷天下
    稀疏圖平方圖的染色數(shù)上界
    基于鄰域競賽的多目標(biāo)優(yōu)化算法
    關(guān)于-型鄰域空間
    改進(jìn)gbest引導(dǎo)的人工蜂群算法
    寧化石壁:客家人的第一中轉(zhuǎn)站
    福建人(2015年11期)2015-02-27 10:23:59
    蜂群夏季高產(chǎn)管理
    基于時(shí)序擴(kuò)展的鄰域保持嵌入算法及其在故障檢測中的應(yīng)用
    午夜免费鲁丝| 色网站视频免费| 国产毛片在线视频| 捣出白浆h1v1| 免费大片黄手机在线观看| 看十八女毛片水多多多| 中文字幕人妻丝袜制服| 欧美激情高清一区二区三区 | 一区二区日韩欧美中文字幕| 少妇精品久久久久久久| 国产精品久久久av美女十八| 在线观看免费视频网站a站| 午夜福利乱码中文字幕| 少妇的逼水好多| 国产午夜精品一二区理论片| 亚洲人成网站在线观看播放| 深夜精品福利| 国产成人一区二区在线| a级片在线免费高清观看视频| 考比视频在线观看| 午夜福利视频在线观看免费| 久久久久网色| 免费黄色在线免费观看| 狂野欧美激情性bbbbbb| 侵犯人妻中文字幕一二三四区| 香蕉精品网在线| 精品99又大又爽又粗少妇毛片| www.自偷自拍.com| 免费黄频网站在线观看国产| 80岁老熟妇乱子伦牲交| 亚洲第一青青草原| 午夜免费观看性视频| 免费在线观看完整版高清| 人妻一区二区av| 日本黄色日本黄色录像| 久久久久久久久久久免费av| 黄色一级大片看看| 人体艺术视频欧美日本| 丝袜美腿诱惑在线| 国产成人精品无人区| 满18在线观看网站| 中文天堂在线官网| 欧美亚洲 丝袜 人妻 在线| 欧美日韩视频精品一区| 69精品国产乱码久久久| 最黄视频免费看| 可以免费在线观看a视频的电影网站 | 狠狠精品人妻久久久久久综合| 国产在线一区二区三区精| 97在线人人人人妻| 久久精品熟女亚洲av麻豆精品| 亚洲视频免费观看视频| 亚洲,欧美精品.| 日韩,欧美,国产一区二区三区| 伦精品一区二区三区| 伦理电影大哥的女人| 黑丝袜美女国产一区| 美女脱内裤让男人舔精品视频| 国产av精品麻豆| 欧美日韩综合久久久久久| 你懂的网址亚洲精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 成人午夜精彩视频在线观看| 美女xxoo啪啪120秒动态图| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级专区第一集| 九色亚洲精品在线播放| xxx大片免费视频| 亚洲精品中文字幕在线视频| 一区二区三区激情视频| 婷婷色av中文字幕| 午夜精品国产一区二区电影| 精品亚洲乱码少妇综合久久| 青春草视频在线免费观看| 国产av一区二区精品久久| 精品人妻熟女毛片av久久网站| 亚洲国产日韩一区二区| 99国产综合亚洲精品| 视频在线观看一区二区三区| 欧美日韩精品成人综合77777| 精品一区二区免费观看| 欧美精品国产亚洲| 国产片内射在线| 日本av手机在线免费观看| 国产免费一区二区三区四区乱码| 国语对白做爰xxxⅹ性视频网站| 最近2019中文字幕mv第一页| 美女视频免费永久观看网站| 精品亚洲成a人片在线观看| 美女午夜性视频免费| 亚洲国产精品999| 黑丝袜美女国产一区| 久久久久久久久久久久大奶| 午夜日本视频在线| 亚洲美女黄色视频免费看| 青草久久国产| 亚洲国产av影院在线观看| 有码 亚洲区| 日本爱情动作片www.在线观看| 热re99久久国产66热| 色婷婷久久久亚洲欧美| 精品国产超薄肉色丝袜足j| 免费黄频网站在线观看国产| 亚洲男人天堂网一区| 久久久久国产网址| 国产成人精品久久二区二区91 | 人妻一区二区av| 免费观看av网站的网址| 亚洲av日韩在线播放| 99热国产这里只有精品6| 国产一区亚洲一区在线观看| 亚洲国产精品一区三区| 国产在视频线精品| 在线观看www视频免费| av在线老鸭窝| 日韩,欧美,国产一区二区三区| 中文天堂在线官网| 精品国产一区二区久久| 国产日韩一区二区三区精品不卡| 丝袜美足系列| 一边摸一边做爽爽视频免费| 国产乱人偷精品视频| 伦精品一区二区三区| 久久精品国产亚洲av高清一级| 国产熟女午夜一区二区三区| 日韩中文字幕视频在线看片| 黄色怎么调成土黄色| 欧美国产精品va在线观看不卡| 欧美日韩亚洲高清精品| 精品人妻在线不人妻| 成人亚洲欧美一区二区av| 日产精品乱码卡一卡2卡三| 日韩中文字幕视频在线看片| 久久青草综合色| 又黄又粗又硬又大视频| 亚洲精品中文字幕在线视频| 国产成人精品一,二区| 国产探花极品一区二区| 天堂中文最新版在线下载| 中文欧美无线码| 少妇的丰满在线观看| 国产一区二区 视频在线| 99久久精品国产国产毛片| 天堂俺去俺来也www色官网| av在线app专区| 国产精品不卡视频一区二区| 一级毛片黄色毛片免费观看视频| 欧美亚洲 丝袜 人妻 在线| 99久久中文字幕三级久久日本| 中国国产av一级| 国产1区2区3区精品| 亚洲精品成人av观看孕妇| 女人高潮潮喷娇喘18禁视频| 久热这里只有精品99| 黑丝袜美女国产一区| 18在线观看网站| 一级毛片 在线播放| 久久青草综合色| 欧美日韩av久久| 18禁动态无遮挡网站| 久久狼人影院| 久久久国产精品麻豆| 电影成人av| 亚洲国产av影院在线观看| 久久热在线av| 免费在线观看视频国产中文字幕亚洲 | 免费不卡的大黄色大毛片视频在线观看| 国产精品一国产av| 99久久中文字幕三级久久日本| 午夜免费鲁丝| 亚洲精品第二区| tube8黄色片| 色94色欧美一区二区| 久久精品亚洲av国产电影网| 亚洲综合色网址| 亚洲婷婷狠狠爱综合网| 99久久中文字幕三级久久日本| 男女高潮啪啪啪动态图| 午夜日韩欧美国产| 999精品在线视频| 久久国产精品男人的天堂亚洲| 国产精品国产av在线观看| 成人亚洲欧美一区二区av| 大香蕉久久成人网| 中文精品一卡2卡3卡4更新| av卡一久久| 香蕉精品网在线| 热99久久久久精品小说推荐| 欧美xxⅹ黑人| 国产麻豆69| 人人澡人人妻人| 9热在线视频观看99| 日日啪夜夜爽| 考比视频在线观看| 你懂的网址亚洲精品在线观看| av免费观看日本| 国产高清国产精品国产三级| a 毛片基地| 亚洲第一av免费看| 人人妻人人澡人人爽人人夜夜| 91国产中文字幕| 亚洲av.av天堂| 久久久久久久精品精品| 亚洲欧美精品综合一区二区三区 | 一级毛片黄色毛片免费观看视频| 亚洲国产av新网站| 欧美少妇被猛烈插入视频| 在线观看国产h片| 黄色视频在线播放观看不卡| 丁香六月天网| 成年人免费黄色播放视频| 只有这里有精品99| 国产视频首页在线观看| a级毛片黄视频| 男女午夜视频在线观看| 成人亚洲精品一区在线观看| 大香蕉久久成人网| 国产人伦9x9x在线观看 | 大话2 男鬼变身卡| 美女国产视频在线观看| 卡戴珊不雅视频在线播放| 伊人久久大香线蕉亚洲五| 亚洲,欧美,日韩| 免费高清在线观看视频在线观看| 亚洲精品国产一区二区精华液| 又黄又粗又硬又大视频| videos熟女内射| 国产男女超爽视频在线观看| 久久精品久久久久久噜噜老黄| 中国国产av一级| 一级a爱视频在线免费观看| 国产精品麻豆人妻色哟哟久久| 五月开心婷婷网| 亚洲熟女精品中文字幕| 综合色丁香网| 国产精品免费视频内射| 婷婷色av中文字幕| 亚洲伊人色综图| 亚洲在久久综合| 日产精品乱码卡一卡2卡三| 欧美日韩视频精品一区| 精品午夜福利在线看| www.精华液| 午夜日韩欧美国产| 中文字幕亚洲精品专区| 亚洲成av片中文字幕在线观看 | 一级毛片电影观看| 欧美 亚洲 国产 日韩一| 国产在线一区二区三区精| 伦理电影大哥的女人| 午夜av观看不卡| 精品少妇久久久久久888优播| 麻豆精品久久久久久蜜桃| 另类亚洲欧美激情| 美女大奶头黄色视频| av网站在线播放免费| 午夜老司机福利剧场| 男女无遮挡免费网站观看| 亚洲精品自拍成人| 午夜免费鲁丝| 亚洲av在线观看美女高潮| 亚洲国产精品一区二区三区在线| 国产男人的电影天堂91| 热re99久久国产66热| 天天躁夜夜躁狠狠躁躁| 美女高潮到喷水免费观看| 免费观看a级毛片全部| 黄色怎么调成土黄色| 飞空精品影院首页| 欧美+日韩+精品| 五月伊人婷婷丁香| 啦啦啦中文免费视频观看日本| 亚洲国产欧美日韩在线播放| www.av在线官网国产| 日韩伦理黄色片| 啦啦啦视频在线资源免费观看| 久久久久久久大尺度免费视频| 看免费成人av毛片| 免费在线观看黄色视频的| 99九九在线精品视频| 97在线人人人人妻| 国产熟女欧美一区二区| 久久久久久伊人网av| 久久午夜福利片| 丝袜喷水一区| 国产女主播在线喷水免费视频网站| 欧美国产精品一级二级三级| 黄色视频在线播放观看不卡| 两个人免费观看高清视频| 亚洲av中文av极速乱| 精品一区在线观看国产| 热re99久久精品国产66热6| 捣出白浆h1v1| 又大又黄又爽视频免费| 黄色毛片三级朝国网站| 乱人伦中国视频| 国产成人精品久久久久久| 黑人猛操日本美女一级片| 国产无遮挡羞羞视频在线观看| av在线app专区| 满18在线观看网站| 伦理电影免费视频| 丰满饥渴人妻一区二区三| 老熟女久久久| 久久精品人人爽人人爽视色| 成人国产av品久久久| 叶爱在线成人免费视频播放| 一本色道久久久久久精品综合| a级毛片黄视频| 久久久精品免费免费高清| 日韩av在线免费看完整版不卡| 色视频在线一区二区三区| 国产在线一区二区三区精| 最近最新中文字幕大全免费视频 | 九草在线视频观看| 亚洲综合色惰| 中文乱码字字幕精品一区二区三区| 亚洲人成网站在线观看播放| 91国产中文字幕| 色播在线永久视频| 天天躁狠狠躁夜夜躁狠狠躁| 男女下面插进去视频免费观看| www.熟女人妻精品国产| 在线观看www视频免费| 欧美激情极品国产一区二区三区| 一区二区av电影网| 久久久久网色| 亚洲精品久久成人aⅴ小说| 一级片免费观看大全| 国产精品一区二区在线不卡| 国产精品久久久久久久久免| 精品酒店卫生间| 一本色道久久久久久精品综合| 狂野欧美激情性bbbbbb| 欧美黄色片欧美黄色片| 免费观看a级毛片全部| av在线播放精品| 极品少妇高潮喷水抽搐| 寂寞人妻少妇视频99o| 91久久精品国产一区二区三区| 另类精品久久| 欧美最新免费一区二区三区| 成年女人毛片免费观看观看9 | 免费观看无遮挡的男女| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 中文字幕制服av| 大话2 男鬼变身卡| 成人漫画全彩无遮挡| 久久久久久久大尺度免费视频| 肉色欧美久久久久久久蜜桃| 国产男女内射视频| 国产精品国产三级国产专区5o| 精品少妇黑人巨大在线播放| 成人国产麻豆网| 亚洲四区av| 国产成人午夜福利电影在线观看| 亚洲色图 男人天堂 中文字幕| 卡戴珊不雅视频在线播放| 亚洲三级黄色毛片| 只有这里有精品99| 美女福利国产在线| 国产精品麻豆人妻色哟哟久久| 天天躁夜夜躁狠狠久久av| 免费观看性生交大片5| 十八禁高潮呻吟视频| 亚洲一区二区三区欧美精品| 在线观看免费高清a一片| 成年动漫av网址| 成人黄色视频免费在线看| 亚洲欧美一区二区三区久久| 考比视频在线观看| 免费播放大片免费观看视频在线观看| 99九九在线精品视频| 国产午夜精品一二区理论片| 免费观看a级毛片全部| 久久精品国产综合久久久| 18+在线观看网站| 免费在线观看黄色视频的| a级毛片黄视频| 老汉色av国产亚洲站长工具| 一级毛片 在线播放| 国产精品国产三级国产专区5o| videosex国产| www.精华液| 狠狠精品人妻久久久久久综合| 成年美女黄网站色视频大全免费| 美女xxoo啪啪120秒动态图| 丁香六月天网| 免费高清在线观看视频在线观看| 亚洲婷婷狠狠爱综合网| 久热这里只有精品99| av福利片在线| 免费日韩欧美在线观看| 深夜精品福利| 日韩大片免费观看网站| 人妻系列 视频| 国产男女超爽视频在线观看| www.精华液| 久久精品国产鲁丝片午夜精品| 18在线观看网站| 美国免费a级毛片| 18+在线观看网站| av片东京热男人的天堂| 精品少妇一区二区三区视频日本电影 | 青青草视频在线视频观看| 欧美精品高潮呻吟av久久| 在线观看免费日韩欧美大片| 久久国内精品自在自线图片| 久久av网站| 欧美精品国产亚洲| 久久人人97超碰香蕉20202| 久久免费观看电影| 欧美日韩国产mv在线观看视频| 另类亚洲欧美激情| 少妇人妻久久综合中文| 麻豆乱淫一区二区| 亚洲成人av在线免费| 国产精品一区二区在线不卡| 午夜福利视频在线观看免费| 两性夫妻黄色片| 国产老妇伦熟女老妇高清| 久久久久久久久久久免费av| 国产一区二区 视频在线| 黄色视频在线播放观看不卡| 人人妻人人澡人人爽人人夜夜| 七月丁香在线播放| 少妇熟女欧美另类| 岛国毛片在线播放| 色哟哟·www| 18禁动态无遮挡网站| 久久久久久久久久久久大奶| 亚洲精品久久午夜乱码| 在线观看一区二区三区激情| 99香蕉大伊视频| 丝袜脚勾引网站| 国产探花极品一区二区| 亚洲内射少妇av| 999久久久国产精品视频| 2022亚洲国产成人精品| 午夜福利影视在线免费观看| 不卡视频在线观看欧美| √禁漫天堂资源中文www| 丁香六月天网| 少妇被粗大猛烈的视频| 麻豆精品久久久久久蜜桃| 欧美精品人与动牲交sv欧美| 国产精品久久久久久av不卡| av在线老鸭窝| 免费观看a级毛片全部| 免费黄色在线免费观看| 精品国产露脸久久av麻豆| 久久97久久精品| 亚洲精品一二三| 啦啦啦在线观看免费高清www| 国产精品蜜桃在线观看| 制服丝袜香蕉在线| 99久久综合免费| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| 日本91视频免费播放| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 欧美精品av麻豆av| 美女午夜性视频免费| 精品国产一区二区久久| 纯流量卡能插随身wifi吗| 亚洲精品,欧美精品| 欧美精品亚洲一区二区| 中文字幕人妻丝袜制服| 国产1区2区3区精品| 99re6热这里在线精品视频| 欧美另类一区| 久久狼人影院| 久久久久视频综合| 伊人久久国产一区二区| 五月伊人婷婷丁香| 天堂俺去俺来也www色官网| 欧美日韩一区二区视频在线观看视频在线| 香蕉丝袜av| 9191精品国产免费久久| 久久鲁丝午夜福利片| 在线观看国产h片| 国产色婷婷99| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| 青草久久国产| 麻豆精品久久久久久蜜桃| 日韩人妻精品一区2区三区| 男人爽女人下面视频在线观看| 十八禁网站网址无遮挡| 日本wwww免费看| 香蕉精品网在线| 国产精品成人在线| 国产精品国产三级国产专区5o| 日韩av不卡免费在线播放| 久久青草综合色| 免费黄频网站在线观看国产| 人人妻人人澡人人看| 777久久人妻少妇嫩草av网站| 99re6热这里在线精品视频| 黑人猛操日本美女一级片| 精品一品国产午夜福利视频| 又粗又硬又长又爽又黄的视频| 国产亚洲av片在线观看秒播厂| 国产成人91sexporn| 亚洲美女搞黄在线观看| 日韩中字成人| 少妇人妻精品综合一区二区| 亚洲精品aⅴ在线观看| www日本在线高清视频| 国产av码专区亚洲av| 久久国产亚洲av麻豆专区| 国产成人一区二区在线| 欧美精品av麻豆av| 少妇的逼水好多| 国产亚洲午夜精品一区二区久久| 欧美在线黄色| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 日韩伦理黄色片| 999久久久国产精品视频| 人人妻人人澡人人看| 久久久久久久国产电影| 久久精品国产自在天天线| 国产精品久久久久久久久免| 国产精品一国产av| 丝袜在线中文字幕| 国产 一区精品| 亚洲精品中文字幕在线视频| 日本av免费视频播放| 性色av一级| 国产精品一国产av| 韩国av在线不卡| 免费观看av网站的网址| 中文字幕av电影在线播放| 国产探花极品一区二区| 欧美精品av麻豆av| 日韩在线高清观看一区二区三区| 欧美精品国产亚洲| 免费少妇av软件| 成人国语在线视频| 丝瓜视频免费看黄片| 免费播放大片免费观看视频在线观看| 国产爽快片一区二区三区| 视频区图区小说| av在线app专区| 成人国产av品久久久| 久久亚洲国产成人精品v| 国产视频首页在线观看| 久久精品熟女亚洲av麻豆精品| 有码 亚洲区| h视频一区二区三区| 国产精品 欧美亚洲| h视频一区二区三区| 多毛熟女@视频| 赤兔流量卡办理| 久久热在线av| 十八禁网站网址无遮挡| 日本av手机在线免费观看| 日韩,欧美,国产一区二区三区| 国产1区2区3区精品| 26uuu在线亚洲综合色| 丰满饥渴人妻一区二区三| 91aial.com中文字幕在线观看| 人妻系列 视频| 汤姆久久久久久久影院中文字幕| 免费在线观看视频国产中文字幕亚洲 | 新久久久久国产一级毛片| 少妇人妻久久综合中文| 免费黄频网站在线观看国产| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 欧美日韩精品成人综合77777| 熟妇人妻不卡中文字幕| 亚洲国产av新网站| 在线观看一区二区三区激情| 免费久久久久久久精品成人欧美视频| 国产午夜精品一二区理论片| 精品第一国产精品| 日本黄色日本黄色录像| 毛片一级片免费看久久久久| 国产高清不卡午夜福利| 亚洲综合精品二区| 国产女主播在线喷水免费视频网站| 97人妻天天添夜夜摸| 国产精品国产三级专区第一集| 999精品在线视频| 国产一区二区在线观看av| videos熟女内射| 国产精品三级大全| 国产成人午夜福利电影在线观看| 国产综合精华液| 日韩熟女老妇一区二区性免费视频| 天天躁日日躁夜夜躁夜夜| 天堂中文最新版在线下载| 欧美日韩av久久| 久久99精品国语久久久| 黄色 视频免费看| 国产成人91sexporn| 欧美在线黄色| 国产人伦9x9x在线观看 | 欧美日韩视频高清一区二区三区二| 国产在线免费精品| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠久久av| 久久精品久久久久久久性| 视频区图区小说| 男人舔女人的私密视频| 性少妇av在线| 巨乳人妻的诱惑在线观看| 欧美国产精品一级二级三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲一区中文字幕在线|