• <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)用
    夜夜躁狠狠躁天天躁| 国产午夜精品久久久久久| 欧美黑人欧美精品刺激| 久久精品成人免费网站| 制服人妻中文乱码| 亚洲av五月六月丁香网| 美女午夜性视频免费| 俄罗斯特黄特色一大片| 国产麻豆成人av免费视频| 国产人伦9x9x在线观看| 欧美乱码精品一区二区三区| 国语自产精品视频在线第100页| √禁漫天堂资源中文www| 日韩欧美一区二区三区在线观看| 淫秽高清视频在线观看| 久久久久久大精品| 欧美日韩中文字幕国产精品一区二区三区 | 很黄的视频免费| 桃色一区二区三区在线观看| 午夜老司机福利片| 一区二区三区精品91| 一区二区三区精品91| 非洲黑人性xxxx精品又粗又长| 免费不卡黄色视频| 亚洲av五月六月丁香网| 老司机靠b影院| 亚洲成人国产一区在线观看| 在线观看午夜福利视频| 最新美女视频免费是黄的| 看免费av毛片| 在线播放国产精品三级| av视频免费观看在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲欧美日韩无卡精品| 18禁观看日本| 一边摸一边抽搐一进一出视频| 久热爱精品视频在线9| 一区二区三区激情视频| 久久婷婷成人综合色麻豆| 久久香蕉精品热| 精品一区二区三区视频在线观看免费| 亚洲精品中文字幕在线视频| 韩国精品一区二区三区| 亚洲av成人不卡在线观看播放网| 啦啦啦免费观看视频1| 午夜免费观看网址| 欧美乱妇无乱码| 午夜久久久久精精品| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 一夜夜www| 在线av久久热| 国产精品乱码一区二三区的特点 | 一级作爱视频免费观看| 亚洲av片天天在线观看| 精品电影一区二区在线| 黑丝袜美女国产一区| 人人妻人人澡人人看| 国产成人免费无遮挡视频| 悠悠久久av| 国产精品一区二区三区四区久久 | 非洲黑人性xxxx精品又粗又长| 黄色毛片三级朝国网站| 色在线成人网| 满18在线观看网站| 级片在线观看| 成人特级黄色片久久久久久久| 精品国产美女av久久久久小说| 亚洲中文av在线| 一个人观看的视频www高清免费观看 | 亚洲成国产人片在线观看| 国产精品永久免费网站| 男女下面插进去视频免费观看| tocl精华| 亚洲性夜色夜夜综合| 亚洲国产精品成人综合色| 嫁个100分男人电影在线观看| 韩国av一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜免费激情av| 18禁美女被吸乳视频| 最近最新免费中文字幕在线| 真人一进一出gif抽搐免费| 精品一区二区三区四区五区乱码| 法律面前人人平等表现在哪些方面| 亚洲国产精品久久男人天堂| 成人三级做爰电影| 50天的宝宝边吃奶边哭怎么回事| 色综合婷婷激情| 久久精品亚洲熟妇少妇任你| 精品少妇一区二区三区视频日本电影| 亚洲五月色婷婷综合| 国产精品影院久久| 大陆偷拍与自拍| 夜夜夜夜夜久久久久| 咕卡用的链子| 久久婷婷人人爽人人干人人爱 | 在线国产一区二区在线| 99re在线观看精品视频| 亚洲人成电影观看| 一级毛片高清免费大全| 97人妻精品一区二区三区麻豆 | 成年人黄色毛片网站| 欧美中文综合在线视频| 国产一区二区三区视频了| 一二三四社区在线视频社区8| 99国产综合亚洲精品| 日韩av在线大香蕉| 久久精品91蜜桃| 亚洲一区二区三区色噜噜| 青草久久国产| 99在线视频只有这里精品首页| 女人精品久久久久毛片| 亚洲第一电影网av| 久久 成人 亚洲| 国产精品亚洲美女久久久| 久久香蕉精品热| 久久人人精品亚洲av| 午夜福利,免费看| 国产成人精品久久二区二区91| 欧美成人性av电影在线观看| 欧美成人一区二区免费高清观看 | 日韩成人在线观看一区二区三区| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 国产欧美日韩一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| 中文字幕色久视频| 女同久久另类99精品国产91| 国产亚洲欧美在线一区二区| 欧美日韩亚洲国产一区二区在线观看| 亚洲五月色婷婷综合| 亚洲av成人av| 中出人妻视频一区二区| 久久久久久免费高清国产稀缺| 变态另类丝袜制服| 欧美丝袜亚洲另类 | 国产精品久久久久久人妻精品电影| 伦理电影免费视频| 91字幕亚洲| 国产精品av久久久久免费| 中出人妻视频一区二区| 亚洲三区欧美一区| av有码第一页| 国产精品 欧美亚洲| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 亚洲av成人不卡在线观看播放网| 成在线人永久免费视频| 国产亚洲欧美精品永久| 亚洲成人久久性| 天天添夜夜摸| 正在播放国产对白刺激| 黑人巨大精品欧美一区二区蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产精品久久男人天堂| 韩国精品一区二区三区| 一边摸一边抽搐一进一小说| 9191精品国产免费久久| 国产精品久久久久久精品电影 | 亚洲 欧美一区二区三区| 波多野结衣高清无吗| 日日干狠狠操夜夜爽| 91大片在线观看| 国产精品综合久久久久久久免费 | 淫妇啪啪啪对白视频| 欧美一区二区精品小视频在线| 久久精品人人爽人人爽视色| 午夜久久久在线观看| 在线观看免费视频日本深夜| or卡值多少钱| 欧美国产日韩亚洲一区| 婷婷精品国产亚洲av在线| 日韩av在线大香蕉| 热re99久久国产66热| 亚洲国产精品999在线| 中亚洲国语对白在线视频| 午夜成年电影在线免费观看| 无限看片的www在线观看| 久久久久精品国产欧美久久久| 亚洲色图 男人天堂 中文字幕| 女人精品久久久久毛片| 国产精品电影一区二区三区| 又紧又爽又黄一区二区| 久久久久久免费高清国产稀缺| 少妇被粗大的猛进出69影院| 国内毛片毛片毛片毛片毛片| 国产成人av教育| 可以在线观看毛片的网站| 999精品在线视频| 99在线人妻在线中文字幕| 国产精品美女特级片免费视频播放器 | 日韩中文字幕欧美一区二区| 久久婷婷人人爽人人干人人爱 | 国产精品免费视频内射| 国产在线观看jvid| 亚洲无线在线观看| 久久久久国内视频| 法律面前人人平等表现在哪些方面| 国产精品乱码一区二三区的特点 | 电影成人av| x7x7x7水蜜桃| 国产高清videossex| 亚洲午夜精品一区,二区,三区| 午夜精品国产一区二区电影| 1024香蕉在线观看| 成人亚洲精品一区在线观看| 日韩大码丰满熟妇| 久久久国产成人免费| 男女之事视频高清在线观看| 极品人妻少妇av视频| 午夜成年电影在线免费观看| 女同久久另类99精品国产91| 90打野战视频偷拍视频| 久久草成人影院| 亚洲一区中文字幕在线| 婷婷六月久久综合丁香| 亚洲精品久久国产高清桃花| 看免费av毛片| 精品久久久久久,| 岛国视频午夜一区免费看| 国产精品二区激情视频| 50天的宝宝边吃奶边哭怎么回事| 在线播放国产精品三级| 在线av久久热| 久久国产亚洲av麻豆专区| 欧美 亚洲 国产 日韩一| 999精品在线视频| 波多野结衣av一区二区av| 国产成年人精品一区二区| 久久中文字幕人妻熟女| 亚洲人成电影观看| 狂野欧美激情性xxxx| 久久伊人香网站| 两个人看的免费小视频| 国产精品综合久久久久久久免费 | 99国产精品一区二区蜜桃av| 国产亚洲精品综合一区在线观看 | 久久人人精品亚洲av| 免费少妇av软件| 色综合站精品国产| 午夜福利视频1000在线观看 | 国产99白浆流出| 精品久久久久久,| 美女高潮到喷水免费观看| 非洲黑人性xxxx精品又粗又长| 亚洲在线自拍视频| 国产精品亚洲一级av第二区| 亚洲男人天堂网一区| 日韩av在线大香蕉| 大码成人一级视频| 精品卡一卡二卡四卡免费| 最好的美女福利视频网| 国产欧美日韩精品亚洲av| 欧美成人一区二区免费高清观看 | 欧美成狂野欧美在线观看| 午夜成年电影在线免费观看| 麻豆一二三区av精品| 国产亚洲精品av在线| 午夜两性在线视频| 丁香欧美五月| 久久久久久久午夜电影| 大香蕉久久成人网| 男女下面进入的视频免费午夜 | 国产成人一区二区三区免费视频网站| 亚洲第一电影网av| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| 国产97色在线日韩免费| 亚洲精品一卡2卡三卡4卡5卡| 亚洲性夜色夜夜综合| 午夜成年电影在线免费观看| 国产一区二区三区视频了| av在线天堂中文字幕| 午夜免费鲁丝| 国产一区在线观看成人免费| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 777久久人妻少妇嫩草av网站| 亚洲欧美日韩高清在线视频| 色综合婷婷激情| 在线观看www视频免费| 久热这里只有精品99| 欧美日韩亚洲综合一区二区三区_| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久免费视频| 亚洲国产欧美网| 欧美激情久久久久久爽电影 | 亚洲一区中文字幕在线| 亚洲国产毛片av蜜桃av| 久久中文字幕人妻熟女| 97超级碰碰碰精品色视频在线观看| 成人精品一区二区免费| 亚洲一区二区三区色噜噜| 一区二区三区高清视频在线| 99国产极品粉嫩在线观看| 99精品在免费线老司机午夜| 色播在线永久视频| 老熟妇乱子伦视频在线观看| 高清在线国产一区| 黑人操中国人逼视频| 看黄色毛片网站| 国产精品野战在线观看| 精品不卡国产一区二区三区| 久久人人精品亚洲av| 国产av一区二区精品久久| 美女扒开内裤让男人捅视频| 香蕉国产在线看| 99精品在免费线老司机午夜| 中文亚洲av片在线观看爽| 18禁观看日本| 露出奶头的视频| 18禁黄网站禁片午夜丰满| 亚洲天堂国产精品一区在线| 国产熟女xx| 999精品在线视频| 午夜福利18| 男人舔女人下体高潮全视频| 成在线人永久免费视频| 别揉我奶头~嗯~啊~动态视频| 久久影院123| 国产精品国产高清国产av| 亚洲精品在线观看二区| 亚洲成人久久性| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 国产国语露脸激情在线看| 精品电影一区二区在线| 久久久国产欧美日韩av| 首页视频小说图片口味搜索| 给我免费播放毛片高清在线观看| 精品熟女少妇八av免费久了| 又大又爽又粗| 成人三级做爰电影| 精品久久久精品久久久| 亚洲全国av大片| 一a级毛片在线观看| 大陆偷拍与自拍| 日韩有码中文字幕| 国产精品一区二区在线不卡| 757午夜福利合集在线观看| 婷婷精品国产亚洲av在线| 亚洲中文av在线| 色综合欧美亚洲国产小说| 久久人人精品亚洲av| 91国产中文字幕| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 99riav亚洲国产免费| 一级毛片高清免费大全| 亚洲精品在线观看二区| 免费看十八禁软件| 国产熟女xx| 免费高清视频大片| 亚洲精品国产色婷婷电影| 亚洲av成人一区二区三| 国产私拍福利视频在线观看| 欧美中文日本在线观看视频| 麻豆久久精品国产亚洲av| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕最新亚洲高清| 一级黄色大片毛片| 大陆偷拍与自拍| 最新美女视频免费是黄的| 国产高清激情床上av| 国产极品粉嫩免费观看在线| av免费在线观看网站| 国产精品久久久久久精品电影 | 在线免费观看的www视频| 国产麻豆成人av免费视频| 老司机在亚洲福利影院| 极品教师在线免费播放| 国产精品电影一区二区三区| 宅男免费午夜| 精品日产1卡2卡| 国产一级毛片七仙女欲春2 | 国产99白浆流出| 麻豆成人av在线观看| 国产精品久久视频播放| 久久久久国产精品人妻aⅴ院| 国产精品久久久人人做人人爽| 亚洲一区高清亚洲精品| 色哟哟哟哟哟哟| 亚洲av电影在线进入| 最近最新中文字幕大全电影3 | www.精华液| 18禁观看日本| 免费看十八禁软件| 国产精品亚洲av一区麻豆| 一二三四在线观看免费中文在| 免费观看人在逋| 国产激情欧美一区二区| 亚洲美女黄片视频| 如日韩欧美国产精品一区二区三区| 不卡av一区二区三区| 在线免费观看的www视频| 精品熟女少妇八av免费久了| 久久人妻av系列| 色婷婷久久久亚洲欧美| 欧美在线黄色| 1024视频免费在线观看| 精品日产1卡2卡| 天天一区二区日本电影三级 | 欧美大码av| 亚洲国产日韩欧美精品在线观看 | 中文亚洲av片在线观看爽| 国产蜜桃级精品一区二区三区| 午夜福利在线观看吧| 国内久久婷婷六月综合欲色啪| 亚洲国产高清在线一区二区三 | 日日摸夜夜添夜夜添小说| 波多野结衣高清无吗| 给我免费播放毛片高清在线观看| 欧美日本视频| 欧美激情 高清一区二区三区| 99久久综合精品五月天人人| 欧美日本亚洲视频在线播放| 国产成人欧美在线观看| 不卡一级毛片| 亚洲五月天丁香| 亚洲自偷自拍图片 自拍| 色综合站精品国产| 黄片播放在线免费| 变态另类丝袜制服| 人人澡人人妻人| 一区二区日韩欧美中文字幕| 在线观看日韩欧美| 午夜免费成人在线视频| 男人操女人黄网站| 黄网站色视频无遮挡免费观看| 国产欧美日韩一区二区精品| 亚洲中文字幕一区二区三区有码在线看 | 欧美精品啪啪一区二区三区| 国产成人精品在线电影| 国产欧美日韩综合在线一区二区| 精品高清国产在线一区| 久久狼人影院| 男人操女人黄网站| 成人永久免费在线观看视频| 欧美老熟妇乱子伦牲交| 好男人在线观看高清免费视频 | 超碰成人久久| 免费少妇av软件| 桃红色精品国产亚洲av| 精品国内亚洲2022精品成人| 天天躁夜夜躁狠狠躁躁| 成人亚洲精品一区在线观看| 午夜免费激情av| 精品国产乱子伦一区二区三区| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 狠狠狠狠99中文字幕| 正在播放国产对白刺激| 国产精品综合久久久久久久免费 | 国产一级毛片七仙女欲春2 | 国产麻豆成人av免费视频| 中文字幕精品免费在线观看视频| 日韩欧美一区视频在线观看| 亚洲五月婷婷丁香| 色综合欧美亚洲国产小说| 最新在线观看一区二区三区| 精品福利观看| 美女高潮喷水抽搐中文字幕| 一进一出抽搐动态| 不卡av一区二区三区| 国产一卡二卡三卡精品| 久久久国产欧美日韩av| 欧美日韩中文字幕国产精品一区二区三区 | 欧美精品啪啪一区二区三区| 国产免费av片在线观看野外av| 18美女黄网站色大片免费观看| 99国产精品99久久久久| 久久久国产欧美日韩av| 高清毛片免费观看视频网站| 亚洲av成人一区二区三| 男女做爰动态图高潮gif福利片 | 成人亚洲精品一区在线观看| 看片在线看免费视频| 欧美日本中文国产一区发布| 色综合欧美亚洲国产小说| 搡老岳熟女国产| 精品欧美一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 国产成人欧美在线观看| 波多野结衣av一区二区av| 中文字幕色久视频| 久久欧美精品欧美久久欧美| 黄片大片在线免费观看| 久久久久国内视频| 黄色毛片三级朝国网站| 国产精品亚洲一级av第二区| 日韩大码丰满熟妇| 亚洲在线自拍视频| 91麻豆精品激情在线观看国产| 久久伊人香网站| 久久久国产欧美日韩av| 久久久久国产一级毛片高清牌| 国产精品影院久久| 色尼玛亚洲综合影院| 日韩欧美一区二区三区在线观看| 伦理电影免费视频| 97人妻精品一区二区三区麻豆 | av电影中文网址| 人妻丰满熟妇av一区二区三区| 纯流量卡能插随身wifi吗| 91av网站免费观看| 久久久久久久久免费视频了| 最好的美女福利视频网| 在线观看免费日韩欧美大片| 法律面前人人平等表现在哪些方面| 精品电影一区二区在线| 1024视频免费在线观看| 咕卡用的链子| 成人三级黄色视频| 亚洲精品久久国产高清桃花| 韩国av一区二区三区四区| 国产一卡二卡三卡精品| x7x7x7水蜜桃| 国产午夜福利久久久久久| 一夜夜www| 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 最近最新免费中文字幕在线| 午夜久久久久精精品| 无遮挡黄片免费观看| 亚洲人成77777在线视频| 亚洲av片天天在线观看| 91精品三级在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产亚洲av香蕉五月| 成人18禁在线播放| 露出奶头的视频| 91九色精品人成在线观看| 在线观看午夜福利视频| 一边摸一边做爽爽视频免费| 一级,二级,三级黄色视频| 亚洲全国av大片| 国产成人av教育| 国产精品免费视频内射| 亚洲五月婷婷丁香| 97人妻天天添夜夜摸| 99国产综合亚洲精品| 99久久综合精品五月天人人| 韩国av一区二区三区四区| 12—13女人毛片做爰片一| 69av精品久久久久久| 欧美乱码精品一区二区三区| 变态另类丝袜制服| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 19禁男女啪啪无遮挡网站| 亚洲精品国产色婷婷电影| 精品福利观看| 不卡一级毛片| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 色综合亚洲欧美另类图片| 男人操女人黄网站| 精品福利观看| 久久久水蜜桃国产精品网| 国产精品一区二区免费欧美| 欧美亚洲日本最大视频资源| 久久久水蜜桃国产精品网| 久久伊人香网站| av天堂在线播放| 亚洲av成人av| 性少妇av在线| 波多野结衣高清无吗| 国产成人欧美在线观看| 国产精品香港三级国产av潘金莲| 搡老岳熟女国产| 午夜福利成人在线免费观看| 1024香蕉在线观看| 一本大道久久a久久精品| 欧美激情 高清一区二区三区| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| av在线天堂中文字幕| 免费搜索国产男女视频| 女人爽到高潮嗷嗷叫在线视频| 日本精品一区二区三区蜜桃| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 亚洲精品久久国产高清桃花| 欧美日韩精品网址| 欧美黄色片欧美黄色片| 可以在线观看的亚洲视频| 成人特级黄色片久久久久久久| av有码第一页| 人妻丰满熟妇av一区二区三区| 男人舔女人下体高潮全视频| 桃色一区二区三区在线观看| 婷婷精品国产亚洲av在线| 深夜精品福利| 欧美激情高清一区二区三区| 国产成人免费无遮挡视频| 啦啦啦 在线观看视频| 免费久久久久久久精品成人欧美视频| 国产亚洲精品久久久久久毛片| 中文字幕人妻丝袜一区二区| 18美女黄网站色大片免费观看| 黄色a级毛片大全视频| 国产蜜桃级精品一区二区三区| 人人妻人人澡人人看| 黄色 视频免费看| 黄片大片在线免费观看| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 啦啦啦观看免费观看视频高清 | 欧美日韩亚洲国产一区二区在线观看| 日本欧美视频一区| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 欧美大码av|