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

    Ergodic stationary distribution of a stochastic rumor propagation model with general incidence function

    2022-06-29 08:52:52YuhuaiZhang張宇槐andJianjunZhu朱建軍
    Chinese Physics B 2022年6期
    關(guān)鍵詞:建軍

    Yuhuai Zhang(張宇槐) and Jianjun Zhu(朱建軍)

    College of Economics and Management,Nanjing University of Aeronautics and Astronautics,Nanjing 211106,China

    Keywords: rumor propagation model,general incidence function,It?o’s formula,ergodic stationary distribution

    1. Introduction

    Recently,due to the pandemic of COVID-19,various rumors related to COVID-19 have been circulated on the internet. The rumors not only easily cause people’s emotional panic, but also cause other social problems. For example,“Novel coronavirus can be eliminated by using the blow dryer to blow the face and hands for a long time”, “Virus can be killed by drinking bleaching disinfectant” and other pseudoscience or malicious“conspiracy theory”. Although some authoritative media and experts have voiced against the rumors,there still are exponential forwarding and diffusion of rumors about the epidemic on social media,just like a virus with great replication ability, which urgently needs to respond with actions matching it. Thus, it is of great significance to study rumor propagation in online social networks for strengthening the supervision of network public opinion and purifying the cyberspace.

    To minimize the harmfulness of rumors on online social network platforms as much as possible, it is necessary to understand the rules of rumor spreading on online social network platforms.In the past period of time,some mathematical models were established to investigate rumor spreading in the complex network.[10,13,23]In fact, the research of classical rumor spreading model started from 1960s. In 1964 and 1965,in order to compare the differences between the rumor spreading and virus spreading, Daley and Kendal[5,6]proposed the typical DK model of rumor spreading. In their model, they first divided people into three categories:ignorants(people who do not hear the rumor),spreaders(people who spread the rumor)and stiflers (people who know but do not spread the rumor).Since then,the original DK model had been modified and extended continuously, which made it possible to interpret the rumors in more detail. By changing the propagation mechanism, Maki and Thomson[26]further modified the DK model and established a new rumor propagation model-MT model.At the beginning of the 21st century, Zanette first extended the problems of rumor spreading to complex networks, who mainly investigated the rules of rumor propagation on a smallworld network and proposed an epidemic-like model.[34,35]Afterwards, some scholars devoted to studying some rumor propagation models in order to better control the rumors, see Refs. [8,18]. Recently, Tuwairqiet al.[1]proposed a rumor propagation model with incubation mechanism. In their work,some important results,such as the basic reproduction number of the model, the local stability and global asymptotic stability of the rumor-free and rumor equilibrium, were obtained under suitable conditions. Additional literature relating to the rumor propagation models in different patterns of groups,see Refs.[22,28,33]and the references therein.

    To ensure these models can provide the realistic description of the rumor propagation, the incidence rate plays an important role. In 1927, Kermack and McKendrick[19]introduced the specific bilinear incidence rateβSIfor some disease spreading models. However, there were many problems when it was used in describing the disease spreading.To fit the phenomena of disease propagation, Capasso and Serio[4]proposed the saturated incidence rateβSI/(1+αI).After the works, many kinds of incidence rate in various biological models were considered by many scholars, see Refs.[2,9,17,24,25,30–32].From these literature,one can find that models with more general incidence can provide more complex dynamics than those with the bilinear or standard incidence rate.[3,20]Therefore, on social networks, one may consider the following deterministic rumor propagation model with general incidence function:

    In sociology,I(t) represents individuals on online social networks who have never heard rumors or been contacted to forward rumors at timet,H(t) indicates those people that have heard rumors but hesitated to spread the rumors at timet,S(t)refers to those people who choose to believe and forward rumors but not actively create rumors at timet,R(t)means the people who have no interest to the rumors or no longer believe in rumors,and hence will not spread the rumors at timet. The functionf(·)represents the general incidence function. More sociological meanings about the parameters appearing in system(1)are explained in the sequel:

    Λ:the number of entering individuals in the whole group;

    μ:leaving rate of the population;

    δ:conversion rate from hesitating individuals to stiflers;

    γ:conversion rate from hesitating individuals to spreaders;

    β:transmission rate from spreaders to stiflers due to lose interest in spreading rumors or forget rumors;

    θ:portion of those hearing the rumor who will spread the it,andθ ∈[0,1].

    However,in the real world,the individuals’psychological state(e.g.,hesitating,forgetting etc.) and individuals’behavior(e.g., forwarding, comment etc.) sometimes are disturbed by the environmental noise through the process of rumor dissemination. Therefore, in the noisy environment, the deterministic models not only have difficulty in fitting the data but also have strong limitations. Thus,it is meaningful to consider the influence of the social randomness in rumor propagation.More precisely,the effect of uncertainty should be considered in the parameters estimations, such as transmission rates between the four groups of people. Recently, some researchers have devoted in studying the rumor propagation models with stochastic perturbations,see e.g.,Refs.[7,12,14–16,37].From the literature,[38]it is known that the deterministic models can be transformed into stochastic models by introducing the time Markov chain,parameter perturbation,and robustness of positive equilibria to the deterministic models or the environmental noise that is proportional to the variables. This gives us the motivation to investigate the stochastic effect to the model(1).

    Adding the environmental noise into the model(1),in this paper,we consider the following stochastic rumor propagation model:Using the method established by Khasminskii,[21]we provide the sufficient condition to guarantee the existence of an ergodic stationary distribution for system(2). The main result is stated as follows.

    2. Preliminaries

    To obtain the ergodic stationary distribution of system(2),we first recall some definitions and known results in the sequel.

    Let (Ω,?,P) be a complete probability space with filtration{?t}t≥0satisfying the usual conditions(i.e.,it is right continuous and increasing while?0contains all P-null sets),and the values of all parameters are positive. Moreover, we denote

    a∨b:=max{a,b}, a∧b:=min{a,b}for alla,b ∈R,

    andC2,1(Rd×(0,∞);R+)by the function space consisting of all functionsV(x,t):Rd×(0,∞)→R+which are continuously twice differentiable inxand once int. Then we consider initial value problem of thed-dimensional stochastic differential equation:

    3. Existence of an ergodic stationary distribution

    Based on the assumptions (A1) and (A2), the right hand side of system (2) satisfies the local Lipschitz condition. Then, given any initial value (I(0),H(0),S(0),R(0))∈R4+, it is known that (2) admits a unique local solution(I(t),H(t),S(t),R(t))ont ∈[0,τe]for someτe>0. Note thatτemeans the explosion time.[27]To derive the result of Theorem 1,we only need to show thatτe=∞a.s.

    Proof of Theorem 1 Letn0be a large number such thatI(0),H(0),S(0),andR(0)all lie in the interval[1/n0,n0]. For each integern ≥n0,we define the following sequence of stopping time:

    where 1Ωnmeans the indicator function ofΩn. Lettingn →∞,we have

    To verify(H2),let us define theC2-nonnegative functionQ(·):R4+→R by

    Summing up the above estimations,we can derive

    Then it follows from theF-inequality of Eq.(43)thatLV ≤-1 for all(I,H,S,R)∈D5.

    By the first inequality of Eq.(42)and similar argument for Case 1, we can also prove thatLV ≤-1 for all (I,H,S,R)∈D2∪∈D3∪∈D4. Similarly, by theG,J,K-inequalities of Eq.(43)and similar argument for Case 2,we haveLV ≤-1 for all (I,H,S,R)∈D6∪∈D7∪∈D8. Therefore, the claim (50)holds,i.e.,the condition(H2)holds.

    By Proposition 1, we conclude that system (2) admits a unique stationary distribution which exhibits the ergodicity.The proof of Theorem 2 is complete.

    4. Numerical simulations

    To better understand the theoretical results, we carry out some numerical simulations to illustrate the dynamic results of the stochastic rumor model under some parameters. Let us choose

    Then, by using the Milstein scheme,[11]the following discretization system is obtained:

    whereτi,k(i= 1,2,3,4) are four independent the Gaussian random variables withN(0,1)and time increment Δt >0. In addition, we assume that the values of parameters are as follows:

    Next,through the numerical simulations of discretization system(58),we discuss its influence on rumor propagation from the following two aspects: the intensities of noise and key parameters.

    4.1. Impact of the intensities of noise

    From their density functions, one can also see that the solution of system (2) has a stationary distribution. In addition,the numerical results reveal the following interesting fact.

    4.2. Impact of different parameters

    In addition to the stationary distribution of model(2),we are also interested in what parameters will affect rumor propagation. Firstly, we analyze the contact transmission rateα,which is of great significance in communication because it is considered to be a key indicator to reflect whether rumors can be further spread. Here, under the same values of other parameters, we select different transmission rateα=0.1, 0.4,and 0.8 to study their influence on rumor propagation.

    Fig.1. Stationary distribution of system(2)with σ1=σ2=σ3=σ4=0.1.Panels(a)and(b)represent the solution and the density function respectively.

    Fig.2. Stationary distribution of system(2)with σ1=σ2=σ3=σ4=0.02.Panels(a)and(b)represent the solution and the density function respectively.

    From Fig.3,we observe that when the value of the transmission rateαis small, the number of rumor spreadersS(t)will decrease rapidly and is close to zero; when the value of the transmission rateαis big,the number of rumor spreadersS(t) will increase in the fluctuation, and the range of the oscillation will also keep increasing. This shows that reasonably cutting off rumor propagation channels is an effective method to curb rumor propagation on social networks. Especially in the case of emergencies,it is necessary to limit the flow,screen the relevant accounts and other methods to avoid serious negative impact.

    Fig.3. Trajectories simulations of rumor spreaders S(t)of stochastic rumor propagation model(2)is graphed for α =0.1,0.4,and 0.8.

    Moreover,we discuss the influence of the conversion rate from rumor-hesitating individuals to spreaders on rumor propagation. Thus, we change the value ofγwhile keeping all other parameter values unchanged. As shown in Fig. 4, one can also see that as the value ofγincreases,the number of rumor spreaders increases,and the amplitude of oscillation also increases significantly under the interference of social environment noise.This shows that netizens express their opinions on rumors,which lead to rumor-hesitating individuals can not correctly judge the truth or falseness of information, so that the number of people fluctuates at every moment.In the meanwhile,some rumor-hesitating individuals believe and forward false information only based on the comments of some netizens,resulting in the continuous spread of rumors on the internet. Therefore, government departments should publicize the truth of emergencies to the public in a timely manner, eliminate the doubts of rumor-hesitant individuals,reduce the conversion rate from rumor-hesitant individuals to rumor spreaders,and effectively curb the spread of rumors.

    Fig.4. Trajectories simulations of rumor spreaders S(t)of stochastic rumor propagation model(2)is graphed for γ =0.1,0.3,and 0.65.

    Besides the above mentioned,we believe that the scale of rumor spreaders is closely related to individual factors (such as the degree of interest in rumor propagation, and the degree of information memory and forgetting). Therefore, we investigate the influence ofβon the rumor spreaders under the premise of other parameter values unchanged.From Fig.5,we can see that with the increase of the value ofβ,the number of rumor spreaders shows a downward trend, and the amplitude of the oscillation also decreased significantly. This illustrates that whenβis small, rumor spreaders are vulnerable to the comments of netizens and constantly forward rumors to their friends around them; With the increase ofβ,the attraction of rumor propagation decreases,which makes the rumor spreaders lose interest in rumor propagation or forget rumors and have no intention to spread, and reduce the influence of rumor propagation. Therefore, government needs to strengthen science education for netizens,and timely release effective information to clarify the facts when an emergency occurs,then the value ofβwill increase correspondingly,so that the rumor will have less influence.

    Fig.5. Trajectories simulations of rumor spreaders S(t)of stochastic rumor propagation model(2)is graphed for β =0.2,0.4,and 0.6.

    5. A brief conclusion

    In this paper, we propose and consider the dynamics of stochastic rumor propagation model(2)with general incidence function. In fact, we extend the classical rumor propagation model by introducing the random interference,and discuss the influence of random factors on the rumor propagation. More precisely,we show the existence of global positive solution of system (2) and provide some sufficient conditions to guarantee the existence of a unique ergodic stationary distribution.Furthermore,the numerical simulation results show that when?0>1, the noise intensity has a great influence on rumor propagation. This means that when rumors diffuse in social media,the intensity of government control will play a key role in the rumor propagation.In addition,the change of parameter value will significantly affect the number of rumor spreaders,so as to provide decision-making departments with effective methods to control rumor propagation. Therefore,this work is helpful to analyze the rumor spreading model in theory and put forward corresponding measures for controlling rumor propagation. In practice, we can change the behavior of netizens by strengthening government intervention, authoritative media reports,publicity and education to wake up the awareness of netizens and other environmental interference factors,so as to reduce the interaction among ignorants, rumor-hesitant individuals and spreaders.

    In addition to Theorems 1 and 2,we further illustrate the following interesting topics and leave them for our future consideration.

    °Examining the existence of positive periodic solutions of system(2).

    °Investigating the strategy to effectively control the dynamics of system(2).

    °Extending our results to the case with L′evy noise, or other non-autonomous systems.

    Acknowledgments

    Project supported by the Funding for Outstanding Doctoral Dissertation in NUAA (Grant No. BCXJ18-09),the National Natural Science Foundation of China (Grant No. 72071106), and Postgraduate Research & Practice Innovation Program of Jiangsu Province, China (Grant No.KYCX180234).

    猜你喜歡
    建軍
    慶祝建軍95周年
    GENERALIZED CES`ARO OPERATORS ON DIRICHLET-TYPE SPACES*
    Spatio-temporal evolution characteristics and pattern formation of a gas–liquid interfacial AC current argon discharge plasma with a deionized water electrode
    古建軍
    無論等多久
    當代人(2020年4期)2020-05-08 08:18:04
    建軍90周年有感
    中華魂(2017年8期)2017-11-22 12:21:09
    建軍90周年
    綠色中國(2017年15期)2017-01-25 08:55:36
    慶祝中國人民解放軍建軍89周年
    Experimental investigation of velocity fluctuations in a radial diffuser pump*
    Totally laparoscopic Billroth Ⅱ gastrectomy without intracorporeal hand-sewn sutures
    大片电影免费在线观看免费| 美女大奶头黄色视频| 国产成人精品一,二区| 国产极品天堂在线| a级毛片黄视频| 国产精品国产三级专区第一集| 黑人欧美特级aaaaaa片| 美女中出高潮动态图| 99热网站在线观看| 国产精品国产三级国产专区5o| 亚洲在久久综合| 老司机影院成人| 久久久久视频综合| 欧美激情国产日韩精品一区| 黑人欧美特级aaaaaa片| 亚洲精品国产av蜜桃| 黑人高潮一二区| 国产高清三级在线| 欧美国产精品一级二级三级| 久久精品国产a三级三级三级| 精品久久蜜臀av无| 成年av动漫网址| 精品国产乱码久久久久久小说| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| 国产毛片在线视频| 啦啦啦中文免费视频观看日本| 视频在线观看一区二区三区| 午夜影院在线不卡| 国产在视频线精品| xxxhd国产人妻xxx| 亚洲综合色惰| 天堂中文最新版在线下载| 在线亚洲精品国产二区图片欧美| 国产精品人妻久久久影院| 精品人妻熟女毛片av久久网站| 蜜桃在线观看..| 欧美日韩一区二区视频在线观看视频在线| 中文字幕人妻熟女乱码| 亚洲欧洲精品一区二区精品久久久 | 亚洲av电影在线观看一区二区三区| av电影中文网址| 纯流量卡能插随身wifi吗| 久久亚洲国产成人精品v| 在线看a的网站| 日本vs欧美在线观看视频| 久久久久久久精品精品| 免费看av在线观看网站| 97超碰精品成人国产| 亚洲国产精品一区二区三区在线| 各种免费的搞黄视频| 美女福利国产在线| 久久久久久久大尺度免费视频| 亚洲精品一二三| 交换朋友夫妻互换小说| xxx大片免费视频| 久久99蜜桃精品久久| 一级爰片在线观看| 色哟哟·www| 日日摸夜夜添夜夜爱| 人妻系列 视频| 午夜免费观看性视频| 免费观看在线日韩| av片东京热男人的天堂| 国产1区2区3区精品| 建设人人有责人人尽责人人享有的| 久久精品熟女亚洲av麻豆精品| 国产精品久久久久久久电影| 99热网站在线观看| 一本久久精品| 99热全是精品| 午夜精品国产一区二区电影| 少妇高潮的动态图| 人人妻人人澡人人看| 啦啦啦视频在线资源免费观看| 婷婷成人精品国产| 久久鲁丝午夜福利片| 亚洲美女黄色视频免费看| 9色porny在线观看| 蜜桃国产av成人99| 精品酒店卫生间| 少妇被粗大猛烈的视频| 卡戴珊不雅视频在线播放| 亚洲国产成人一精品久久久| 久久精品人人爽人人爽视色| 极品人妻少妇av视频| 18禁在线无遮挡免费观看视频| 亚洲美女搞黄在线观看| 国产高清不卡午夜福利| 五月玫瑰六月丁香| 下体分泌物呈黄色| 中国美白少妇内射xxxbb| 亚洲精品,欧美精品| 在线观看国产h片| 国产精品国产三级国产专区5o| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av网站免费在线观看视频| 777米奇影视久久| 熟女人妻精品中文字幕| 久久亚洲国产成人精品v| 91aial.com中文字幕在线观看| 欧美日韩视频高清一区二区三区二| 国产无遮挡羞羞视频在线观看| 高清视频免费观看一区二区| 99香蕉大伊视频| 丝袜美足系列| 丝袜喷水一区| 纵有疾风起免费观看全集完整版| 激情五月婷婷亚洲| 另类精品久久| 看十八女毛片水多多多| 我的女老师完整版在线观看| 曰老女人黄片| 97超碰精品成人国产| 男人操女人黄网站| 热re99久久国产66热| 成年av动漫网址| 宅男免费午夜| 久久久精品免费免费高清| www.熟女人妻精品国产 | 精品亚洲乱码少妇综合久久| 久久精品aⅴ一区二区三区四区 | 久久婷婷青草| 国产精品国产av在线观看| 亚洲久久久国产精品| 天天躁夜夜躁狠狠躁躁| 在线观看人妻少妇| 在线天堂中文资源库| 视频在线观看一区二区三区| 最新中文字幕久久久久| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久噜噜老黄| 搡老乐熟女国产| 最近的中文字幕免费完整| 国产在线一区二区三区精| 亚洲精品国产av蜜桃| 内地一区二区视频在线| 国产深夜福利视频在线观看| 涩涩av久久男人的天堂| 精品一品国产午夜福利视频| 国产免费一级a男人的天堂| av片东京热男人的天堂| 亚洲 欧美一区二区三区| av卡一久久| 在线观看www视频免费| 久热这里只有精品99| 免费观看在线日韩| 亚洲国产欧美在线一区| 午夜91福利影院| 国产精品久久久久久精品古装| 国产日韩欧美在线精品| 久久精品人人爽人人爽视色| 亚洲精品aⅴ在线观看| 天天躁夜夜躁狠狠久久av| 你懂的网址亚洲精品在线观看| 大香蕉久久成人网| 欧美bdsm另类| 午夜日本视频在线| 国产极品粉嫩免费观看在线| 一区二区三区精品91| 男女啪啪激烈高潮av片| 观看美女的网站| 自线自在国产av| 免费观看在线日韩| 黑人欧美特级aaaaaa片| 伊人久久国产一区二区| 国产精品 国内视频| 欧美日韩视频高清一区二区三区二| 高清毛片免费看| 寂寞人妻少妇视频99o| 亚洲精品国产色婷婷电影| 欧美 日韩 精品 国产| 在线天堂最新版资源| 黑人猛操日本美女一级片| 丝袜在线中文字幕| 寂寞人妻少妇视频99o| 桃花免费在线播放| 国产一区二区三区av在线| 国产69精品久久久久777片| 下体分泌物呈黄色| 国产不卡av网站在线观看| 亚洲熟女精品中文字幕| 国产精品嫩草影院av在线观看| 亚洲精品一二三| 亚洲伊人久久精品综合| 少妇人妻精品综合一区二区| 国产男女超爽视频在线观看| av一本久久久久| 18禁观看日本| 欧美精品人与动牲交sv欧美| 国产成人a∨麻豆精品| 男女无遮挡免费网站观看| av免费观看日本| 成人国语在线视频| 亚洲精品国产av成人精品| 国产精品不卡视频一区二区| 搡女人真爽免费视频火全软件| 水蜜桃什么品种好| 看十八女毛片水多多多| 国产片内射在线| 久久精品国产亚洲av天美| 一级,二级,三级黄色视频| 亚洲欧美日韩卡通动漫| 久久精品熟女亚洲av麻豆精品| 欧美日韩精品成人综合77777| 老司机影院成人| 国产在线视频一区二区| 国产激情久久老熟女| 国产精品久久久久久久久免| 永久免费av网站大全| 一区二区三区四区激情视频| 精品视频人人做人人爽| 免费播放大片免费观看视频在线观看| 免费黄色在线免费观看| 国产高清不卡午夜福利| 视频中文字幕在线观看| 巨乳人妻的诱惑在线观看| 亚洲美女黄色视频免费看| av不卡在线播放| 人妻一区二区av| 中国三级夫妇交换| 中文字幕精品免费在线观看视频 | 日本午夜av视频| av女优亚洲男人天堂| 日本黄色日本黄色录像| 日韩中文字幕视频在线看片| 国产精品蜜桃在线观看| 母亲3免费完整高清在线观看 | 97人妻天天添夜夜摸| 久久久国产精品麻豆| 亚洲精品国产色婷婷电影| av在线app专区| 蜜臀久久99精品久久宅男| 亚洲精品国产av成人精品| 久久久久久久久久成人| 久久人人爽人人爽人人片va| 熟妇人妻不卡中文字幕| 丝袜人妻中文字幕| 欧美人与性动交α欧美软件 | 日本欧美国产在线视频| 黄色一级大片看看| 伦理电影免费视频| 精品亚洲成a人片在线观看| 国产精品无大码| 亚洲国产色片| 国产男女超爽视频在线观看| 观看美女的网站| 精品一品国产午夜福利视频| 七月丁香在线播放| 国产69精品久久久久777片| 久久99一区二区三区| 亚洲精品av麻豆狂野| 午夜免费观看性视频| av片东京热男人的天堂| 纯流量卡能插随身wifi吗| 亚洲精品自拍成人| 99久久精品国产国产毛片| 捣出白浆h1v1| 亚洲精品国产色婷婷电影| 少妇精品久久久久久久| 国产一区二区在线观看日韩| 秋霞在线观看毛片| 大话2 男鬼变身卡| 丰满少妇做爰视频| 亚洲精品美女久久久久99蜜臀 | 成年女人在线观看亚洲视频| 男人添女人高潮全过程视频| 黄色毛片三级朝国网站| 另类精品久久| 欧美最新免费一区二区三区| 熟女电影av网| 桃花免费在线播放| 亚洲精品成人av观看孕妇| 久久精品aⅴ一区二区三区四区 | 97人妻天天添夜夜摸| 黑人巨大精品欧美一区二区蜜桃 | av国产精品久久久久影院| 日韩制服丝袜自拍偷拍| 观看av在线不卡| 欧美成人午夜精品| 伊人亚洲综合成人网| av国产久精品久网站免费入址| 国产又色又爽无遮挡免| 大码成人一级视频| 日韩制服丝袜自拍偷拍| 中文字幕人妻熟女乱码| 99热全是精品| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 观看美女的网站| 日韩视频在线欧美| 日韩中字成人| 最近的中文字幕免费完整| 久久国产精品大桥未久av| xxxhd国产人妻xxx| 人人妻人人爽人人添夜夜欢视频| 欧美国产精品一级二级三级| 精品亚洲成a人片在线观看| 午夜av观看不卡| 美女国产高潮福利片在线看| 黑丝袜美女国产一区| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 亚洲美女视频黄频| 女的被弄到高潮叫床怎么办| 黄片无遮挡物在线观看| 少妇熟女欧美另类| 美女大奶头黄色视频| 亚洲欧美清纯卡通| 国产av码专区亚洲av| 久久97久久精品| 亚洲精品一二三| √禁漫天堂资源中文www| 一个人免费看片子| 精品国产一区二区久久| 亚洲精品国产av蜜桃| 国产精品一国产av| 狂野欧美激情性bbbbbb| 在线看a的网站| 最近手机中文字幕大全| 久久久久久久久久人人人人人人| 大片免费播放器 马上看| 精品人妻熟女毛片av久久网站| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 精品久久国产蜜桃| 亚洲丝袜综合中文字幕| 欧美精品一区二区免费开放| 亚洲,一卡二卡三卡| www.色视频.com| 一级,二级,三级黄色视频| 女人精品久久久久毛片| 卡戴珊不雅视频在线播放| 国产在线一区二区三区精| 亚洲欧美精品自产自拍| 天堂8中文在线网| 久久久久久人人人人人| 亚洲成av片中文字幕在线观看 | 日韩熟女老妇一区二区性免费视频| 搡女人真爽免费视频火全软件| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 九九爱精品视频在线观看| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 精品一区二区三区视频在线| 亚洲国产日韩一区二区| 国产日韩一区二区三区精品不卡| 亚洲欧美成人精品一区二区| 久久ye,这里只有精品| 亚洲国产看品久久| 亚洲国产欧美在线一区| 国产xxxxx性猛交| 亚洲精品乱久久久久久| 欧美精品高潮呻吟av久久| 七月丁香在线播放| 亚洲图色成人| 人成视频在线观看免费观看| 男人舔女人的私密视频| 亚洲丝袜综合中文字幕| 91午夜精品亚洲一区二区三区| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 国产成人精品无人区| 亚洲 欧美一区二区三区| 在线观看免费视频网站a站| 99热网站在线观看| 国产成人a∨麻豆精品| 性色av一级| 午夜免费观看性视频| 亚洲精品久久久久久婷婷小说| 又粗又硬又长又爽又黄的视频| 在线免费观看不下载黄p国产| 女人久久www免费人成看片| 亚洲精品久久久久久婷婷小说| 美女福利国产在线| 精品视频人人做人人爽| av播播在线观看一区| 性色avwww在线观看| 日韩在线高清观看一区二区三区| 成人国产麻豆网| 日本欧美国产在线视频| 不卡视频在线观看欧美| 亚洲综合色网址| 另类精品久久| 午夜福利网站1000一区二区三区| 天美传媒精品一区二区| 欧美精品高潮呻吟av久久| av福利片在线| 国产精品一区www在线观看| 不卡视频在线观看欧美| 国产极品天堂在线| 九色亚洲精品在线播放| 日产精品乱码卡一卡2卡三| 国产69精品久久久久777片| 国产精品成人在线| 精品99又大又爽又粗少妇毛片| 国产极品粉嫩免费观看在线| 乱码一卡2卡4卡精品| 99久久中文字幕三级久久日本| 少妇精品久久久久久久| 一边亲一边摸免费视频| 80岁老熟妇乱子伦牲交| 人人妻人人添人人爽欧美一区卜| 精品酒店卫生间| 亚洲高清免费不卡视频| 免费女性裸体啪啪无遮挡网站| 国产老妇伦熟女老妇高清| 最近最新中文字幕大全免费视频 | videosex国产| 性色avwww在线观看| 永久网站在线| 亚洲国产精品成人久久小说| 大香蕉久久成人网| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 在线观看一区二区三区激情| 国产精品欧美亚洲77777| 亚洲欧美成人精品一区二区| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说| 国产亚洲精品久久久com| 最近最新中文字幕大全免费视频 | 亚洲天堂av无毛| 在线观看www视频免费| 波多野结衣一区麻豆| 久久人人97超碰香蕉20202| 日韩制服骚丝袜av| 国产色爽女视频免费观看| 中文精品一卡2卡3卡4更新| 亚洲国产毛片av蜜桃av| 777米奇影视久久| 赤兔流量卡办理| 久久久久人妻精品一区果冻| 中国美白少妇内射xxxbb| 亚洲色图综合在线观看| 亚洲成色77777| 精品卡一卡二卡四卡免费| 成年美女黄网站色视频大全免费| 国产精品嫩草影院av在线观看| 美女大奶头黄色视频| 韩国高清视频一区二区三区| 欧美日韩精品成人综合77777| 亚洲av电影在线观看一区二区三区| 精品一品国产午夜福利视频| 中文字幕最新亚洲高清| 日韩,欧美,国产一区二区三区| 91在线精品国自产拍蜜月| 国产一区二区三区综合在线观看 | 欧美xxⅹ黑人| 欧美xxxx性猛交bbbb| 1024视频免费在线观看| 亚洲国产成人一精品久久久| av片东京热男人的天堂| 亚洲丝袜综合中文字幕| 校园人妻丝袜中文字幕| 国产高清不卡午夜福利| 亚洲国产色片| 成人黄色视频免费在线看| 晚上一个人看的免费电影| 日韩不卡一区二区三区视频在线| 国产色婷婷99| 国产高清国产精品国产三级| 九草在线视频观看| 春色校园在线视频观看| 2022亚洲国产成人精品| 精品99又大又爽又粗少妇毛片| 久久久久视频综合| 久久毛片免费看一区二区三区| 国产成人av激情在线播放| 捣出白浆h1v1| 人妻人人澡人人爽人人| 五月开心婷婷网| 国产男人的电影天堂91| xxx大片免费视频| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| 两个人看的免费小视频| 久久久久精品人妻al黑| 国产成人a∨麻豆精品| 狠狠婷婷综合久久久久久88av| 五月天丁香电影| 日韩av在线免费看完整版不卡| 三级国产精品片| 国产麻豆69| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 尾随美女入室| 亚洲,一卡二卡三卡| 久久女婷五月综合色啪小说| 水蜜桃什么品种好| 制服诱惑二区| 十八禁高潮呻吟视频| 97在线视频观看| 亚洲,欧美精品.| 成人18禁高潮啪啪吃奶动态图| 亚洲中文av在线| 2021少妇久久久久久久久久久| 亚洲av免费高清在线观看| 久久久亚洲精品成人影院| 妹子高潮喷水视频| 国产成人免费观看mmmm| 欧美97在线视频| 九色成人免费人妻av| 捣出白浆h1v1| 国内精品宾馆在线| 亚洲色图 男人天堂 中文字幕 | 日本免费在线观看一区| 欧美成人午夜精品| 国产国语露脸激情在线看| av在线老鸭窝| 国产成人av激情在线播放| 成人亚洲精品一区在线观看| 九色亚洲精品在线播放| 国产熟女午夜一区二区三区| 在线观看人妻少妇| 黑人欧美特级aaaaaa片| 免费高清在线观看日韩| 少妇熟女欧美另类| 日韩精品有码人妻一区| 免费人成在线观看视频色| 日韩成人av中文字幕在线观看| 国产精品一区二区在线观看99| 国产黄色视频一区二区在线观看| 午夜av观看不卡| 午夜影院在线不卡| 激情五月婷婷亚洲| 99re6热这里在线精品视频| 精品久久久精品久久久| 亚洲精品久久久久久婷婷小说| 亚洲av在线观看美女高潮| 在线观看三级黄色| 在线天堂中文资源库| 女人精品久久久久毛片| 新久久久久国产一级毛片| 99久久中文字幕三级久久日本| 男女午夜视频在线观看 | 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 一级片免费观看大全| 亚洲欧美中文字幕日韩二区| 国产黄频视频在线观看| 飞空精品影院首页| 日韩制服骚丝袜av| 九九爱精品视频在线观看| 久久精品久久久久久噜噜老黄| 久久久久久人妻| 哪个播放器可以免费观看大片| 免费看av在线观看网站| 成人手机av| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻熟女乱码| 天天躁夜夜躁狠狠躁躁| 极品人妻少妇av视频| 91成人精品电影| 国产精品人妻久久久影院| 久久精品aⅴ一区二区三区四区 | av福利片在线| 国产亚洲精品第一综合不卡 | 欧美精品高潮呻吟av久久| 国产国拍精品亚洲av在线观看| 巨乳人妻的诱惑在线观看| 九色成人免费人妻av| 亚洲欧美成人精品一区二区| 人人澡人人妻人| www.av在线官网国产| videossex国产| 国产免费福利视频在线观看| 丝袜脚勾引网站| 国产极品天堂在线| 欧美 日韩 精品 国产| 日本wwww免费看| 日本91视频免费播放| 欧美精品国产亚洲| 18+在线观看网站| 看免费av毛片| 国产精品.久久久| 成人无遮挡网站| 男女下面插进去视频免费观看 | 天天躁夜夜躁狠狠久久av| 国产日韩欧美视频二区| 国产极品粉嫩免费观看在线| 热re99久久精品国产66热6| 国产免费又黄又爽又色| 国产欧美另类精品又又久久亚洲欧美| 韩国精品一区二区三区 | 色婷婷久久久亚洲欧美| 全区人妻精品视频| av国产久精品久网站免费入址| 欧美激情极品国产一区二区三区 | 亚洲国产精品999| 久久97久久精品| 久久女婷五月综合色啪小说| 精品久久久精品久久久| 亚洲国产精品一区二区三区在线| 在现免费观看毛片| 男女啪啪激烈高潮av片| 超碰97精品在线观看| 国产熟女午夜一区二区三区| 午夜精品国产一区二区电影| 日韩不卡一区二区三区视频在线| 亚洲欧美色中文字幕在线| 国产激情久久老熟女| 亚洲国产精品999| 啦啦啦中文免费视频观看日本| 国产av精品麻豆| 日韩免费高清中文字幕av| 欧美xxⅹ黑人| 国产 精品1| 日本与韩国留学比较| 国产爽快片一区二区三区| 18禁裸乳无遮挡动漫免费视频| 欧美+日韩+精品| a 毛片基地|