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

    開放式帶時間窗車輛路徑問題及變鄰域搜索算法

    2021-11-10 04:33:14陳久梅李英娟
    計算機集成制造系統(tǒng) 2021年10期
    關鍵詞:搜索算法算例鄰域

    陳久梅,李英娟,胡 婷,但 斌,李 俊

    (1.重慶工商大學 長江上游經(jīng)濟研究中心,重慶 400067;2.重慶工商大學 管理科學與工程學院,重慶 400067;3.重慶工商大學 工商管理學院,重慶 400067;4.重慶大學 經(jīng)濟與工商管理學院,重慶 400044)

    0 引言

    近年來,隨著O2O平臺的快速發(fā)展,物流配送需求大幅度增長,人們的消費模式從線下向線上轉變。相比線下而言,線上平臺面臨更多挑戰(zhàn),如配送時效性與準確性、配送成本等。以美團外賣為例,在消費者下單、商家接單和騎手配送的幾大環(huán)節(jié)中,騎手配送是關鍵,面臨最大的挑戰(zhàn)便是車輛調度問題。與傳統(tǒng)車輛路徑問題(Vehicle Routing Problem, VRP)不同,騎手配送環(huán)節(jié)往往由第三方平臺負責,配送結束后車輛不必返回商家,這屬于開放式車輛路徑問題。此外,顧客往往對接受配送有最早和最晚時間的要求。該車輛路徑具有受時間窗約束的特征,這類問題屬于學術界研究的開放式帶時間窗車輛路徑問題(Open Vehicle Routing Problem with Time Windows, OVRPTW)。

    OVRPTW在現(xiàn)實中應用廣泛,如報紙配送、牛奶配送、快遞服務、校園班車、大型生產(chǎn)制造業(yè)原材料的運輸?shù)取R孕@班車為例,校車下午從學校出發(fā),沿著預定路線,依次按照預定時間到達學生下車點,直到最后一位學生下車,早上接學生的情況正好相反。又如,當某些公司不擁有車隊,或其車隊不合適或不足以滿足顧客需求時,這家公司需要將全部或部分貨物交付第三方或借助外部載體進行運輸,這種情況下車輛不必返回倉庫,費用主要取決于行駛距離,即開放路線的長度。類似地,對于危險品、易耗品、易腐食品、建筑材料、原材料的運輸?shù)龋约拌F路部門或者航空公司在運送貨物時,車輛在完成一批貨物運輸后,為節(jié)省物流成本無需返回原出發(fā)點。上述物流配送過程中,貨物的到達時間是影響消費者服務滿意度的關鍵因素之一,因此車輛路徑通常是執(zhí)行具有時間窗約束的開放路線。

    雖然OVRPTW的應用廣泛,但目前的研究還相對較少。REPOUSSIS等[1]針對OVRPTW設計了一種貪婪的前瞻性路徑構造啟發(fā)式求解算法,該算法通過組合顧客選擇和路徑插入來利用時間窗相關信息;FAIZ等[2]針對OVRPTW分別建立了兩個整數(shù)線性規(guī)劃模型:基于弧的混合整數(shù)線性規(guī)劃模型并采用通用求解器求解和基于路徑公式設計一個列生成框架,最后經(jīng)過數(shù)值實驗比較了兩種模型的性能;RUSSELL[3]提出一種利用建模語言和CP優(yōu)化器約束編程的方法,解決多產(chǎn)品報紙生產(chǎn)和投遞的協(xié)調問題;CHEN等[4]基于緊急程度采用插入啟發(fā)式算法構造初始解,運用帶強化學習的變鄰域搜索算法求解具有時間窗和開放路徑的周期車輛路徑問題;REPOUSSIS等[5]結合新達爾文主義進化模型的基本原理和短時記憶的禁忌搜索算法,提出一種進化算法求解OVRPTW;BRANDO[6]針對OVRPTW設計了迭代局部搜索算法,并采用Solomon、Homberger和Gehring等算例進行求解;李三彬等[7]提出一種多開始禁忌搜索算法求解OVRPTW,同樣使用Solomon算例進行驗證;RUSSELL等[8]在構建初始路徑的基礎上,使用禁忌搜索對路徑改進;孫國華[9]設計了改進的自適應遺傳算法求解開放車輛路徑。

    變鄰域搜索(Variable Neighbourhood Search, VNS)算法作為一種新元啟發(fā)式算法,最早由MLADENOVI等[10]提出,該算法通過改變鄰域結構使搜索方向多元化,從而增強搜索能力,優(yōu)化計算性能[11]。鑒于VNS算法的優(yōu)越性,其已經(jīng)被廣泛應用于解決VRP的各種擴展問題。XU等[11]使用VNS算法求解連續(xù)性車輛路徑問題;SEVKLI等[12]提出一種求解OVRP的變鄰域搜索算法,并在各種標準算例和大規(guī)模實際問題上進行了測試;REZGUI等[13]將VNS算法應用于解決具有異車型和時間窗等約束條件的問題;姚冠新等[14]將改進VNS算法用于研究多隔室車輛路徑優(yōu)化問題。

    本文針對OVRPTW,在建立數(shù)學模型的基礎上,借鑒文獻[1]的貪婪插入法生成初始解,設計VNS算法對其進行求解。算法框架中融入GLOVER[15]的路徑重連和KINDERVATER等[16]的鄰域搜索算子,并將鄰域搜索嵌套在VNS算法的抖動過程中。通過求解兩組算例,并將求解結果與已有文獻進行比較,驗證了VNS算法的有效性。

    1 問題描述與模型構建

    OVRPTW是經(jīng)典VRPTW的一個擴展問題,可考慮為配送中心安排車輛將物品配送給一系列顧客,每輛車完成自己所負責配送的最后一位顧客后,服務結束,無需返回配送中心。其中,車輛有容量限制,在配送過程中不得超載;顧客在地理位置上是分散的,且各自有不同的需求量及接受配送服務的時間窗。OVRPTW的目標是構建一組滿足所有顧客點需求的總行駛成本最小的配送路徑。

    OVRPTW可用有向圖D=(V,A)來表示,其中V表示點集,A表示邊集。點集V=Vc∪{0},其中Vc表示顧客點集,0表示配送中心。邊集A=Ac∪E(0),其中Ac={(i,j)|i,j∈Vc,i≠j}表示顧客點i到顧客點j的邊集,E(0)表示離開配送中心的邊。給定點子集G?Vc,δ(G)表示到達G的邊,E(G)表示離開G的邊。d(G)表示G中所有顧客的需求量之和,Q表示車輛容量。車輛集合為K={1,2,…,k}。

    借鑒KALLEHAUGE[17]的方法,OVRPTW的解記為R={r1,r2,…,rk},其中rk=(0,vk1,vk2,…),vkj∈Vc。路徑rk上所有顧客的需求量之和記為d(rk),車輛行駛成本記為c(rk),路徑rk上顧客vkj的實際開始服務時間記為ssvkj,服務時間記為svkj,車輛從點i到點j的行駛時間記為tij,節(jié)點i的開始服務時間時間窗為[ai,bi],其中a0=0。若路徑rk上任意顧客vkj都滿足avkj≤ssvkj≤bvkj且d(rk)≤Q,則路徑rk為可行路徑。若任意rk∈R均為可行路徑,則R為OVRPTW的一個可行解,X(?)表示“?”中邊的數(shù)量。由此,可構建OVRPTW的模型。

    目標函數(shù):

    (1)

    s.t.

    X(δ(i))=1,?i∈Vc;

    (2)

    X(E(i))=1,?i∈V,i≠vk(|rk|-1),?k∈K;

    (3)

    X(E(i))=0,?i=vk(|rk|-1),?k∈K;

    (4)

    X(E(G))≤|G|-1,G?Vc,|G|≥2;

    (5)

    ssvkj=max{ssvk(j-1)+svk(j-1)+tvk(j-1)vkj,avkj},

    ?j∈Vc,?k∈K;

    (6)

    avkj≤ssvkj≤bvkj,?j∈Vc,?k∈K;

    (7)

    d(rk)≤Q,?k∈K。

    (8)

    目標函數(shù)(1)表示車輛總行駛成本最小。約束(2)表示有且僅有一輛車到達顧客點為其提供服務;約束(3)表示車輛到達任意路徑rk上除最后一個顧客點外,必須離開;約束(2)和約束(3)保證了路徑的連續(xù)性;約束(4)表示任意路徑rk上車輛服務完最后一個顧客即結束;約束(5)用于消除子回路徑。約束(6)和約束(7)表示任意路徑上任意顧客vkj的實際開始服務時間必須在其時間窗內(nèi);約束(8)表示車輛容量約束,即任意路徑rk上顧客需求量之和不得超過車輛容量。

    OVRPTW中,車輛不需要像VRPTW那樣返回配送中心。從圖論角度看,求解VRPTW是要找到一組哈密頓回路(如圖1),而求解OVRPTW是要找到一組哈密頓路徑(如圖2)。SYSLO等[18]認為哈密頓路徑問題是NP-hard問題,因為它可以轉換成等價的哈密頓回路問題。因此,OVRPTW也是NP-hard問題??紤]到變鄰域搜索算法的優(yōu)越性以及該算法在車輛路徑問題方面的成功應用[18-19],在此運用變鄰域搜索算法對OVRPTW進行求解。

    2 算法設計

    VNS算法的基本思想是通過系統(tǒng)改變鄰域,從一個解跳到另一個解,不斷改善現(xiàn)有解,以得到更好的解。與禁忌搜索、模擬退火等傳統(tǒng)的元啟發(fā)式算法相比,VNS算法不但結構簡單,容易實現(xiàn),而且與求解問題無關,適用于各種優(yōu)化問題[20]。為提高求解效果,本文在基本VNS算法的框架下,在其抖動階段,采用當前解向個體歷史最優(yōu)解和種群歷史最優(yōu)解靠近的路徑重連來實現(xiàn);在其鄰域搜索階段,采用路徑內(nèi)和路徑間的交換、插入、2-opt算子來實現(xiàn),并將鄰域搜索嵌套在抖動過程中。

    定義1抖動算子記為SNk,k=1,2,其中SN1指當前解向個體最優(yōu)解靠近的路徑重連,SN2指當前解向種群最優(yōu)解靠近的路徑重連;鄰域搜索記為LNl,l=1,2,3,其中LN1指交換算子,LN2指插入算子,LN3指2-opt算子。求解OVRPTW的VNS算法偽代碼如下:

    Generate the initial solution x

    Definition:a set of neighborhood structures SNkfor k=1,2 for shaking

    a set of neighborhood structures LNlfor l=1,2,3 for local search

    while(stopping criterion is not met)do

    k=1

    l=1

    while k≤2 do //shaking

    x′←SNk(x)

    while l≤3 do // neighborhood search

    x″←LNl(x′)

    If Z(x″)

    x′←x″

    else

    l=l+1

    End if

    End while

    If Z(x′)

    x←x′

    Else

    k=k+1

    End if

    End while

    End while

    Output x

    2.1 解的表示

    OVRPTW的路徑集合R={r1,r2,…,rk},其中路徑rk表示車輛k服務的顧客點及其服務順序,rk=(0,vk1,vk2,…)。為保證每輛車從配送中心出發(fā),服務完最后一個顧客后不再返回,每條路徑第一個元素均為0,其余元素為非0。每個非0數(shù)字對應一個顧客,數(shù)字不重復。所有路徑的集合構成一個解,如路徑r1=(0,1,3),r2=(0,2,5,6),r3=(0,4)表示由3輛車服務6個顧客點的一個解,如圖3所示。

    圖3中,車輛1從配送中心出發(fā),依次服務顧客1和3;車輛2從配送中心出發(fā),依次服務顧客2、5和6;車輛3從配送中心出發(fā),服務顧客4。這個解可表示為{0,1,3,0,2,5,6,0,4}。該種表示具有如下優(yōu)點:①每條路徑可以表示車輛提供服務的顧客點及其服務的先后順序;②每個顧客點的服務時間可以根據(jù)式(6)計算得到,然后根據(jù)式(7)判斷是否滿足時間窗約束;③可以很方便地計算整條路徑上所有顧客的需求總量,從而根據(jù)約束條件(8)判斷是否超過車輛容量。

    2.2 初始解的生成

    初始解的生成采用貪婪插入法,將顧客點依次插入最優(yōu)插入位置,形成一條條行駛路徑,直到完成所有顧客點的插入,即生成初始解。

    定義2最優(yōu)插入位置。給定路徑rk=(0,vk1,vk2,…),在vk1之前和vkj之后共有|rk|個插入候選位置。如果插入后新路徑增加的車輛行駛成本最小,則該位置稱為最優(yōu)插入位置。

    生成初始解的步驟如下:

    步驟1隨機選擇一個不在任意一條路徑上的顧客點,將其與配送中心相連接。

    步驟2以該顧客點和配送中心的連線為軸,圍繞配送中心順時針進行掃描,以時間窗和車輛容量為約束條件,將滿足條件的顧客點按照貪婪插入的方法,插入到最優(yōu)插入位置,循環(huán)該步驟,直到不滿足約束為止,便形成一條車輛行駛路徑。

    步驟3若系統(tǒng)中存在不在任意一條路徑上的顧客點,則轉步驟1。

    步驟4若系統(tǒng)中所有顧客均在某條路徑上,則形成一系列行駛路徑,得到初始解。

    步驟5算法結束。

    2.3 抖動

    抖動是變鄰域搜索的關鍵階段,它能夠改變解的搜索方向,實現(xiàn)搜索空間多樣化,避免陷入局部最優(yōu)。本文算法中的抖動階段采用路徑重連實現(xiàn)當前解到另一個解的轉換。

    定義3可行的路徑重連。路徑重連(path relinking)是GLOVER[15]在1997年引入的元啟發(fā)式算法,該方法通過在當前解與導向解之間建立聯(lián)系來生成新解,若生成的新解可行,則該操作稱為可行的路徑重連。本文分別以個體歷史最優(yōu)解和種群歷史最優(yōu)解作為導向解進行路徑重連。

    抖動的步驟如下:

    步驟1選擇導向解中某條路徑rk。

    步驟2將路徑rk包含的顧客點從當前解中刪除。

    步驟3當前解中,將有顧客點被刪除的路徑中剩下的顧客點,在滿足式(7)和式(8)有關時間窗和車輛容量約束條件下,添加到?jīng)]有刪除顧客點的其他路徑上。

    步驟4將路徑rk直接加入到當前解中,便得到了一個新解。

    步驟5算法結束。

    假定隨機選擇導向解中的一條路徑為(0,2,4)。當前解中含有顧客點2和顧客點4的路徑分別為r1和r2。首先分別從r1和r2中刪除顧客點2和顧客點4;r1和r2中剩余顧客點,在滿足時間窗和車輛容量約束條件下,依次添加到?jīng)]有刪除顧客點的其他路徑上;不滿足上述約束的顧客點,則保留在原路徑上。將導向解中路徑(0,2,4)加入到當前解中,由此生成新解。

    2.4 鄰域搜索

    鄰域搜索采用交換、插入和2-opt三種算子[16]。三種算子同時考慮選擇的顧客點或邊在同一路徑和不同路徑2種情況。

    定義4可行交換?;诳尚薪庵型粭l路徑,交換2個顧客點產(chǎn)生新路徑。若交換后,新路徑上交換位之后的顧客點滿足式(7)顧客時間窗約束,則該交換為可行交換?;诳尚薪庵胁煌窂剑粨Q2個顧客點產(chǎn)生兩條新路徑。若形成的新路徑均滿足式(7)顧客時間窗約束和式(8)車輛容量約束,則該交換為可行交換。

    交換步驟如下:

    步驟1在路徑集中隨機選擇2個顧客點。

    步驟4算法結束。

    定義5可行插入。在路徑中隨機選擇一個顧客,若將該顧客點從當前位置刪除,然后插入到同一路徑或不同路徑中,形成新的路徑依然是可行路徑,則這個插入操作被稱為可行插入。

    插入步驟如下:

    步驟1從路徑rk=(0,vk1,…,vkl,…,vkn,…)中隨機選擇一個顧客(記為vkl)并將其從rk中刪除。

    步驟4算法結束。

    定義6可行2-opt。在同一路徑rk或不同路徑rk和rk′中,將兩條不相鄰的邊斷開重新連接,若形成路徑是可行路徑,則稱該2-opt操作為可行2-opt。

    2-opt步驟如下:

    步驟1在路徑集中隨機選擇兩條不相鄰的邊,將其斷開。

    步驟3若斷開的兩條邊位于不同路徑上,分別記為rk和rk′。

    步驟5算法結束。

    由圖4可知,當斷開同一條路徑上不相連的兩條邊(1,3)和(4,2)后,將被斷開的路徑段(3,4)反轉,得到新的路徑(0,1,4,3,2,5)。

    由圖5可知,當斷開兩條路徑上兩條邊(3,5)和(4,7)后,將被斷開的兩個路徑段進行交叉連接,得到新的路徑(0,1,3,7,8)和(0,2,4,5,6)。

    需要注意的是,鄰域搜索時,由于顧客點位置發(fā)生變化,需要判斷路徑是否可行。由于在同一路徑上執(zhí)行3種操作算子時,該路徑包含的顧客點沒有發(fā)生變化,必定滿足車輛容量的約束,只需考慮是否滿足時間窗。而在不同路徑上執(zhí)行3種操作算子時,每條路徑上包含的顧客點發(fā)生變化,此時除了考慮是否滿足時間窗的約束以外,還要考慮是否滿足車輛容量的約束。

    2.5 算法復雜度分析

    VNS算法時間復雜度分析如下:對于最大迭代次數(shù)為Tmax、問題規(guī)模為I的OVRPTW,生成初始解解的時間復雜度為O(1);每個解的操作中,顧客移動、顧客交換、2-opt和路徑重連的時間復雜度均為O(I2);是否接受新解的時間復雜度為O(1);每一代的時間復雜度為a×O(I2)+b×O(1),其中a,b為正整數(shù);Tmax代的時間復雜度為Tmax×(a×O(I2)+b×O(1))。因此,整個算法的時間復雜度為O(1)+Tmax×(a×O(I2)+b×O(1))。可見,所提算法的時間復雜度與算法迭代次數(shù)和問題規(guī)模有關。

    2.6 算法證明

    引理1最優(yōu)解存在性。OVRPTW存在最優(yōu)解x*。

    證明OVRPTW為線性規(guī)劃問題,屬于凸優(yōu)化問題[21],因此存在最優(yōu)解x*[22],證畢。

    定理1鄰域存在性。模型存在最優(yōu)解,因此在可行解集F內(nèi),給定一個解x(x≠x*),則至少存在一個鄰域結構N,滿足N=|x→x*|,并稱解x*在解x的鄰域結構N內(nèi),即解x與最優(yōu)解x*之間至少存在一個廣義距離,廣義距離可以是歐式距離、漢明距離或者k-opt算子等[23]。

    證明反證法。若模型對于任意給定可行解x(x∈F),不存在一個鄰域結構N=|x→x*|,即解x*不在可行解集F中,則解x*不是模型的最優(yōu)解,與已知矛盾。證畢。

    定理2局部最優(yōu)解可窮性。模型為組合優(yōu)化問題且存在最優(yōu)解x*,任意給定一個可行解x,則在解x的鄰域結構N(x)范圍內(nèi)可以找到所有局部最優(yōu)解。

    證明因為模型為組合優(yōu)化問題,則N(x)鄰域內(nèi)解空間有限,所以可以通過枚舉算法獲得所有局部最優(yōu)解,即局部最優(yōu)解是可窮的。證畢。

    引理2路徑可行性。使用VNS算法得到的開放式帶時間窗的車輛路徑是可行的。

    證明可行解集F中所有解x均滿足定義的等式約束(2)~約束(8),變鄰域也是在可行解集F內(nèi)尋找不同的鄰域結構Ng,從而尋找局部最優(yōu)解xg,即xg∈F,因此局部最優(yōu)解xg滿足定義的等式約束(2)~約束(8),局部最優(yōu)解中的路徑也滿足模型條件約束。證畢。

    基于以上定理和引理,證明使用VNS算法求解OVRPTW具有可行性。

    3 實驗結果與分析

    下面采用兩組算例測試本文VNS算法求解OVRPTW的性能。算例1來自文獻[24];算例2采用100個顧客點的Solomon標準測試算例(http://web.cba.neu.edu/~msolomon/problems.htm)。

    算法采用Microsoft Visual Studio 2017編程,在Intel(R)Core(TM)i7-8750H CPU@2.20 GHz 2.21 GHz環(huán)境下運行。本文參數(shù)設置如下:種群規(guī)模為N=50,最大迭代次數(shù)Tmax=1 000,隨機運行10次。

    3.1 算例1結果對比分析

    算例1中顧客坐標是在[1,150]×[1,150]的區(qū)域內(nèi)隨機生成,其編號為1,2,…,19,20,配送中心坐標為(70,70),其中顧客具體信息(坐標、貨物需求量及時間窗)參見文獻[24]。本文將VNS算法與文獻[24]中的簡單遺傳算法、基于or-opt的模擬退火算法、混合遺傳算法,文獻[25]中的基本蟻群算法、單點單親遺傳混合蟻群算法、多點單親遺傳混合蟻群算法,文獻[26]中的基本蝙蝠、精英遺傳混合蝙蝠、多點重組精英混合蝙蝠、單點重組精英混合蝙蝠共11種算法的求解結果進行比較,具體結果如表1所示。

    表1 11種算法結果比較

    由表1可知,使用同一組算例求解VRPTW問題時,本文VNS算法的平均路徑長度以及最優(yōu)路徑長度均最短,而使用的車輛數(shù)比簡單遺傳算法、基本蟻群算法、單點單親遺傳混合蟻群算法、多點單親遺傳混合蟻群算法、基本蝙蝠以及精英遺傳混合蝙蝠都要少。這表明,在求解VRPTW時,本文VNS算法求解質量相比其他算法更優(yōu)。

    盡管表1所呈現(xiàn)的結果能初步驗證本文的VNS算法在求解質量上具有優(yōu)越性,但還需進一步通過統(tǒng)計學上的顯著性檢驗,比較排序以上11種算法在同一個數(shù)據(jù)集上的性能優(yōu)劣。Friedman檢驗[27]作為非參數(shù)檢驗的方法,能對多種算法在每個數(shù)據(jù)集上的計算結果進行比較排序,并檢驗算法在置信區(qū)間內(nèi)是否存在顯著性差異。詳細比較結果和統(tǒng)計結果分別如表2和表3所示。

    表2 11種算法Friedman檢驗排序及秩平均值

    表3 Friedman統(tǒng)計結果

    Friedman檢驗秩平均值越低,算法性能越優(yōu)。因此,由表2的結果可知,11種算法性能排序為:VNS算法(1.67)>單點(多點)重組精英混合蝙蝠(3.00)>混合遺傳算法(3.67)>精英遺傳混合蝙蝠(5.50)>單點單親遺傳混合蟻群算法(6.17)>基于or-opt的模擬退火算法(6.33)>多點單親遺傳混合蟻群算法(7.17)>基本蟻群算法(9.50)>簡單遺傳算法(9.83)>基本蝙蝠(10.17)。

    表3顯示了在自由度為10、置信區(qū)間為99%的卡方分布下,F(xiàn)riedman統(tǒng)計量為25.262,p值為0.005<0.01,由此可得出11種算法之間存在1%水平下顯著差異的結論。綜上所述,根據(jù)表1的求解結果,如表2和表3的統(tǒng)計結果可知,本文VNS算法在各方面均較優(yōu),11種算法之間存在顯著性差異。因此,初步驗證了本文VNS算法在求解VRPTW時具有可行性和有效性。

    3.2 算例2結果對比分析

    算例2采用經(jīng)典的Solomon Benchmark算例集,共包括R1,C1,RC1,R2,C2,RC2六類。每類算例均是在100×100的單位區(qū)域內(nèi)生成了1個配送中心和100個顧客點,其中R類、C類和RC類客戶的地理位置分別由隨機分布、聚類分布以及混合隨機聚類分布生成;R1、C1、RC1類算例的時間窗緊,而R2、C2、RC2類算例的時間窗寬,因此在每條路徑上能服務更多的客戶。Solomon Benchmark的算例,原本是VRPTW的算例。在此,借鑒文獻[5-6]的作法,即采用該算例的基本數(shù)據(jù),車輛服務完最后一個顧客即結束,不用返回配送中心。本文將VNS算法與文獻[5]的改進的禁忌搜索啟發(fā)式算法(Improved Tabu Search Heuristic,TS)、貪婪的前瞻性路線構造啟發(fā)式算法(Greedy Look-Ahead Route Construction heuristic Algorithm,GLARCA)、進化算法(Evolutionary Algorithm, EA)和貪婪的隨機多啟動軌跡局部搜索算法(Greedy Randomized Multi-Start Trajectory local Search algorithm, GRMSTS)以及文獻[6]的統(tǒng)一的變鄰域搜索算法(Unified Variable Neighborhood Search,UVNS)、統(tǒng)一的混合遺傳搜索算法(Unified Hybrid Genetic Search,UHGS)、彈射鏈迭代局部搜索算法(Iterated Local Search Algorithm,ILSA)和求解OVRPTW的當前已知最好解(Best Known Solutions,BKS)的求解結果進行比較。

    (1)解比較

    通過與已有文獻最好解,以及其他算法得到的解進行比較,分析本文VNS算法的求解性能,如表4和表5所示。表中:TD表示總行駛成本;GAP表示VNS算法的解比當前已知最好解之間的差距,GAP=(已有文獻最好解-已知最好解)/已知最好解×100%。

    表4 Solomon中R1、C1和RC1數(shù)據(jù)集的詳細結果

    表5 Solomon中R2、C2和RC2數(shù)據(jù)集的詳細結果

    由表4可知,本文VNS算法在求解R1、C1和RC1類算例中,有16個算例的VNS解比當前已知最好解更優(yōu),5個算例的VNS解與當前已知最好解相同,其行駛成本平均節(jié)約了0.60%。由表5可知,在求解R2、C2、RC2類算例中,有19個算例的VNS解比當前已知最好解更優(yōu),3個算例的VNS解與當前已知最好解相同,其行程成本平均節(jié)約了7.80%。因此,由表4和表5的詳細結果及比較可知,VNS算法能表現(xiàn)出良好的求解質量和性能,由此進一步驗證了本文VNS在求解OVRPTW上的可行性和有效性。

    (2)收斂性分析

    為分析VNS算法的收斂性,本文以算例2中R101~R104、RC101~RC104、R201~R204、和RC201~RC204為例,根據(jù)VNS算法求解OVRPTW隨機運行10次,迭代1 000次的過程數(shù)據(jù),繪制兩組VNS算法收斂圖,如圖6和圖7所示。

    由圖6可知,算法在300代之前收斂速度快,300~900代中收斂速度較慢,在900代左右向各自的最好解收斂。由圖7可知,算法在200代之前收斂速度較快,在200~900代中收斂速度緩慢,在900代后收斂于各自的最好解。因此,根據(jù)兩組收斂圖可知,VNS算法不管是在求解時間窗寬的算例還是時間窗窄的算例,最終都能收斂于各自的最好解。由此可知,本文VNS算法具有較好的收斂性。

    (3)穩(wěn)定性分析

    為了分析本文VNS算法的穩(wěn)定性,采用箱形圖來觀測VNS算法求解R101~R104、RC101~RC104、R201~R204、和RC201~RC204每組算例10次后得到解的分布情況,如圖8和圖9所示。其中,黑色粗線表示中位數(shù),該箱形圖分別由上四分位數(shù)和下四分位數(shù)組成,豎線表示該組結果的平均偏差,豎線兩端的橫線分別為上限(最大值)和下限(最小值),空心的圓圈表示異常值。

    由圖8可知,R101~R104和RC101~RC104的箱盒長度短,上四分位數(shù)和下四分位數(shù)間距較小,且部分算例其上下四分位數(shù)接近于0。相比圖8來說,圖9中R201~R204和RC201~RC204的箱盒長度雖較長,但異常值較少,其中R201~R204的解不存在異常值。從兩組箱形圖整體來看,平均偏差均較小。由此可知,本文利用VNS算法求解OVRPTW 10次后解的離散程度較小,從而表明本文VNS算法具有穩(wěn)定性。

    4 結束語

    在外賣配送過程中,通常需要騎手滿足配送時效性、準確性以及開放性等要求。本文基于此背景,對開放式帶時間窗車輛路徑問題進行探討,考慮配送車輛無需返回配送中心和時間窗等影響,構建配送行駛總成本最小化為目標函數(shù)的OVRPTW模型,并設計了一種變鄰域搜索的啟發(fā)式算法進行問題求解。本文通過使用兩組算例,利用VNS算法分別對VRPTW和OVRPTW進行求解,并與已有文獻和現(xiàn)有已知最好解進行對比分析。實驗結果表明:①本文設計的VNS算法在求解VRPTW時,能獲得更短的運輸里程,且與眾多種算法相比,該算法優(yōu)勢明顯;②在求解OVRPTW時,大多數(shù)算例的VNS算法最好解都優(yōu)于已有文獻最好解,通過該算法求解,明顯降低了車輛行駛成本,由此表現(xiàn)出VNS算法良好的求解質量;③通過收斂圖和箱形圖可知,本文VNS算法在求解OVRPTW時,具有較好的收斂性和穩(wěn)定性,最終驗證了該算法良好的求解性能。

    未來,將進一步考慮天氣、交通、實時路況和突發(fā)事件等對路徑優(yōu)化的影響,建立相應的數(shù)學模型并設計算法求解。

    猜你喜歡
    搜索算法算例鄰域
    改進的和聲搜索算法求解凸二次規(guī)劃及線性規(guī)劃
    稀疏圖平方圖的染色數(shù)上界
    基于鄰域競賽的多目標優(yōu)化算法
    自動化學報(2018年7期)2018-08-20 02:59:04
    關于-型鄰域空間
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補問題算例分析
    基于汽車接力的潮流轉移快速搜索算法
    基于逐維改進的自適應步長布谷鳥搜索算法
    基于CYMDIST的配電網(wǎng)運行優(yōu)化技術及算例分析
    基于跳點搜索算法的網(wǎng)格地圖尋路
    亚洲高清免费不卡视频| 亚洲最大成人中文| av天堂在线播放| 97超视频在线观看视频| 亚洲精品久久国产高清桃花| 一卡2卡三卡四卡精品乱码亚洲| 91精品一卡2卡3卡4卡| 色哟哟哟哟哟哟| av免费在线看不卡| 日韩高清综合在线| 国产精品麻豆人妻色哟哟久久 | 精品人妻一区二区三区麻豆| 久久精品国产亚洲av涩爱 | 国产一区二区在线av高清观看| 99在线人妻在线中文字幕| 九色成人免费人妻av| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 三级男女做爰猛烈吃奶摸视频| 久久人人精品亚洲av| 亚洲18禁久久av| 亚洲va在线va天堂va国产| 亚洲国产精品久久男人天堂| 女同久久另类99精品国产91| 久久99热这里只有精品18| av在线播放精品| 婷婷色av中文字幕| 国内精品一区二区在线观看| 黄片wwwwww| 综合色丁香网| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人中文| 久久精品国产亚洲av香蕉五月| 亚洲精品国产成人久久av| 亚洲精品成人久久久久久| 内射极品少妇av片p| www.av在线官网国产| 欧美日韩国产亚洲二区| 久久鲁丝午夜福利片| 不卡视频在线观看欧美| 亚洲av二区三区四区| 永久网站在线| 国产av在哪里看| 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说 | 99久久精品一区二区三区| 日韩欧美在线乱码| 少妇熟女aⅴ在线视频| 在线观看免费视频日本深夜| av天堂在线播放| 亚洲精品色激情综合| 最近手机中文字幕大全| 亚洲人成网站高清观看| 天天躁夜夜躁狠狠久久av| 又粗又硬又长又爽又黄的视频 | 网址你懂的国产日韩在线| 久久精品夜夜夜夜夜久久蜜豆| 久久6这里有精品| 欧美成人免费av一区二区三区| 欧美+日韩+精品| 国产真实伦视频高清在线观看| 国产黄色小视频在线观看| 免费看日本二区| 波野结衣二区三区在线| 国产精品无大码| 18禁在线播放成人免费| 又爽又黄无遮挡网站| 麻豆成人av视频| 国产一区二区三区av在线 | 国产老妇女一区| 色综合亚洲欧美另类图片| 亚洲婷婷狠狠爱综合网| 国产大屁股一区二区在线视频| 在线观看66精品国产| 亚洲综合色惰| 国产精品久久久久久精品电影小说 | 国产在线精品亚洲第一网站| 久久精品人妻少妇| 啦啦啦观看免费观看视频高清| 色哟哟哟哟哟哟| 国产成人福利小说| 午夜激情欧美在线| 国产毛片a区久久久久| 久久精品久久久久久噜噜老黄 | 色尼玛亚洲综合影院| .国产精品久久| 免费av毛片视频| 国模一区二区三区四区视频| 久久这里有精品视频免费| 亚洲精品自拍成人| 免费搜索国产男女视频| 高清在线视频一区二区三区 | 亚洲av中文av极速乱| 我要看日韩黄色一级片| 蜜桃亚洲精品一区二区三区| 亚洲av二区三区四区| 国产私拍福利视频在线观看| 久久久色成人| 欧美日本视频| 丰满的人妻完整版| 在线免费十八禁| 午夜免费男女啪啪视频观看| 少妇被粗大猛烈的视频| 一级毛片我不卡| 久久精品国产亚洲av天美| 尾随美女入室| 人妻夜夜爽99麻豆av| 色哟哟·www| 欧美色欧美亚洲另类二区| 亚洲,欧美,日韩| 看非洲黑人一级黄片| 嫩草影院入口| 国产午夜精品一二区理论片| 男女下面进入的视频免费午夜| 成人午夜高清在线视频| 欧美xxxx性猛交bbbb| 内地一区二区视频在线| 亚洲av中文字字幕乱码综合| 九九在线视频观看精品| 搡女人真爽免费视频火全软件| 国产黄色小视频在线观看| 国产亚洲欧美98| 中文资源天堂在线| 午夜福利在线在线| 亚洲精品国产成人久久av| 插阴视频在线观看视频| 性欧美人与动物交配| 啦啦啦韩国在线观看视频| 夫妻性生交免费视频一级片| 久久久午夜欧美精品| 国产精品一及| 亚洲国产精品成人综合色| 欧美成人a在线观看| 国产一区二区在线观看日韩| 97超视频在线观看视频| 久久久精品欧美日韩精品| 一进一出抽搐动态| 欧美不卡视频在线免费观看| 99久久中文字幕三级久久日本| 久久久国产成人免费| 亚洲成av人片在线播放无| 精品久久久久久久久亚洲| 在线观看午夜福利视频| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久亚洲| 亚洲va在线va天堂va国产| 久久这里有精品视频免费| 国产真实伦视频高清在线观看| 天天躁日日操中文字幕| 日韩一区二区视频免费看| 黄片无遮挡物在线观看| videossex国产| 一级毛片aaaaaa免费看小| 自拍偷自拍亚洲精品老妇| 亚洲av熟女| 丝袜喷水一区| 亚洲精品日韩av片在线观看| 日韩亚洲欧美综合| 欧美+日韩+精品| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 亚洲成av人片在线播放无| 最新中文字幕久久久久| 国产麻豆成人av免费视频| 两性午夜刺激爽爽歪歪视频在线观看| av又黄又爽大尺度在线免费看 | 久久午夜亚洲精品久久| 一个人看视频在线观看www免费| 久久精品国产亚洲av香蕉五月| 麻豆成人午夜福利视频| 看十八女毛片水多多多| 乱人视频在线观看| 亚洲在线自拍视频| 日本色播在线视频| 国产亚洲精品久久久com| 在线播放无遮挡| 国产高清激情床上av| 国产淫片久久久久久久久| 欧美激情国产日韩精品一区| 国产国拍精品亚洲av在线观看| 日本黄色片子视频| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 国产av一区在线观看免费| 久久热精品热| 插阴视频在线观看视频| 国产成人freesex在线| 高清日韩中文字幕在线| 特大巨黑吊av在线直播| 中文字幕av成人在线电影| 国产高清有码在线观看视频| 少妇的逼水好多| 国产不卡一卡二| 国内精品宾馆在线| 婷婷亚洲欧美| 嫩草影院精品99| 久久精品国产鲁丝片午夜精品| 大型黄色视频在线免费观看| 久久久久久九九精品二区国产| 中国国产av一级| 尤物成人国产欧美一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久久久亚洲| 久久6这里有精品| 国产精品久久电影中文字幕| 一区二区三区高清视频在线| 岛国在线免费视频观看| 18+在线观看网站| 国产精品无大码| 久久6这里有精品| 国产精品一区www在线观看| 激情 狠狠 欧美| 晚上一个人看的免费电影| 亚洲成人av在线免费| 又粗又硬又长又爽又黄的视频 | 欧美最黄视频在线播放免费| 欧美三级亚洲精品| 天堂√8在线中文| 99在线视频只有这里精品首页| 啦啦啦观看免费观看视频高清| 国产乱人视频| 嫩草影院新地址| 又爽又黄无遮挡网站| 免费无遮挡裸体视频| 五月玫瑰六月丁香| 人妻系列 视频| 国产伦理片在线播放av一区 | 日本爱情动作片www.在线观看| 国产单亲对白刺激| 国产免费男女视频| 只有这里有精品99| 免费看美女性在线毛片视频| 亚洲精品456在线播放app| 天天一区二区日本电影三级| www日本黄色视频网| 国产免费视频播放在线视频| 亚洲人成网站在线播| 精品酒店卫生间| 九九久久精品国产亚洲av麻豆| 丝袜美足系列| 国产色爽女视频免费观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲av.av天堂| 一区二区三区乱码不卡18| 国产精品麻豆人妻色哟哟久久| 国产免费现黄频在线看| kizo精华| 男女边摸边吃奶| 亚洲国产精品一区三区| 五月天丁香电影| 老女人水多毛片| 在线播放无遮挡| 国产 一区精品| 国产在视频线精品| 大又大粗又爽又黄少妇毛片口| 国产黄色免费在线视频| 国产免费视频播放在线视频| av在线老鸭窝| 免费av不卡在线播放| 久久久久视频综合| 国产成人aa在线观看| 人人妻人人爽人人添夜夜欢视频| 国产极品粉嫩免费观看在线 | 亚洲成人一二三区av| 黄色一级大片看看| 国产成人免费观看mmmm| 男男h啪啪无遮挡| 夜夜爽夜夜爽视频| 22中文网久久字幕| 成人漫画全彩无遮挡| 免费黄网站久久成人精品| 边亲边吃奶的免费视频| 久久ye,这里只有精品| 99热全是精品| 日韩成人伦理影院| 蜜臀久久99精品久久宅男| 青青草视频在线视频观看| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 日韩av不卡免费在线播放| 少妇的逼好多水| 性色avwww在线观看| kizo精华| 日韩制服骚丝袜av| 2021少妇久久久久久久久久久| 伦精品一区二区三区| 亚洲一级一片aⅴ在线观看| 少妇熟女欧美另类| 免费高清在线观看日韩| 欧美另类一区| 三上悠亚av全集在线观看| 国产精品.久久久| 亚洲av福利一区| 欧美亚洲日本最大视频资源| 熟女电影av网| av一本久久久久| 丝瓜视频免费看黄片| 亚洲精品aⅴ在线观看| 中文字幕最新亚洲高清| 久久精品夜色国产| 十八禁网站网址无遮挡| 大码成人一级视频| 日韩一区二区三区影片| 免费高清在线观看日韩| 精品久久国产蜜桃| 亚洲精品456在线播放app| 日本免费在线观看一区| 国产精品欧美亚洲77777| 妹子高潮喷水视频| 日本黄大片高清| 丰满迷人的少妇在线观看| 99久久人妻综合| 久久99一区二区三区| 久久综合国产亚洲精品| 少妇 在线观看| 成人综合一区亚洲| 久久热精品热| 日日摸夜夜添夜夜爱| 国产色爽女视频免费观看| 大香蕉久久成人网| 婷婷色综合www| 国产伦理片在线播放av一区| 国产精品一国产av| av电影中文网址| 18+在线观看网站| 亚洲av男天堂| 国产精品久久久久久精品电影小说| 蜜桃在线观看..| 少妇丰满av| 中文乱码字字幕精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻熟女av久视频| 在线精品无人区一区二区三| 亚洲精品日韩av片在线观看| 精品少妇内射三级| 中文字幕久久专区| 交换朋友夫妻互换小说| 综合色丁香网| 免费黄频网站在线观看国产| 精品人妻偷拍中文字幕| 国产免费视频播放在线视频| 一级毛片aaaaaa免费看小| 美女主播在线视频| 国产亚洲精品久久久com| 国产伦理片在线播放av一区| 国产成人免费观看mmmm| 一边摸一边做爽爽视频免费| 国产白丝娇喘喷水9色精品| 日韩av免费高清视频| 欧美日韩亚洲高清精品| 男女免费视频国产| 国产精品女同一区二区软件| 国产成人精品一,二区| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩一区二区| 欧美精品国产亚洲| 丰满饥渴人妻一区二区三| www.av在线官网国产| 日日撸夜夜添| 久久人人爽人人片av| 熟女电影av网| av天堂久久9| 亚洲国产成人一精品久久久| 交换朋友夫妻互换小说| 在线观看免费高清a一片| 成人综合一区亚洲| 欧美人与性动交α欧美精品济南到 | 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| 国产成人精品久久久久久| 热99国产精品久久久久久7| 亚洲欧美一区二区三区黑人 | 嫩草影院入口| 精品国产露脸久久av麻豆| 午夜老司机福利剧场| 精品国产露脸久久av麻豆| 婷婷色综合www| 免费观看av网站的网址| a级毛片黄视频| 母亲3免费完整高清在线观看 | 免费av中文字幕在线| av专区在线播放| 男女边吃奶边做爰视频| 尾随美女入室| 丰满迷人的少妇在线观看| 伦精品一区二区三区| 色视频在线一区二区三区| 十八禁网站网址无遮挡| 亚洲成人手机| 在线天堂最新版资源| 人妻人人澡人人爽人人| 永久免费av网站大全| 国产成人aa在线观看| 婷婷色麻豆天堂久久| 国产精品一二三区在线看| 如日韩欧美国产精品一区二区三区 | 亚洲精华国产精华液的使用体验| 免费大片18禁| 成人免费观看视频高清| 亚洲国产最新在线播放| 成人毛片60女人毛片免费| 丝袜脚勾引网站| 精品人妻偷拍中文字幕| 妹子高潮喷水视频| 亚洲av成人精品一二三区| 极品少妇高潮喷水抽搐| 少妇的逼好多水| 热99久久久久精品小说推荐| 日本爱情动作片www.在线观看| 大香蕉久久网| 一本大道久久a久久精品| 高清黄色对白视频在线免费看| 简卡轻食公司| 欧美亚洲 丝袜 人妻 在线| 亚洲av综合色区一区| 精品国产乱码久久久久久小说| 亚洲图色成人| 亚洲av福利一区| 日韩欧美一区视频在线观看| 成人漫画全彩无遮挡| 欧美日本中文国产一区发布| 能在线免费看毛片的网站| 国产成人精品久久久久久| 午夜激情av网站| av在线观看视频网站免费| 嘟嘟电影网在线观看| 亚洲av成人精品一二三区| 在线观看人妻少妇| √禁漫天堂资源中文www| 中文字幕免费在线视频6| 久久99精品国语久久久| 日韩av在线免费看完整版不卡| 免费av中文字幕在线| av在线app专区| 亚洲精品久久久久久婷婷小说| 麻豆精品久久久久久蜜桃| 欧美日韩视频精品一区| 久久精品国产亚洲av天美| 爱豆传媒免费全集在线观看| 另类亚洲欧美激情| 亚洲精品av麻豆狂野| 九九在线视频观看精品| 国产黄色视频一区二区在线观看| 国产 一区精品| 久久久久久久久久久丰满| 国产一区亚洲一区在线观看| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 国产熟女欧美一区二区| 少妇人妻精品综合一区二区| 一区在线观看完整版| 亚洲精品一区蜜桃| 最新中文字幕久久久久| 一边摸一边做爽爽视频免费| 日韩av不卡免费在线播放| 国产片特级美女逼逼视频| 精品人妻偷拍中文字幕| 香蕉精品网在线| 日日撸夜夜添| 赤兔流量卡办理| 久久精品久久精品一区二区三区| 欧美最新免费一区二区三区| 中文乱码字字幕精品一区二区三区| 赤兔流量卡办理| 亚洲欧美色中文字幕在线| 制服诱惑二区| 国产视频内射| 亚洲av中文av极速乱| 黄色视频在线播放观看不卡| 成年美女黄网站色视频大全免费 | 男人添女人高潮全过程视频| 久久久久视频综合| 亚洲美女搞黄在线观看| 国产一区二区在线观看日韩| 国产无遮挡羞羞视频在线观看| 最新的欧美精品一区二区| 国产黄色视频一区二区在线观看| 国产成人精品无人区| 亚洲精品日本国产第一区| 蜜臀久久99精品久久宅男| 中文字幕免费在线视频6| 一级毛片 在线播放| 赤兔流量卡办理| 中文字幕av电影在线播放| 插逼视频在线观看| 97精品久久久久久久久久精品| 久久精品国产自在天天线| 一级爰片在线观看| 天堂8中文在线网| 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 日韩精品有码人妻一区| 亚洲丝袜综合中文字幕| 日日爽夜夜爽网站| 人妻夜夜爽99麻豆av| 亚洲av综合色区一区| 成人国产av品久久久| 欧美精品国产亚洲| 99国产精品免费福利视频| 午夜老司机福利剧场| 99精国产麻豆久久婷婷| 最近2019中文字幕mv第一页| 视频区图区小说| 男女无遮挡免费网站观看| 校园人妻丝袜中文字幕| 免费观看在线日韩| 性色av一级| 成人国产av品久久久| 日韩三级伦理在线观看| 麻豆成人av视频| 插阴视频在线观看视频| 久久99精品国语久久久| 精品亚洲成国产av| 男女国产视频网站| 免费av不卡在线播放| 狂野欧美白嫩少妇大欣赏| 亚洲欧洲国产日韩| 一区二区日韩欧美中文字幕 | 精品一区二区三区视频在线| 哪个播放器可以免费观看大片| 午夜影院在线不卡| 日本爱情动作片www.在线观看| 午夜福利在线观看免费完整高清在| 美女cb高潮喷水在线观看| 国产精品偷伦视频观看了| 极品少妇高潮喷水抽搐| 亚洲av电影在线观看一区二区三区| 国产白丝娇喘喷水9色精品| 精品人妻熟女av久视频| 人人妻人人澡人人看| 午夜影院在线不卡| 欧美精品国产亚洲| av专区在线播放| 久久影院123| www.av在线官网国产| 亚洲国产最新在线播放| 欧美日韩成人在线一区二区| 丝袜喷水一区| 天美传媒精品一区二区| 人妻夜夜爽99麻豆av| 亚洲国产成人一精品久久久| 色哟哟·www| 久久精品人人爽人人爽视色| 51国产日韩欧美| 精品久久久久久电影网| 欧美老熟妇乱子伦牲交| 美女cb高潮喷水在线观看| 如日韩欧美国产精品一区二区三区 | 久久久久久久久久人人人人人人| 最近的中文字幕免费完整| 日本wwww免费看| 少妇被粗大的猛进出69影院 | 久久久久国产精品人妻一区二区| 精品一区二区免费观看| 国产伦理片在线播放av一区| 国产亚洲精品久久久com| 亚洲色图 男人天堂 中文字幕 | 飞空精品影院首页| 亚洲怡红院男人天堂| 黄色欧美视频在线观看| 中文字幕av电影在线播放| 国产精品 国内视频| 久久精品国产鲁丝片午夜精品| 在线观看国产h片| 美女脱内裤让男人舔精品视频| 三级国产精品片| 日韩三级伦理在线观看| 国产欧美亚洲国产| 黄色视频在线播放观看不卡| 极品人妻少妇av视频| 国产毛片在线视频| 色婷婷久久久亚洲欧美| 成人国产麻豆网| kizo精华| 日本欧美国产在线视频| 最新中文字幕久久久久| 亚洲精华国产精华液的使用体验| 高清欧美精品videossex| 欧美人与性动交α欧美精品济南到 | 欧美性感艳星| 三级国产精品欧美在线观看| 欧美成人精品欧美一级黄| √禁漫天堂资源中文www| 国产av精品麻豆| 婷婷成人精品国产| 色哟哟·www| 99精国产麻豆久久婷婷| 国产一级毛片在线| 国国产精品蜜臀av免费| 99视频精品全部免费 在线| 999精品在线视频| 夫妻午夜视频| 熟女电影av网| 久久久久久久大尺度免费视频| 有码 亚洲区| 美女脱内裤让男人舔精品视频| 999精品在线视频| 国国产精品蜜臀av免费| 在线天堂最新版资源| 欧美日韩亚洲高清精品| 亚洲欧洲精品一区二区精品久久久 | 色吧在线观看| 久久鲁丝午夜福利片| 午夜激情福利司机影院| 成人国产av品久久久| 久久青草综合色| 亚洲精品乱码久久久久久按摩| 国产高清三级在线| 在线观看免费高清a一片| 在线观看免费日韩欧美大片 | 国产高清有码在线观看视频| 男女边摸边吃奶|