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

    兩層星型網(wǎng)絡(luò)上的傳染病建模和控制

    2017-07-07 02:15:16鄭國慶唐清干祝光湖
    關(guān)鍵詞:星型邊界點中心點

    鄭國慶,唐清干,祝光湖

    (桂林電子科技大學(xué)數(shù)學(xué)與計算科學(xué)學(xué)院,桂林 541004)

    0 引言

    傳染病數(shù)學(xué)模型可以追溯到1927年Kermack與McKendrick創(chuàng)立的“倉室模型”,至今仍在廣泛應(yīng)用,并被不斷發(fā)展和推廣[1-2]。傳染病模型能夠描繪疾病的傳播過程,解釋疾病的流行因素,并為疾病的防治工作給予理論指導(dǎo)[1]。模型中倉室結(jié)構(gòu)以及個體在倉室的轉(zhuǎn)移由具體疾病決定,考慮不同的疾病和傳播方式,現(xiàn)已建立了非常多的傳染病模型[1-2]。經(jīng)典的有SIS和SIR模型,如在SIS模型中,所有個體將會被分成兩類:S-易感者(Susceptibles)和I-染病者(Infected)。易感者被感染者感染后可以恢復(fù),但不具備免疫力。模型用常微分方程描述如下:

    其中β=pC(N)表示傳染率,γ表示恢復(fù)率,C(N)是接觸率,p表示每次接觸的傳染概率。傳統(tǒng)模型假設(shè)群體是均勻混合,個體間的接觸機會相等,如上式中C和β,這種假設(shè)往往無法體現(xiàn)個體的異質(zhì)性,相應(yīng)的模型往往不夠準(zhǔn)確。而新興的復(fù)雜網(wǎng)絡(luò)理論為傳染病模型的構(gòu)建提供了很好的工具,能夠更精確地描述傳播過程[3-4]。

    從20世紀(jì)末開始,復(fù)雜網(wǎng)絡(luò)研究逐步滲透到數(shù)理、生命和工程等眾多不同的領(lǐng)域,并獲得了很多開創(chuàng)性的成果[3-16]。在眾多的網(wǎng)絡(luò)結(jié)構(gòu)中,星型網(wǎng)絡(luò)是一種典型的規(guī)則網(wǎng)絡(luò),由一個中心節(jié)點和多個邊界點構(gòu)成,每一個邊界點都與中心節(jié)點相連,但邊界點之間不連接,如圖1a所示。如局域網(wǎng),內(nèi)部所有個人計算機與一個總服務(wù)器相連,形成以總服務(wù)器為中心的一個星型網(wǎng)絡(luò)。星型網(wǎng)絡(luò)也可推廣到多個中心的情形,如在社會網(wǎng)絡(luò)中以少數(shù)人為中心的社區(qū)或單位[4]。最近研究表明系統(tǒng)都不是孤立存在的,因而由不同子網(wǎng)相互作用構(gòu)成的耦合網(wǎng)絡(luò)更能準(zhǔn)確地刻畫真實系統(tǒng),如交通網(wǎng)絡(luò)和計算機網(wǎng)絡(luò)等[4-6]。本文將把單層星型網(wǎng)絡(luò)推廣到兩層星型耦合網(wǎng)絡(luò),如圖1b所示。這是一類特色的社團(tuán)結(jié)構(gòu),特別的,在中心城市及其周邊縣城構(gòu)成的交通系統(tǒng)中,如果把兩個中心城看成是中心節(jié)點,周邊縣城看出邊緣節(jié)點,而忽略交通量較少的縣城連接,則此系統(tǒng)可用兩層星型耦合網(wǎng)絡(luò)表示。研究此網(wǎng)絡(luò)上的傳播能夠揭示疾病在中心和周邊地區(qū)傳播擴散的規(guī)律并獲得有效的防控策略?,F(xiàn)實生活也有類似案例,如2003年非典先從中國經(jīng)濟(jì)繁榮、人口流動性大的廣州和北京這兩個中心城市開始蔓延至周邊地區(qū),之后引起大范圍的擴散。

    圖1 兩種星型網(wǎng)絡(luò)兩層星際網(wǎng)絡(luò)由兩個單層星型模型構(gòu)成,兩個中心節(jié)相連,每層中的邊界點除了與該層的中心節(jié)點相連還與另一層中對應(yīng)的邊界點相連。

    傳染病的防控策略主要有3種:控制感染源、切斷傳播途徑和保護(hù)感染者,具體有隔離和接種等。例如在非典全面爆發(fā)時期,政府將被感染者和疑似感染者進(jìn)行隔離,對所有單位和學(xué)校等公共場所定期消毒,并號召人們出行帶著口罩。這些措施有效減輕了SARS病毒擴散。另外,疫苗可以刺激機體產(chǎn)生抗體,通過提高對病毒的免疫力來預(yù)防疾病,進(jìn)而減少疾病傳播的范圍。世界衛(wèi)生組織的數(shù)據(jù)顯示,通過接種全球麻疹死亡人數(shù)從2000年的54.4萬例下降到2013年的14.5萬例,降幅約75%。最近針對復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)特點,人們提出了幾類有效的免疫策略,如熟人免疫、目標(biāo)免疫和主動免疫等[7-9]。研究發(fā)現(xiàn)網(wǎng)絡(luò)的異質(zhì)性和個體行為的自適應(yīng)性很大程度上左右著疾病的傳播和控制[3,7-10],因此研究特定網(wǎng)絡(luò)上的傳染病擴散和防控有著實際意義。

    本文利用Jennifer等人提出的有效度模型[11-12],研究兩層星型網(wǎng)絡(luò)上疾病的傳播,探索有效的疾病控制策略。下一節(jié)我們給出刻畫模型的常微分方程組。第3節(jié)我們采用The next generation matrix方法分析疾病的基本再生數(shù)。第4節(jié)我們介紹兩種防控疾病的方法:截斷傳播途徑和對易感者免疫。第5節(jié)中,我們利用Matlab分別對兩種防控方法做數(shù)值模擬,通過依次改變傳染率,群體規(guī)模和易感者的初值比較兩種方法的優(yōu)劣性。最后我們對結(jié)果進(jìn)行分析與討論。

    1 有效度模型

    有效度模型是由Jennifer等人于2010年提出的一種新的網(wǎng)絡(luò)傳染病的建模思想[11-13]。與傳統(tǒng)的平均場模型相比,它不僅將所有節(jié)點根據(jù)狀態(tài)分為易感者或感染者,而且還記錄了相連的鄰居個數(shù)和狀態(tài),例如用Ssi表示具有s個易感鄰居和i個感染鄰居的易感節(jié)點數(shù)量;如果一個Ssi類節(jié)點被它的i個感染鄰居中的某個感染,則從Ssi類變?yōu)镮si類節(jié)點,節(jié)點狀態(tài)發(fā)生了變化;如果一個Ssi類節(jié)點的i個感染鄰居中某個恢復(fù),節(jié)點就會從Ssi類變?yōu)镾s+1,i-1類;如果某個Ssi類節(jié)點的s個易感鄰居中某個被感染,節(jié)點就從Ssi類變?yōu)镾s-1,i+1類。

    由此靜態(tài)網(wǎng)絡(luò)上SIS有效度模型可表示為[13]:

    (1)

    (2)

    根據(jù)上述模型,本文分析兩層星型網(wǎng)絡(luò)傳染病的有效度模型的常微分系統(tǒng)。假設(shè)每層邊界點的個數(shù)為n,中心點個數(shù)為1,恢復(fù)率γ=1。對于邊界點,它的度是2,即k=s+i=2。對應(yīng)的有效度模型為

    而對于中心節(jié)點,它的度是n+1,即k=s+i=n+1,有效度模型為

    其中,0≤t≤n+1,G=Δ/Λ,H=Φ/Ψ。

    2 基本再生數(shù)

    在傳染病動力學(xué)研究中,基本再生數(shù)是一個十分重要的概念。假設(shè)群體中的所有個體均是易感者,此時一個感染者在平均患病期內(nèi)所感染的人數(shù)就定義為基本再生數(shù)[1,2,14],一般用R0表示。如果R0<1,則一個染病者在其患病期間感染人數(shù)不足一個,疾病將會迅速消亡;相反如果R0>1,則疾病會始終存在并形成地方病。因此為了防止疾病的流行,必須采取相應(yīng)的方法盡量減少R0。通常用R0=1作為決定疾病是否消亡的閾值。

    有效度模型(1)和(2)滿足平衡條件,即所有易感節(jié)點的感染鄰居總數(shù)和所有感染節(jié)點的易感鄰居總數(shù)均等于所有連接易感節(jié)點和感染節(jié)點邊的總數(shù)[11]:

    (3)

    根據(jù)式(3),方程(2)可以改寫成:

    (4)

    在無病平衡點附近,感染節(jié)點的數(shù)量很少,一個易感節(jié)點的感染鄰居數(shù)量多于1的情況不大可能,即對一切i>1,Ssi≈0。則有:

    因此在無病平衡點附近有:

    (5)

    下面采用The next generation matrix方法來計算模型的基本再生數(shù)R0。將(1)和(5)在無病平衡點處線性化,其雅克比矩陣寫成F-V的形式,其中F矩陣表示疾病生成項。V矩陣表示轉(zhuǎn)移項[14]。對于一個固定的度k(1≤k≤M),(1)和(5)可以寫成2(k+1)個方程,但由于Sk0是全易感類,沒有機會接觸染病者,這類方程可以忽略,這樣需要考慮的方程共有2k+1個。

    根據(jù)Jennifer等人的研究[11,13],矩陣F的表達(dá)式為

    (6)

    其中,uk和vk都是2(k+1)×1的向量。uk的第一個元素是kSk0,其它元素都是0。而vk的前k-1個元素分別是k-1,2(k-2),…,t(k-t),…,k-1,其余元素均為0。這里F是一個M(M+2)階矩陣,對它進(jìn)行初等變換可知道矩陣F的秩是1。

    根據(jù)上面的方程可寫出矩陣V2:

    O.Diekmann等人用譜半徑ρ來表示基本再生數(shù)R0=ρ(FV-1),ρ(A)表示矩陣A的譜半徑,即矩陣A的最大特征值的絕對值[16]。有效度模型的基本再生數(shù)可以表示為[13]:

    (7)

    根據(jù)式(7),兩層星型網(wǎng)絡(luò)傳染病的有效度模型的基本再生數(shù)可表示為

    3 防控策略

    為了控制疾病擴散,我們應(yīng)該采取一些防控策略來降低基本再生數(shù)。這里我們介紹兩種最常見的方法。第一種是控制感染源,切斷傳播途徑,如禽流感病毒爆發(fā)時,應(yīng)盡可能的減少與家禽的直接接觸;第二種是保護(hù)易感者,對易感者接種疫苗,使其對病毒具有免疫能力。

    3.1 切斷連邊

    假設(shè)兩層星型網(wǎng)絡(luò)中每層邊界點的個數(shù)為n,中心點個數(shù)為1。采取每種方法都會有3種不同的方案。我們先考慮第一種方法——切斷傳播路徑,即斷邊,有下面3類實施方案。

    1)切斷兩個邊界點之間的邊,如圖2a所示。假設(shè)選擇10%的邊斷開,此時模型中有3類節(jié)點,分別是度為1的節(jié)點,度為2的節(jié)點以及度為n+1的中心節(jié)點。根據(jù)(7)式,基本再生數(shù)可表示為

    圖2 切斷傳播路徑的3種方案

    2)切斷邊界點與中心點之間的邊,如圖2b所示。假設(shè)上下各選擇10%的邊斷開,此時模型中有3類節(jié)點,分別是度為1和度為2的的邊界點以及度為(0.9n+1)的中心點。根據(jù)式(7),基本再生數(shù)為

    3)切斷兩個中心點之間的邊,如圖2c所示。此時模型中有兩類節(jié)點,分別是度為2的邊界點以及度為n的中心點。那么基本再生數(shù)為:

    3.2 免疫節(jié)點

    考慮第二種方法防控方法,即免疫節(jié)點,同樣也有三種實施方案。

    圖3 免疫節(jié)點的3種實施方案

    1)隨機選擇一些易感的邊界點,對其進(jìn)行免疫,如圖3a所示。假設(shè)上下各選擇10%的易感邊界點進(jìn)行免疫,則此時模型中有3類節(jié)點,分別是度為1和度為2的邊界點以及度為(0.9n+1)的中心點。此時的基本再生數(shù)為:

    2)免疫其中一個中心點,如圖3b所示。則此時模型中有3類節(jié)點,分別是度為1點和度為2的邊界點以及度為n的中心點。此時的基本再生數(shù):

    3)對兩個中心點都免疫,如圖3c所示。則此時模型中只有度為1的節(jié)點。此時的基本再生數(shù)為:

    4 數(shù)值分析

    我們已經(jīng)計算了各種控制方案的基本再生數(shù),為了揭示模型參數(shù)對傳播過程的影響,從而出最佳的疾病控制方案,下面我們利用Matlab軟件模擬傳播模型與斷邊和免疫的3種方案的基本再生數(shù)R0隨感染率β、S20的初值、邊界點規(guī)模n變化的曲線圖。下面假設(shè)γ=1,初始時刻有一個中心節(jié)點被感染。

    4.1 切斷連邊的效果分析

    圖4表示在不同的參數(shù)條件下,切斷連邊對基本再生數(shù)的影響。左圖表明基本再生數(shù)與感染率是成正比,因此降低感染率能有效地減少疾病的傳播;中圖表明基本再生數(shù)與S20的初值變化是成反比的關(guān)系,說明網(wǎng)絡(luò)中易感節(jié)點越多疾病越不容易傳播;右圖表明基本再生數(shù)與邊界點規(guī)模是成正比,因此減少邊界點的規(guī)模能有效降低疾病的傳播??偠灾袛嗤緩降?種方案的曲線都能降低基本再生數(shù),但第1種方案(切斷10%邊界連邊)效果不理想,而若只切斷中心節(jié)點之間的一條邊,效果也一般,第2種方案(切斷10%邊界與中心的連邊)效果最佳。

    圖4 在不同的情況下,切斷連邊對基本再生數(shù)R0的影響Fig.4 The effects of disconnecting links on the basic reproduction number R0in case of different situations

    4.2 免疫節(jié)點的效果分析

    圖4表示在不同的參數(shù)條件下,免疫節(jié)點對基本再生數(shù)的影響。結(jié)果發(fā)現(xiàn)免疫的3種方案都能有效地降低基本再生數(shù)。其中第2種方案(免疫一個中心節(jié)點)效果最差,第3種方案(免疫兩個中心節(jié)點)效果最佳,基本能達(dá)到將基本再生數(shù)控制在0的附近,使疾病徹底消除。

    圖5 在不同的情況下,免疫節(jié)點對基本再生數(shù)R0的影響Fig.5 The effects of vaccinating nodes on the basic reproduction numberR0in case of different situations

    從圖5也可以看出免疫10%的邊界點比免疫其中一個中心節(jié)點更有效。下面我們通過改變接受免疫的邊界點數(shù)量,分析免疫邊界點和一個中心節(jié)點的等效性,結(jié)果如圖6所示。我們發(fā)現(xiàn)將接受免疫的邊界點數(shù)量設(shè)定為1%時,兩種方案的基本再生數(shù)曲線基本重合;當(dāng)免疫的邊界點越多,控制效果就會越好。由此可知,當(dāng)網(wǎng)絡(luò)規(guī)模很大時,免疫中心節(jié)點即有效又經(jīng)濟(jì),而且規(guī)模越大,效果越明顯。

    β=0.03,S10(0)=0.01n,S20(0)=0.5n圖6 免疫邊界點與中心點的比較

    5 結(jié)論

    本文運用復(fù)雜網(wǎng)絡(luò)和傳染病動力學(xué)理論研究了傳染病在兩層星型網(wǎng)絡(luò)上傳播動力性。首先利用SIS有效度模型將個體根據(jù)其鄰居的狀態(tài)進(jìn)行分類,從而細(xì)化了傳播過程。接著根據(jù)The next generation matrix方法分析了模型的基本再生數(shù)。為了降低基本再生數(shù)從而控制疾病的傳播,我們提出了切斷傳播途徑和免疫個體兩種方法,根據(jù)不同情況又將其各自細(xì)分為3種不同的實施方案,并求出相應(yīng)的基本再生數(shù)。最后我們利用Matlab數(shù)值分析,獲得了兩層星型網(wǎng)絡(luò)上參數(shù)對疾病的傳播影響以及有效的控制方案。

    研究表明基本再生數(shù)與感染率和邊界點規(guī)模成正比,降低感染率或邊界點規(guī)模都能有效的減少疾病的傳播;基本再生數(shù)與易感者的初值規(guī)模成反比,網(wǎng)絡(luò)中初始時刻感染者越多會使疾病迅速擴散。這說明兩層星型網(wǎng)絡(luò)的特殊結(jié)構(gòu)對疫情傳播有直接影響,那些與感染者有直接接觸的易感者會很容易被感染,而且網(wǎng)絡(luò)節(jié)點數(shù)越多疾病傳播越迅猛。我們還發(fā)現(xiàn)在兩層星型模型中,我們提出的六種防控方案都能降低基本再生數(shù),但是通過比較可以看出切斷部分邊界點與中心點的接觸途徑以及免疫兩個中心點對疾病控制最有效。因此我們在疾病防控時要重點關(guān)注那些活動頻繁和度大的節(jié)點,將這些點隔離或者免疫能大大降低傳染病對社會財產(chǎn)和人們生命安全的威脅。

    本文結(jié)果是基于靜態(tài)網(wǎng)絡(luò)上的疾病傳播,適用于傳播速度相對于群體結(jié)構(gòu)變化快得多的情況。我們設(shè)定的兩層星型模型中心節(jié)點上下各一個,但在現(xiàn)實中,一個群體內(nèi)的中心節(jié)點或許有多個,這樣網(wǎng)絡(luò)結(jié)構(gòu)就會變得復(fù)雜多變。如何將我們的理論研究貼近實際情況,解決現(xiàn)實問題是我們今后研究的動力和方向。

    [1]馬知恩,周義倉,王穩(wěn)地,等.傳染病動力學(xué)的數(shù)學(xué)建模與研究[M].北京:科學(xué)出版社,2004.

    [2]Diekmann O, Heesterbeek J A P. Mathematical epidemiology of infectious diseases: model building, analysis, and interpretation [M]. New York: Wiley, 2000.

    [3]Read J M, Keeling M J. Disease evolution on networks: the role of contact structure [J]. Proc R Soc Lond B, 2003, 270: 699-708.

    [4]汪小帆,李翔,陳關(guān)榮.網(wǎng)絡(luò)科學(xué)導(dǎo)論[M].北京:高等教育出版社,2012.

    [5]Newmann M E J. The structure and function of complex networks [J]. SIAM Rev, 2003, 45: 167-256.

    [6]Gao J, Buldyrev S V, Stanley H E, Havlin S. Networks formed from interdependent networks [J]. Nature Phys, 2012, 8: 40-48.

    [7]Cohen R, Havlin S, et al. Efficient immunization strategies for computer networks and populations [J]. Phys Rev Lett, 2003, 91(24): 247901.

    [8]王琴,祝光湖,傅新楚.有向網(wǎng)絡(luò)上流行病閾值比較和免疫分析[J].復(fù)雜系統(tǒng)與復(fù)雜性科學(xué), 2012, 9(4): 26-33.

    Wang Qin, Zhu Guanghu, Fu Xinchu. Comparison of epidemic thresholds on directed networks and immunization analysis [J]. Complex Systems and Complexity Science, 2002, 9(4): 26-33.

    [9]Fu X C, Michael S, et al. Epidemic dynamics on scale-free networks with piecewise linear infectivity and immunization [J]. Phys Rev E, 2008, 77: 036113.

    [10] 王陽陽,闞佳倩,張海峰.個體的多策略行為不利于控制疾病傳播[J].復(fù)雜系統(tǒng)與復(fù)雜性科學(xué), 2015, 12(3): 14-18.

    Wang Yangyang, Kan Jiaqian, Zhang Haifeng. Individuals' multi-strategy behavior is not conducive to control the spread of epidemic [J]. Complex Systems and Complexity Science, 2015, 12 (3):14-18.

    [11] Jennifer L, Ma J, P van den D, et al. Effective degree network disease models [J]. Mathematical Biology, 2011, 62:143-164.

    [12] Ma J, Pvan den D, Frederick H W. Effective degree household network disease model [J]. Mathematical Biology, 2013, 66: 75-94.

    [13] 鄭國慶,唐清干,祝光湖.帶接種免疫的網(wǎng)絡(luò)傳染病的有效度模型[J].數(shù)學(xué)的實踐與認(rèn)識,2015, 45(15): 315-322.

    Zheng Guoqing, Tang Qingan, Zhu Guanghu. An effective degree epidemic model with vaccination on networks [J]. Mathematics in Practice and Theory, 2015, 45(15): 315-322.

    [14] Pvan den D, James W. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission [J]. Mathematical Biosciences, 2002, 180: 29-48.

    猜你喜歡
    星型邊界點中心點
    道路空間特征與測量距離相結(jié)合的LiDAR道路邊界點提取算法
    層次化點云邊界快速精確提取方法研究
    增加斷電連鎖 減少絞傷風(fēng)險
    Scratch 3.9更新了什么?
    電腦報(2020年12期)2020-06-30 19:56:42
    如何設(shè)置造型中心點?
    電腦報(2019年4期)2019-09-10 07:22:44
    金銀點綴
    漢字藝術(shù)結(jié)構(gòu)解析(二)中心點處筆畫應(yīng)緊奏
    D-π-A星型分子的合成及非線性光學(xué)性質(zhì)
    尋找視覺中心點
    大眾攝影(2015年9期)2015-09-06 17:05:41
    一種去除掛網(wǎng)圖像鋸齒的方法及裝置
    電腦與電信(2014年6期)2014-03-22 13:21:06
    亚洲精品色激情综合| 亚洲精品第二区| 国产亚洲午夜精品一区二区久久| 日韩制服骚丝袜av| 国产精品一区二区性色av| 日韩人妻高清精品专区| 蜜臀久久99精品久久宅男| 欧美bdsm另类| 中国三级夫妇交换| 亚洲av福利一区| 国产视频首页在线观看| 国产亚洲91精品色在线| 女人精品久久久久毛片| xxx大片免费视频| 中国三级夫妇交换| 亚洲欧美成人综合另类久久久| 欧美bdsm另类| 亚洲内射少妇av| 久久6这里有精品| 国产淫片久久久久久久久| 韩国高清视频一区二区三区| 久久久久久久亚洲中文字幕| 久久久精品免费免费高清| 亚洲国产日韩一区二区| 国产成人精品福利久久| 中文字幕人妻丝袜制服| 国产黄频视频在线观看| 极品人妻少妇av视频| 黑人猛操日本美女一级片| 一级a做视频免费观看| 久久综合国产亚洲精品| 美女福利国产在线| 亚洲高清免费不卡视频| a 毛片基地| 婷婷色av中文字幕| 一个人看视频在线观看www免费| 街头女战士在线观看网站| 亚洲av电影在线观看一区二区三区| 亚洲欧美中文字幕日韩二区| 日本欧美视频一区| 国产午夜精品一二区理论片| 欧美亚洲 丝袜 人妻 在线| √禁漫天堂资源中文www| 久久久久网色| 97超碰精品成人国产| 亚洲av国产av综合av卡| 午夜久久久在线观看| 简卡轻食公司| 久久这里有精品视频免费| 永久免费av网站大全| 国产成人精品婷婷| 日韩亚洲欧美综合| 免费看av在线观看网站| 亚洲av综合色区一区| 极品少妇高潮喷水抽搐| 中文字幕免费在线视频6| 久久久久久久国产电影| 人妻系列 视频| 久久精品国产鲁丝片午夜精品| 欧美一级a爱片免费观看看| 丝袜在线中文字幕| 欧美成人午夜免费资源| 中文在线观看免费www的网站| 国产精品国产av在线观看| 久久久久久久久久成人| 内射极品少妇av片p| 日韩,欧美,国产一区二区三区| 在线观看人妻少妇| 亚洲,欧美,日韩| 高清黄色对白视频在线免费看 | 国产av码专区亚洲av| 亚洲激情五月婷婷啪啪| 亚洲av男天堂| 日日撸夜夜添| 人人澡人人妻人| 亚洲国产最新在线播放| a级毛色黄片| 久久精品国产a三级三级三级| 久久这里有精品视频免费| 精品酒店卫生间| 全区人妻精品视频| 久久99蜜桃精品久久| 亚洲美女黄色视频免费看| 女人精品久久久久毛片| 国产精品女同一区二区软件| 国产综合精华液| 国产成人免费无遮挡视频| 国产精品久久久久久av不卡| 亚洲欧美日韩卡通动漫| 热re99久久国产66热| 99久国产av精品国产电影| 一本大道久久a久久精品| 熟女人妻精品中文字幕| 国产精品一二三区在线看| 乱系列少妇在线播放| 国产淫片久久久久久久久| 国产精品久久久久久精品电影小说| 久久精品夜色国产| 少妇熟女欧美另类| av女优亚洲男人天堂| 亚洲欧美日韩另类电影网站| 国产亚洲一区二区精品| 秋霞在线观看毛片| 极品教师在线视频| 欧美97在线视频| 国产精品国产av在线观看| 国产成人午夜福利电影在线观看| 国产精品人妻久久久久久| av在线老鸭窝| 欧美老熟妇乱子伦牲交| 亚洲av二区三区四区| 成人毛片a级毛片在线播放| 男女边摸边吃奶| 五月伊人婷婷丁香| 七月丁香在线播放| 99久久中文字幕三级久久日本| 99久久精品热视频| 色吧在线观看| 国产精品免费大片| 免费观看无遮挡的男女| 好男人视频免费观看在线| 亚洲色图综合在线观看| 亚洲国产色片| 亚洲av综合色区一区| 亚洲国产精品国产精品| 亚洲丝袜综合中文字幕| 久热这里只有精品99| 少妇的逼水好多| av天堂中文字幕网| 国产精品一区二区性色av| 草草在线视频免费看| 九色成人免费人妻av| 啦啦啦中文免费视频观看日本| 国产在视频线精品| 亚洲成人一二三区av| 国产成人精品一,二区| www.av在线官网国产| 免费人妻精品一区二区三区视频| 最近最新中文字幕免费大全7| 国产综合精华液| 国产视频内射| 80岁老熟妇乱子伦牲交| 黄色一级大片看看| 亚洲第一区二区三区不卡| 国产极品粉嫩免费观看在线 | 能在线免费看毛片的网站| 一级av片app| 久久精品国产亚洲av涩爱| 日本与韩国留学比较| 国产亚洲91精品色在线| 一级毛片久久久久久久久女| 在现免费观看毛片| 亚洲欧美一区二区三区国产| 国产亚洲一区二区精品| 精品视频人人做人人爽| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美精品专区久久| 亚洲av男天堂| 免费不卡的大黄色大毛片视频在线观看| 日韩成人av中文字幕在线观看| 久久免费观看电影| 人人妻人人澡人人看| 99久久综合免费| 亚洲精品一区蜜桃| 十八禁网站网址无遮挡 | 建设人人有责人人尽责人人享有的| 成人综合一区亚洲| 三级国产精品欧美在线观看| 精品久久久久久久久av| 性色av一级| 在线观看免费视频网站a站| 亚洲av欧美aⅴ国产| 国产精品欧美亚洲77777| 国产黄色视频一区二区在线观看| 日日摸夜夜添夜夜添av毛片| 日韩欧美 国产精品| av一本久久久久| 国产日韩欧美视频二区| 久久久欧美国产精品| 男人舔奶头视频| 国产亚洲午夜精品一区二区久久| 午夜福利影视在线免费观看| 精品亚洲成国产av| 国产精品一区二区在线观看99| 国产黄色免费在线视频| 免费播放大片免费观看视频在线观看| 插阴视频在线观看视频| 91成人精品电影| 老司机影院毛片| 国产亚洲91精品色在线| 国产片特级美女逼逼视频| 美女中出高潮动态图| 观看免费一级毛片| 亚洲av成人精品一区久久| 秋霞伦理黄片| av.在线天堂| 国产精品秋霞免费鲁丝片| 精品久久久精品久久久| 欧美精品高潮呻吟av久久| 国内少妇人妻偷人精品xxx网站| 精品午夜福利在线看| 亚洲欧美中文字幕日韩二区| h视频一区二区三区| 亚洲av国产av综合av卡| 午夜激情久久久久久久| 春色校园在线视频观看| 久久国产精品男人的天堂亚洲 | 一区二区三区四区激情视频| 少妇人妻精品综合一区二区| 日本欧美视频一区| 亚洲av日韩在线播放| 国产亚洲欧美精品永久| 国产免费一区二区三区四区乱码| 国产亚洲av片在线观看秒播厂| 国产av码专区亚洲av| 777米奇影视久久| 欧美精品亚洲一区二区| 欧美老熟妇乱子伦牲交| 老司机影院毛片| 特大巨黑吊av在线直播| 久久精品熟女亚洲av麻豆精品| 亚洲av日韩在线播放| 国产欧美日韩综合在线一区二区 | 国产一区亚洲一区在线观看| 国产成人a∨麻豆精品| 亚洲精品成人av观看孕妇| 国产毛片在线视频| 日日啪夜夜撸| 9色porny在线观看| 亚洲精品456在线播放app| 欧美高清成人免费视频www| 婷婷色麻豆天堂久久| 一级毛片电影观看| 高清欧美精品videossex| 亚洲成人一二三区av| 秋霞伦理黄片| 国产精品国产三级专区第一集| 一级二级三级毛片免费看| 我要看黄色一级片免费的| 国产精品蜜桃在线观看| a级片在线免费高清观看视频| 九九在线视频观看精品| 亚洲av中文av极速乱| 成人18禁高潮啪啪吃奶动态图 | 日韩精品有码人妻一区| 18禁在线无遮挡免费观看视频| 国产色爽女视频免费观看| 日本免费在线观看一区| 少妇的逼好多水| 成人二区视频| 麻豆成人av视频| 两个人的视频大全免费| 国产男女内射视频| 国产高清有码在线观看视频| 免费黄频网站在线观看国产| 高清黄色对白视频在线免费看 | 最近的中文字幕免费完整| 91成人精品电影| 久久国产精品大桥未久av | 三上悠亚av全集在线观看 | 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 我的老师免费观看完整版| 少妇人妻精品综合一区二区| 人体艺术视频欧美日本| 国产亚洲5aaaaa淫片| 久久6这里有精品| 国产精品久久久久久精品电影小说| 婷婷色av中文字幕| 不卡视频在线观看欧美| freevideosex欧美| 多毛熟女@视频| 最近的中文字幕免费完整| 女性生殖器流出的白浆| 少妇的逼好多水| 国产精品伦人一区二区| av福利片在线观看| 国产男女内射视频| 国产精品久久久久久久久免| 九九在线视频观看精品| 一级a做视频免费观看| 免费看日本二区| 插逼视频在线观看| 男女国产视频网站| 亚洲美女搞黄在线观看| 国内少妇人妻偷人精品xxx网站| 国产高清有码在线观看视频| videos熟女内射| 2022亚洲国产成人精品| 亚洲va在线va天堂va国产| 久久国产精品男人的天堂亚洲 | 三上悠亚av全集在线观看 | 一二三四中文在线观看免费高清| 免费在线观看成人毛片| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| a级毛片免费高清观看在线播放| 国产无遮挡羞羞视频在线观看| 亚洲无线观看免费| 日本wwww免费看| 黄色视频在线播放观看不卡| videos熟女内射| 国产成人精品无人区| 成人影院久久| 超碰97精品在线观看| 十八禁网站网址无遮挡 | 中文乱码字字幕精品一区二区三区| 欧美成人午夜免费资源| 国产成人精品福利久久| 国产精品秋霞免费鲁丝片| 亚洲国产精品成人久久小说| 国产色婷婷99| 国产有黄有色有爽视频| 亚洲av福利一区| 久久精品久久久久久久性| 美女脱内裤让男人舔精品视频| 久久午夜综合久久蜜桃| 亚洲精品一区蜜桃| 亚洲av中文av极速乱| 一本久久精品| 人妻一区二区av| 日日撸夜夜添| 久久久久久久国产电影| 国产av精品麻豆| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区国产| 成人亚洲精品一区在线观看| 91久久精品国产一区二区三区| 欧美 日韩 精品 国产| 国产男女超爽视频在线观看| 免费久久久久久久精品成人欧美视频 | 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 亚洲人成网站在线观看播放| 日本wwww免费看| 2021少妇久久久久久久久久久| kizo精华| 少妇丰满av| 精品国产国语对白av| 国产一区有黄有色的免费视频| 一区二区三区精品91| 国模一区二区三区四区视频| a级一级毛片免费在线观看| 午夜av观看不卡| 国产免费视频播放在线视频| 久久精品国产鲁丝片午夜精品| 纵有疾风起免费观看全集完整版| 2018国产大陆天天弄谢| 成人二区视频| 91久久精品电影网| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 久久久久久久久大av| 九九爱精品视频在线观看| 婷婷色综合大香蕉| 在线亚洲精品国产二区图片欧美 | 久久影院123| 欧美高清成人免费视频www| av一本久久久久| 国产精品国产av在线观看| 成人特级av手机在线观看| 80岁老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡 | 色网站视频免费| 国产av国产精品国产| 国产av一区二区精品久久| 亚洲美女视频黄频| 中文字幕精品免费在线观看视频 | 国产欧美另类精品又又久久亚洲欧美| 乱系列少妇在线播放| 99热这里只有是精品50| 中文欧美无线码| 久久6这里有精品| 国产伦精品一区二区三区四那| 街头女战士在线观看网站| 一本大道久久a久久精品| 国产在线视频一区二区| 一本—道久久a久久精品蜜桃钙片| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 欧美 亚洲 国产 日韩一| 免费观看av网站的网址| 国产精品久久久久成人av| 在线天堂最新版资源| 春色校园在线视频观看| 交换朋友夫妻互换小说| 麻豆乱淫一区二区| 五月伊人婷婷丁香| 亚洲精品一区蜜桃| 五月玫瑰六月丁香| av视频免费观看在线观看| 成人特级av手机在线观看| 赤兔流量卡办理| 国产精品一区二区在线不卡| 国产在线视频一区二区| 王馨瑶露胸无遮挡在线观看| 亚洲欧美精品自产自拍| a 毛片基地| 亚洲精品日韩在线中文字幕| a级毛色黄片| 少妇人妻一区二区三区视频| 蜜桃久久精品国产亚洲av| 久久久亚洲精品成人影院| 色哟哟·www| 亚洲内射少妇av| av线在线观看网站| 国产亚洲5aaaaa淫片| 大片免费播放器 马上看| 乱码一卡2卡4卡精品| 久久99一区二区三区| 欧美日韩国产mv在线观看视频| 国产精品人妻久久久久久| 一区在线观看完整版| 久久国产亚洲av麻豆专区| 亚洲精品第二区| 高清视频免费观看一区二区| 九色成人免费人妻av| 99热这里只有精品一区| 91久久精品电影网| 久久6这里有精品| 国产成人精品一,二区| 日韩精品有码人妻一区| 国产精品伦人一区二区| 国产精品偷伦视频观看了| 久久久久网色| 久久久国产精品麻豆| 91在线精品国自产拍蜜月| 97超碰精品成人国产| 伦理电影免费视频| 久热这里只有精品99| 国产成人精品婷婷| 日日啪夜夜爽| 国产熟女午夜一区二区三区 | 国产中年淑女户外野战色| 久久久久久久久久久免费av| 国产精品三级大全| 日日啪夜夜爽| 欧美亚洲 丝袜 人妻 在线| 插阴视频在线观看视频| 成年人免费黄色播放视频 | 国产美女午夜福利| 99精国产麻豆久久婷婷| 在线观看免费高清a一片| 不卡视频在线观看欧美| 男女国产视频网站| a级一级毛片免费在线观看| 国产男人的电影天堂91| 久久久久久久久久久丰满| 最新中文字幕久久久久| 国产精品蜜桃在线观看| 精品国产乱码久久久久久小说| 欧美 日韩 精品 国产| 亚洲精品亚洲一区二区| 春色校园在线视频观看| 国产男女超爽视频在线观看| 嫩草影院新地址| kizo精华| 91久久精品国产一区二区三区| 少妇的逼水好多| 色哟哟·www| 大又大粗又爽又黄少妇毛片口| 色婷婷av一区二区三区视频| 国产在线一区二区三区精| 一区二区av电影网| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 日日啪夜夜撸| 丰满乱子伦码专区| 欧美97在线视频| 午夜激情久久久久久久| av黄色大香蕉| 九九久久精品国产亚洲av麻豆| 亚洲美女搞黄在线观看| 观看av在线不卡| 国产精品久久久久久精品古装| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 久久久午夜欧美精品| 国产精品国产三级国产专区5o| 男女边摸边吃奶| 免费黄色在线免费观看| 伊人亚洲综合成人网| 老司机亚洲免费影院| .国产精品久久| 啦啦啦视频在线资源免费观看| 亚洲在久久综合| 久久综合国产亚洲精品| 日本黄色日本黄色录像| 我要看日韩黄色一级片| 日韩成人伦理影院| 狂野欧美激情性xxxx在线观看| 国产成人精品一,二区| 亚洲国产av新网站| 亚洲精品国产av蜜桃| 成人国产麻豆网| 国产深夜福利视频在线观看| 视频区图区小说| 色视频www国产| 一级,二级,三级黄色视频| 美女cb高潮喷水在线观看| 精品酒店卫生间| 最近中文字幕高清免费大全6| 亚洲在久久综合| 大码成人一级视频| 久久狼人影院| 精品久久久久久久久av| videos熟女内射| 亚州av有码| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 精品人妻一区二区三区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一区蜜桃| 一级毛片 在线播放| 国产高清三级在线| 一级毛片 在线播放| 乱码一卡2卡4卡精品| 大码成人一级视频| 在线观看免费视频网站a站| 欧美国产精品一级二级三级 | 久久久久久久久大av| 日本91视频免费播放| 亚洲一级一片aⅴ在线观看| 老女人水多毛片| 精品一区二区免费观看| 免费看av在线观看网站| 亚洲国产欧美日韩在线播放 | 日本与韩国留学比较| 女性被躁到高潮视频| 人人妻人人澡人人爽人人夜夜| 高清不卡的av网站| 欧美人与善性xxx| 亚洲精品色激情综合| 欧美成人精品欧美一级黄| 少妇丰满av| av在线app专区| av免费在线看不卡| 又爽又黄a免费视频| 精品午夜福利在线看| 久久狼人影院| 国产视频首页在线观看| 午夜福利,免费看| 黄色怎么调成土黄色| 欧美一级a爱片免费观看看| 91成人精品电影| 色吧在线观看| 久久久国产欧美日韩av| 十分钟在线观看高清视频www | 精品少妇黑人巨大在线播放| 久久97久久精品| 国产高清不卡午夜福利| 国产亚洲精品久久久com| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 国产伦在线观看视频一区| 一本—道久久a久久精品蜜桃钙片| 亚洲第一av免费看| 国产亚洲最大av| 欧美一级a爱片免费观看看| 高清av免费在线| 热re99久久精品国产66热6| 美女大奶头黄色视频| 下体分泌物呈黄色| 成人毛片60女人毛片免费| 99热全是精品| 亚洲av中文av极速乱| 丝袜在线中文字幕| 哪个播放器可以免费观看大片| av黄色大香蕉| 精品久久国产蜜桃| 久久亚洲国产成人精品v| 亚洲,一卡二卡三卡| 最近中文字幕2019免费版| 国产精品熟女久久久久浪| 久久久亚洲精品成人影院| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 精品视频人人做人人爽| 欧美性感艳星| 老司机影院毛片| 老司机影院成人| 伊人亚洲综合成人网| 国产美女午夜福利| 一级,二级,三级黄色视频| 亚洲成人一二三区av| 色婷婷av一区二区三区视频| a级一级毛片免费在线观看| 久久久久视频综合| 免费看不卡的av| 午夜福利视频精品| 在线天堂最新版资源| 国产乱来视频区| 国产黄片美女视频| 久久久久人妻精品一区果冻| 久久鲁丝午夜福利片| 亚洲情色 制服丝袜| 国产在线视频一区二区| 插阴视频在线观看视频| 天天躁夜夜躁狠狠久久av| 免费看日本二区| 最新中文字幕久久久久| 国产一区有黄有色的免费视频| 一边亲一边摸免费视频| 精品一品国产午夜福利视频| av播播在线观看一区| 国产黄色视频一区二区在线观看| 自拍偷自拍亚洲精品老妇| 精品99又大又爽又粗少妇毛片| 91久久精品电影网| 国产精品久久久久久精品古装| 亚洲国产av新网站| 国产淫语在线视频| 校园人妻丝袜中文字幕| 插逼视频在线观看| 久久国产亚洲av麻豆专区|