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

    基于自行車(chē)共享系統(tǒng)靜態(tài)再平衡問(wèn)題的分支定界算法

    2018-12-25 13:10:16王天宇夏曉梅
    物流科技 2018年11期
    關(guān)鍵詞:定界下界站點(diǎn)

    王天宇,韓 印,夏曉梅

    (上海理工大學(xué) 管理學(xué)院,上海 200093)

    0 引言

    隨著全球油價(jià)的上漲、交通擁擠和噪音、空氣等污染,世界上許多城市鼓勵(lì)使用公共交通代替私家車(chē)出行[1]。最有趣的解決方案之一是自行車(chē)共享系統(tǒng)(BSS),它是一個(gè)日益流行的系統(tǒng)。一位顧客使用自行車(chē)從一個(gè)車(chē)站到另一個(gè)車(chē)站。只要有一個(gè)空的上鎖泊位,一輛自行車(chē)可以從任何一個(gè)車(chē)站取出,并可在任何時(shí)候返回任何站點(diǎn)。BSS需要的不僅僅是自行車(chē)和站點(diǎn),各種其他設(shè)備和系統(tǒng)(例如:自行車(chē)車(chē)隊(duì),停車(chē)和鎖定機(jī)制,用戶(hù)界面和結(jié)賬退出協(xié)議和站點(diǎn)網(wǎng)絡(luò)),以及需要維護(hù)和管理需求(例如:車(chē)隊(duì)和站點(diǎn)維修、狀態(tài)信息系統(tǒng)和自行車(chē)再分配系統(tǒng))來(lái)維持自行車(chē)和站點(diǎn)功能保持在足夠的服務(wù)水平。

    大多數(shù)有關(guān)自行車(chē)共享系統(tǒng)的研究都集中在其歷史和發(fā)展上[1],促銷(xiāo)政策和安全問(wèn)題[2]。在過(guò)去的幾年中,很多科學(xué)文章正在研究這些系統(tǒng)中的各種問(wèn)題,包括戰(zhàn)略和運(yùn)營(yíng)設(shè)計(jì)。現(xiàn)有的關(guān)于公共自行車(chē)的研究主要集中在宏觀層面分析其發(fā)展所遇到的政策性問(wèn)題[3]、立法監(jiān)管[4]、共享經(jīng)濟(jì)發(fā)展[5]以及共享單車(chē)調(diào)度優(yōu)化[6]等問(wèn)題。總體而言,較為成熟的研究方向可以歸納為兩類(lèi):戰(zhàn)略設(shè)計(jì)和運(yùn)營(yíng)設(shè)計(jì)。

    在戰(zhàn)略設(shè)計(jì)層面,一些工程研究BSS的戰(zhàn)略規(guī)劃,如選址、網(wǎng)絡(luò)設(shè)計(jì)和戰(zhàn)略規(guī)劃目標(biāo)的站點(diǎn)的能力等。國(guó)內(nèi)相關(guān)專(zhuān)家,例如吳瑤、龔翔等,都是選取我國(guó)某一個(gè)城市,將理論研究投入到設(shè)計(jì)案例中,進(jìn)行分析解剖,根據(jù)公共自行車(chē)分擔(dān)率和居民出行需求測(cè)算出公共自行車(chē)的租借需求規(guī)模,再結(jié)合自行車(chē)和停車(chē)樁的周轉(zhuǎn)率預(yù)測(cè)公共自行車(chē)的總需求和停車(chē)樁的總規(guī)模大小,提出為強(qiáng)化不同的分區(qū)內(nèi)公共自行車(chē)的特色功能,各區(qū)站點(diǎn)布局應(yīng)各有所側(cè)重,有所不同,做到了針對(duì)性的研究和設(shè)計(jì)[7-8]。而國(guó)外專(zhuān)家研究則與國(guó)內(nèi)略有不同,例如Vogel等,提出了一個(gè)合適的模型來(lái)測(cè)試動(dòng)態(tài)定位策略對(duì)系統(tǒng)性能的影響,但是模型對(duì)于戰(zhàn)略規(guī)劃來(lái)說(shuō)是有效的,對(duì)于平衡運(yùn)營(yíng)卻沒(méi)有足夠的細(xì)節(jié)來(lái)支持[9]。

    此外,在運(yùn)營(yíng)設(shè)計(jì)層面而言,系統(tǒng)的運(yùn)營(yíng)必須確保至少具有最低服務(wù)質(zhì)量,以滿足用戶(hù)的最大需求并降低運(yùn)營(yíng)成本。國(guó)外專(zhuān)家,Dell'Olio等調(diào)查了站點(diǎn)的選址問(wèn)題,并為自行車(chē)共享系統(tǒng)的設(shè)計(jì)和實(shí)現(xiàn)開(kāi)發(fā)了一種方法。核心是基于對(duì)一個(gè)給定時(shí)期的需求的估算,考慮到站點(diǎn)和租金的價(jià)格[10]。而Kadri,et al[11-13],則是預(yù)測(cè)用戶(hù)未來(lái)的需求、自行車(chē)的可用性。我國(guó)學(xué)者則指出存在公共自行車(chē)系統(tǒng)不兼容、站點(diǎn)分布不合理、站點(diǎn)自行車(chē)保護(hù)措施欠佳等問(wèn)題,認(rèn)為要促進(jìn)公共自行車(chē)的發(fā)展,不僅要解決以上問(wèn)題,還應(yīng)尋求長(zhǎng)期可靠的資金來(lái)源以維持系統(tǒng)的可持續(xù)發(fā)展[14]。

    綜上所述,系統(tǒng)平衡已經(jīng)成為一個(gè)重要的研究方向了,平衡運(yùn)營(yíng)可以在不同的模式下進(jìn)行,有些操作人員使用靜態(tài)平衡,有些則使用動(dòng)態(tài)平衡,也可以通過(guò)考慮這兩種平衡模式來(lái)實(shí)現(xiàn)[15-16]。在現(xiàn)有文獻(xiàn)中,靜態(tài)平衡問(wèn)題(SBP)比動(dòng)態(tài)平衡問(wèn)題有更多的解決辦法。大多數(shù)研究采用混合整數(shù)規(guī)劃(MIP)技術(shù)。然而,對(duì)這些工作的比較是比較困難的,因?yàn)榇蠖鄶?shù)文獻(xiàn)都考慮了不同的問(wèn)題特征[17-19]。

    針對(duì)這類(lèi)問(wèn)題的不足,本文將重點(diǎn)放在單一車(chē)輛路徑問(wèn)題上,目標(biāo)是在不平衡狀態(tài)下最小化站點(diǎn)的總等待時(shí)間,以達(dá)到最大程度使用現(xiàn)有資源的目的。本文結(jié)構(gòu)如下:在第1節(jié)中,本文提出問(wèn)題并對(duì)問(wèn)題進(jìn)行建模。在第2節(jié)中,本文描述了所提出的分支定界方法以及相關(guān)的下界和上界。在第3節(jié)中,本文案例仿真并討論了實(shí)驗(yàn)結(jié)果。最后,以一些評(píng)論和總結(jié)來(lái)結(jié)束全文。

    1 問(wèn)題描述和數(shù)學(xué)建模

    圖1 各區(qū)間示意圖

    本文的問(wèn)題可以正式定義如下。問(wèn)題的輸入是一組站點(diǎn),用一個(gè)指定的倉(cāng)庫(kù),每個(gè)站(i)的當(dāng)前狀態(tài)(Ei),一個(gè)已知的閾值水平和)表示對(duì)于每個(gè)車(chē)站(i)所需的最小和最大的自行車(chē)數(shù)量,最后,在車(chē)站之間的時(shí)間擴(kuò)展矩陣(dij)表示從站點(diǎn)i到站點(diǎn)j的在分配車(chē)輛的位移時(shí)間。平衡車(chē)輛的路線是由一個(gè)給定的序列表示在訪問(wèn)順序上的站點(diǎn)。其目的是將在不平衡狀態(tài)下等待時(shí)間總和的加權(quán)目標(biāo)最小化,直到車(chē)輛到達(dá)。一個(gè)站點(diǎn)(i)的不平衡(w)表示在平衡操作前的(i)車(chē)站的自行車(chē)數(shù)量(E)和之前那個(gè)站的自行車(chē)[R,R]的ii信用等級(jí)之間的差距。因此,目標(biāo)函數(shù)值取決于站點(diǎn)的不平衡和車(chē)輛的到達(dá)站點(diǎn)的時(shí)間(ti)。各區(qū)間示意圖如圖1所示:

    1.1 問(wèn)題描述

    問(wèn)題是找到再平衡車(chē)輛的一個(gè)最優(yōu)調(diào)度,使其加權(quán)時(shí)間Σiwi(ti- t0)的總和最小,這里 (ti- t0)是站點(diǎn)仍留在站點(diǎn)(i)的不平衡(wi)的時(shí)間。符合條件的問(wèn)題解決方案的數(shù)量,即使是小的也相當(dāng)可觀,因?yàn)榻Y(jié)構(gòu)類(lèi)似于出行商的問(wèn)題。因此,可行的Hamiltonian circuits(哈密頓回路) 的數(shù)量是平常的(n-1)!,這里n是站點(diǎn)的數(shù)量)。在多項(xiàng)式時(shí)間里精確地解決這個(gè)問(wèn)題太難了。問(wèn)題的復(fù)雜性變得越來(lái)越重要,尤其是當(dāng)n很大時(shí)。這一問(wèn)題屬于強(qiáng)烈的NP—困難問(wèn)題。

    1.2 模型建立

    在這一節(jié)中,本文提出SBP的數(shù)學(xué)公式。為了解決這個(gè)問(wèn)題,首先引入以下符號(hào)、變量和參數(shù)。

    符號(hào):

    n——站點(diǎn)的總數(shù);

    N——表示網(wǎng)絡(luò)站點(diǎn)的節(jié)點(diǎn)的集合,i=1,2,…,n;

    N0——節(jié)點(diǎn)集合,包括站點(diǎn)和倉(cāng)庫(kù)(0表示),i=0,1,…,n;

    Ei——在重新定位操作開(kāi)始之前,在節(jié)點(diǎn)i的自行車(chē)的當(dāng)前數(shù)量;

    Ci——站點(diǎn)i的容量(i=1,2,…,n);

    di,j——從站點(diǎn)i到j(luò)的出行時(shí)間;

    ti——到達(dá)站點(diǎn)i的時(shí)間,i=1,2,…,n;

    wi——站點(diǎn)i的權(quán)重,表示Ei和Ri的自信水平之間的差;

    決策變量:xi,j——二元變量,如果站點(diǎn)i在j之前被訪問(wèn)就為1,否則為0,i=0,1,…,n;j=1,…,n且i≠j。

    SBP問(wèn)題可以由以下模型計(jì)算:

    目標(biāo)函數(shù)式(1)使加權(quán)次數(shù)的總和最小。約束式(2)迫使從倉(cāng)庫(kù)離開(kāi)的車(chē)輛。約束式(3)表示車(chē)輛從倉(cāng)庫(kù)出發(fā)的時(shí)間。約束式(4) 確保從(i)到(j)站的車(chē)輛的位移所需的最少時(shí)間。約束條件式(5) (或式(6)) 是車(chē)輛的流量守恒方程(即一個(gè)車(chē)站只能一次被一輛車(chē)訪問(wèn))。最后,式(7)是二元約束。

    2 分支定界方法

    由于問(wèn)題的復(fù)雜性和分支定界(B&B)方法的搜索策略,本文開(kāi)發(fā)了這種類(lèi)型的算法來(lái)解決問(wèn)題。這種方法的目的是通過(guò)減少搜索空間找到一個(gè)最優(yōu)解。該方法基于不同的工具(下界、支配規(guī)則、上界)。在此部分,給出了建議分支和定界過(guò)程的特征。

    搜索樹(shù)最初由一個(gè)根節(jié)點(diǎn)組成,在過(guò)程中以動(dòng)態(tài)的方式開(kāi)發(fā)。上界由最佳可行解的值初始化,這是由一個(gè)或多個(gè)元啟發(fā)式算法預(yù)先計(jì)算的。選擇的方法是,通過(guò)使用選擇策略,從對(duì)應(yīng)于未開(kāi)發(fā)的可行子問(wèn)題的活節(jié)點(diǎn)組中選擇一個(gè)節(jié)點(diǎn)探究。選擇基本上是基于節(jié)點(diǎn)下界的值,它代表了車(chē)輛的部分調(diào)度。

    分支方案包括在部分調(diào)度之后分配一個(gè)新任務(wù)。然后研究選定的節(jié)點(diǎn),并構(gòu)造節(jié)點(diǎn)的子節(jié)點(diǎn)(序列中未調(diào)度的站點(diǎn))。這樣,子空間就被分成了其他子空間。然后,為每個(gè)子空間計(jì)算一個(gè)下界。如果這樣一個(gè)下界的值等于或大于上界,則會(huì)消除子問(wèn)題(因此,子問(wèn)題的可行解不可能比最優(yōu)解更好)。使用深度優(yōu)先和最佳優(yōu)先策略來(lái)探索搜索空間,它允許在可能的情況下快速更新上界,因此當(dāng)它們的下界值等于或大于上界時(shí),就會(huì)消除節(jié)點(diǎn)。

    2.1 下 界

    邊界函數(shù)是B&B算法效率的主要組成部分,不能被節(jié)點(diǎn)選擇和分支的應(yīng)用策略所抵消。在最好的情況下,給定子問(wèn)題的邊界函數(shù)的值應(yīng)該等于最優(yōu)解的值,然而,在大多數(shù)情況下,實(shí)現(xiàn)這樣的目標(biāo)是NP—困難的。

    由于目標(biāo)的目的是使函數(shù)最小化,因此有必要開(kāi)發(fā)低邊界的良好質(zhì)量。一般來(lái)說(shuō),計(jì)算下界包括放松一些約束,然后解決一個(gè)更簡(jiǎn)單的問(wèn)題。

    這個(gè)從模型SBP的約束式(5)和式(6)放松的下界結(jié)果和加權(quán)最短處理時(shí)間的應(yīng)用規(guī)則(WSPT和被稱(chēng)為史密斯的規(guī)則) 是最優(yōu)的單機(jī)調(diào)度問(wèn)題 (1| Σ wjCj),目的在于加權(quán)總完成時(shí)間的最小化。這個(gè)規(guī)則包括排序工作,以增加它們的比率pj/wj(pj表示工作j的持續(xù)時(shí)間,wj表示那項(xiàng)工作的相關(guān)權(quán)重,wj≥0,pj>0)。

    在構(gòu)建路徑(序列)時(shí),放松授權(quán)重新使用相同的弧。下界的值是基于到達(dá)節(jié)點(diǎn)j的最小的出行時(shí)間di,j,對(duì)于給定的節(jié)點(diǎn),本文分兩步來(lái)計(jì)算這個(gè)下界:

    (1)對(duì)于已執(zhí)行的調(diào)度,本文將累積次數(shù)的總和加上相應(yīng)的值wi。

    (2)對(duì)于未執(zhí)行的(未調(diào)度的)作業(yè),本文將關(guān)聯(lián)人工調(diào)度問(wèn)題1|ΣwjCj,其中x(x是未執(zhí)行作業(yè)之一)的處理時(shí)間是到達(dá)相關(guān)節(jié)點(diǎn)的最短運(yùn)行時(shí)間,并且權(quán)重等于wi。人工問(wèn)題由史密斯的規(guī)則解決,加權(quán)完成時(shí)間被添加到步驟(1)中發(fā)現(xiàn)的值中。

    2.2 上界

    利用上界的質(zhì)量是提高B&B算法效率的關(guān)鍵因素?;谶@個(gè)原因,本文利用遺傳算法來(lái)定義上界。

    (1)解決方法編碼。使用遺傳算法解決問(wèn)題的編碼過(guò)程通常是最困難的方面。將它們應(yīng)用于特定的問(wèn)題時(shí),通常很難找到解決方案的適當(dāng)表示,從而易于執(zhí)行交叉過(guò)程。SBP解決方案的良好編碼必須確定被訪問(wèn)站點(diǎn)的順序。一個(gè)合適的表示是一個(gè)染色體組成的站點(diǎn)的路線。應(yīng)該按照它們?cè)谌旧w中出現(xiàn)的順序來(lái)訪問(wèn)每個(gè)站點(diǎn)。圖2給出了一個(gè)9站點(diǎn)的網(wǎng)絡(luò)表示的例子。

    (2)交叉。主要的遺傳算子是交叉,它模擬了兩個(gè)個(gè)體(或者同卵)之間的繁殖。 這是通過(guò)從父母的其中一個(gè)基因的某些部分和另一人的其他部分進(jìn)行的,并且在同一個(gè)孩子中進(jìn)行。

    (3)停止條件。在遺傳算法中使用的停止條件是給定的迭代次數(shù)。即使本文沒(méi)有找到局部或全局最優(yōu),也不可能收斂,但在給定的時(shí)間之后,程序就停止了。停止遺傳算法的其他標(biāo)準(zhǔn)也可以使用(例如,如果在固定的迭代次數(shù)下沒(méi)有改進(jìn)最好的記錄解決方案)。

    3 案例仿真與分析

    在本節(jié)中,本文將進(jìn)行數(shù)值實(shí)驗(yàn),以評(píng)估和比較算法。在區(qū)間[0,+10]內(nèi)選取代表各站不平衡的權(quán)重值。對(duì)于n={10,20 },本文定義了8組實(shí)例,具體結(jié)果見(jiàn)表1和表2。對(duì)于n=30,本文簡(jiǎn)單定義了5個(gè)組,具體結(jié)果見(jiàn)表3。因?yàn)锽&B算法的時(shí)間要求更大。本文注意到所有的CPU時(shí)間都是秒。為了評(píng)估B&B的性能,計(jì)算實(shí)驗(yàn)是在相同的實(shí)例上進(jìn)行的,用于評(píng)估下界和上界。計(jì)算時(shí)間限制在7 200s。B&B的性能表現(xiàn)在下面兩張表中,本文使用“BB-NxGy”這個(gè)符號(hào),其中x是n的值,y表示組號(hào)。報(bào)告的參數(shù)如下:

    圖2 9站點(diǎn)網(wǎng)絡(luò)示意圖

    OPT——最優(yōu)值的平均值;

    UB*——最好上界值的平均值;

    CPU——在幾秒內(nèi)分支定界計(jì)算時(shí)間的平均值;

    Nbr——探索節(jié)點(diǎn)的平均數(shù)量;

    Nd——消除節(jié)點(diǎn)的平均數(shù)量;

    Np——在初步消除過(guò)程中消除節(jié)點(diǎn)的平均數(shù)量;

    Gap(%)——最優(yōu)上限與最優(yōu)值之間的平均差值。

    表1 10個(gè)站點(diǎn)實(shí)驗(yàn)結(jié)果匯總

    表2 20個(gè)站點(diǎn)實(shí)驗(yàn)結(jié)果匯總

    由上表可以得到以下幾點(diǎn)結(jié)論:

    B&B能夠解決多達(dá)30個(gè)站點(diǎn)的實(shí)例;

    (1)當(dāng)n更大的時(shí)候,初步排除是更有效的,因?yàn)橥ㄟ^(guò)消除一個(gè)節(jié)點(diǎn)包含x剩余任務(wù)(x<n)本文避免x!可能性和結(jié)果本文減少了B&B的計(jì)算時(shí)間;

    (2)當(dāng)n較大時(shí),所探索的節(jié)點(diǎn)的計(jì)算時(shí)間和平均數(shù)量都較高。這是由于搜索空間的增加引起的。因此,搜索樹(shù)越大,問(wèn)題就越難解決。

    表3 30個(gè)站點(diǎn)實(shí)驗(yàn)結(jié)果匯總

    得到的結(jié)果表明,在合理的時(shí)間內(nèi)可以計(jì)算出(n≤30)的最優(yōu)解。此外,開(kāi)發(fā)的B&B算法可以轉(zhuǎn)化為一種貪婪的搜索算法,它能夠在非常短的計(jì)算時(shí)間內(nèi)提供良好的解決方案。

    4 結(jié)束語(yǔ)

    本文研究了考慮靜態(tài)模式的自行車(chē)共享系統(tǒng)中的單個(gè)車(chē)輛調(diào)度問(wèn)題。本文定義并制定了一個(gè)數(shù)學(xué)模型,它的動(dòng)機(jī)是需要將不平衡狀態(tài)下的車(chē)站等待時(shí)間降到最低。在分支定界的算法中,提出并使用了下界和上界。為了評(píng)估和比較本文的算法,在大量的實(shí)例上進(jìn)行了數(shù)值實(shí)驗(yàn)。結(jié)果表明,B&B算法可以解決多達(dá)30個(gè)節(jié)點(diǎn)的實(shí)例。但是,此方法在短計(jì)算時(shí)間內(nèi)解決更大數(shù)量問(wèn)題仍然有些力不從心。

    在未來(lái)的工作中,本文的目標(biāo)是在其他類(lèi)型的放松的基礎(chǔ)上開(kāi)發(fā)新的下界,以及一個(gè)分支切割的方法來(lái)加速本文的算法。此外,新的元啟發(fā)式算法和表示方法的研究似乎很有趣,以提高上界質(zhì)量,從而提高B&B的性能。

    猜你喜歡
    定界下界站點(diǎn)
    RTK技術(shù)在土地勘測(cè)定界中的應(yīng)用研究
    一類(lèi)DC規(guī)劃問(wèn)題的分支定界算法
    基于Web站點(diǎn)的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    2017~2018年冬季西北地區(qū)某站點(diǎn)流感流行特征分析
    Lower bound estimation of the maximum allowable initial error and its numerical calculation
    基于外定界橢球集員估計(jì)的純方位目標(biāo)跟蹤
    首屆歐洲自行車(chē)共享站點(diǎn)協(xié)商會(huì)召開(kāi)
    怕被人認(rèn)出
    矩陣Hadamard積的上下界序列
    最大度為10的邊染色臨界圖邊數(shù)的新下界
    99久久精品一区二区三区| 中文亚洲av片在线观看爽| 亚洲图色成人| 亚洲成a人片在线一区二区| 久久草成人影院| 一卡2卡三卡四卡精品乱码亚洲| 免费黄网站久久成人精品| 免费黄网站久久成人精品| 麻豆国产av国片精品| 亚洲欧美精品自产自拍| 色5月婷婷丁香| 联通29元200g的流量卡| 级片在线观看| 国产精品久久电影中文字幕| 中文字幕免费在线视频6| 日本一二三区视频观看| 18禁在线无遮挡免费观看视频| 一边摸一边抽搐一进一小说| 一本久久中文字幕| 欧美成人一区二区免费高清观看| 卡戴珊不雅视频在线播放| 91久久精品国产一区二区成人| 一级毛片我不卡| 男人和女人高潮做爰伦理| 午夜激情欧美在线| 欧美成人免费av一区二区三区| 99久久成人亚洲精品观看| 亚洲精华国产精华液的使用体验 | 日日干狠狠操夜夜爽| 久久99热6这里只有精品| 99久久精品一区二区三区| 日韩欧美 国产精品| 亚洲国产欧美在线一区| 99久国产av精品国产电影| 成人亚洲精品av一区二区| 欧美人与善性xxx| 国产亚洲5aaaaa淫片| 男女边吃奶边做爰视频| 亚洲成a人片在线一区二区| 91精品一卡2卡3卡4卡| 又粗又爽又猛毛片免费看| 99在线人妻在线中文字幕| 亚洲四区av| 老师上课跳d突然被开到最大视频| 国产精品一区二区在线观看99 | www.av在线官网国产| 亚洲乱码一区二区免费版| 免费电影在线观看免费观看| 麻豆国产97在线/欧美| 久久久色成人| 国产在线男女| 国产三级在线视频| 偷拍熟女少妇极品色| 三级毛片av免费| 国产探花在线观看一区二区| 又粗又爽又猛毛片免费看| 亚洲无线在线观看| 啦啦啦啦在线视频资源| 国产精品一区二区三区四区免费观看| 国产精品一及| 亚洲国产色片| 在线播放无遮挡| 国产亚洲5aaaaa淫片| 18禁在线播放成人免费| 校园春色视频在线观看| 高清午夜精品一区二区三区 | 午夜视频国产福利| 亚洲乱码一区二区免费版| 在线观看美女被高潮喷水网站| 亚洲成a人片在线一区二区| 日韩在线高清观看一区二区三区| 3wmmmm亚洲av在线观看| 能在线免费看毛片的网站| 午夜福利成人在线免费观看| 国产老妇女一区| www.色视频.com| 一级毛片我不卡| 又黄又爽又刺激的免费视频.| 人体艺术视频欧美日本| 国产一区亚洲一区在线观看| 午夜激情福利司机影院| 久久久久久久久久成人| www.色视频.com| 国产精华一区二区三区| 青春草视频在线免费观看| 精品熟女少妇av免费看| or卡值多少钱| 精品久久久久久成人av| 不卡一级毛片| 国产一区二区三区在线臀色熟女| 亚洲av电影不卡..在线观看| 亚洲国产高清在线一区二区三| 97超视频在线观看视频| 国产精品一区www在线观看| 中文字幕av在线有码专区| 国产极品天堂在线| 91久久精品国产一区二区成人| 亚洲欧美精品综合久久99| 国内精品一区二区在线观看| 91麻豆精品激情在线观看国产| 在线观看午夜福利视频| 国产成人精品一,二区 | 在线播放国产精品三级| 人妻夜夜爽99麻豆av| 又粗又硬又长又爽又黄的视频 | 网址你懂的国产日韩在线| 26uuu在线亚洲综合色| 69av精品久久久久久| 九九爱精品视频在线观看| 白带黄色成豆腐渣| 爱豆传媒免费全集在线观看| 精品欧美国产一区二区三| 免费人成在线观看视频色| 欧美一区二区精品小视频在线| 亚洲精品粉嫩美女一区| 国产69精品久久久久777片| 伦精品一区二区三区| 亚洲精品日韩av片在线观看| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕| 日韩一本色道免费dvd| or卡值多少钱| 中文字幕免费在线视频6| 免费一级毛片在线播放高清视频| 色哟哟哟哟哟哟| 色综合亚洲欧美另类图片| 国国产精品蜜臀av免费| 麻豆av噜噜一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久久久丰满| 日本一本二区三区精品| 免费黄网站久久成人精品| 成人亚洲精品av一区二区| 久久久久久久久久久丰满| av黄色大香蕉| 成人高潮视频无遮挡免费网站| 亚洲精品乱码久久久久久按摩| 国产成人91sexporn| 日本爱情动作片www.在线观看| 精品一区二区三区视频在线| 亚洲美女搞黄在线观看| 日日摸夜夜添夜夜添av毛片| 熟女电影av网| 免费av毛片视频| 免费观看人在逋| 国产精品久久久久久精品电影小说 | 亚洲欧美成人精品一区二区| 欧美bdsm另类| 99久久精品国产国产毛片| 男女那种视频在线观看| 久久热精品热| 国产黄片视频在线免费观看| 欧美极品一区二区三区四区| 午夜福利在线观看免费完整高清在 | 老司机福利观看| 三级男女做爰猛烈吃奶摸视频| 美女黄网站色视频| 亚洲成人久久爱视频| 欧美日本视频| 3wmmmm亚洲av在线观看| av女优亚洲男人天堂| 97超碰精品成人国产| 亚洲精品日韩av片在线观看| 日韩亚洲欧美综合| 免费看光身美女| 国产精品久久久久久精品电影| 亚洲av免费高清在线观看| 日韩欧美在线乱码| 高清日韩中文字幕在线| 99热这里只有是精品在线观看| av在线播放精品| 中文字幕免费在线视频6| 国产成人精品久久久久久| 综合色av麻豆| 亚洲国产精品国产精品| 国产亚洲5aaaaa淫片| 精品久久久久久久末码| 亚洲人成网站在线观看播放| 91麻豆精品激情在线观看国产| 狂野欧美白嫩少妇大欣赏| 看免费成人av毛片| a级毛片免费高清观看在线播放| 麻豆精品久久久久久蜜桃| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 内地一区二区视频在线| av天堂在线播放| 卡戴珊不雅视频在线播放| 欧美最黄视频在线播放免费| 精品国产三级普通话版| 简卡轻食公司| 高清午夜精品一区二区三区 | 三级经典国产精品| 欧美精品国产亚洲| 免费黄网站久久成人精品| 寂寞人妻少妇视频99o| 国产视频首页在线观看| 亚洲在线观看片| 国产av在哪里看| 丝袜美腿在线中文| 亚洲成av人片在线播放无| 国产高清激情床上av| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添av毛片| 少妇的逼好多水| 亚洲av电影不卡..在线观看| 久久99精品国语久久久| 美女cb高潮喷水在线观看| 欧美潮喷喷水| 男女做爰动态图高潮gif福利片| 国产人妻一区二区三区在| 岛国毛片在线播放| 激情 狠狠 欧美| 99久久成人亚洲精品观看| 久久久午夜欧美精品| 亚洲婷婷狠狠爱综合网| 亚洲,欧美,日韩| 免费不卡的大黄色大毛片视频在线观看 | 变态另类丝袜制服| 亚洲第一电影网av| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 国产探花极品一区二区| 国产黄a三级三级三级人| 中文字幕精品亚洲无线码一区| 91精品国产九色| 成人漫画全彩无遮挡| 国产老妇伦熟女老妇高清| 亚洲三级黄色毛片| 三级国产精品欧美在线观看| 亚洲av成人av| 国产日本99.免费观看| 少妇熟女aⅴ在线视频| 91麻豆精品激情在线观看国产| 亚洲欧美清纯卡通| 校园春色视频在线观看| 熟妇人妻久久中文字幕3abv| 亚洲自拍偷在线| 国产精品精品国产色婷婷| 国产单亲对白刺激| 夫妻性生交免费视频一级片| 九九久久精品国产亚洲av麻豆| 日韩高清综合在线| 啦啦啦观看免费观看视频高清| 日韩一本色道免费dvd| av在线播放精品| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看 | 欧美bdsm另类| 男插女下体视频免费在线播放| 啦啦啦观看免费观看视频高清| 日韩人妻高清精品专区| 免费大片18禁| 成年av动漫网址| 亚洲中文字幕日韩| 亚洲18禁久久av| 国产成人aa在线观看| 91精品一卡2卡3卡4卡| 哪个播放器可以免费观看大片| 男女视频在线观看网站免费| 国产在视频线在精品| 国产高潮美女av| 深爱激情五月婷婷| 性欧美人与动物交配| 亚洲国产精品合色在线| 免费av毛片视频| 国产真实乱freesex| 色5月婷婷丁香| www.色视频.com| 日本欧美国产在线视频| 国产高清三级在线| 免费观看在线日韩| 国产高清视频在线观看网站| 美女 人体艺术 gogo| 欧美激情国产日韩精品一区| 一区福利在线观看| 狠狠狠狠99中文字幕| 午夜久久久久精精品| 国产高清不卡午夜福利| 亚洲欧美成人综合另类久久久 | av在线播放精品| 内地一区二区视频在线| av天堂中文字幕网| 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产| 国产伦一二天堂av在线观看| 国产不卡一卡二| 禁无遮挡网站| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 色尼玛亚洲综合影院| 夜夜看夜夜爽夜夜摸| 国产视频首页在线观看| 干丝袜人妻中文字幕| .国产精品久久| 床上黄色一级片| 成人鲁丝片一二三区免费| 91久久精品电影网| 久久久久网色| 国产免费一级a男人的天堂| 性色avwww在线观看| 伦理电影大哥的女人| 亚洲av一区综合| 18+在线观看网站| 国产av不卡久久| 国产在视频线在精品| 男人舔女人下体高潮全视频| 欧美极品一区二区三区四区| 26uuu在线亚洲综合色| 国产亚洲欧美98| 日韩一本色道免费dvd| 亚洲欧美日韩高清专用| 国内精品美女久久久久久| 久久午夜福利片| 日韩av在线大香蕉| 天堂网av新在线| 99久久精品热视频| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 热99re8久久精品国产| 不卡一级毛片| 久久久久网色| 国产精品综合久久久久久久免费| 久久九九热精品免费| 美女脱内裤让男人舔精品视频 | 国产伦在线观看视频一区| 国产一区二区激情短视频| 日本一本二区三区精品| av又黄又爽大尺度在线免费看 | 久久久久久久久中文| 少妇的逼好多水| 国产av一区在线观看免费| 国内揄拍国产精品人妻在线| 精品久久久久久成人av| 婷婷六月久久综合丁香| 又爽又黄a免费视频| 婷婷六月久久综合丁香| 级片在线观看| 爱豆传媒免费全集在线观看| 国内精品久久久久精免费| 美女黄网站色视频| 久久久精品大字幕| 久久人人爽人人爽人人片va| 久久这里有精品视频免费| 亚洲av成人av| 免费一级毛片在线播放高清视频| 一级毛片电影观看 | 久久精品国产亚洲av涩爱 | 精品不卡国产一区二区三区| 一边摸一边抽搐一进一小说| 在线天堂最新版资源| 一级毛片我不卡| 欧美日韩综合久久久久久| 亚洲高清免费不卡视频| 欧美日韩综合久久久久久| 国模一区二区三区四区视频| 亚洲国产高清在线一区二区三| 国产成人影院久久av| 欧美日韩综合久久久久久| 欧美xxxx黑人xx丫x性爽| 一个人看的www免费观看视频| 69av精品久久久久久| 国产精品一区二区三区四区免费观看| 性插视频无遮挡在线免费观看| 国产精品.久久久| 一区二区三区四区激情视频 | 午夜福利视频1000在线观看| 亚洲中文字幕一区二区三区有码在线看| 99久久精品国产国产毛片| 日本黄色视频三级网站网址| 日韩亚洲欧美综合| 中文字幕av成人在线电影| 国产精品一区www在线观看| 欧美精品一区二区大全| 最近的中文字幕免费完整| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 哪个播放器可以免费观看大片| 亚洲欧美精品综合久久99| eeuss影院久久| 如何舔出高潮| 久久精品国产亚洲av香蕉五月| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄 | 午夜免费男女啪啪视频观看| 久久综合国产亚洲精品| 黄色配什么色好看| 午夜福利在线观看吧| 国产熟女欧美一区二区| 久久久国产成人免费| 两个人的视频大全免费| 国产精品久久久久久精品电影小说 | 国产成年人精品一区二区| 国内精品宾馆在线| 久久精品国产亚洲av涩爱 | 成人综合一区亚洲| 精品久久久久久久久亚洲| 成人性生交大片免费视频hd| 国产精品.久久久| 我要看日韩黄色一级片| 又黄又爽又刺激的免费视频.| av国产免费在线观看| 亚洲成人久久爱视频| 国产精品一区二区在线观看99 | 午夜免费男女啪啪视频观看| 美女黄网站色视频| 精品日产1卡2卡| 少妇被粗大猛烈的视频| 午夜久久久久精精品| 99热网站在线观看| 中文字幕制服av| 99国产极品粉嫩在线观看| 久久精品久久久久久久性| 国产三级中文精品| 欧美色视频一区免费| 久久精品91蜜桃| 久久久国产成人免费| 男人舔奶头视频| 久久精品久久久久久久性| 精品久久久久久久久亚洲| 日本在线视频免费播放| 亚州av有码| a级一级毛片免费在线观看| 我要看日韩黄色一级片| 黄色一级大片看看| 男插女下体视频免费在线播放| 国产免费一级a男人的天堂| 美女 人体艺术 gogo| 欧美日韩在线观看h| 国产91av在线免费观看| 黑人高潮一二区| 国产伦精品一区二区三区四那| 听说在线观看完整版免费高清| 不卡视频在线观看欧美| 最后的刺客免费高清国语| av在线蜜桃| 秋霞在线观看毛片| 婷婷色综合大香蕉| 国产大屁股一区二区在线视频| 哪里可以看免费的av片| 少妇人妻一区二区三区视频| 精品日产1卡2卡| 国产精品久久久久久精品电影| 精品久久久噜噜| a级毛片a级免费在线| 国产伦在线观看视频一区| 色播亚洲综合网| 在线天堂最新版资源| 乱系列少妇在线播放| 村上凉子中文字幕在线| 日本在线视频免费播放| 99在线人妻在线中文字幕| 久久亚洲国产成人精品v| 久久6这里有精品| 高清午夜精品一区二区三区 | 婷婷色综合大香蕉| 国产色爽女视频免费观看| 九色成人免费人妻av| 亚洲无线观看免费| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区性色av| 三级毛片av免费| 97人妻精品一区二区三区麻豆| 国内久久婷婷六月综合欲色啪| 中文精品一卡2卡3卡4更新| 亚洲成av人片在线播放无| 女人十人毛片免费观看3o分钟| 国产老妇女一区| 精品一区二区三区视频在线| 亚洲自拍偷在线| 久久精品国产鲁丝片午夜精品| 久久人人爽人人爽人人片va| 熟妇人妻久久中文字幕3abv| 插阴视频在线观看视频| 神马国产精品三级电影在线观看| 人妻久久中文字幕网| 亚洲一级一片aⅴ在线观看| 可以在线观看毛片的网站| 午夜亚洲福利在线播放| 熟女电影av网| 国产日本99.免费观看| 久久精品夜夜夜夜夜久久蜜豆| 男女那种视频在线观看| 美女xxoo啪啪120秒动态图| 一个人看的www免费观看视频| 好男人在线观看高清免费视频| 久久久午夜欧美精品| 国产精品久久电影中文字幕| 床上黄色一级片| 日韩国内少妇激情av| 99riav亚洲国产免费| 精品久久久久久久久久久久久| av又黄又爽大尺度在线免费看 | 欧美日韩综合久久久久久| 精品人妻熟女av久视频| 26uuu在线亚洲综合色| 天堂网av新在线| 一级二级三级毛片免费看| 三级男女做爰猛烈吃奶摸视频| 久99久视频精品免费| 婷婷色综合大香蕉| 天堂影院成人在线观看| a级毛片a级免费在线| 国产黄片美女视频| 在线免费十八禁| 亚洲国产精品sss在线观看| 最近最新中文字幕大全电影3| 九九热线精品视视频播放| 成年版毛片免费区| 可以在线观看的亚洲视频| 97人妻精品一区二区三区麻豆| av在线老鸭窝| 青青草视频在线视频观看| 精品不卡国产一区二区三区| 人妻系列 视频| 美女国产视频在线观看| 国产精品人妻久久久久久| 亚洲高清免费不卡视频| 亚洲精品亚洲一区二区| 色尼玛亚洲综合影院| 一级毛片电影观看 | 一进一出抽搐gif免费好疼| 国产v大片淫在线免费观看| av又黄又爽大尺度在线免费看 | 久久精品91蜜桃| 成年免费大片在线观看| 久久久久久久久久黄片| 91精品一卡2卡3卡4卡| 日韩欧美精品免费久久| 免费看a级黄色片| 国产亚洲av片在线观看秒播厂 | 精华霜和精华液先用哪个| 亚洲在线观看片| 嫩草影院新地址| 日本黄色视频三级网站网址| 久久久国产成人精品二区| 久久中文看片网| 亚洲性久久影院| 最近手机中文字幕大全| 日韩欧美精品v在线| 国产 一区 欧美 日韩| 国产精品女同一区二区软件| 国产午夜福利久久久久久| 色视频www国产| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱 | 久久国产乱子免费精品| 啦啦啦观看免费观看视频高清| 日日干狠狠操夜夜爽| 国产高潮美女av| 99热这里只有精品一区| 91久久精品国产一区二区三区| h日本视频在线播放| 国产中年淑女户外野战色| 色综合亚洲欧美另类图片| 欧美日本视频| 久久国内精品自在自线图片| 18禁黄网站禁片免费观看直播| 啦啦啦韩国在线观看视频| 卡戴珊不雅视频在线播放| 亚洲精品亚洲一区二区| 一本精品99久久精品77| 天美传媒精品一区二区| 悠悠久久av| 免费黄网站久久成人精品| 91久久精品国产一区二区三区| 精品不卡国产一区二区三区| 欧美丝袜亚洲另类| 深夜精品福利| 黄色一级大片看看| videossex国产| 悠悠久久av| a级毛色黄片| av免费观看日本| 99久久人妻综合| 国产精品国产三级国产av玫瑰| 午夜精品一区二区三区免费看| 久久精品影院6| kizo精华| 嫩草影院新地址| 人体艺术视频欧美日本| 精品一区二区三区视频在线| 久久久午夜欧美精品| 欧美成人精品欧美一级黄| 日韩欧美一区二区三区在线观看| 久久久久久久亚洲中文字幕| 国产精品爽爽va在线观看网站| 国产成人freesex在线| 97热精品久久久久久| 我的老师免费观看完整版| 韩国av在线不卡| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美精品专区久久| 午夜久久久久精精品| 国产一区二区三区在线臀色熟女| 真实男女啪啪啪动态图| 亚洲内射少妇av| 99热精品在线国产| 狂野欧美激情性xxxx在线观看| 成年版毛片免费区| 男人和女人高潮做爰伦理| 日韩精品青青久久久久久| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 成人欧美大片| 91久久精品国产一区二区成人| 亚洲一区二区三区色噜噜| 国产美女午夜福利| 26uuu在线亚洲综合色| 午夜精品一区二区三区免费看| 九九热线精品视视频播放| 亚洲欧美成人精品一区二区| av在线播放精品|