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

    基于多旅行商問題的應(yīng)急設(shè)施服務(wù)區(qū)劃分模型

    2020-10-31 03:30:22星,吉康,申
    關(guān)鍵詞:集散中心服務(wù)區(qū)旅行

    趙 星,吉 康,申 珂

    (河海大學(xué)土木與交通學(xué)院,南京210098)

    0 引 言

    近年來,全球范圍內(nèi)各類突發(fā)公共衛(wèi)生事件頻繁發(fā)生,例如,2002年SARS 事件,2019年新型冠狀病毒事件.面對此類區(qū)域性的突發(fā)公共衛(wèi)生事件,疫區(qū)的應(yīng)急設(shè)施服務(wù)區(qū)劃分,例如臨時醫(yī)院,物資周轉(zhuǎn)中心,避難所等,可以幫助物資配送及緊急救援,是控制疫情發(fā)展,保障民眾生活需要的基本前提.

    應(yīng)急設(shè)施位置選擇對服務(wù)區(qū)域劃分具有重要影響,較多研究將應(yīng)急設(shè)施選址問題與服務(wù)區(qū)域劃分問題聯(lián)合求解.典型選址—分配(Location-Allocation,L-A)模型根據(jù)其優(yōu)化目標(biāo)不同可以分為:P-中值模型,P-中心模型,覆蓋模型,均已有成熟的研究應(yīng)用[1].Pan[2]以各個疏散集結(jié)點到避難所的加權(quán)距離和最小為優(yōu)化目標(biāo),構(gòu)建了一種基于P-中值問題的避難所選址—分配模型,采用遺傳算法進(jìn)行求解,可同步獲得避難所最優(yōu)選址與服務(wù)區(qū)劃分方案.郭鵬輝等[3]研究了災(zāi)后救援物資選址—路徑—配給的集成優(yōu)化問題,并基于進(jìn)化算法對問題進(jìn)行求解.

    在應(yīng)急設(shè)施位置已知的情況下,設(shè)施服務(wù)區(qū)劃分可轉(zhuǎn)化為多旅行商問題求解.多旅行商問題(Multiple Traveling Salesman Problem,MTSP)是在傳統(tǒng)的旅行商問題(Travelling Salesman Problem,TSP)基礎(chǔ)上引入多個旅行商和對應(yīng)多個出發(fā)城市的變種問題.事實上,存在多個應(yīng)急設(shè)施的服務(wù)區(qū)劃分問題,與多起點的MTSP(旅行商位于不同的出發(fā)城市)存在共通之處.通常可將應(yīng)急設(shè)施視為虛擬的旅行商,應(yīng)急設(shè)施服務(wù)容量為旅行商攜帶的商品數(shù)量,M為應(yīng)急設(shè)施個數(shù),則虛擬旅行商的旅行路徑即可匹配為應(yīng)急設(shè)施的服務(wù)區(qū)域.通過應(yīng)急設(shè)施服務(wù)區(qū)劃分,可以將多倉庫多約束的應(yīng)急調(diào)度問題轉(zhuǎn)換為多個單倉庫的應(yīng)急調(diào)度問題,并預(yù)處理相關(guān)約束條件,再利用傳統(tǒng)的VRP 或TSP 方法進(jìn)行應(yīng)急路徑規(guī)劃.針對此類問題,Steven等[4]采用聚集聚類方法對多起點的MTSP進(jìn)行求解.Liu等[5]利用圖論將MTSP分解為TSP,再進(jìn)行模型最優(yōu)解的獲取.此外,Alencar等[6]具體研究MTSP在物資配送、路徑規(guī)劃方面的實際應(yīng)用,并提出了相應(yīng)算法對問題求解,基本思路為先通過旅行商旅行路徑劃分服務(wù)區(qū)域,再進(jìn)行路徑規(guī)劃.

    基于MTSP 的應(yīng)急設(shè)施服務(wù)區(qū)劃分研究較為不足,可以看出:MTSP 契合多設(shè)施多集散中心的應(yīng)急設(shè)施服務(wù)區(qū)劃分問題特征,具有重要的研究應(yīng)用價值.本文構(gòu)建一種基于MTSP的應(yīng)急設(shè)施服務(wù)區(qū)劃分模型,并利用針對P-中值選址模型的整數(shù)規(guī)劃方法及禁忌搜索算法,結(jié)合LKH 求解器提出一種復(fù)合算法對模型進(jìn)行求解,從而為物資配送及各類緊急救援方案的制定提供依據(jù).

    1 模型構(gòu)建

    基于MTSP 構(gòu)建一種應(yīng)急設(shè)施服務(wù)區(qū)劃分模型,將各應(yīng)急設(shè)施視為虛擬旅行商,在網(wǎng)絡(luò)需求,資源供應(yīng),流量平衡等約束條件下,以最小化遍歷區(qū)域內(nèi)全部集散中心的綜合旅行時間成本為優(yōu)化目標(biāo),獲取模型最優(yōu)的區(qū)域應(yīng)急設(shè)施服務(wù)區(qū)劃分方案.

    1.1 基本假設(shè)

    (1)具備先進(jìn)的應(yīng)急調(diào)度系統(tǒng)及信息更新方式,所有應(yīng)急設(shè)施信息共享,均服從系統(tǒng)的聯(lián)合管理,以綜合效益最大化進(jìn)行服務(wù)區(qū)劃分.各應(yīng)急設(shè)施的位置及資源儲存情況已知,且能夠滿足區(qū)域內(nèi)的全部資源需求.

    (2)各個集散中心的位置及資源需求已知,且基于公平原則進(jìn)行資源獲取.

    (3)應(yīng)急車輛的每次運輸均從其所屬的應(yīng)急設(shè)施出發(fā),前往服務(wù)區(qū)內(nèi)的集散中心發(fā)放物資,并最終返回該應(yīng)急設(shè)施.

    1.2 參數(shù)及變量

    (1)參 數(shù).

    N——節(jié)點集合,N=E∪D;

    Z——應(yīng)急設(shè)施服務(wù)區(qū)域集合;

    D——應(yīng)急設(shè)施集合;

    E——集散中心集合;

    nz——應(yīng)急設(shè)施服務(wù)區(qū)數(shù)量;

    nd——應(yīng)急設(shè)施數(shù)量;

    ne——集散中心數(shù)量;

    ——集散中心i的資源需求數(shù)量,i∈E;

    rj——應(yīng)急設(shè)施j的資源供應(yīng)數(shù)量,j∈D;

    dxy——鄰接矩陣中節(jié)點x到節(jié)點y的應(yīng)急車輛行程時間,x,y∈N;

    (2)變 量.

    決策變量:

    ——若,表示集散中心i被分配至應(yīng)急設(shè)施j;若,表示集散中心i未被分配至應(yīng)急設(shè)施j.

    中間變量:

    μxy——旅行商從節(jié)點x出發(fā)前往節(jié)點y的次數(shù),若μxy=0,表示旅行商未從節(jié)點x出發(fā)前往節(jié)點y;

    Ej——應(yīng)急設(shè)施j所服務(wù)的集散中心集合;

    Pj——應(yīng)急設(shè)施j所在的設(shè)施服務(wù)區(qū)的全部節(jié)點集合,Pj=Ej∪{j} .

    1.3 模型結(jié)構(gòu)

    (1)目標(biāo)函數(shù).

    (2)約束條件.

    式(1)為模型優(yōu)化目標(biāo)函數(shù),要求最小化遍歷區(qū)域內(nèi)全部集散中心的綜合旅行時間成本.式(2)確保每個服務(wù)區(qū)均圍繞一個應(yīng)急設(shè)施進(jìn)行劃分.式(3)確保在同一個應(yīng)急設(shè)施服務(wù)區(qū)內(nèi),任意集散中心必然與該服務(wù)區(qū)內(nèi)的應(yīng)急設(shè)施或其他集散中心鄰接,換而言之,保證設(shè)施服務(wù)區(qū)的聚團(tuán)劃分.式(4)確保任意一個應(yīng)急設(shè)施服務(wù)區(qū)內(nèi)的集散中心物資需求總量不能超過該服務(wù)區(qū)應(yīng)急設(shè)施的資源供應(yīng)數(shù)量.式(5)確保各個應(yīng)急設(shè)施被分配的集散中心數(shù)量之和等于集散中心總數(shù).式(6)要求每個集散中心僅被分配至一個應(yīng)急設(shè)施.因此,式(5)和式(6)共同構(gòu)成集散中心全覆蓋約束.式(7)確保網(wǎng)絡(luò)中的應(yīng)急設(shè)施資源總量可以滿足全部資源集散需求.式(8)為流量平衡約束,保證任意節(jié)點的旅行商抵達(dá)次數(shù)等于該節(jié)點的旅行商出發(fā)次數(shù).

    2 求解算法

    MTSP 是TSP 的變種形式,屬于組合優(yōu)化問題,在計算復(fù)雜度上屬于NP 困難問題,難以精確求解.因此,本文基于P-中值選址模型的優(yōu)化理念及禁忌搜索(Tabu Search,TS)算法,結(jié)合LKH求解器(Lin-Kernighan-Helsgaun TSP solver)對所構(gòu)建的應(yīng)急設(shè)施服務(wù)區(qū)劃分模型求解.

    禁忌搜索算法是一種全局逐步尋優(yōu)的亞啟發(fā)式算法,可有效求解NP困難問題[7].LKH求解器是當(dāng)前求解TSP 及其變種問題的有效工具,其求解速度快,解的優(yōu)越性高,可以為目前所有能夠精確求解的TSP實例找到最優(yōu)解.

    首先參考P-中值選址模型的優(yōu)化理念,利用整數(shù)規(guī)劃方法形成初步可行的設(shè)施服務(wù)區(qū)劃分方案;然后基于禁忌搜索算法的理念,通過迭代方式,不斷選取集散中心執(zhí)行移動、交換操作;調(diào)用LKH 求解器計算目標(biāo)函數(shù)值的變化,使初始方案向目標(biāo)函數(shù)優(yōu)化方向搜索,最終獲得模型的局部最優(yōu)解.算法流程如圖1所示,詳細(xì)步驟如下.

    Step 0初始化.輸入應(yīng)急設(shè)施集合D,設(shè)施資源供應(yīng)數(shù)量rj,集散中心集合E,各個集散中心的資源需求,網(wǎng)絡(luò)中全部鄰接節(jié)點之間的行程時間dxy.

    Step 1參考P-中值選址模型,在約束條件下,以最小化總資源旅行成本,即集散中心到應(yīng)急設(shè)施的旅行時間成本與集散中心資源需求數(shù)量的乘積和為優(yōu)化目標(biāo),初步劃分應(yīng)急設(shè)施服務(wù)區(qū).具體而言,可以通過Dijkstra算法獲取各個集散中心到各個應(yīng)急設(shè)施的最短路徑,然后將集散中心資源需求數(shù)量作為權(quán)重系數(shù),以歸屬關(guān)系作為0-1變量直接求解,形成服務(wù)區(qū)劃分初始可行解.

    Step 2將應(yīng)急設(shè)施j作為虛擬旅行商,以總旅行時間成本最小為目標(biāo)對其服務(wù)區(qū)域內(nèi)的集散中心進(jìn)行遍歷,即將問題轉(zhuǎn)化為TSP問題.在此基礎(chǔ)上,通過LKH 求解器快速求解其最小的TSP 時間成本.以此類推,獲得各應(yīng)急設(shè)施服務(wù)區(qū)的最小TSP時間成本,得到初始可行解的目標(biāo)函數(shù)值Ttsp.

    Step 3將當(dāng)前模型最優(yōu)解作為禁忌對象,執(zhí)行禁忌搜索.設(shè)置備選搜索步長為{0,1,2,3},禁忌長度為4,最大迭代次數(shù)為nmax,T0=Ttsp,n=0,l=0.

    Step 4從當(dāng)前模型最優(yōu)解中隨機選擇應(yīng)急設(shè)施服務(wù)區(qū)域z1,從剩余區(qū)域集合中隨機選擇一個與z1鄰接的區(qū)域z2.執(zhí)行z1→z2,根據(jù)鄰接矩陣篩選出區(qū)域z1中與區(qū)域z2直接相連的集散中心,并將這些集散中心納入集合r1.隨機選擇搜索步長Mmove1,r1中隨機選擇Mmove1個集散中心,并將這些集散中心從區(qū)域z1中移到區(qū)域z2中.

    Step 5更新區(qū)域z1,z2.執(zhí)行z2→z1,根據(jù)鄰接矩陣篩選出區(qū)域z2中與區(qū)域z1直接相連的集散中心,并將這些集散中心納入集合r2.隨機選擇搜索步長Mmove2,r2中隨機選擇Mmove2個集散中心,并將這些集散中心從區(qū)域z2中移到區(qū)域z1中.

    圖1 算法流程Fig.1 Algorithm flowchart

    Step 6基于式(3),檢驗集散中心的連續(xù)性,確保服務(wù)區(qū)域聚團(tuán)分布. 如 果?x∈Pj,j∈D,?y∈Pj使得dxy≠∞,則進(jìn)入Step 7;反之,重新進(jìn)入Step 4.

    Step 7基于式(4),檢驗應(yīng)急設(shè)施資源容量約束.如果,則進(jìn)入Step 8;反之,重新進(jìn)入Step 4.

    Step 8計算當(dāng)前候選解的目標(biāo)函數(shù)值Ttsp,檢驗是否滿足特赦規(guī)則.若Ttsp<T0,則直接將該候選解替換為模型最優(yōu)解,更新禁忌列表,令T0=Ttsp,n=n+1,l=1,并重新進(jìn)入Step 4;反之,令l=l+1,進(jìn)入Step 9.

    Step 9檢驗是否已達(dá)到禁忌長度. 如果l≥4,則進(jìn)入Step 10;反之,重新進(jìn)入Step 4.

    Step 10以模型目標(biāo)函數(shù)作為評價標(biāo)準(zhǔn),選擇候選集合中的最優(yōu)解及其目標(biāo)函數(shù)值Ttsp,替換該候選最優(yōu)解為模型最優(yōu)解,更新禁忌列表,令T0=Ttsp,n=n+1,l=0,并重新進(jìn)入Step 4.

    終止條件:迭代次數(shù)n≥nmax.此時,對應(yīng)于當(dāng)前T0的應(yīng)急設(shè)施服務(wù)區(qū)劃分方案即為模型最優(yōu)解.

    3 案例分析

    基于浙江省寧波市北侖區(qū)實際拓?fù)渚W(wǎng)絡(luò)進(jìn)行案例分析.網(wǎng)絡(luò)布局如圖2 所示,該區(qū)域內(nèi)共有39條路段,23個交叉口,以及分別位于節(jié)點24,25,26的3 處應(yīng)急設(shè)施,具體道路參數(shù)如表1 所示.現(xiàn)假設(shè)該區(qū)域內(nèi)發(fā)生突發(fā)公共衛(wèi)生事件,并產(chǎn)生應(yīng)急資源調(diào)度需求.為確保資源集散的可達(dá)性,集散中心設(shè)置于各路段中點,網(wǎng)絡(luò)中各集散中心的資源需求數(shù)量及應(yīng)急設(shè)施的資源供應(yīng)數(shù)量如表2和表3所示(以防護(hù)服為例).

    圖2 北侖區(qū)道路網(wǎng)絡(luò)Fig.2 Road network of Beilun District

    表1 網(wǎng)絡(luò)道路屬性表Table 1 Road property in network

    表2 集散中心資源需求數(shù)量Table 2 Resource demand of distribution centers

    表3 應(yīng)急設(shè)施資源供應(yīng)數(shù)量Table 3 Resource supply of emergency facilities

    假設(shè)網(wǎng)絡(luò)上設(shè)有應(yīng)急專用車道,應(yīng)急車輛路段通行速度為60 km/h,獲得模型最優(yōu)解如圖3 所示.該方案滿足模型的全部約束條件,劃分出的3個服務(wù)區(qū)域Ⅰ、Ⅱ、Ⅲ的TSP旅行時間成本分別為632,895,825 s,服務(wù)區(qū)資源需求分別為1 229,964,1 671 件,方案遍歷區(qū)域內(nèi)全部集散中心的綜合旅行時間成本Ttsp=2 352 s.

    圖3 模型最優(yōu)服務(wù)區(qū)劃分方案Fig.3 Model optimal service area division scheme

    為證明所提模型與算法的優(yōu)越性,基于參考文獻(xiàn)[4],以應(yīng)急設(shè)施到所服務(wù)集散中心的總旅行時間成本最小為優(yōu)化目標(biāo),在相同的模型約束條件下,采用聚類蟻群算法進(jìn)行服務(wù)區(qū)劃分對照,獲得對照方案如圖4所示.由圖4可知,劃分出的3個服務(wù)區(qū)域Ⅰ、Ⅱ、Ⅲ的TSP 旅行時間成本分別為591,1 088,873 s,服務(wù)區(qū)資源需求分別為1 267,936,1 661 件,遍歷區(qū)域內(nèi)全部集散中心的綜合旅行時間成本Ttsp=2 552 s.

    圖4 服務(wù)區(qū)劃分參照方案Fig.4 Service area division control scheme

    對比本文模型最優(yōu)方案與參照方案可以發(fā)現(xiàn):在模型最優(yōu)方案中,盡管部分集散中心被分配到較遠(yuǎn)的應(yīng)急設(shè)施(如集散中心52距離應(yīng)急設(shè)施25、26的時間成本分別為368 s、270 s,在優(yōu)化過程中被從應(yīng)急設(shè)施26分配至應(yīng)急設(shè)施25),但是模型最優(yōu)解的目標(biāo)函數(shù)得到一定程度的優(yōu)化,Ttsp降低了7.84%.這意味著,在各個應(yīng)急設(shè)施服務(wù)區(qū)內(nèi),車輛從應(yīng)急設(shè)施出發(fā)遍歷該區(qū)域內(nèi)全部集散中心的綜合時間成本得到降低,有利于減少運輸時間成本,從整體上提高物資配送或緊急救援方案的效率.

    在此基礎(chǔ)上,本文進(jìn)一步假設(shè)應(yīng)急設(shè)施24、25、26 處各有1 輛應(yīng)急運輸車輛,其車輛滿載為400件防護(hù)服,基于簡單插入算法[8],獲得綜合物資配送方案如表4 所示,其中,車輛每次從應(yīng)急設(shè)施出發(fā)時更新為滿載,不足滿載時則更新為應(yīng)急設(shè)施剩余物資容量,配送完畢后返回應(yīng)急設(shè)施.

    表4 綜合物資配送方案Table 4 Resource distribution scheme

    4 結(jié) 論

    本文構(gòu)建了一種基于多旅行商問題的應(yīng)急設(shè)施服務(wù)區(qū)劃分模型,將各應(yīng)急設(shè)施視為虛擬旅行商,以最小化遍歷區(qū)域內(nèi)全部集散中心的綜合旅行時間成本為優(yōu)化目標(biāo)劃分應(yīng)急設(shè)施服務(wù)區(qū)域,并設(shè)計一種復(fù)合算法對模型進(jìn)行求解.經(jīng)過實例驗證,該模型可以有效解決應(yīng)急設(shè)施服務(wù)區(qū)劃分問題,提高物資配送或緊急救援方案的效率.下一步將在本文研究基礎(chǔ)上,探究如何采用更合理的方法進(jìn)行各應(yīng)急設(shè)施物資配送或緊急救援的路徑規(guī)劃,以及應(yīng)急設(shè)施位置未知情況下的應(yīng)急設(shè)施選址—路徑問題.

    猜你喜歡
    集散中心服務(wù)區(qū)旅行
    基于AIoT+GIS的智慧服務(wù)區(qū)構(gòu)建
    高速公路服務(wù)區(qū)信息技術(shù)的應(yīng)用
    農(nóng)產(chǎn)品集散中心選址問題研究
    卷宗(2019年29期)2019-11-11 12:18:11
    包裹如山
    不可能旅行
    小黑的旅行
    建言高速公路服務(wù)區(qū)實現(xiàn)“雙提升”
    中國公路(2017年5期)2017-06-01 12:10:10
    小黑去旅行
    夏日旅行
    淺談威海旅游集散中心導(dǎo)游人員流失的原因及對策分析
    科技視界(2014年1期)2014-07-29 06:39:35
    国国产精品蜜臀av免费| 99热6这里只有精品| 亚洲天堂av无毛| 一级毛片aaaaaa免费看小| 国产欧美日韩一区二区三区在线 | 我的老师免费观看完整版| 国产在线视频一区二区| 国产高清三级在线| 韩国av在线不卡| 韩国高清视频一区二区三区| 久久久色成人| 嫩草影院入口| 一区二区三区四区激情视频| 美女xxoo啪啪120秒动态图| av一本久久久久| 国产国拍精品亚洲av在线观看| 国产成人91sexporn| 亚洲天堂av无毛| 国产黄频视频在线观看| 在线观看av片永久免费下载| 三级国产精品欧美在线观看| 国产女主播在线喷水免费视频网站| 欧美xxⅹ黑人| 性高湖久久久久久久久免费观看| 亚洲欧美日韩东京热| 欧美+日韩+精品| 亚洲性久久影院| 亚洲精华国产精华液的使用体验| 亚洲精品成人av观看孕妇| 国产精品女同一区二区软件| 欧美一级a爱片免费观看看| 2022亚洲国产成人精品| 欧美xxⅹ黑人| 亚洲av国产av综合av卡| 国产黄片美女视频| 日韩成人伦理影院| 国产伦在线观看视频一区| 另类亚洲欧美激情| 日韩电影二区| 一级爰片在线观看| 免费黄频网站在线观看国产| 欧美 日韩 精品 国产| 亚洲精品日本国产第一区| 亚洲国产最新在线播放| 欧美成人a在线观看| 丰满少妇做爰视频| 午夜福利视频精品| av在线播放精品| 91在线精品国自产拍蜜月| 综合色丁香网| 少妇裸体淫交视频免费看高清| 麻豆国产97在线/欧美| 多毛熟女@视频| 国产精品爽爽va在线观看网站| 亚洲精品国产色婷婷电影| 男女国产视频网站| 狂野欧美激情性xxxx在线观看| 视频中文字幕在线观看| 欧美性感艳星| 国产成人免费无遮挡视频| av国产久精品久网站免费入址| 99久久精品国产国产毛片| 日本欧美国产在线视频| 夜夜骑夜夜射夜夜干| 亚洲va在线va天堂va国产| 精品久久久久久久末码| 卡戴珊不雅视频在线播放| 91久久精品国产一区二区成人| 精品久久久精品久久久| 日韩免费高清中文字幕av| 国产片特级美女逼逼视频| 免费在线观看成人毛片| 亚洲精品中文字幕在线视频 | 亚洲国产色片| 国产美女午夜福利| 亚洲av综合色区一区| 久久精品久久久久久久性| 国产69精品久久久久777片| 久久久色成人| 99久久综合免费| 久久久久久久久久人人人人人人| 18+在线观看网站| 男人添女人高潮全过程视频| 高清在线视频一区二区三区| 99久久人妻综合| 狠狠精品人妻久久久久久综合| 99热这里只有精品一区| 青春草国产在线视频| 亚洲四区av| 亚洲国产色片| 大香蕉97超碰在线| 欧美成人一区二区免费高清观看| 女性被躁到高潮视频| 好男人视频免费观看在线| 久久 成人 亚洲| 97在线视频观看| 日韩中字成人| 99热网站在线观看| 一区二区av电影网| 一区在线观看完整版| 永久免费av网站大全| 网址你懂的国产日韩在线| 午夜福利高清视频| 亚洲av欧美aⅴ国产| 美女脱内裤让男人舔精品视频| 九九在线视频观看精品| 国产成人一区二区在线| 久久久久久久亚洲中文字幕| 香蕉精品网在线| 国产亚洲精品久久久com| 欧美成人a在线观看| 免费观看a级毛片全部| 午夜福利在线观看免费完整高清在| 国产免费一级a男人的天堂| 成人一区二区视频在线观看| 内地一区二区视频在线| 国产成人a区在线观看| 久久久久久久久久久丰满| 女的被弄到高潮叫床怎么办| 久久精品国产a三级三级三级| 麻豆成人午夜福利视频| 亚洲精品成人av观看孕妇| 日韩一区二区三区影片| 最近的中文字幕免费完整| 国产69精品久久久久777片| 一级av片app| 麻豆乱淫一区二区| 国产精品一二三区在线看| 最近2019中文字幕mv第一页| 春色校园在线视频观看| 春色校园在线视频观看| 欧美xxxx黑人xx丫x性爽| 日韩亚洲欧美综合| 午夜福利视频精品| 亚洲,一卡二卡三卡| 久久久久久久精品精品| 中国国产av一级| 久久久久人妻精品一区果冻| 国产精品欧美亚洲77777| 欧美少妇被猛烈插入视频| 亚洲内射少妇av| 久久青草综合色| 高清不卡的av网站| 国精品久久久久久国模美| 日韩成人av中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 热99国产精品久久久久久7| 国产精品熟女久久久久浪| 一区二区三区乱码不卡18| 国产免费一区二区三区四区乱码| 精品少妇久久久久久888优播| 国产黄片视频在线免费观看| 久久韩国三级中文字幕| 大片免费播放器 马上看| 91精品国产国语对白视频| 最后的刺客免费高清国语| 亚洲成人手机| 欧美成人a在线观看| 毛片一级片免费看久久久久| 国产免费一区二区三区四区乱码| av国产精品久久久久影院| av在线播放精品| 一区二区三区乱码不卡18| 亚洲欧美清纯卡通| 嘟嘟电影网在线观看| 亚洲国产成人一精品久久久| 国产色爽女视频免费观看| 国产成人一区二区在线| 成人高潮视频无遮挡免费网站| av福利片在线观看| 看十八女毛片水多多多| 熟女电影av网| 熟女电影av网| 日本-黄色视频高清免费观看| 少妇丰满av| 精品一区二区免费观看| 日韩av免费高清视频| 丰满乱子伦码专区| 婷婷色综合www| 毛片女人毛片| 亚洲四区av| 亚洲欧美一区二区三区国产| 日本黄色片子视频| 男人爽女人下面视频在线观看| 国产精品久久久久久久久免| 久久精品国产鲁丝片午夜精品| 涩涩av久久男人的天堂| 在线观看一区二区三区激情| 黄色欧美视频在线观看| 精品一区二区三卡| 久久久a久久爽久久v久久| 国产日韩欧美亚洲二区| 国产日韩欧美在线精品| 国产精品人妻久久久久久| 中国美白少妇内射xxxbb| 人人妻人人爽人人添夜夜欢视频 | av免费观看日本| 搡女人真爽免费视频火全软件| 一本—道久久a久久精品蜜桃钙片| 久久精品久久久久久噜噜老黄| 色吧在线观看| 高清不卡的av网站| 夜夜骑夜夜射夜夜干| 熟女人妻精品中文字幕| 国产乱来视频区| 国产一级毛片在线| 亚洲精品成人av观看孕妇| 国产黄片视频在线免费观看| 久久青草综合色| 老司机影院成人| 麻豆成人午夜福利视频| 国产精品国产三级国产专区5o| 欧美成人精品欧美一级黄| 亚洲av中文字字幕乱码综合| 久久亚洲国产成人精品v| 亚洲性久久影院| 九草在线视频观看| 黄色配什么色好看| 91午夜精品亚洲一区二区三区| 日日啪夜夜撸| 午夜日本视频在线| 亚洲成人中文字幕在线播放| 欧美bdsm另类| 日本vs欧美在线观看视频 | 久久久久久久久久久免费av| 不卡视频在线观看欧美| 精品国产一区二区三区久久久樱花 | 婷婷色综合大香蕉| 少妇人妻一区二区三区视频| 少妇人妻久久综合中文| 2018国产大陆天天弄谢| 国产高清国产精品国产三级 | av专区在线播放| 国产av一区二区精品久久 | 国产黄片美女视频| 日韩av不卡免费在线播放| 最黄视频免费看| 久久久久久伊人网av| av国产久精品久网站免费入址| 亚洲人与动物交配视频| 特大巨黑吊av在线直播| 亚洲图色成人| 精品一品国产午夜福利视频| 亚洲无线观看免费| 中文字幕av成人在线电影| 婷婷色综合大香蕉| 高清日韩中文字幕在线| 九色成人免费人妻av| 国产精品爽爽va在线观看网站| 五月天丁香电影| 午夜精品国产一区二区电影| 国产精品麻豆人妻色哟哟久久| 一级毛片我不卡| 妹子高潮喷水视频| 亚洲人成网站在线观看播放| 国产精品一区二区性色av| 一本色道久久久久久精品综合| 少妇猛男粗大的猛烈进出视频| 免费看av在线观看网站| 日日摸夜夜添夜夜添av毛片| 亚洲精品视频女| 精品久久国产蜜桃| av网站免费在线观看视频| 亚洲欧美日韩另类电影网站 | 全区人妻精品视频| 日韩精品有码人妻一区| 大片免费播放器 马上看| 久久久久久久大尺度免费视频| 久久久久久久久久成人| 亚洲欧洲国产日韩| 久久6这里有精品| 如何舔出高潮| 国产伦在线观看视频一区| 亚洲欧美日韩无卡精品| 五月玫瑰六月丁香| 最近最新中文字幕大全电影3| 国产成人免费观看mmmm| av国产精品久久久久影院| 亚洲一级一片aⅴ在线观看| 欧美性感艳星| 97超视频在线观看视频| 一区二区三区精品91| 日本-黄色视频高清免费观看| 成人午夜精彩视频在线观看| 97在线人人人人妻| 久热久热在线精品观看| 三级国产精品欧美在线观看| h视频一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲精品成人av观看孕妇| 最后的刺客免费高清国语| 超碰av人人做人人爽久久| 欧美精品国产亚洲| 欧美日韩精品成人综合77777| 五月开心婷婷网| 日韩,欧美,国产一区二区三区| 欧美97在线视频| 五月天丁香电影| 国产亚洲欧美精品永久| 我要看日韩黄色一级片| 欧美性感艳星| 夫妻午夜视频| 日韩中字成人| 肉色欧美久久久久久久蜜桃| 久久精品熟女亚洲av麻豆精品| av专区在线播放| 亚洲精品中文字幕在线视频 | 久久久久国产网址| 欧美+日韩+精品| 精品午夜福利在线看| 日本午夜av视频| 精品人妻偷拍中文字幕| 欧美三级亚洲精品| 美女脱内裤让男人舔精品视频| 国产精品一区二区三区四区免费观看| 免费黄频网站在线观看国产| 久久精品国产亚洲av天美| 国产亚洲午夜精品一区二区久久| 国产在线免费精品| 国产精品伦人一区二区| 97超视频在线观看视频| 最黄视频免费看| 交换朋友夫妻互换小说| 亚洲欧美清纯卡通| 国产精品.久久久| 国产综合精华液| 午夜福利在线在线| 亚洲图色成人| 日本一二三区视频观看| 免费不卡的大黄色大毛片视频在线观看| 日韩制服骚丝袜av| 在线观看三级黄色| 97热精品久久久久久| 亚洲精品乱码久久久v下载方式| 亚洲图色成人| 日日摸夜夜添夜夜添av毛片| 国产精品.久久久| 极品少妇高潮喷水抽搐| 王馨瑶露胸无遮挡在线观看| 亚洲熟女精品中文字幕| 中文字幕制服av| 少妇人妻久久综合中文| videos熟女内射| 亚洲内射少妇av| 日韩在线高清观看一区二区三区| 久久久久人妻精品一区果冻| 内射极品少妇av片p| 国产黄片视频在线免费观看| 校园人妻丝袜中文字幕| 欧美极品一区二区三区四区| 亚洲第一区二区三区不卡| 黄色配什么色好看| 麻豆成人av视频| 国产色爽女视频免费观看| 久久97久久精品| av国产精品久久久久影院| av国产精品久久久久影院| 亚洲精品成人av观看孕妇| 干丝袜人妻中文字幕| 高清毛片免费看| 日本猛色少妇xxxxx猛交久久| 尾随美女入室| 天美传媒精品一区二区| 人人妻人人澡人人爽人人夜夜| 热99国产精品久久久久久7| 日韩亚洲欧美综合| 大香蕉97超碰在线| 秋霞在线观看毛片| 国产欧美日韩一区二区三区在线 | 色婷婷av一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 国产成人a区在线观看| 免费人成在线观看视频色| 国产淫语在线视频| a 毛片基地| 日韩 亚洲 欧美在线| 亚洲图色成人| 国产黄片美女视频| 高清午夜精品一区二区三区| av福利片在线观看| 男人舔奶头视频| 亚洲精品色激情综合| 91久久精品国产一区二区三区| 欧美成人一区二区免费高清观看| 免费看av在线观看网站| 国产av一区二区精品久久 | 亚洲av成人精品一区久久| 五月天丁香电影| 嘟嘟电影网在线观看| 晚上一个人看的免费电影| 老司机影院成人| 三级国产精品片| 青春草国产在线视频| 日韩免费高清中文字幕av| 一区二区三区乱码不卡18| 久久午夜福利片| 欧美成人一区二区免费高清观看| 大香蕉97超碰在线| 亚洲精品成人av观看孕妇| 久久精品国产a三级三级三级| 最近最新中文字幕大全电影3| 亚洲人成网站在线观看播放| 黄色配什么色好看| 夫妻性生交免费视频一级片| 国产成人精品久久久久久| 婷婷色综合大香蕉| 国产高清有码在线观看视频| 免费黄色在线免费观看| 丝瓜视频免费看黄片| 久久精品夜色国产| 韩国高清视频一区二区三区| 精品99又大又爽又粗少妇毛片| 久热久热在线精品观看| 国产 一区 欧美 日韩| 人妻 亚洲 视频| 日本欧美国产在线视频| 狂野欧美白嫩少妇大欣赏| 老司机影院成人| 日韩亚洲欧美综合| 国产淫片久久久久久久久| 偷拍熟女少妇极品色| 亚洲欧美成人综合另类久久久| 黑丝袜美女国产一区| 亚洲av二区三区四区| 国产精品国产三级国产av玫瑰| 午夜老司机福利剧场| 精品亚洲乱码少妇综合久久| 久久久久国产网址| 午夜日本视频在线| av卡一久久| 精品人妻一区二区三区麻豆| 夜夜看夜夜爽夜夜摸| 五月玫瑰六月丁香| 免费看日本二区| 中文资源天堂在线| 亚洲精品日韩av片在线观看| av卡一久久| 久久毛片免费看一区二区三区| 人体艺术视频欧美日本| 国产 一区精品| 深爱激情五月婷婷| 国产精品久久久久成人av| 中文字幕制服av| 国产一区二区三区综合在线观看 | 国产成人午夜福利电影在线观看| 国产乱人视频| 七月丁香在线播放| 男人和女人高潮做爰伦理| 97精品久久久久久久久久精品| 国产成人a区在线观看| 熟女电影av网| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 欧美97在线视频| 色5月婷婷丁香| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 丝袜喷水一区| 男女免费视频国产| 午夜激情久久久久久久| 最近手机中文字幕大全| 内射极品少妇av片p| 三级经典国产精品| 少妇丰满av| 日本爱情动作片www.在线观看| 日韩亚洲欧美综合| 亚洲国产毛片av蜜桃av| 高清视频免费观看一区二区| 精品国产乱码久久久久久小说| 欧美成人一区二区免费高清观看| 亚洲欧美日韩东京热| 国产精品无大码| 成人漫画全彩无遮挡| 秋霞在线观看毛片| 国产精品国产av在线观看| 国产精品国产三级专区第一集| 各种免费的搞黄视频| 亚洲av免费高清在线观看| 少妇精品久久久久久久| 亚洲精品日韩av片在线观看| 国语对白做爰xxxⅹ性视频网站| 国产深夜福利视频在线观看| 国产有黄有色有爽视频| 97超视频在线观看视频| 欧美成人a在线观看| 国产久久久一区二区三区| 成人亚洲精品一区在线观看 | 国产欧美另类精品又又久久亚洲欧美| 国产深夜福利视频在线观看| 看非洲黑人一级黄片| 免费看av在线观看网站| 欧美最新免费一区二区三区| 大片电影免费在线观看免费| 国产精品蜜桃在线观看| 亚洲精品视频女| 国产av国产精品国产| 丝瓜视频免费看黄片| 波野结衣二区三区在线| 久久99热这里只有精品18| 日本欧美国产在线视频| 联通29元200g的流量卡| 久久久亚洲精品成人影院| 超碰av人人做人人爽久久| 九九爱精品视频在线观看| 网址你懂的国产日韩在线| 一边亲一边摸免费视频| 国产精品一区www在线观看| 午夜视频国产福利| 国产精品偷伦视频观看了| 少妇裸体淫交视频免费看高清| 欧美成人午夜免费资源| 国产乱人视频| 毛片一级片免费看久久久久| 亚洲精品乱码久久久v下载方式| 亚洲精品视频女| 一区二区三区乱码不卡18| 91狼人影院| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 成人亚洲精品一区在线观看 | 亚洲av电影在线观看一区二区三区| 亚洲av不卡在线观看| 精品国产一区二区三区久久久樱花 | 久久精品久久精品一区二区三区| 男人添女人高潮全过程视频| 日本色播在线视频| 一级毛片我不卡| 80岁老熟妇乱子伦牲交| 日本色播在线视频| 大片免费播放器 马上看| 国产淫语在线视频| 身体一侧抽搐| 亚洲,一卡二卡三卡| av视频免费观看在线观看| 亚洲国产精品专区欧美| 人妻少妇偷人精品九色| 秋霞伦理黄片| 亚洲精品视频女| 国产av国产精品国产| 国产成人精品一,二区| 日韩伦理黄色片| 国产亚洲一区二区精品| 身体一侧抽搐| 男人和女人高潮做爰伦理| 国产白丝娇喘喷水9色精品| 欧美国产精品一级二级三级 | 欧美国产精品一级二级三级 | 日产精品乱码卡一卡2卡三| a级毛色黄片| 成年美女黄网站色视频大全免费 | 伦理电影大哥的女人| 嫩草影院新地址| kizo精华| 久久久精品94久久精品| 男女无遮挡免费网站观看| 一级a做视频免费观看| 国产av一区二区精品久久 | 中文字幕亚洲精品专区| 国产伦精品一区二区三区视频9| 国产大屁股一区二区在线视频| 亚洲国产欧美在线一区| h视频一区二区三区| 夫妻午夜视频| 久久亚洲国产成人精品v| 久久久久久久精品精品| 狠狠精品人妻久久久久久综合| 国产黄色免费在线视频| 99久久精品国产国产毛片| av在线蜜桃| 噜噜噜噜噜久久久久久91| 久久久久精品性色| 肉色欧美久久久久久久蜜桃| 亚洲精品国产成人久久av| 乱码一卡2卡4卡精品| av又黄又爽大尺度在线免费看| 久久久久性生活片| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 亚洲精品,欧美精品| 极品教师在线视频| 亚洲怡红院男人天堂| 国产日韩欧美在线精品| 亚洲精品国产成人久久av| 亚洲国产日韩一区二区| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 性色av一级| 日本wwww免费看| 日本欧美视频一区| 激情五月婷婷亚洲| 久久久久国产网址| 高清欧美精品videossex| 国产爱豆传媒在线观看| 色视频www国产| 一区二区av电影网| 国产淫语在线视频| 观看美女的网站| 成人免费观看视频高清| videossex国产| 成人毛片60女人毛片免费| 久久精品国产亚洲网站| 3wmmmm亚洲av在线观看| 老师上课跳d突然被开到最大视频| 日韩av免费高清视频| 日韩国内少妇激情av| av国产久精品久网站免费入址| 高清午夜精品一区二区三区| 久久久久精品性色| 成人黄色视频免费在线看| 熟妇人妻不卡中文字幕| 亚洲自偷自拍三级| 人人妻人人看人人澡| 91久久精品国产一区二区成人| 插逼视频在线观看| 99热全是精品|