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

    分層PageRank算法和逆序系數(shù)在鐵路網(wǎng)絡(luò)中的應(yīng)用

    2021-12-16 08:21:14劉芝秀黃小杰
    關(guān)鍵詞:逆序結(jié)點(diǎn)排序

    劉芝秀,熊 可,黃小杰*

    (1.南昌工程學(xué)院理學(xué)院,江西 南昌 330099;2.江西財(cái)經(jīng)大學(xué)統(tǒng)計(jì)學(xué)院,江西 南昌 330077)

    鐵路交通是最為重要的公共基礎(chǔ)設(shè)施之一,是社會(huì)經(jīng)濟(jì)發(fā)展的重要支撐。隨著我國鐵路和高鐵建設(shè)的不斷推進(jìn),其建設(shè)規(guī)模與技術(shù)都取得了巨大的成就,國內(nèi)各個(gè)城市和區(qū)域之間的交通往來更為便利。通常經(jīng)濟(jì)發(fā)達(dá)的地方其城際鐵路交通也較發(fā)達(dá),顯然各城市列車站點(diǎn)的重要性是一個(gè)值得探討的問題,至少有助于分析鐵路網(wǎng)絡(luò)的安全性與效能。關(guān)于城市網(wǎng)絡(luò)結(jié)點(diǎn)重要性的計(jì)算,文獻(xiàn)[1]提出了一種適應(yīng)PageRank算法,此算法考慮了來源數(shù)據(jù)本身的重要程度,使城市網(wǎng)絡(luò)中來自相鄰結(jié)點(diǎn)的數(shù)據(jù)比來自較遠(yuǎn)結(jié)點(diǎn)的數(shù)據(jù)更有可能被考慮,并據(jù)此計(jì)算每個(gè)城市結(jié)點(diǎn)的重要性。

    其中PageRank方法最初是用于評價(jià)互聯(lián)網(wǎng)上網(wǎng)站的知名度和重要性,它是由Page及其同事于1998年提出的[2]。經(jīng)過不斷的發(fā)展和改進(jìn),PageRank方法已被廣泛的應(yīng)用于社會(huì)生活的各個(gè)方面,例如近期的文獻(xiàn)[3]將其應(yīng)用于挖掘各學(xué)科的優(yōu)秀和潛在人才,文獻(xiàn)[4]將其用于對負(fù)責(zé)生物狀態(tài)動(dòng)態(tài)變化的轉(zhuǎn)錄因子進(jìn)行排序,文獻(xiàn)[5]采用PageRank方法識(shí)別了光伏并網(wǎng)系統(tǒng)中的關(guān)鍵線路,文獻(xiàn)[6]將PageRank方法應(yīng)用于分析陜西省COVID-19患者動(dòng)態(tài)接觸網(wǎng)絡(luò),文獻(xiàn)[7]也基于PageRank算法對號(hào)稱全球貿(mào)易先行“晴雨表”的中間品貿(mào)易的關(guān)鍵節(jié)點(diǎn)進(jìn)行了識(shí)別和分析。一般的,如果所研究的事物能夠用網(wǎng)絡(luò)連接圖模型來描述,PageRank方法就有其潛在的應(yīng)用價(jià)值,而城市間的鐵路交通網(wǎng)絡(luò)正是一個(gè)典型的網(wǎng)絡(luò)連接圖,本文基于這一事實(shí),并考慮到我國直轄和省會(huì)城市在城市群中的特殊地位,而提出了便于并行計(jì)算的分層PageRank算法,結(jié)合各城市鐵路站點(diǎn)間往返列車班次數(shù)據(jù),將其用于計(jì)算各城市鐵路站點(diǎn)的重要性及排名;為說明所計(jì)算的PageRank值能夠體現(xiàn)城市列車站點(diǎn)的重要性,本文還將城市列車站點(diǎn)的PageRank排名與2020年城市GDP排名進(jìn)行了比較,比較過程中引入了逆序系數(shù)的概念,逆序系數(shù)越大,兩種排序差異越大,而逆序系數(shù)越小,兩種排序越相近。數(shù)據(jù)實(shí)驗(yàn)的結(jié)果顯示城市列車站點(diǎn)的PageRank排名與城市GDP排名的逆序系數(shù)較小,兩者排名較近似,從而佐證了采用分層PageRank算法計(jì)算城市列車站點(diǎn)的重要性是有效的。

    1 數(shù)據(jù)來源

    本文所使用的鐵路列車班次數(shù)據(jù)來自12306網(wǎng)站。收集數(shù)據(jù)時(shí),在12306允許查詢的時(shí)間范圍內(nèi),隨機(jī)選擇了對2021年4月13日的列車班次進(jìn)行了查詢,收集了當(dāng)日國內(nèi)31個(gè)直轄和省會(huì)城市之間往返列車班次數(shù)據(jù)和江西、江蘇兩省省內(nèi)地級市之間往返列車班次數(shù)據(jù)。同時(shí),通過百度查詢采集了2020年相應(yīng)各城市的GDP數(shù)據(jù)。

    2 主要算法和概念

    下面介紹本文的主要算法和概念,它包括用于計(jì)算城市列車站點(diǎn)重要性的分層PageRank算法和用于度量城市列車站點(diǎn)PageRank排名與城市GDP排名差異的逆序系數(shù)的概念,其基本步驟和原理分別如下:

    2.1 分層PageRank算法

    可以認(rèn)為對于一個(gè)城市,如果開往其站點(diǎn)的列車越多、開往的列車班次越多,則該城市的列車站點(diǎn)就越重要,這一點(diǎn)與互聯(lián)網(wǎng)網(wǎng)頁排名的基本假設(shè)完全類似,所以可以采用網(wǎng)頁的PageRank算法對各城市列車站點(diǎn)的重要性進(jìn)行計(jì)算,它仍然是分層PageRank算法的主體思想。不過在運(yùn)用PageRank算法對城市列車站點(diǎn)的重要性進(jìn)行計(jì)算時(shí),可以充分利用已有信息,考慮城市行政級別的層次性和特殊地位,通常直轄和省會(huì)城市能在很大程度上反映本省級行政區(qū)域與區(qū)域外的鐵路連通情況。利用這一點(diǎn)可將我國城市鐵路網(wǎng)絡(luò)進(jìn)行分層,將整個(gè)城市鐵路網(wǎng)絡(luò)劃分為一些子網(wǎng)絡(luò)(以省級行政區(qū)為子網(wǎng)絡(luò)進(jìn)行劃分),再在每個(gè)子網(wǎng)絡(luò)中選擇一個(gè)重要結(jié)點(diǎn)(以直轄和省會(huì)城市為該結(jié)點(diǎn)),這些重要結(jié)點(diǎn)構(gòu)成第1層次的鐵路站點(diǎn)網(wǎng)絡(luò),而這些重要結(jié)點(diǎn)所在的子網(wǎng)絡(luò)均為第2層次的鐵路站點(diǎn)網(wǎng)絡(luò)(有多個(gè)第2層次子網(wǎng)絡(luò)),即把直轄和省會(huì)城市放在第一個(gè)網(wǎng)絡(luò)層次,把各省省會(huì)和省內(nèi)地級市放在第二個(gè)網(wǎng)絡(luò)層次。直觀表述如下圖1所示,細(xì)線圓圈結(jié)點(diǎn)構(gòu)成第1層次的網(wǎng)絡(luò),細(xì)線圓圈表示的結(jié)點(diǎn)就是細(xì)線圓圈所圍的粗線圓圈所表示的結(jié)點(diǎn)中的某個(gè)重要的結(jié)點(diǎn),而細(xì)線圓圈所圍繞的所有粗線圓圈所表示的結(jié)點(diǎn)則構(gòu)成了第2層次的網(wǎng)絡(luò)。

    圖1 分層網(wǎng)絡(luò)示意圖

    分層將一個(gè)大的網(wǎng)絡(luò)分成了若干小網(wǎng)絡(luò),其優(yōu)點(diǎn)是可以并行計(jì)算第1層網(wǎng)絡(luò)和第2層各個(gè)子網(wǎng)絡(luò)結(jié)點(diǎn)的PageRank值,再綜合計(jì)算出整個(gè)綜合網(wǎng)絡(luò)的PageRank值,即未分層網(wǎng)絡(luò)的PageRank值。

    分層PageRank算法步驟:

    (a) 輸入第1層次n個(gè)列車站點(diǎn)間的往返列車班次數(shù)據(jù)矩陣P={pij}n×n,其中pij是從站點(diǎn)j到開往站點(diǎn)i的班車次數(shù)。

    (b) 計(jì)算第1層次n個(gè)站點(diǎn)間的轉(zhuǎn)移概率矩陣M={mij}n×n,其中mij表示旅客從站點(diǎn)j到站點(diǎn)i的概率。根據(jù)往返列車班次數(shù)據(jù)矩陣P={pij}n×n計(jì)算轉(zhuǎn)移概率的公式如下:

    這里假設(shè)了旅客以等可能性乘坐從站點(diǎn)j開往站點(diǎn)i的任意一趟班次的列車。

    圖2 站點(diǎn)連接與往返列車班次示意圖

    為直觀表述,見上圖2,設(shè)有A,B,C三地站點(diǎn),它們間往返列車班次數(shù)據(jù)矩陣P={pij}為

    則轉(zhuǎn)移概率矩陣M={mij}為

    (c) 輸入阻尼因子d和初始向量R0,其中d∈(0,1),此處取d=0.85,而初始向量R0=(R0(1),R0(2),…,R0(n))T設(shè)定了各個(gè)站點(diǎn)的初始重要性,滿足

    (d) 對初始向量關(guān)于轉(zhuǎn)移概率矩陣做帶平滑項(xiàng)的迭代,即

    其中I表示所有分量均為1的n維向量,t=0,1,2,…。

    (e) 如果Rt+1與Rt充分接近,令R=Rt+1,停止迭代。

    (f) 否則,令t=t+1,繼續(xù)執(zhí)行(d)步。

    最后得到的R=(R(1),R(2),…,R(n))T即是第1層次各個(gè)站點(diǎn)的PageRank值,它表示第1層次各個(gè)列車站點(diǎn)的重要性。

    (g) 輸入第1層次站點(diǎn)所對應(yīng)的第2層次子網(wǎng)絡(luò)中各列車站點(diǎn)間的往返列車班次數(shù)據(jù)矩陣,類似步驟(b)~(f),計(jì)算第2層次子網(wǎng)絡(luò)各列車站點(diǎn)的PageRank值,設(shè)第2層次的第k個(gè)子網(wǎng)絡(luò)中含有的結(jié)點(diǎn)數(shù)為mk,k=1,2,…,n,記它們的PageRank值為

    rk=(rk(1),rk(2),…,rk(mk))Tk=1,2,…,n

    注意:此步驟可以與(a)~(f)同步進(jìn)行。

    R(k)rk(mk))T

    (1)

    k=1,2,…,n。綜合PageRank值的計(jì)算亦可如下表1所示。

    表1 綜合PageRank值的計(jì)算

    2.2 逆序系數(shù)—兩種排序差異的度量指標(biāo)

    利用上述分層PageRank算法計(jì)算各城市列車站點(diǎn)的PageRank值并排序,通過對比以GDP為標(biāo)準(zhǔn)而進(jìn)行的城市排名,可以佐證PageRank值是否較好的反映了城市站點(diǎn)的重要性,這是因?yàn)槌鞘蠫DP反映了城市的經(jīng)濟(jì)發(fā)展情況,從某種角度可以認(rèn)為經(jīng)濟(jì)越發(fā)達(dá)的城市越重要,該城市列車站點(diǎn)自然也越重要,那么用城市GDP排序佐證城市列車站點(diǎn)的PageRank排序的關(guān)鍵的問題就是如何度量兩種排序的差異。下面采用逆序系數(shù)來進(jìn)行衡量,其思想來自定義方陣行列式時(shí)所采用的序列逆序數(shù)的概念[8],它的具體概念如下:

    設(shè)有n個(gè)城市,根據(jù)其GDP排序依次為A1,A2,…,An,而根據(jù)其列車站點(diǎn)的PageRank值排序?yàn)锳k1,Ak2,…,Akn,其中Aki∈{A1,A2,…,An}。

    考慮自然數(shù)序列k1,k2,…,kn,如果kp>kq(p

    可以看出逆序系數(shù)取值在[0,1]之內(nèi),逆序系數(shù)越小,兩種排序越接近,逆序系數(shù)越大,兩種排序差異越大。設(shè)有A、B兩組城市,按城市的GDP排名分別為A1,A2,A3,A4和B1,B2,B3,B4,B5,如果它們列車站點(diǎn)的PageRank排名分別為A3,A2,A1,A4和B4,B1,B2,B3,B5,則A組城市有(3,2)、(3,1)、(2,1)這3個(gè)逆序,逆序系數(shù)為τA=3/[4(4-1)/2]=1/2;而B組城市有(4,1)、(4,2)、(4,3)這3個(gè)逆序,逆序系數(shù)為τB=3/[5(5-1)/2]=3/10,因而認(rèn)為B組城市列車站點(diǎn)的PageRank排名比A組更接近其城市GDP排名。

    3 數(shù)據(jù)實(shí)驗(yàn)

    根據(jù)所采集的列車班次數(shù)據(jù),并將直轄和省會(huì)城市作為第1層網(wǎng)絡(luò)結(jié)點(diǎn),各省會(huì)城市和本省內(nèi)地級市作為第2層子網(wǎng)絡(luò)結(jié)點(diǎn)。由2.1節(jié)分層PageRank算法的(a)~(f)步,計(jì)算各直轄和省會(huì)城市的列車站點(diǎn)的第1層PageRank值如下表2的第2列,按PageRank值從大到小向下排列;表2的第3列城市名前面的序號(hào)是根據(jù)第4列GDP(單位為:億元)從大到小排序的序號(hào)。

    以城市GDP的排序序號(hào)為準(zhǔn),表2的城市根據(jù)列車站點(diǎn)的第1層PageRank值調(diào)整后,城市對應(yīng)的序號(hào)列表如下,如北京對應(yīng)2,鄭州對應(yīng)11,上海對應(yīng)1等等:

    [2,11,1,8,7,10,19,6,9,14,12,3,18,13,23,5,20,4,26,17,21,16,24,15,22,30,28,27,25,31,29]

    表2 直轄和省會(huì)城市的第1層PageRank值

    由2.1節(jié)分層PageRank算法的(g)步,計(jì)算江西和江蘇兩省省內(nèi)地級市列車站點(diǎn)的第2層PageRank值,如下表3和表4的第2列。

    表3 江西城市第2層PageRank值

    表4 江蘇城市第2層PageRank值

    同上,由表3與表4可知,江西和江蘇兩省城市站點(diǎn)的第2層PageRank值排序和城市GDP排序的逆序系數(shù)分別為:

    由2.1節(jié)分層PageRank算法的(h)步,計(jì)算江西和江蘇兩省省內(nèi)地級市列車站點(diǎn)的綜合PageRank值如下表5的第2列

    表5 江西、江蘇兩省城市綜合PageRank值

    由表5可知,江西江蘇兩省城市站點(diǎn)的綜合PageRank值排序和城市GDP排序的逆序系數(shù)為:

    對于通常PageRank的計(jì)算有許多不同的數(shù)值計(jì)算方法,例如文獻(xiàn)[9]和[10]分別提出了一種基于子空間的啟發(fā)式搜索算法和采用多冪次、多分裂的內(nèi)外迭代算法。而采用2.1節(jié)分層PageRank算法對上述表2至表5中的PageRank值進(jìn)行計(jì)算時(shí),可以合并步驟(c)(d)(e)(f)。因?yàn)橥ㄟ^下述公式

    迭代時(shí),如果Rt+1與Rt能充分接近,說明極限存在,將上述迭代格式兩邊取極限(Rt→R),那么迭代計(jì)算R則變?yōu)榍蠼饩€性方程組

    (2)

    對于線性方程組的求解,文獻(xiàn)[11-12]歸納了多種求解算法,可以使用已有且成熟的數(shù)學(xué)工具包進(jìn)行實(shí)驗(yàn),其中R為PageRank值,d取0.85,n表示列車站點(diǎn)的個(gè)數(shù),M表示列車站點(diǎn)間的轉(zhuǎn)移概率矩陣,I是所有分量均為1的n維向量。

    表2至表5中PageRank值計(jì)算流程圖如下

    圖3 表2至表5中的PageRank值的計(jì)算流程圖

    4 分析與結(jié)果

    在表2中可以看到,首都北京的第1層PageRank值排名第一,而表3和表4又可以看到江西和江蘇兩省的省會(huì)城市南昌和南京的第2層PageRank值都排第一,而綜合PageRank值排序仍然保持子網(wǎng)絡(luò)的原順序,見第5節(jié)補(bǔ)充說明部分的圖4解釋.直轄和省會(huì)城市的PageRank值排名第一與鐵路布局政策的傾向是吻合的。此外,在表2中還可以看到鄭州鐵路站點(diǎn)的PageRank值很大,這與鄭州站是國內(nèi)鐵路網(wǎng)絡(luò)的一個(gè)樞紐事實(shí)相吻合,同樣的,由表3和表4可以發(fā)現(xiàn)上饒和蘇州是江西和江蘇兩省省內(nèi)的重要鐵路交匯點(diǎn),這些PageRank數(shù)據(jù)與事實(shí)的相符性在一定程度上說明了分層PageRank算法用于計(jì)算我國城市鐵路站點(diǎn)的重要性是合理的。

    更進(jìn)一步的對比城市GDP排序,可以看出國內(nèi)31個(gè)直轄和省會(huì)城市構(gòu)成的第1層城市網(wǎng)絡(luò)的PageRank排序和GDP排序是比較接近的,逆序系數(shù)僅為0.2,與0更接近而與1相距更遠(yuǎn),小于0.5;第2層PageRank值(0.4、0.205128)和綜合PageRank值(0.205128)也均小于0.5。這說明如果城市GDP可以反映城市的重要性(至少在經(jīng)濟(jì)發(fā)展角度上可以承認(rèn)),那么用分層PageRank算法計(jì)算的城市站點(diǎn)的重要性可以從城市GDP方面得到佐證。繼而可看到,江西城市站點(diǎn)的第2層網(wǎng)絡(luò)逆序系數(shù)0.4比江蘇城市站點(diǎn)的第2層網(wǎng)絡(luò)逆序系數(shù)0.205128要大;另外,江西和江蘇兩省城市網(wǎng)絡(luò)的綜合PageRank排序與GDP排序的逆序系數(shù)為0.286232,恰好在江西第2層城市網(wǎng)絡(luò)逆序系數(shù)0.4與江蘇第2層城市網(wǎng)絡(luò)逆序系數(shù)0.205128之間,說明如果把江西和江蘇作為一個(gè)綜合城市網(wǎng)絡(luò)來看,兩省城市站點(diǎn)綜合網(wǎng)絡(luò)比江蘇城市站點(diǎn)網(wǎng)絡(luò)變的偏離了GDP排序,但比江西城市站點(diǎn)網(wǎng)絡(luò)變得貼近了GDP排序,這反映了江西鐵路交通相比江蘇而言有進(jìn)一步提升加強(qiáng)的空間。這些計(jì)算結(jié)果更進(jìn)一步的說明了由分層PageRank算法計(jì)算的鐵路站點(diǎn)的PageRank值有效的刻畫了列車站點(diǎn)的重要性。

    5 補(bǔ)充說明

    最后需要補(bǔ)充說明的是,由兩個(gè)第2層子網(wǎng)絡(luò)構(gòu)成的綜合網(wǎng)絡(luò)的逆序系數(shù)不一定總在第2層子網(wǎng)絡(luò)的兩個(gè)逆序系數(shù)之間,但是已知第2層子網(wǎng)絡(luò)的逆序系數(shù),是可以給出綜合網(wǎng)絡(luò)逆序系數(shù)的上下界估計(jì)的。推而廣之,由d個(gè)第2層子網(wǎng)絡(luò)構(gòu)成的綜合網(wǎng)絡(luò),其逆序系數(shù)的估計(jì)也可以給出,其下界估計(jì)是

    (3)

    而上界估計(jì)是

    (4)

    其中τi,ni分別表示第2層網(wǎng)絡(luò)中第i個(gè)子網(wǎng)絡(luò)的逆序系數(shù)和第i個(gè)子網(wǎng)絡(luò)的結(jié)點(diǎn)個(gè)數(shù),i=1,2,…,d。

    為便于說明,下面以第2層三個(gè)結(jié)點(diǎn)的A網(wǎng)絡(luò)和四個(gè)結(jié)點(diǎn)的B網(wǎng)絡(luò)構(gòu)成的綜合網(wǎng)絡(luò)逆序系數(shù)估計(jì)為例進(jìn)行分析,根據(jù)2.1節(jié)分層PageRank算法(h)步的(1)式可知,子網(wǎng)絡(luò)結(jié)點(diǎn)在第2層網(wǎng)絡(luò)的PageRank排序與在綜合網(wǎng)絡(luò)中的綜合PageRank排序是保持一致的,綜合后可能發(fā)生兩類種情況,一類是兩組結(jié)點(diǎn)按原順序交叉排序,一類是兩組結(jié)點(diǎn)按原順序分離排序,如下圖4所示。

    無論哪種情況,逆序個(gè)數(shù)都不會(huì)減少,而新增逆序個(gè)數(shù)要么為0,要么最多為3×4,原有的逆序數(shù)為

    圖4 子網(wǎng)絡(luò)結(jié)點(diǎn)的綜合排序示意圖

    τA3(3-1)/2和τB4(4-1)/2,所以綜合網(wǎng)絡(luò)的逆序個(gè)數(shù)在

    之間,而綜合網(wǎng)絡(luò)的總結(jié)點(diǎn)數(shù)變?yōu)榱?+4,所以逆序系數(shù)在

    之間,其中τA,τB分別是第2層子網(wǎng)絡(luò)A,B的逆序系數(shù)。用完全類似的推理可得,一般綜合網(wǎng)絡(luò)的逆序系數(shù)上下界估計(jì)式(3)和(4),用(3)式和(4)式計(jì)算江西和江蘇兩省城市綜合網(wǎng)絡(luò)的逆序系數(shù)在

    之間,這與江西江蘇兩省綜合逆序系數(shù)為0.286 232吻合。

    猜你喜歡
    逆序結(jié)點(diǎn)排序
    排序不等式
    恐怖排序
    有界線性算子的Drazin逆的逆序律
    關(guān)于矩陣廣義BottDuffin逆的逆序律
    新中國70年漢語逆序詞研究(1949—2019)
    節(jié)日排序
    刻舟求劍
    兒童繪本(2018年5期)2018-04-12 16:45:32
    Ladyzhenskaya流體力學(xué)方程組的確定模與確定結(jié)點(diǎn)個(gè)數(shù)估計(jì)
    對外漢語教學(xué)中AB-BA式逆序詞教學(xué)分析
    基于Raspberry PI為結(jié)點(diǎn)的天氣云測量網(wǎng)絡(luò)實(shí)現(xiàn)
    91大片在线观看| 日韩人妻精品一区2区三区| 多毛熟女@视频| 亚洲激情五月婷婷啪啪| 日韩,欧美,国产一区二区三区| 日韩大片免费观看网站| 在线精品无人区一区二区三| 国产精品久久久久久精品电影小说| 99久久国产精品久久久| 欧美激情极品国产一区二区三区| 久久人人97超碰香蕉20202| 9191精品国产免费久久| 色综合欧美亚洲国产小说| 五月开心婷婷网| 亚洲熟女毛片儿| 人人妻人人添人人爽欧美一区卜| 国产麻豆69| 国产精品 国内视频| 国产三级黄色录像| 久久狼人影院| 亚洲成人手机| 桃花免费在线播放| 久久人妻福利社区极品人妻图片| 国产不卡av网站在线观看| 99久久人妻综合| 日韩 亚洲 欧美在线| 久久久久久亚洲精品国产蜜桃av| 男女之事视频高清在线观看| 国产精品久久久久久精品古装| 50天的宝宝边吃奶边哭怎么回事| tocl精华| 国产有黄有色有爽视频| 亚洲国产精品999| 新久久久久国产一级毛片| 建设人人有责人人尽责人人享有的| 亚洲欧美日韩高清在线视频 | 建设人人有责人人尽责人人享有的| www日本在线高清视频| 欧美另类亚洲清纯唯美| 免费不卡黄色视频| 亚洲激情五月婷婷啪啪| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 国产福利在线免费观看视频| 黄色视频不卡| 国产精品熟女久久久久浪| 操美女的视频在线观看| 老鸭窝网址在线观看| 久久亚洲国产成人精品v| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产中文字幕在线视频| 波多野结衣一区麻豆| 老司机在亚洲福利影院| 久久精品亚洲av国产电影网| 国产精品免费视频内射| 久久九九热精品免费| 午夜福利影视在线免费观看| 亚洲色图综合在线观看| 久久热在线av| 黄色怎么调成土黄色| 国产1区2区3区精品| 久久ye,这里只有精品| 欧美少妇被猛烈插入视频| √禁漫天堂资源中文www| 欧美精品av麻豆av| 久久久精品94久久精品| 爱豆传媒免费全集在线观看| 69精品国产乱码久久久| 91麻豆精品激情在线观看国产 | 国产区一区二久久| 久久香蕉激情| 一区二区三区精品91| 十八禁网站免费在线| av视频免费观看在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 伊人久久大香线蕉亚洲五| 婷婷成人精品国产| 夜夜骑夜夜射夜夜干| 美女国产高潮福利片在线看| a 毛片基地| 黑人操中国人逼视频| 色综合欧美亚洲国产小说| 欧美另类一区| 亚洲久久久国产精品| 久久精品国产亚洲av香蕉五月 | av电影中文网址| 久久国产亚洲av麻豆专区| 亚洲国产精品成人久久小说| 不卡一级毛片| 一本综合久久免费| 法律面前人人平等表现在哪些方面 | 秋霞在线观看毛片| 老司机影院成人| 亚洲精品国产av成人精品| 一区二区三区四区激情视频| 一级毛片电影观看| 久久久久国产一级毛片高清牌| 日本一区二区免费在线视频| 汤姆久久久久久久影院中文字幕| 欧美xxⅹ黑人| 丝袜美腿诱惑在线| 黄色视频不卡| 一二三四社区在线视频社区8| 老汉色av国产亚洲站长工具| 久热这里只有精品99| 国产精品麻豆人妻色哟哟久久| 在线亚洲精品国产二区图片欧美| 精品少妇一区二区三区视频日本电影| 99国产精品一区二区三区| 国产成人精品久久二区二区免费| 久久热在线av| 美女扒开内裤让男人捅视频| 日韩大片免费观看网站| 三级毛片av免费| 午夜精品久久久久久毛片777| 免费观看a级毛片全部| 欧美成狂野欧美在线观看| 大香蕉久久成人网| av网站免费在线观看视频| 亚洲七黄色美女视频| 午夜福利免费观看在线| 成人免费观看视频高清| 亚洲精品国产精品久久久不卡| 人妻一区二区av| kizo精华| 国产一卡二卡三卡精品| 在线永久观看黄色视频| 免费少妇av软件| 国产淫语在线视频| 欧美97在线视频| 国产成人系列免费观看| 午夜福利视频在线观看免费| 久久精品国产a三级三级三级| 嫩草影视91久久| 亚洲专区中文字幕在线| 亚洲专区中文字幕在线| 国产在线一区二区三区精| www.自偷自拍.com| 日韩制服丝袜自拍偷拍| 丝瓜视频免费看黄片| 午夜精品国产一区二区电影| 飞空精品影院首页| 婷婷色av中文字幕| 国产精品一二三区在线看| 欧美成人午夜精品| 久久人人爽人人片av| 成年人午夜在线观看视频| 国产精品99久久99久久久不卡| 免费av中文字幕在线| 成人三级做爰电影| 精品欧美一区二区三区在线| 色综合欧美亚洲国产小说| 老司机福利观看| 国产激情久久老熟女| 欧美另类亚洲清纯唯美| 日本猛色少妇xxxxx猛交久久| 亚洲欧洲精品一区二区精品久久久| 50天的宝宝边吃奶边哭怎么回事| 国产成人啪精品午夜网站| 操美女的视频在线观看| 精品欧美一区二区三区在线| 91麻豆av在线| 国产成人啪精品午夜网站| 中国国产av一级| 制服人妻中文乱码| 欧美久久黑人一区二区| 在线观看免费高清a一片| 黄色毛片三级朝国网站| 热re99久久国产66热| 丰满少妇做爰视频| 一本综合久久免费| 五月天丁香电影| 午夜激情av网站| 一边摸一边做爽爽视频免费| 一边摸一边做爽爽视频免费| 成人av一区二区三区在线看 | 国产成人精品久久二区二区91| 色精品久久人妻99蜜桃| 男女边摸边吃奶| 中亚洲国语对白在线视频| 国产精品偷伦视频观看了| 久久香蕉激情| 国产成人av教育| 少妇人妻久久综合中文| 国产区一区二久久| 人妻久久中文字幕网| 中文字幕最新亚洲高清| 久久青草综合色| 王馨瑶露胸无遮挡在线观看| 国产精品熟女久久久久浪| 一二三四在线观看免费中文在| 建设人人有责人人尽责人人享有的| 国产成人欧美| 好男人电影高清在线观看| 91成人精品电影| 国产真人三级小视频在线观看| 久热这里只有精品99| 亚洲男人天堂网一区| 亚洲男人天堂网一区| 19禁男女啪啪无遮挡网站| 99国产综合亚洲精品| 国产成人精品久久二区二区91| 午夜福利在线观看吧| 亚洲美女黄色视频免费看| 国产国语露脸激情在线看| 一本—道久久a久久精品蜜桃钙片| 亚洲国产中文字幕在线视频| 亚洲情色 制服丝袜| 亚洲精品av麻豆狂野| 午夜福利影视在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久大尺度免费视频| 国产深夜福利视频在线观看| 国产日韩欧美在线精品| 欧美老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 免费观看av网站的网址| 午夜福利在线观看吧| tocl精华| 最新在线观看一区二区三区| av网站在线播放免费| 亚洲情色 制服丝袜| netflix在线观看网站| 亚洲午夜精品一区,二区,三区| 国产亚洲一区二区精品| 丝瓜视频免费看黄片| h视频一区二区三区| 国产一区二区三区综合在线观看| 成年人免费黄色播放视频| 国产有黄有色有爽视频| 又大又爽又粗| 久久中文字幕一级| 亚洲av日韩在线播放| 久久精品人人爽人人爽视色| 视频区欧美日本亚洲| 99国产精品一区二区蜜桃av | 秋霞在线观看毛片| 亚洲伊人色综图| 久久天堂一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区 | 日本a在线网址| 中文欧美无线码| 亚洲精品美女久久久久99蜜臀| 女人久久www免费人成看片| 香蕉丝袜av| 又大又爽又粗| 免费不卡黄色视频| 久久国产精品人妻蜜桃| 女人爽到高潮嗷嗷叫在线视频| 激情视频va一区二区三区| 在线观看人妻少妇| a级片在线免费高清观看视频| 亚洲精品中文字幕一二三四区 | 大香蕉久久网| 国产精品 国内视频| 亚洲国产av新网站| 国产高清国产精品国产三级| 欧美日韩国产mv在线观看视频| 亚洲欧美精品综合一区二区三区| 国产深夜福利视频在线观看| 少妇被粗大的猛进出69影院| 视频区欧美日本亚洲| 电影成人av| 国产麻豆69| 色视频在线一区二区三区| 午夜福利影视在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品久久成人aⅴ小说| 亚洲av电影在线进入| 汤姆久久久久久久影院中文字幕| 性少妇av在线| 精品一区二区三区av网在线观看 | 精品久久久久久电影网| 正在播放国产对白刺激| 丰满少妇做爰视频| 中亚洲国语对白在线视频| 91av网站免费观看| 99香蕉大伊视频| 搡老乐熟女国产| 男人爽女人下面视频在线观看| 国产精品影院久久| 老汉色av国产亚洲站长工具| 美女大奶头黄色视频| 亚洲国产av新网站| 国产精品一区二区免费欧美 | 欧美激情 高清一区二区三区| 免费观看人在逋| 在线观看免费午夜福利视频| 国产精品免费大片| 欧美乱码精品一区二区三区| 热re99久久精品国产66热6| 国产xxxxx性猛交| 91成年电影在线观看| 中文欧美无线码| 亚洲国产中文字幕在线视频| 天天躁日日躁夜夜躁夜夜| 不卡一级毛片| 久久久精品免费免费高清| 欧美日韩亚洲综合一区二区三区_| 大香蕉久久成人网| 亚洲精品国产一区二区精华液| 午夜日韩欧美国产| 每晚都被弄得嗷嗷叫到高潮| 国产又色又爽无遮挡免| 天堂中文最新版在线下载| 老司机午夜十八禁免费视频| 精品国内亚洲2022精品成人 | avwww免费| 天堂中文最新版在线下载| 国产成人免费观看mmmm| 亚洲成人手机| 婷婷成人精品国产| 1024香蕉在线观看| cao死你这个sao货| 成人国语在线视频| 色播在线永久视频| 欧美 日韩 精品 国产| 久久av网站| 汤姆久久久久久久影院中文字幕| 高清av免费在线| 五月天丁香电影| 免费黄频网站在线观看国产| 十八禁网站免费在线| 免费高清在线观看视频在线观看| 免费久久久久久久精品成人欧美视频| av电影中文网址| 国产在线一区二区三区精| 亚洲国产日韩一区二区| a级毛片在线看网站| 一本色道久久久久久精品综合| 久久国产精品人妻蜜桃| 天天躁日日躁夜夜躁夜夜| 啦啦啦在线免费观看视频4| 国产精品国产三级国产专区5o| 日韩视频在线欧美| 国产成人精品无人区| 国产国语露脸激情在线看| 黄片播放在线免费| 一二三四在线观看免费中文在| 欧美在线黄色| 最近最新中文字幕大全免费视频| av又黄又爽大尺度在线免费看| 欧美 亚洲 国产 日韩一| 日本91视频免费播放| 最近中文字幕2019免费版| 亚洲伊人久久精品综合| 一级毛片精品| 成人av一区二区三区在线看 | 777米奇影视久久| 搡老熟女国产l中国老女人| 亚洲专区中文字幕在线| 五月开心婷婷网| bbb黄色大片| 777米奇影视久久| 淫妇啪啪啪对白视频 | 国产精品影院久久| 日韩三级视频一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 国产不卡av网站在线观看| a级毛片在线看网站| 男女国产视频网站| 又紧又爽又黄一区二区| 国产精品久久久人人做人人爽| 一级,二级,三级黄色视频| 法律面前人人平等表现在哪些方面 | 91精品国产国语对白视频| 国产成+人综合+亚洲专区| 国产一级毛片在线| 1024香蕉在线观看| 久久久水蜜桃国产精品网| 欧美精品一区二区大全| 最新在线观看一区二区三区| 少妇 在线观看| 91麻豆精品激情在线观看国产 | 亚洲黑人精品在线| 一级片'在线观看视频| av电影中文网址| 久久精品国产a三级三级三级| 国产一区二区激情短视频 | 美女国产高潮福利片在线看| 久久国产精品男人的天堂亚洲| 黄频高清免费视频| 可以免费在线观看a视频的电影网站| 色播在线永久视频| 9热在线视频观看99| 亚洲国产精品999| 亚洲人成电影观看| 中国国产av一级| 亚洲三区欧美一区| 久久亚洲国产成人精品v| 亚洲人成电影观看| 久久国产精品人妻蜜桃| 久热爱精品视频在线9| 一区二区日韩欧美中文字幕| 欧美精品一区二区大全| 免费av中文字幕在线| 久久久久国产精品人妻一区二区| 午夜免费观看性视频| 亚洲精品日韩在线中文字幕| 国产成人欧美在线观看 | 亚洲成av片中文字幕在线观看| 青草久久国产| 黄色怎么调成土黄色| 久久久久久免费高清国产稀缺| 亚洲精品一二三| 1024视频免费在线观看| 中国国产av一级| 亚洲精品在线美女| 老汉色av国产亚洲站长工具| 中文字幕人妻熟女乱码| 久久久久久久精品精品| 纯流量卡能插随身wifi吗| 成人免费观看视频高清| 极品人妻少妇av视频| 午夜福利影视在线免费观看| 国产精品麻豆人妻色哟哟久久| 最近最新中文字幕大全免费视频| 高清av免费在线| 少妇 在线观看| 深夜精品福利| 国产熟女午夜一区二区三区| 国产在视频线精品| 极品少妇高潮喷水抽搐| 人妻 亚洲 视频| 国产av一区二区精品久久| 麻豆av在线久日| 亚洲精品久久午夜乱码| 纯流量卡能插随身wifi吗| 免费观看人在逋| 中国国产av一级| 在线亚洲精品国产二区图片欧美| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频| 国产无遮挡羞羞视频在线观看| 亚洲九九香蕉| 日日爽夜夜爽网站| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区 | 午夜激情久久久久久久| www日本在线高清视频| 一本一本久久a久久精品综合妖精| 狂野欧美激情性xxxx| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区精品| 欧美国产精品va在线观看不卡| 另类亚洲欧美激情| 久久久国产一区二区| 在线观看人妻少妇| 少妇人妻久久综合中文| 日日爽夜夜爽网站| 午夜两性在线视频| 亚洲人成电影观看| 狂野欧美激情性bbbbbb| cao死你这个sao货| 老熟妇乱子伦视频在线观看 | 麻豆av在线久日| 国产精品一区二区精品视频观看| 亚洲精品中文字幕一二三四区 | 亚洲五月色婷婷综合| 亚洲精品中文字幕在线视频| 亚洲一区二区三区欧美精品| 一级a爱视频在线免费观看| 国产福利在线免费观看视频| 色老头精品视频在线观看| 免费观看av网站的网址| 不卡av一区二区三区| 久久女婷五月综合色啪小说| 国产99久久九九免费精品| 亚洲av美国av| 中国国产av一级| 他把我摸到了高潮在线观看 | 看免费av毛片| 深夜精品福利| 91大片在线观看| 在线永久观看黄色视频| 多毛熟女@视频| 美女午夜性视频免费| 成人国语在线视频| 精品人妻1区二区| 香蕉国产在线看| 咕卡用的链子| 精品国产一区二区三区久久久樱花| 99九九在线精品视频| 香蕉丝袜av| 色播在线永久视频| 叶爱在线成人免费视频播放| 又黄又粗又硬又大视频| 国内毛片毛片毛片毛片毛片| 亚洲专区字幕在线| 男女免费视频国产| 91精品国产国语对白视频| 亚洲人成77777在线视频| 欧美日韩国产mv在线观看视频| 久久精品人人爽人人爽视色| 91精品国产国语对白视频| 国产免费现黄频在线看| 亚洲欧美一区二区三区久久| 女人久久www免费人成看片| 99久久99久久久精品蜜桃| 国产有黄有色有爽视频| 亚洲精品国产av成人精品| 欧美日韩亚洲综合一区二区三区_| 王馨瑶露胸无遮挡在线观看| 99国产精品免费福利视频| 中国国产av一级| 久久av网站| 亚洲欧洲精品一区二区精品久久久| 亚洲精品久久成人aⅴ小说| 国产成人啪精品午夜网站| 日本欧美视频一区| 精品欧美一区二区三区在线| 美女午夜性视频免费| 亚洲avbb在线观看| 久久国产精品男人的天堂亚洲| 欧美另类亚洲清纯唯美| 日韩欧美一区视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 天天影视国产精品| 欧美久久黑人一区二区| 亚洲av美国av| 窝窝影院91人妻| 最近中文字幕2019免费版| 亚洲专区字幕在线| 天堂中文最新版在线下载| 一本大道久久a久久精品| 少妇裸体淫交视频免费看高清 | 免费高清在线观看视频在线观看| 激情视频va一区二区三区| 男人操女人黄网站| 久久久久久久大尺度免费视频| 精品久久久久久电影网| 天天躁狠狠躁夜夜躁狠狠躁| 成人国语在线视频| 亚洲国产看品久久| 国产成人啪精品午夜网站| 久久亚洲精品不卡| 国产高清videossex| 宅男免费午夜| 啦啦啦中文免费视频观看日本| 视频区欧美日本亚洲| 中文字幕高清在线视频| 欧美黑人精品巨大| 国产成人啪精品午夜网站| 免费在线观看完整版高清| 国产99久久九九免费精品| 五月开心婷婷网| 亚洲中文日韩欧美视频| 黄片小视频在线播放| 国内毛片毛片毛片毛片毛片| 在线观看舔阴道视频| 精品国产乱码久久久久久男人| 超碰成人久久| a级毛片在线看网站| 97人妻天天添夜夜摸| 麻豆av在线久日| 亚洲成人免费电影在线观看| 欧美精品一区二区免费开放| 老汉色∧v一级毛片| 欧美成狂野欧美在线观看| 亚洲中文字幕日韩| 国产精品偷伦视频观看了| 韩国高清视频一区二区三区| 中国美女看黄片| 亚洲精品第二区| 亚洲精品自拍成人| 亚洲精品一卡2卡三卡4卡5卡 | 性少妇av在线| 国产精品麻豆人妻色哟哟久久| 欧美另类亚洲清纯唯美| 久久天堂一区二区三区四区| 国产成人av激情在线播放| videosex国产| 国产成人a∨麻豆精品| 免费高清在线观看日韩| 亚洲欧美清纯卡通| 日韩欧美免费精品| 又黄又粗又硬又大视频| 中文字幕另类日韩欧美亚洲嫩草| 女人爽到高潮嗷嗷叫在线视频| 国内毛片毛片毛片毛片毛片| av不卡在线播放| 中文字幕人妻丝袜制服| 国产免费一区二区三区四区乱码| 啦啦啦 在线观看视频| 高潮久久久久久久久久久不卡| 欧美日本中文国产一区发布| 亚洲av片天天在线观看| 狠狠婷婷综合久久久久久88av| 国产精品熟女久久久久浪| 亚洲国产欧美一区二区综合| 欧美成人午夜精品| 中文字幕人妻熟女乱码| 女人高潮潮喷娇喘18禁视频| av又黄又爽大尺度在线免费看| 午夜福利乱码中文字幕| 91大片在线观看| 日本黄色日本黄色录像| 少妇猛男粗大的猛烈进出视频| videosex国产| 亚洲av片天天在线观看| 99精国产麻豆久久婷婷| 一本大道久久a久久精品| 人成视频在线观看免费观看| 色老头精品视频在线观看| 亚洲人成电影观看| 国产欧美日韩一区二区三区在线| 免费不卡黄色视频| 两个人免费观看高清视频| xxxhd国产人妻xxx| 国产精品一区二区在线观看99| 欧美乱码精品一区二区三区| 免费观看av网站的网址| 9热在线视频观看99| 国产成人影院久久av| 飞空精品影院首页|