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

    帶中式用餐約束的乘務(wù)調(diào)度問題

    2013-07-31 21:35:43陳仕軍沈吟東陳賀命
    關(guān)鍵詞:方法

    陳仕軍,沈吟東,蘇 璇,陳賀命

    (華中科技大學(xué)控制科學(xué)與工程系,武漢430074)

    帶中式用餐約束的乘務(wù)調(diào)度問題

    陳仕軍,沈吟東*,蘇 璇,陳賀命

    (華中科技大學(xué)控制科學(xué)與工程系,武漢430074)

    有效的乘務(wù)調(diào)度能夠為公交企業(yè)帶來巨大的成本節(jié)約,但是,公交乘務(wù)調(diào)度問題因受制于一系列勞動法規(guī)的約束變得十分復(fù)雜.我國公交普遍存在“中式用餐”約束,進(jìn)一步加大了問題的復(fù)雜性,使西方主流調(diào)度系統(tǒng)在國內(nèi)實(shí)施面臨困難.本文基于“生成與選擇”方法解決乘務(wù)調(diào)度問題,關(guān)鍵在于“生成”階段處理“中式用餐”難題;利用“中式用餐”約束和乘務(wù)問題特點(diǎn),設(shè)計一種基于啟發(fā)式規(guī)則的換班機(jī)會篩選方法;在所選換班機(jī)會集合的基礎(chǔ)上構(gòu)造能滿足“中式用餐”約束的潛在乘務(wù)班次集合.對實(shí)際公交乘務(wù)調(diào)度問題中的12組實(shí)例進(jìn)行測試,表明本文方法不僅能處理 “中式用餐”約束,而且能極大減少所求問題的規(guī)模,因此適用于解決大規(guī)模的帶有“中式用餐”約束的乘務(wù)調(diào)度問題.

    城市交通;乘務(wù)調(diào)度;中式用餐;啟發(fā)式;班次生成

    1 引 言

    乘務(wù)調(diào)度問題是公共交通領(lǐng)域中研究的重要難題之一,屬于世界公認(rèn)的NP難題[1].編制優(yōu)化的乘務(wù)調(diào)度方案,可以幫助公交企業(yè)節(jié)約成本、提高乘務(wù)資源效益.因此,自上世紀(jì)六十年代起,西方發(fā)達(dá)國家就開始利用計算機(jī)進(jìn)行公共交通乘務(wù)調(diào)度問題的研究,至今已有許多解決方法和技術(shù)問世[2],并取得成功應(yīng)用.其中,基于整數(shù)線性規(guī)劃(ILP)的方法仍占主導(dǎo)地位,其應(yīng)用也最為廣泛[3].

    經(jīng)典的ILP方法包括一個“生成班次”和一個“選擇班次”的過程,因此,又被歸為“生成與選擇”方法[4].該類方法首先要根據(jù)勞動法規(guī)等約束生成一個大的乘務(wù)班次(簡稱“候選班次”)集合;然后基于ILP方法從中選擇一個最優(yōu)的班次子集覆蓋全部車輛工作.盡管該方法已有許多成功應(yīng)用[4],但面對一些特殊問題時有如下困難:一方面,當(dāng)乘務(wù)調(diào)度問題具有復(fù)雜的勞動法規(guī)和約束(如本文的“中式用餐”)時,生成滿足勞動法規(guī)和約束的全部候選班次會很耗時;另一方面,當(dāng)乘務(wù)調(diào)度問題規(guī)模較大(現(xiàn)實(shí)問題基本都屬于大規(guī)模問題)時,全部候選班次的數(shù)量通常很大(可能是天文數(shù)字),使得“選擇”階段超出ILP的求解能力.目前,大多數(shù)研究側(cè)重于解決后一問題,即擴(kuò)大ILP的求解能力,例如典型的方法有列生成法[5]、元啟發(fā)式方法[6]等.然而,對于如何生成潛在班次集合,一般描述比較簡單或假定其已經(jīng)存在.實(shí)際上,潛在班次的構(gòu)造過程依賴于具體的勞動法規(guī)和約束,特別是在求解具有復(fù)雜法規(guī)和約束的乘務(wù)調(diào)度問題時,生成全部潛在班次往往非常費(fèi)時,有時甚至占據(jù)總求解時間的85%以上[7].因此,研究有效構(gòu)造潛在班次集合的方法,使其滿足各種勞動法規(guī)和約束(尤其是復(fù)雜約束),是應(yīng)用ILP方法解決乘務(wù)調(diào)度問題的首要工作和基礎(chǔ).

    我國的乘務(wù)調(diào)度問題除了具有一般乘務(wù)調(diào)度問題的常見勞動法規(guī)和約束外(例如總工時、連續(xù)工時等的限制),還有一些具有普遍性的特色約束.其中,要求有“中國式”的固定時間段用餐(簡稱“中式用餐”)這一特色約束[8],使其比已有的乘務(wù)調(diào)度問題更加復(fù)雜,進(jìn)而,使得在西方國家得到成功應(yīng)用的各種乘務(wù)調(diào)度方法在國內(nèi)實(shí)施面臨求解困難.本文將采取“生成與選擇”的方法解決我國的具有“中式用餐”約束的乘務(wù)調(diào)度問題,主要是在“生成班次”階段來處理復(fù)雜的“中式用餐”約束.先通過分析“中式用餐”的特點(diǎn),結(jié)合相關(guān)的勞動法規(guī)等約束和乘務(wù)問題特征,設(shè)計出基于啟發(fā)式規(guī)則的換班機(jī)會篩選方法;再利用篩選出的換班機(jī)會集合生成所有潛在班次.本文方法不僅能滿足“中式用餐”的約束要求,而且能極大減小所求問題的規(guī)模,有利于解決大規(guī)模的“中式用餐”乘務(wù)調(diào)度問題.

    2 乘務(wù)調(diào)度問題

    公共交通乘務(wù)調(diào)度問題就是合理安排一組乘務(wù)班次來完成一日內(nèi)的全部車輛運(yùn)營任務(wù).單個車輛的運(yùn)營任務(wù)通常是從車場出發(fā),歷經(jīng)若干趟在線路上的運(yùn)營,最終返回車場,稱為一個車輛運(yùn)營任務(wù)塊(Block).在每輛車的運(yùn)營線路上一般會有一個或多個指定的可供乘務(wù)員換班的地點(diǎn),稱為“換班點(diǎn)”;換班點(diǎn)與車輛到達(dá)該地點(diǎn)的時間結(jié)合起來就構(gòu)成一個“換班機(jī)會(RO)”;在任意一個換班機(jī)會,乘務(wù)員都可以(但不一定)換班;在兩個連續(xù)的“換班機(jī)會”之間,乘務(wù)員通常是不允許換班的,這段連續(xù)的工作被稱為“最小值乘段(Piece)”;一個或幾個連續(xù)的“最小值乘段”可以構(gòu)成“連續(xù)值乘段(Spell)”.一個或者多個“連續(xù)值乘段”進(jìn)行合理組合,可構(gòu)成一個乘務(wù)班次,簡稱“班次(Shift)”,也就是乘務(wù)員一天的工作安排,通常從場站簽到開始到場站簽退結(jié)束[4,9].根據(jù)“班次”的特點(diǎn),一般可以將其劃分成整班、大單班、單班3種類型.整班(又稱連續(xù)班),期間有一個短暫的時間供乘務(wù)員用餐和休息;大單班的乘務(wù)員上班時間跨度(從簽到至簽退)比較長,兩個連續(xù)值乘段間有一段較長休息時間;單班只包含一個連續(xù)乘務(wù)段的班次.圖1描述了部分與乘務(wù)調(diào)度相關(guān)的概念.

    乘務(wù)調(diào)度是一個極為復(fù)雜的組合優(yōu)化問題,其目標(biāo)是采用最少數(shù)量和最小運(yùn)營成本的乘務(wù)班次集合去執(zhí)行完所有車輛任務(wù),且所有乘務(wù)班次需要同時滿足運(yùn)營企業(yè)和國家規(guī)定的相關(guān)勞動法規(guī)和約束.對于一般的乘務(wù)調(diào)度問題,其乘務(wù)班次需要滿足如下約束[9-11]:

    (1)乘務(wù)員連續(xù)值乘時長不可超過規(guī)定的最長時長Tmax和最短時長Tmin;

    (2)班次的總工時不可超過相應(yīng)班型所規(guī)定的最長總工時TSmax與最短總工時TSmin;

    (3)班次的有效工時不可超過相應(yīng)班型所規(guī)定的最長有效工時TWmax與最短有效工時TWmin;

    (4)整班型班次的用餐時間不得超出規(guī)定的用餐最長時間TBmax與最短時間TBmin;

    (5)大單班型班次中兩個連續(xù)值乘段間的休息時間不得少于規(guī)定的最短休息時間TBSmin.

    圖1 車輛任務(wù)和乘務(wù)班次Fig.1 Vehicle work and crew shift

    式中 xj為決策變量,cj是班次j的成本;式(1)是目標(biāo)函數(shù);式(2)表示任意最小值乘段至少被一個被選擇班次所覆蓋.

    由模型可以看出,乘務(wù)調(diào)度問題的規(guī)模與n(即候選班次的數(shù)量)的大小直接相關(guān),而生成候選班次的數(shù)量直接由RO數(shù)量決定.假設(shè)有p輛車的運(yùn)營計劃,平均每輛車的運(yùn)營計劃中含有r個RO(含有r-1個最小值乘段).若暫不考慮班次的合法性約束,則所有車輛運(yùn)行計劃可以生成p×r× ( r-1 )/2個連續(xù)值乘段.當(dāng)考慮生成最多包含k個連續(xù)值乘段的班次時,理論上我們可以算出可能的候選班次數(shù)量

    當(dāng)k較大時,候選班次的數(shù)量n以RO的數(shù)量r呈高次多項式級增長,給ILP模型求解帶來極大困難.因此,在不損失最優(yōu)解前提下盡可能減小候選RO數(shù)量r是解決該問題的關(guān)鍵.下文將基于“中式用餐”約束和乘務(wù)問題特點(diǎn)設(shè)計一些啟發(fā)式策略,篩選出一些“好”的RO作為候選換班機(jī)會集合,從而縮小候選乘務(wù)班次集J的規(guī)模.

    3 “中式用餐”約束下的乘務(wù)班次集生成

    3.1 換班機(jī)會的選擇

    由于“中式用餐”約束要求用餐時間在規(guī)定時間段[tstart,tend]內(nèi)進(jìn)行,故規(guī)定[tstart,tend]內(nèi)的任何RO都可能被用作換班,所以應(yīng)保留這些RO;另一方面,因為用餐必須要換班,且在[tstart,tend]內(nèi)進(jìn)行,故用[tstart,tend]外的大部分換班機(jī)會可能不會被使用,因此可以剔除其中一些“不好”或“不合理”的RO,以縮小問題的規(guī)模.本文在用餐規(guī)定的時間段[tstart,tend]外,根據(jù)乘務(wù)調(diào)度問題的本身特點(diǎn)選擇出一些“好”或“合理”的RO,剔除剩余的RO,使得生成的班次既能滿足“中式用餐”約束,且不會損失關(guān)鍵的RO.選擇一些“好”的RO時,主要基于乘務(wù)問題的如下特點(diǎn):乘務(wù)人員在執(zhí)行車輛運(yùn)營任務(wù)時,一般應(yīng)盡可能持續(xù)較長的工作,以免造成多次換班,從而影響值乘效率;特別,車輛出車場后的工作,一般由同一乘務(wù)員執(zhí)行盡可能長的任務(wù)后,才考慮換班;同理,返回車場的最后一段值乘任務(wù)也盡可能長,以提高工作效率.

    首先,對給定時間區(qū)間[t1,t2]內(nèi)的車輛任務(wù),給出一個簡單啟發(fā)式RO選擇方法(SHA),使所生成連續(xù)值乘段盡可能長.記[t1,t2]內(nèi)最早的換班機(jī)會為rc(對應(yīng)時間tc),SHA如下:

    Step 1 若t2-tc≤Tmax,則停止;

    Step 2 從rc開始,向前尋找新的換班機(jī)會rf(對應(yīng)時間tr),使得rc與rf連接構(gòu)成的連續(xù)值乘段盡可能長,且滿足tr-tc≤Tmax;

    Step 3 選擇rf,記當(dāng)前換班機(jī)會(rf)為rc,返回Step1.

    現(xiàn)對各個車輛任務(wù)分別選擇出一些RO,以用作生成潛在班次集J.記車輛從車場出發(fā)時間ts,返回車場的時間t,中餐時間段[t1,t1],晚餐時

    estartend間段[t2start,t2end],具體的篩選RO策略如下.

    (1)選擇車輛從車場出發(fā)和返回車場的RO (必用換班機(jī)會,即班次開始和結(jié)束RO),如圖2所示.圖2中,黑色圓點(diǎn),表示選擇的RO;白色圓點(diǎn)表示沒有被選擇的RO;虛線表示未畫出的車輛運(yùn)營任務(wù)段.

    圖2 選擇車輛運(yùn)營任務(wù)的始末ROFig.2 Selecting the start RO and the end RO

    (2)選擇規(guī)定用中餐時間段[t1,t1]和晚餐

    startend時間段[t2,t2]內(nèi)的所有RO,如圖3所示.圖3

    startend中的黑圓點(diǎn)表示被選擇的RO,白圓點(diǎn)表示未被選擇RO.

    圖3 選擇用餐時間段內(nèi)的所有ROFig.3 Selecting all the ROs in meal-break time

    由于只依賴策略(1)、策略(2)所選擇的RO集合,無法保證所有車輛任務(wù)被分割成合法的連續(xù)值乘段,因此還需根據(jù)策略(3)選擇一些其他RO.

    (3)選擇用餐時間段[t1,t1]、[t2,t2]以

    startendstartend外的其他RO,確保所有車輛任務(wù)能形成連續(xù)值乘段(根據(jù)參數(shù)Tmax).

    不妨假定車輛運(yùn)營任務(wù)包含中餐和晚餐時間段(只包含1個中餐或晚餐時間段的方法類似),記[t1start,t1

    end]內(nèi)最早和最晚RO的時間分別為t1f, t1

    l;[t2start,t2

    end]內(nèi)最早和最晚RO的時間分別為t2f, t2.分以下三種情況,依次選擇RO

    l

    ① 當(dāng)t1f-ts>Tmax時,則從 t1f向后逐個尋找新RO(對應(yīng)時間tr),直到tr-ts≤Tmax;選擇時間段[tr,t1start]內(nèi)的所有RO,如圖4所示.

    圖4 向后逐個尋找選擇ROFig.4 Selecting the ROs backward

    ② 當(dāng)te-t2l>Tmax時,則從t2l向前逐個尋找新RO(對應(yīng)時間tr),直到te-tr≤Tmax;選擇時間段[t2end,tr]內(nèi)的所有RO,如圖5所示.

    ③ 當(dāng)t2f-t1l>Tmax時,則對[t1l,t2f]內(nèi)的車輛任務(wù)利用SHA選擇RO,如圖6所示.

    圖6 利用SHA選擇ROFig.6 Selecting the ROs by using SHA

    (4)若車輛任務(wù)計劃不包含用餐時間段,則在[ts,te]時間段內(nèi)應(yīng)用SHA選擇RO.

    通過以上RO的選擇過程,形成可能用到的RO集合R.

    3.2 候選班次集J的生成

    基于3.1節(jié)篩選出的RO集合R,根據(jù)如下步驟形成候選班次集J0:

    Step 1 根據(jù)R中的換班機(jī)會,對所有車輛任務(wù)集進(jìn)行切割劃分,以形成最小值乘段集(Piece集);

    Step 2 對多個相鄰的最小值乘段進(jìn)行連接(要求相連的總長度不超過Tmax),形成所有可行連續(xù)值乘段集(Spell集);

    Step 3 對所有連續(xù)值乘段按起始時間排序;

    Step 4 對排序后的連續(xù)值乘段集,通過完全枚舉法進(jìn)行連接,并判斷是否滿足第2節(jié)中的乘務(wù)約束(2)-(5),得到候選班次集J0.

    Step 5 對J0剔除不滿足中式用餐約束的班次,即刪除不滿足約束(6)的班次,得到班次集J.

    4 案例計算

    生成候選班次集J之后,若|J|不大(例如小于2 000),則可以利用著名數(shù)學(xué)優(yōu)化軟件ILOG Cplex求解式(1)、式(2)對應(yīng)的ILP模型[12];但當(dāng)候選班次集規(guī)模|J|較大時,還需要利用列生成法[10]、元啟發(fā)式方法[6]等大規(guī)模優(yōu)化方法求解.主要測試在生成班次滿足“中式用餐”約束條件下,篩選換班機(jī)會方法對生成乘務(wù)班次集規(guī)模|J|的影響.

    通過測試12組實(shí)際數(shù)據(jù)來說明本文方法的有效性,測試數(shù)據(jù)主要來自海口市、武漢市和十堰市的部分公交線路行車方案.相關(guān)乘務(wù)約束所采用的參數(shù)設(shè)置如下:乘務(wù)員的連續(xù)值乘時間不得低于2 h且不超過5 h;乘務(wù)員的實(shí)際值乘時間不得少于6.5 h且不超過8 h,其中單班的值乘時間不低于2 h且不超過5 h;中餐時間段:11:00-13:00,晚餐時間段:18:00-20:00;整班用餐時間不得低于30 min,大單班的休息時間不得低于4 h.

    主要考慮和比較兩類潛在班次集生成方法:

    (1)傳統(tǒng)的考慮所有換班機(jī)會的方法(一般方法);

    (2)本文篩選換班機(jī)會的方法(本文方法).

    計算和測試兩類不同方法下 Piece、Spell及Shift的數(shù)量變化.特別是Shift數(shù)量直接影響到求解乘務(wù)調(diào)度集覆蓋模型的求解難度.相關(guān)的實(shí)驗結(jié)果如表1所示.在表1中,12個測試案例按Block數(shù)量從小到大排列;ro、piece、spell和shift,各自表示該方法下ro、piece、spell和最終生成潛在班次的數(shù)量.RPD表示利用本文方法后潛在班次數(shù)量相對于一般方法班次數(shù)量減少的百分比,計算為

    表1 潛在班次集生成實(shí)驗結(jié)果Table 1 Experimental results for generation of potential shift set

    從表1可知,應(yīng)用本文所提篩選換班點(diǎn)方法后,很大程度地減少了ro、piece、spell及候選班次的數(shù)量.特別是,篩選換班機(jī)會后,候選班次數(shù)量平均減少了63.9%.因此,有利于求解更大規(guī)模的問題.

    5 研究結(jié)論

    本文提出了一種能夠有效處理“中式用餐”約束且能極大縮小問題規(guī)模的乘務(wù)候選班次集的生成方法.“中式用餐”約束與西方國家的用餐約束存在很大差異,因此,使得現(xiàn)有的基于ILP的方法無法直接求解.本文通過在“生成”階段預(yù)先生成滿足“中式用餐”約束的候選班次集合,不僅使得傳統(tǒng)的ILP方法可以用以求解具有中國特色的乘務(wù)調(diào)度問題,而且,還能有效利用“中式用餐”約束,達(dá)到降低候選班次集規(guī)模的目的.總之,本文所提的乘務(wù)班次生成方法不僅有助于解決帶有“中式用餐”的乘務(wù)調(diào)度問題,而且有助于擴(kuò)大基于ILP的乘務(wù)調(diào)度方法的求解規(guī)模.

    雖然本文方法能夠有效降低問題的規(guī)模,但對于大規(guī)模的乘務(wù)調(diào)度問題(如本文所測試的多數(shù)案例),潛在班次數(shù)量仍然很大,直接利用整數(shù)規(guī)劃軟件(如ILOG Cplex)求解ILP模型仍然存在時間上的困難.因此,需要進(jìn)一步研究智能優(yōu)化方法或列生成技術(shù)以解決更大規(guī)模的乘務(wù)調(diào)度問題.

    [1]Wren A,Fores S,Kwan A S K,et al.A flexible system for scheduling drivers[J].Journal of Scheduling, 2003,6(5):437-55.

    [2]Hickman M,Mirchandani P.Computer-aided systems in public transport[C]//Lecture Notes in Economics and Mathematical Systems.Berlin:Springer,2008.

    [3]JütteS, Thonemann U W. Divide-and-price: A decomposition algorithm for solving large railway crew scheduling problems[J]. European Journal of Operational Research,2012,219(2):214-223.

    [4]Shen Y,Kwan R S K.Tabu search for driver scheduling [C].Lecture Notes in Economics and Mathematical Systems,Berlin:Springer,2001(505):121-135.

    [5]Fores S.Column generation approaches to bus driver scheduling[D].University of Leeds,1996.

    [6]Li J,Kwan R S K.A fuzzy genetic algorithm for driver scheduling[J]. European JournalofOperational Research,2003,147(2):334-344.

    [7]Shen Y.Tabu search for bus and train driver scheduling with time windows[D].University of Leeds,2001.

    [8]Shen Y,Xia J.Integrated bus transit scheduling for the Beijing bus group based on a unified mode of operation[J].International Transactions in Operational Research,2009,16(2):227-242.

    [9]沈吟東,倪郁東.含時間窗的司售員調(diào)度問題建模及多鄰域結(jié)構(gòu)設(shè)計[J].華中科技大學(xué)學(xué)報,2008, 36(12):31-34.[SHEN Y D,NI Y D.Model for crew scheduling with time windows and multineighborhood structures[J].Journal of Huazhong University of Science and Technology(Natural Science Edition),2008,36(12):31-34.]

    [10]Mesquita M,Paias A.Set partitioning and coveringbased approaches for the integrated vehicle and crew scheduling problem[J].Computers& Operations Research,2008,35(5):1562-1575.

    [11]Shen Yindong.Two neighbourhood search approaches: 2-opt heuristics and tabu search for bus and train driver scheduling[M].Beijing:Science Press,2003.

    [12]Shen Y,Chen S,Su X.Rail crew scheduling based on a pooling mode for high speed passenger lines[C]// Proceedingsof2010 InternationalConference on Logistics Engineering and Intelligent Transportation Systems,2010.

    A Crew Scheduling with Chinese Meal Break Rules

    CHEN Shi-jun,SHEN Yin-dong,SU Xuan,CHEN He-ming
    (Department of Control Science and Engineering,Huazhong University of Science and Technology,Wuhan 430074,China)

    An efficient crew scheduling can significantly reduce the operational cost for transit enterprises. However,the crew scheduling problem,known to be NP-hard,is complicated by the fact that there are many restrictions on the shift generation.Moreover,there are some special requirements in China,for example, meal break is normally required to be taken during the conventional time ranges for lunch or dinner,which is called a Chinese meal break rule to distinguish from the western ones.It also makes the existing crew scheduling approaches encountering difficulties.On the basis of the“generate and select”approach to solve the crew scheduling problem,this paper proposes an approach to handling the Chinese meal break rule in the phase of“generate”.Taking advantages of the characteristics of Chinese meal break rule and problem domain knowledge,a heuristic-based approach is proposed to select some promising relief opportunities (ROs).A shift generation approach is then devised to generate a large set of potential shifts that satisfy the Chinese meal break rule.Experimental results from 12 groups of real-world problem instances demonstrate the success of the proposed approach,which can greatly reduce the number of potential shifts generated. Therefore,it is suggested that the proposed approach be used to solve the large scale crew scheduling problems with Chinese meal break rule.

    U268.6

    A

    U268.6

    A

    1009-6744(2013)02-0090-06

    2012-10-11

    2012-11-20錄用日期:2012-12-07

    國家自然科學(xué)基金(70971044,71171087).

    陳仕軍(1980-),男,湖北襄陽人,博士生.

    *通訊作者:yindong@hust.edu.cn.

    Key words:urban transportation;crew scheduling;Chinese meal break;heuristics;shift generation

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數(shù)學(xué)教學(xué)改革的方法
    河北畫報(2021年2期)2021-05-25 02:07:46
    化學(xué)反應(yīng)多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學(xué)習(xí)方法
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    亚洲最大成人中文| 天天躁夜夜躁狠狠久久av| 日本 av在线| 久久久久久大精品| 熟女人妻精品中文字幕| 国产成人a区在线观看| 国产黄片美女视频| 亚洲美女视频黄频| 91在线精品国自产拍蜜月| 国产黄色视频一区二区在线观看 | 少妇熟女欧美另类| 中文字幕久久专区| 午夜日韩欧美国产| 成年版毛片免费区| 男人舔奶头视频| 免费看a级黄色片| 精品久久久久久久久av| 简卡轻食公司| 亚洲av中文字字幕乱码综合| 99久久久亚洲精品蜜臀av| 亚洲av.av天堂| а√天堂www在线а√下载| 午夜激情欧美在线| 波多野结衣巨乳人妻| 色哟哟哟哟哟哟| 99热网站在线观看| 亚洲欧美清纯卡通| 男女那种视频在线观看| 国产午夜福利久久久久久| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 秋霞在线观看毛片| 色5月婷婷丁香| 久久热精品热| 亚洲无线在线观看| 国产麻豆成人av免费视频| 成人美女网站在线观看视频| 精品一区二区三区人妻视频| 床上黄色一级片| 给我免费播放毛片高清在线观看| 国产精品一区二区免费欧美| 身体一侧抽搐| 国产一区二区激情短视频| 99久国产av精品国产电影| 精品久久久久久久久久久久久| 黄色一级大片看看| 国产高清激情床上av| 日本a在线网址| 麻豆精品久久久久久蜜桃| 麻豆一二三区av精品| 99久久久亚洲精品蜜臀av| 天堂动漫精品| 国产欧美日韩一区二区精品| 麻豆一二三区av精品| 91久久精品电影网| 欧美一区二区国产精品久久精品| 可以在线观看的亚洲视频| 青春草视频在线免费观看| 你懂的网址亚洲精品在线观看 | 看十八女毛片水多多多| 久久精品国产亚洲av涩爱 | 少妇人妻一区二区三区视频| 91在线精品国自产拍蜜月| 国产精品伦人一区二区| 国产精品乱码一区二三区的特点| 草草在线视频免费看| 女人十人毛片免费观看3o分钟| 国产单亲对白刺激| 国产白丝娇喘喷水9色精品| 精华霜和精华液先用哪个| 国产一区二区三区在线臀色熟女| 免费看av在线观看网站| av在线亚洲专区| 最好的美女福利视频网| 亚洲av熟女| 99久国产av精品| 精品少妇黑人巨大在线播放 | 国产爱豆传媒在线观看| 我要看日韩黄色一级片| 99热这里只有是精品在线观看| 成人永久免费在线观看视频| 一a级毛片在线观看| 两个人视频免费观看高清| 晚上一个人看的免费电影| 两性午夜刺激爽爽歪歪视频在线观看| 男女下面进入的视频免费午夜| 欧美zozozo另类| 97碰自拍视频| 中国国产av一级| 最新在线观看一区二区三区| 亚洲最大成人中文| 日韩欧美三级三区| 国产精品美女特级片免费视频播放器| 国产午夜精品论理片| 亚洲美女黄片视频| 久久天躁狠狠躁夜夜2o2o| 级片在线观看| 人妻夜夜爽99麻豆av| 有码 亚洲区| 能在线免费观看的黄片| 看片在线看免费视频| ponron亚洲| 欧美一区二区国产精品久久精品| 国产精品久久久久久久久免| 亚洲人成网站在线播放欧美日韩| 变态另类丝袜制服| 久久99热6这里只有精品| 欧美激情国产日韩精品一区| 亚洲精品成人久久久久久| 日本爱情动作片www.在线观看 | 亚洲av美国av| 在线免费观看不下载黄p国产| 亚洲国产日韩欧美精品在线观看| 久久精品国产自在天天线| 亚洲aⅴ乱码一区二区在线播放| 中文字幕免费在线视频6| 白带黄色成豆腐渣| 香蕉av资源在线| 亚洲精品国产av成人精品 | 日本免费a在线| 免费看日本二区| 国产毛片a区久久久久| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利在线观看吧| 老师上课跳d突然被开到最大视频| 日本黄色视频三级网站网址| 成人特级av手机在线观看| 少妇熟女欧美另类| .国产精品久久| 日韩欧美三级三区| 久久99热6这里只有精品| 亚洲最大成人中文| 亚洲国产欧洲综合997久久,| 久久99热6这里只有精品| 日韩国内少妇激情av| 一a级毛片在线观看| 亚洲精品久久国产高清桃花| 看片在线看免费视频| 国产精品福利在线免费观看| 午夜福利在线观看吧| 午夜精品一区二区三区免费看| 国产爱豆传媒在线观看| 亚洲国产精品成人综合色| 亚洲美女搞黄在线观看 | 最近的中文字幕免费完整| 亚洲成av人片在线播放无| 欧美日韩在线观看h| 91久久精品国产一区二区成人| 亚洲不卡免费看| av.在线天堂| 床上黄色一级片| 亚洲成a人片在线一区二区| a级一级毛片免费在线观看| 国产一级毛片七仙女欲春2| 亚洲精品日韩在线中文字幕 | 有码 亚洲区| 精品一区二区三区视频在线| 亚洲国产精品国产精品| 亚洲经典国产精华液单| 亚洲av美国av| 性欧美人与动物交配| 99久久九九国产精品国产免费| av免费在线看不卡| 成人av在线播放网站| 午夜免费男女啪啪视频观看 | 日本一本二区三区精品| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| 成人欧美大片| 无遮挡黄片免费观看| 成人无遮挡网站| 人妻丰满熟妇av一区二区三区| 人人妻,人人澡人人爽秒播| 最近中文字幕高清免费大全6| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 成年女人毛片免费观看观看9| 亚洲av二区三区四区| 亚洲熟妇中文字幕五十中出| 中文字幕免费在线视频6| 日本-黄色视频高清免费观看| 亚洲精品国产成人久久av| 成人亚洲精品av一区二区| 亚洲国产高清在线一区二区三| 国产一区二区激情短视频| 国产v大片淫在线免费观看| 国产精品久久久久久久电影| 久久综合国产亚洲精品| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 少妇人妻精品综合一区二区 | 亚洲中文字幕一区二区三区有码在线看| 国产精品爽爽va在线观看网站| 极品教师在线视频| 悠悠久久av| 校园人妻丝袜中文字幕| 99热精品在线国产| 国产蜜桃级精品一区二区三区| 欧美色视频一区免费| 欧美+日韩+精品| 免费av毛片视频| 一进一出抽搐动态| 九九在线视频观看精品| 免费观看人在逋| 久久久a久久爽久久v久久| 日本五十路高清| 亚洲五月天丁香| 国产精品野战在线观看| 日日撸夜夜添| 免费看av在线观看网站| 国产爱豆传媒在线观看| 国产成人影院久久av| 午夜日韩欧美国产| 国产精品,欧美在线| 久久精品国产亚洲网站| 寂寞人妻少妇视频99o| 俄罗斯特黄特色一大片| 成人综合一区亚洲| 最新中文字幕久久久久| 最近的中文字幕免费完整| h日本视频在线播放| 男女边吃奶边做爰视频| 国产成人a区在线观看| 国产一区二区三区av在线 | 国内精品久久久久精免费| 欧美日韩在线观看h| 精品无人区乱码1区二区| 免费看光身美女| 十八禁网站免费在线| 啦啦啦啦在线视频资源| 12—13女人毛片做爰片一| 亚洲av熟女| av专区在线播放| 国产不卡一卡二| 精品一区二区三区视频在线| 女同久久另类99精品国产91| 嫩草影院新地址| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 99在线视频只有这里精品首页| 久久精品国产亚洲av涩爱 | 能在线免费观看的黄片| 国产老妇女一区| 综合色av麻豆| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 成熟少妇高潮喷水视频| 亚洲精品亚洲一区二区| 男人和女人高潮做爰伦理| 成人高潮视频无遮挡免费网站| 变态另类丝袜制服| 国产色爽女视频免费观看| 亚洲成人精品中文字幕电影| 国产成人aa在线观看| 国产视频一区二区在线看| 黄色视频,在线免费观看| 精品久久久久久久久久久久久| 最好的美女福利视频网| av女优亚洲男人天堂| 色综合色国产| 日韩强制内射视频| 免费人成在线观看视频色| 最新中文字幕久久久久| 又爽又黄无遮挡网站| 色哟哟哟哟哟哟| 免费搜索国产男女视频| 国产伦在线观看视频一区| 国产一区亚洲一区在线观看| 国产久久久一区二区三区| 亚洲精品一区av在线观看| 午夜精品在线福利| 成人亚洲欧美一区二区av| 精品久久久久久久久av| 少妇熟女欧美另类| 91狼人影院| 亚洲一级一片aⅴ在线观看| 国产69精品久久久久777片| 国产成人精品久久久久久| 日本a在线网址| a级一级毛片免费在线观看| 91麻豆精品激情在线观看国产| 久久精品国产99精品国产亚洲性色| 国产高潮美女av| 毛片一级片免费看久久久久| 天天一区二区日本电影三级| 国产精品女同一区二区软件| 国产精品人妻久久久久久| 久久久久久伊人网av| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站| 97热精品久久久久久| 特级一级黄色大片| 激情 狠狠 欧美| 国产高清三级在线| 久久99热这里只有精品18| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 一区二区三区免费毛片| 成年版毛片免费区| 成人午夜高清在线视频| 成人永久免费在线观看视频| 九九热线精品视视频播放| 十八禁国产超污无遮挡网站| 国产蜜桃级精品一区二区三区| 一级av片app| 国产三级中文精品| 午夜爱爱视频在线播放| 寂寞人妻少妇视频99o| 久久天躁狠狠躁夜夜2o2o| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 精品久久久噜噜| 国内揄拍国产精品人妻在线| 亚洲av电影不卡..在线观看| 美女免费视频网站| 亚洲一区高清亚洲精品| 欧美日本亚洲视频在线播放| 中文字幕精品亚洲无线码一区| 久久精品91蜜桃| 欧美成人a在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲精品国产成人久久av| 亚洲人成网站高清观看| 色吧在线观看| 日韩欧美国产在线观看| 香蕉av资源在线| 精品久久久久久成人av| 精华霜和精华液先用哪个| 99在线视频只有这里精品首页| av天堂中文字幕网| 男女视频在线观看网站免费| 特级一级黄色大片| 深爱激情五月婷婷| 国产精品亚洲美女久久久| a级毛片免费高清观看在线播放| 国产精品一区二区免费欧美| 亚洲av免费在线观看| 欧美高清性xxxxhd video| 日日啪夜夜撸| 国产精品美女特级片免费视频播放器| 日本黄色视频三级网站网址| 禁无遮挡网站| 成人亚洲欧美一区二区av| 99国产精品一区二区蜜桃av| 亚洲最大成人中文| 久久6这里有精品| 精品少妇黑人巨大在线播放 | 高清日韩中文字幕在线| 在线播放国产精品三级| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 色在线成人网| 老司机午夜福利在线观看视频| 最好的美女福利视频网| 干丝袜人妻中文字幕| 亚洲真实伦在线观看| 熟女人妻精品中文字幕| 亚洲激情五月婷婷啪啪| 日本撒尿小便嘘嘘汇集6| 亚洲精品亚洲一区二区| 一级黄片播放器| 亚洲激情五月婷婷啪啪| 三级毛片av免费| 亚洲成人久久性| 美女cb高潮喷水在线观看| 免费av不卡在线播放| 男女视频在线观看网站免费| 亚洲专区国产一区二区| 插阴视频在线观看视频| 在线天堂最新版资源| videossex国产| .国产精品久久| 日韩成人av中文字幕在线观看 | 最后的刺客免费高清国语| 一个人看的www免费观看视频| 亚洲av熟女| 国内精品一区二区在线观看| 六月丁香七月| 国产熟女欧美一区二区| 亚洲欧美清纯卡通| 麻豆久久精品国产亚洲av| 午夜福利视频1000在线观看| 国内精品美女久久久久久| 午夜精品在线福利| 午夜福利高清视频| 美女 人体艺术 gogo| 在线观看av片永久免费下载| 最近的中文字幕免费完整| 亚洲精品亚洲一区二区| 97超级碰碰碰精品色视频在线观看| 久久欧美精品欧美久久欧美| 天天躁日日操中文字幕| 成人精品一区二区免费| 国内少妇人妻偷人精品xxx网站| 亚洲成人av在线免费| 国产亚洲91精品色在线| 麻豆成人午夜福利视频| 国产免费男女视频| 免费av毛片视频| 午夜精品在线福利| 欧美人与善性xxx| 日韩欧美一区二区三区在线观看| 人人妻人人澡欧美一区二区| 亚洲国产高清在线一区二区三| 免费观看人在逋| 亚洲三级黄色毛片| 欧美日本亚洲视频在线播放| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 国产精品国产三级国产av玫瑰| 91精品国产九色| 99久久久亚洲精品蜜臀av| 一区二区三区免费毛片| 午夜免费男女啪啪视频观看 | 色播亚洲综合网| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 国产成人a∨麻豆精品| 婷婷精品国产亚洲av| 在线观看午夜福利视频| 久久精品91蜜桃| 国产片特级美女逼逼视频| 国产爱豆传媒在线观看| 色视频www国产| 五月伊人婷婷丁香| 亚洲第一电影网av| 天堂动漫精品| 看黄色毛片网站| av免费在线看不卡| 午夜福利成人在线免费观看| 国产高清不卡午夜福利| 日本在线视频免费播放| 深夜a级毛片| 日本撒尿小便嘘嘘汇集6| 在线观看av片永久免费下载| 性色avwww在线观看| 一个人看视频在线观看www免费| АⅤ资源中文在线天堂| 亚洲国产精品久久男人天堂| 又粗又爽又猛毛片免费看| 亚洲五月天丁香| 欧美激情国产日韩精品一区| 九色成人免费人妻av| 国产精品久久电影中文字幕| 人人妻,人人澡人人爽秒播| 好男人在线观看高清免费视频| 黄色日韩在线| 亚洲人成网站在线观看播放| 亚洲成人中文字幕在线播放| 啦啦啦啦在线视频资源| av国产免费在线观看| 亚洲人成网站在线播| 一进一出好大好爽视频| 老女人水多毛片| 精品乱码久久久久久99久播| 真人做人爱边吃奶动态| 精品免费久久久久久久清纯| 国产黄片美女视频| 在现免费观看毛片| 男人舔奶头视频| 老司机午夜福利在线观看视频| 久久这里只有精品中国| 尤物成人国产欧美一区二区三区| 哪里可以看免费的av片| 日韩在线高清观看一区二区三区| 校园人妻丝袜中文字幕| 国产老妇女一区| av女优亚洲男人天堂| av中文乱码字幕在线| 亚洲av不卡在线观看| 日韩欧美一区二区三区在线观看| 国产探花在线观看一区二区| 成人漫画全彩无遮挡| 黄色一级大片看看| 国产亚洲精品久久久com| 日韩一本色道免费dvd| 欧美丝袜亚洲另类| 国产一区二区在线观看日韩| 国产高清三级在线| 嫩草影院精品99| 国产毛片a区久久久久| 久久久久九九精品影院| 97碰自拍视频| 国产91av在线免费观看| 亚洲国产欧洲综合997久久,| 中文资源天堂在线| 卡戴珊不雅视频在线播放| 亚洲电影在线观看av| 午夜激情福利司机影院| 精华霜和精华液先用哪个| 午夜久久久久精精品| 国产国拍精品亚洲av在线观看| 久久久久免费精品人妻一区二区| 免费人成视频x8x8入口观看| 日本-黄色视频高清免费观看| 在线免费观看不下载黄p国产| 五月伊人婷婷丁香| 国产成人a区在线观看| 中文字幕熟女人妻在线| h日本视频在线播放| 国产欧美日韩一区二区精品| 久久人人精品亚洲av| 日本在线视频免费播放| 国产一区二区在线av高清观看| 一个人看视频在线观看www免费| 日韩欧美精品免费久久| 97碰自拍视频| 亚洲不卡免费看| 亚洲内射少妇av| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 91午夜精品亚洲一区二区三区| 欧美成人免费av一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲中文字幕一区二区三区有码在线看| 最近的中文字幕免费完整| 3wmmmm亚洲av在线观看| 午夜精品一区二区三区免费看| 18禁在线播放成人免费| 在线观看免费视频日本深夜| 亚洲人与动物交配视频| 一区福利在线观看| 国产欧美日韩精品一区二区| 嫩草影院新地址| 日本-黄色视频高清免费观看| 国产视频内射| 高清日韩中文字幕在线| 在线播放无遮挡| 六月丁香七月| 又粗又爽又猛毛片免费看| 午夜精品国产一区二区电影 | 岛国在线免费视频观看| 亚洲一区二区三区色噜噜| 成人鲁丝片一二三区免费| 久久99热这里只有精品18| 真人做人爱边吃奶动态| 国产精品久久久久久av不卡| 在线观看av片永久免费下载| 色尼玛亚洲综合影院| 狂野欧美激情性xxxx在线观看| 一级黄色大片毛片| 国产精品一二三区在线看| 久久中文看片网| 国产伦在线观看视频一区| 久久国内精品自在自线图片| 97人妻精品一区二区三区麻豆| 国产精品一区二区免费欧美| 搡老妇女老女人老熟妇| 国产高潮美女av| 亚洲久久久久久中文字幕| 久久久精品94久久精品| 免费av不卡在线播放| aaaaa片日本免费| 日日撸夜夜添| 国产免费男女视频| 在线免费十八禁| 男女做爰动态图高潮gif福利片| 精品久久久久久久人妻蜜臀av| 国产精品久久视频播放| 久久久午夜欧美精品| 精品久久久久久久久久久久久| 偷拍熟女少妇极品色| 身体一侧抽搐| 高清毛片免费看| 激情 狠狠 欧美| 免费观看在线日韩| 人妻久久中文字幕网| 久久久久久久午夜电影| 欧美最新免费一区二区三区| 美女高潮的动态| АⅤ资源中文在线天堂| 国产精品1区2区在线观看.| a级毛色黄片| 波多野结衣高清作品| 少妇高潮的动态图| 精品无人区乱码1区二区| 在线免费观看的www视频| 搞女人的毛片| 国产精品一区www在线观看| 精品国产三级普通话版| 久久久久久久久中文| 日本一二三区视频观看| 国产一区二区在线观看日韩| 看黄色毛片网站| 免费观看的影片在线观看| 久久精品国产自在天天线| 在线国产一区二区在线| 给我免费播放毛片高清在线观看| 国产成人91sexporn| 国产伦一二天堂av在线观看| 麻豆成人午夜福利视频| 日韩欧美一区二区三区在线观看| 特级一级黄色大片| 亚洲国产色片| eeuss影院久久| 成人一区二区视频在线观看| 欧美成人一区二区免费高清观看| 亚洲美女黄片视频| 一级毛片aaaaaa免费看小| 亚洲精品影视一区二区三区av| 人人妻人人澡欧美一区二区| 精品国内亚洲2022精品成人| 中国美女看黄片| 一级a爱片免费观看的视频| 美女被艹到高潮喷水动态| 日韩高清综合在线| www日本黄色视频网| 欧美高清成人免费视频www| 国产成人a∨麻豆精品| 日日啪夜夜撸| 亚洲在线自拍视频|