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

    Dynamics and near-optimal control in a stochastic rumor propagation model incorporating media coverage and L′evy noise

    2022-03-12 07:47:04LianganHuo霍良安andYafangDong董雅芳
    Chinese Physics B 2022年3期
    關(guān)鍵詞:雅芳

    Liang’an Huo(霍良安) and Yafang Dong(董雅芳)

    Business School,University of Shanghai for Science and Technology,Shanghai 200093,China

    Keywords: rumor propagation,near-optimal control,popular science education,media coverage,L′evy noise

    1. Introduction

    In recent years, with the rapid development of the Internet, information exchange between individuals has become more rapid and convenient. As a kind of information with negative influence, the spreading speed and influence range of rumors has grown exponentially. In serious cases, rumors also harm social stability, which should cause all society to be highly vigilant.[1-3]With the outbreak of COVID-19,there have been many rumors, such as “You can prevent the virus with multiple layers of masks”, “A patient flees from Peking Union Medical College Hospital” and so on. These rumors spread wildly on the Internet, resulting in the inability to allocate and use material resources such as masks, and also caused strong negative effects and social panic,which largely hindered the timely prevention and control of the epidemic.Therefore, it is necessary to minimize the negative impact of rumor propagation by studying the optimal control theory of the rumor propagation process.

    As for the propagation process,many scholars have done relevant studies.In 1965,Daley and Kendal proposed the classic rumor propagation model (DK model),[4]which brought a breakthrough for further research on rumor propagation.Later,Maki and Thomson further improved the DK model into the MK model.[5]On this basis, scholars have proposed and studied many extended propagation models.[6-8]With the development of science and technology, the Internet and other complex networks have penetrated into people’s lives and become an important platform for rumor propagation. Therefore, many scholars began to study the influence of complex network topology on the rumor propagation process. Zanette studied the dynamics of an infectious disease-like model of rumor propagation in a small-world network,[9]while Albertet al.studied the tolerance of complex networks to errors.[10]Doeret al.analyzed how rumors spread in social networks and simulated the propagation process of rumors in different network topologies.[11]Liuet al.proposed a Susceptible-Hesitated-Infected-Removed (SHIR) model to study the dynamics of competitive dual information diffusion.[12]Donget al.proposed a new model for the rumor spreading process in online social networks.[13]Zhaoet al.extended the SHIR model as in epidemic disease modeling to investigate the propagation dynamics of information online.[14]Zhuet al.proposed a rumor propagation model with nonlinear functions and time delay in social networks.[15]

    The models involved in the above articles are all deterministic, but uncertainty is inherent in most real-world systems,such as changes in national policies and the involvement of media coverage(that is,the numbers of ignorant and spreaders may change suddenly after the involvement of media coverage). The L′evy jumping process is a discontinuous function with right continuity and left limit, which can better describe this mutation. Therefore, it is of great value to study the impact of the L′evy jumping process on the dynamics of rumor propagation. The stochastic model constructed by the combination of a deterministic model and L′evy noise can better describe the propagation process of rumors in the real world. In recent years, the stochastic propagation model of rumors has also increasingly attracted scholars’ attention. Some scholars have considered L′evy noise, but they mainly studied the dynamic propagation behavior of rumors.[16,17]Dauhooet al.converted the deterministic rumor model into a stochastic rumor model by using the stochastic coefficient representing the influence of the system,proving the existence and uniqueness of global and local solutions.[18]Jiaet al.studied a random rumor propagation model and obtained that the threshold of the random model under the influence of white noise was less than that of the deterministic model.[19]Lvet al.studied a random delayed rumor propagation model.[20]Ankuret al.introduced a random rumor model into a social network containing expert intervention,and drew the conclusion that noise might be one of the reasons for the persistence of rumors in a stochastic environment.[21]

    The purpose of studying the dynamic process of rumor transmission is to better control rumors and reduce their negative impact. Therefore, how to control the spread of rumors has become an issue of concern to an increasing number of scholars.Amineet al.studied a new rumor propagation model and theoretically proved the existence of optimal control.[22]Chenet al.discussed the control of rumor propagation in an emergency and analyzed the stability of the model.[23]The literature contains some research in this area. For example,Huiet al.aimed to study the rumor spreading mechanism on social network platforms during the transmission of COVID-19,and to take education as a control measure for rumor spreading.[24]Wanget al.studied the extent to which media reports before the crisis effectively refuted rumors,and evaluated the impact of media reports before the crisis on the effectiveness of rumor refutation in the experiment.[25]Based on the nonlinear dynamics model, Huoet al.analyzed the rumor propagation process and studied the influence of media coverage and government emergency strategies on the information propagation dynamics after emergencies.[26]

    There are two possible reasons for the spread of rumors.First of all, as false information, rumors are extremely misleading and provocative. Some people do not have enough knowledge to distinguish between true and false information,which leads people to listen to rumors and propagate them into facts. Secondly,when a rumor is generated and propagated in an emergency, there is no corresponding real information released in time to refute the rumor,which leads to the constant spread of the rumor. In this case, popular science education and media coverage are two effective methods to prevent the crisis of rumor spreading. Popular science education can improve people’s knowledge level and the ability to identify rumors,and directly transform the ignorant into the immune. In this way, the source of rumor propagation can be reasonably controlled and the path of rumor propagation can be further cut off to weaken the influence of rumors. Media coverage can refute rumors in a targeted way and reduce the probability of rumors propagating,which can not only let people know the truth in time,but also guide people to choose the right behavior and prevent overreaction. In addition,the exact solution of the equation of state and the adjoint equation is difficult to obtain and the exact control is also difficult when the stochastic factors are considered.Therefore,it is of great practical significance to study the approximate optimal control of the rumor propagation model.

    Based on the above analysis and actual situation,a rumor spreading model that joins the L′evy process can describe the dynamic process of rumor propagation more accurately. Most researchers use deterministic models when conducting rumor spreading research,and rarely consider random factors. In addition,some researchers only conduct research on the process,and rarely mention control strategies. Therefore, this paper proposes a stochastic model combining media coverage and L′evy noise to describe the dynamic process of rumors spread in the real world. By introducing two control strategies,popular science education and media coverage, an optimal control problem is proposed to minimize the influence of rumors spread and control costs.

    The structure of this paper is as follows. Section 2 introduces the preparatory knowledge and establishes the stochastic ignorant-spreader-recovered(ISR)model;Section 3 and 4 give the necessary and sufficient conditions for near-optimal control; Section 6 carries on the numerical simulation; and Section 7 summarizes the whole paper.

    2. Model formulation

    In the classical ISR model of rumor propagation,[27]nodes in social networks are divided into three categories: ignorant(I)is a node without any information about the rumor,spreader(S)is the node that accepts the rumor and propagates it, and recovered (R) is the node that fnially decides to stop propagating rumors after contact with rumors. In social networks,rumors propagate among people in accordance with the rules of mass action. The ignorant node has not contacted the rumor and will not propagate it,but it may become a spreader by contacting the spreader nodes.

    The propagation process of rumors can be described as follows:

    (i) The main propagation model of rumors is communication after the interaction of individual opinions, that is,the information exchange between individuals through telephone, instant messaging software, face-to-face conversation and other means.

    (ii)Assume that a spreader propagates rumors with a certain probability when they meet an ignorant person. Media coverage can specifically refute rumors by releasing facts,thus gradually reducing the probability of rumor propagation and finally reaching a stable value.

    (iii) When the spreader meets the recovered or another spreader,they lose the motivation to spread rumors.

    (iv) Rumors spread within a certain community, and the total number of people is variable. Those who initially enter the community are ignorant.

    The total population isN=I+S+R. Due to the contact between different types of individuals,information is transmitted and diffused,and then the types of individuals are changed,as shown in Fig.1.

    Fig.1. Schematic diagram of a rumor under the influence of popular science education and media reports.

    We develop this problem on the basis of a deterministic ISR model[28]as follows:

    Generally, a L′evy process ~L(t) consists of the sum of a linear drift, a Brownian motionB(t) and a superposition of centered independent Poisson processes with different jump sizesλ(dv) which is the rate of arrival intensity of the Poisson process with a jump of sizev. According to the L′evy-Ito decomposition theorem,[29]we get

    In the above model, when rumors spread to a certain scale, they may endanger social harmony and stability, and bring huge losses to the related people. Therefore,it is necessary to study how to control the propagation of rumors. The reason why rumors are widely spread may be that some people are unable to distinguish between true and false information due to their insufficient knowledge reserve, so they take rumors as facts to propagate.Another reason may be that there is no corresponding real information about the breaking event to refute the rumor,so the rumor keeps propagating.

    Under such circumstances,popular science education and media reports are the most effective ways to curb the spread of rumors. On the one hand, popular science education can improve a person’s cognitive level so as to directly transform the ignorant into the immune. Popular science education stops rumors from spreading in the first place and reduces the number of the ignorant while increasing the number of the immune.On the other hand,media coverage can also be a targeted refutation of rumors. By increasing media coverage,the truth can be publicized in a timely manner and public opinion can be actively guided so as to reduce the number of spreaders and the negative impact of rumor propagation. Reasonable coverage and time for media coverage and popular science education can reduce the number of spreaders and ignorants to a lower level. Therefore, we introduce education and media controlu(t)=(uθ(t),um(t))into the stochastic ISR model(3)as follows:

    whereumanduθrepresent the control variables for the influence of science popularization education and media coverage,respectively. LetUad?Rbe a bounded nonempty closed set.uθ,um:[0,1]×Ω →Uis called admissible, if it is anFt-adapted process with values inUad.

    The purpose of optimal control is to minimize the negative effects caused by the propagation of rumors and maximize the positive social utility. This paper chooses the following aspects to consider the social utility of the propagation of untrue information:

    (I) Minimize the number of ignorant individualsIand spreadersSof rumors.

    (II)Minimize the cost of variable control in the process.

    We define the associated objective function as

    3. Analysis of near-optimal control for stochastic ISR model

    Before we study near-optimal control,we give some preliminary knowledge as follows; let (Ω,F,(Ft)0≤t≤T,P) be a complete filtered probability space on a finite time horizon[0,T], and we assume that (Ft)0≤t≤Tis the natural filtration andF=Ft. If the control processu:[0,T]×Ω →Uis anFt-adapted process with values inU,then it is called admissible.We useUadto denote the set that consists of all admissible controls;|·|denotes the form of an Euclidean space;fxis the partial derivative offwith respect tox;χsdenotes the indicator function of a setS;andX+Yis the set{x+y:x ∈X,y ∈Y}for any setsXandY.

    Definition 1[32]Ifu*(·) attains the minimum ofJ(0,x0;u(·)), then an admissible controlu*(·) is called optimal control.

    Definition 2[32]For any givenε >0, if|J(0,x0;uε(·))-V(0,x0)|≤ε, an admissible controluε(·) is calledε-optimal control.

    Definition 3[32]A family of admissible controls{uε(·)}parameterized byε >0 and any elementuε(·) in the family are called near-optimal control,if|J(0,x0;uε(·))-V(0,x0)|≤δ(ε)holds for a sufficiently smallε >0,whereδis a function ofεsatisfyingδ(ε)→0 asε →0. The estimateδ(ε)is called an error bound. Ifδ(ε)=cεkfor somec,k >0,thenuε(·)is called near-optimal with orderεk.

    Lemma 1 Ekeland’s principle.[33]Let (Q,d) be a complete metric space andF(·):Q →Rbe a lower-semicontinuous and bounded from below. For anyε >0, we assume thatuε(·)∈Qsatisfies

    (H3)The control setUis convex.

    Under hypotheses(H1)and(H2),the adjoint equation admits a unique adapted solution. By Mazur’s theorem,[34]we can then obtain the existence of optimal control under hypothesis(H3).

    3.1. Sufficient conditions of near-optimal controls

    In order to obtain sufficient conditions of near-optimal controls, we give the following lemmas for the moment estimation of the state and adjoint equation. For convenience,letx(t) = (x1(t),x2(t),x3(t))T= (I(t),S(t),R(t))T,u(t) =(um(t),uθ(t));then we get

    3.1.1. Some priori estimates of the ignorant, spreader and recovered

    Then we give some prior estimates of the ignorant,spreader and recovered. The following lemma illustrates the boundedness ofεmoment estimates for the ignorant,spreader and recovered.

    The second moment estimate of the adjoint Eq. (10) is given in the following lemma.

    Lemma 3 For anyum(t),uθ(t),u′m(t),u′θ(t)∈Uad, we have

    Using Gronwall’s inequality,we can derive Eq.(11)from the above inequality. Then the proof is complete.

    3.1.2. Sufficient conditions for near-optimal controls

    In this section,we define the Hamiltonian functionH(t,x(t),u(t),p(t),q(t),r(v)),H:[0,T]×R3+×Uad×R3+×R2+×R2+→R,

    According to Eqs.(31)and(23),we can obtain

    Therefore,from Eq.(38)and hypothesis(H2),we get

    According to lemma 3 and the definition of ~d, the expected conclusion can be obtained from Eqs. (30) and (39) and Holder’s inequality.

    3.2. Necessary conditions for near-optimal controls

    In this section, we discuss the necessary conditions of near-optimal control. First, we give the following priori estimates to the ignorant, spread and recovered persons of the system to get our results.

    In order to give the following lemma, we first define a metric on the admissible control fieldUad[0,T]as follows:

    where mes stands for the Lebesgue measure. SinceUadis closed,it can be proven thatUadis a complete metric space underd.[34]The problem of moment estimation of the distance between two state variables is discussed by using the distance of the defined control variables.

    Lemma 4 For allk ≥0 and 0<l <1 satisfyingkl <1, andu(t),^u(t)∈Uad, along with the corresponding trajectoriesx(t),x′(t), there exists a positive imprecise parameterC=C(k,l)such that

    The proof is complete,which is the result of the estimation of the state Eq.(6).

    For the result of moment estimation of the adjoint equation,we have the following lemma.

    Lemma 5 Let hypotheses (H2) and (H3) hold. For any 0<l <1 and 1<k <2 satisfying (1+l)k <2, there exists a positive imprecise parameterC=C(l,k) such that for anyu(t),^u(t)∈Ualong with the corresponding trajectoryx(t),^x(t)and the solutionp(t),q(t),r(v),(^p(t),^q(t),^r(v))of the adjoint equation,we have

    4. Numerical simulations

    In this section, we use MATLAB to perform numerical simulations and use these numerical examples to illustrate the impact of popular science education, media reports and L′evy noise on the dynamics of rumor spreading.We set the parameter values as follows:Λ=0.25,μ=0.15,β1=0.6,β2=0.4,λ=0.2,θ=0.2,σ=0.009.

    According to Milstein’s method,[37]we can get the corresponding discretization equation of the state Eqs.(4)and(10)as follows:

    The intensity of L′evy noise can affect the rumor propagation process. Figure 2 depicts the influence of L′evy noise on the rumor propagation process.r(v)stands for response intensities. As shown in the figure,the trend chart of the number of spreaders over time is different under different response intensitiesr(v). The greater the value ofr(v),the faster the value ofS(t)decreases. Therefore, we can reduce the ultimate spread scale of rumors by strengthening external intervention. For example, when rumors spread wildly and even cause bad social impact,the government should promptly popularize relevant knowledge about rumors to guide people to correctly view the content of rumors,and relevant media should also increase their reporting efforts. In this way,the number of spreaders of rumors can be reduced by strengthening external intervention.

    We also simulated the time evolution diagram of the number of people in different states during rumor propagation under four different control conditions: (1)um=0,uθ=0 represents the uncontrolled situation;(2)um=0,uθ/=0 represents the situation where only popular science education control is added;(3)um/=0,uθ=0 represents the situation where only media coverage control is added; and(4)um/=0,uθ/=0 represents the situation where both media coverage and popular science education control are added. As can be seen from Fig. 4, when media coverage and popular science education control are added at the same time,the peak value and number of spreaders drop sharply, and the influence scope of rumors becomes smaller. Compared with the other three conditions,the rumor propagation process is easier to control, indicating that this situation is nearoptimal control. For example, in the early days of the COVID-19 outbreak,many people were not familiar with the specific case of novel coronavirus. It was rumored that “symptoms such as fever and cough were signs of infection with COVID-19”. The rampant propagation of this rumor was bound to cause public panic. In order to prevent this crisis,the government provided the latest research on COVID-19 in official news and educated the public through multiple media channels. Popular science education allowed people to learn more about the symptoms of COVID-19 infection, which directly reduced the number of the ignorant and enabled them to convert directly to the immune. Media coverage directly reduced the number of communicators through targeted refutation of the rumor.

    In addition, it can be seen from Figs. 3-5 that the dynamic optimal control is effective. Both the control of popular science education and the control of media coverage decline faster than the situation without control. Under the control of government and media intervention,the number of spreaders declines faster than in the absence of control. For example,during the COVID-19 outbreak,the public listened to the rumor that “Shuanhuanglian could prevent the novel coronavirus”,so they rushed to buy Shuanhuanglian,causing chaos in the market. The government quickly educated the public about the novel coronavirus treatment, and many media outlets debunked the rumor. The double control increased the number of people immune to the rumor. Taking a closer analysis,we can see that when the media coverage control strategy is added, the control effect is good but not great. When the government science education control strategy is added, the control effect is significantly improved. This simulation result shows that the government science education control strategy is more effective than the media coverage control strategy.

    Fig.2. Path of S(t)under different response intensities r(v).

    Fig. 3. Path of I(t) for the stochastic ISR model with initial (I0,S0,R0)=(0.8,0.2,0.05).

    Fig. 4. Path of S(t) for the stochastic ISR model with initial (I0,S0,R0)=(0.8,0.2,0.05).

    Fig. 5. Path of R(t) for the stochastic ISR model with initial (I0,S0,R0)=(0.8,0.2,0.05).

    As for the near-optimal control,it can be seen from Fig.6 that solutionsp1,p2, andp3of adjoint Eq. (10) converge to zero in the end, indicating that the objective function reaches the minimum value.

    Fig. 6. Path of p1, p2, p3 for the stochastic ISR model with initial(I0,S0,R0)=(0.8,0.2,0.05).

    5. Conclusion

    This paper studies the problem of near-optimal control of rumor propagation. We first construct a stochastic model that includes media reports and L′evy noise. Then,in order to minimize the impact of rumor propagation and control costs,we introduce two control strategies: popular science education and media coverage. Using the Hamiltonian function and Pontryagin maximum principle, the necessary and sufficient conditions for the model to be close to the optimal control are obtained.Finally,we verify the effects of popular science education,media report control and noise control through numerical simulations. The results show that popular science education and media coverage can effectively control the spread of rumors, and these two control strategies can reduce the cost of rumor dispelling to the lowest level. At the same time,L′evy noise will also affect the process of rumor propagation.The louder the noise, the faster the attenuation speed of the spreader. In this paper,we assume that the parameters are precisely known,but may not be true due to various uncertainties.The influence of uncertain parameters and time delay on the rumor control process is worthy of further study.

    Acknowledgments

    Project supported by the Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning,and the Project for the Natural Science Foundation of Shanghai,China(Grant No.21ZR1444100),and the Project for the National Natural Science Foundation of China(Grant Nos.72174121,71774111,71871144,and 71804047).

    猜你喜歡
    雅芳
    Dynamics of information diffusion and disease transmission in time-varying multiplex networks with asymmetric activity levels
    An improved ISR-WV rumor propagation model based on multichannels with time delay and pulse vaccination
    小情緒
    雅芳 老樹(shù)能否開(kāi)新花
    新雅芳能否浴火重生?
    雅芳再變革
    雅芳 空降新帥
    亚洲精品成人av观看孕妇| 国产一区二区三区av在线| 亚洲国产毛片av蜜桃av| 国产黄色视频一区二区在线观看| 国产日韩一区二区三区精品不卡| 亚洲伊人久久精品综合| 亚洲国产欧美一区二区综合| 久久精品人人爽人人爽视色| 热99国产精品久久久久久7| 又粗又硬又长又爽又黄的视频| 女人久久www免费人成看片| 欧美激情 高清一区二区三区| 亚洲图色成人| 欧美激情极品国产一区二区三区| 亚洲七黄色美女视频| 性色av乱码一区二区三区2| bbb黄色大片| 婷婷丁香在线五月| 亚洲免费av在线视频| 亚洲,一卡二卡三卡| 无限看片的www在线观看| 韩国精品一区二区三区| 日韩欧美一区视频在线观看| 久久av网站| 国产成人精品久久二区二区91| 国产黄频视频在线观看| 国产成人一区二区三区免费视频网站 | 美女大奶头黄色视频| 欧美成狂野欧美在线观看| 老司机午夜十八禁免费视频| 99久久99久久久精品蜜桃| 午夜激情av网站| 女人被躁到高潮嗷嗷叫费观| 女性生殖器流出的白浆| 免费高清在线观看视频在线观看| 亚洲中文av在线| 天堂中文最新版在线下载| 精品一区二区三区四区五区乱码 | 色播在线永久视频| 大片免费播放器 马上看| 欧美日韩亚洲高清精品| 久久毛片免费看一区二区三区| 日本黄色日本黄色录像| 亚洲欧美成人综合另类久久久| 成人三级做爰电影| 精品卡一卡二卡四卡免费| 午夜精品国产一区二区电影| 黄片小视频在线播放| 亚洲国产日韩一区二区| 亚洲一码二码三码区别大吗| 亚洲熟女精品中文字幕| 美女福利国产在线| 日韩中文字幕欧美一区二区 | 后天国语完整版免费观看| 免费不卡黄色视频| 最新在线观看一区二区三区 | 丝袜脚勾引网站| 国产精品一二三区在线看| 蜜桃在线观看..| 99国产综合亚洲精品| 丝袜人妻中文字幕| 天堂8中文在线网| 无限看片的www在线观看| 婷婷色av中文字幕| 无遮挡黄片免费观看| 国产有黄有色有爽视频| 精品国产乱码久久久久久小说| 成人三级做爰电影| 美女福利国产在线| 国产老妇伦熟女老妇高清| 国产av国产精品国产| 欧美大码av| 一级片'在线观看视频| 日韩制服骚丝袜av| 如日韩欧美国产精品一区二区三区| 老司机影院毛片| videos熟女内射| 黑人巨大精品欧美一区二区蜜桃| 99热全是精品| 亚洲三区欧美一区| 黑人巨大精品欧美一区二区蜜桃| 搡老岳熟女国产| 国产淫语在线视频| 三上悠亚av全集在线观看| 欧美亚洲 丝袜 人妻 在线| 国产三级黄色录像| 亚洲男人天堂网一区| avwww免费| 国产午夜精品一二区理论片| 天天影视国产精品| 午夜精品国产一区二区电影| 国产亚洲av高清不卡| 日韩 亚洲 欧美在线| 国产精品久久久久成人av| 精品一区二区三区av网在线观看 | 亚洲国产av影院在线观看| 七月丁香在线播放| 国产精品人妻久久久影院| 欧美亚洲 丝袜 人妻 在线| 建设人人有责人人尽责人人享有的| 国产麻豆69| 青草久久国产| 国产精品国产三级专区第一集| 国产成人精品在线电影| 国产欧美亚洲国产| 欧美人与善性xxx| 日韩欧美一区视频在线观看| 久久人人97超碰香蕉20202| 一区二区三区激情视频| 中文乱码字字幕精品一区二区三区| 99九九在线精品视频| 97人妻天天添夜夜摸| 老汉色∧v一级毛片| 午夜影院在线不卡| 90打野战视频偷拍视频| 高清av免费在线| 一级毛片我不卡| 在线观看一区二区三区激情| 成在线人永久免费视频| 国产一区二区在线观看av| 一本综合久久免费| 巨乳人妻的诱惑在线观看| 女性生殖器流出的白浆| 观看av在线不卡| 国产av国产精品国产| 国产日韩欧美视频二区| 精品人妻熟女毛片av久久网站| 国产精品 国内视频| 女警被强在线播放| 女性被躁到高潮视频| 搡老岳熟女国产| 老司机影院毛片| 一级片'在线观看视频| 人妻人人澡人人爽人人| 女人高潮潮喷娇喘18禁视频| 老鸭窝网址在线观看| 少妇的丰满在线观看| 真人做人爱边吃奶动态| 伊人久久大香线蕉亚洲五| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费| 久久久久久亚洲精品国产蜜桃av| 中文乱码字字幕精品一区二区三区| 国产亚洲精品第一综合不卡| 亚洲欧美精品自产自拍| svipshipincom国产片| av不卡在线播放| 老熟女久久久| 欧美激情极品国产一区二区三区| 精品久久久久久电影网| 丰满少妇做爰视频| 久久久欧美国产精品| 最近最新中文字幕大全免费视频 | 美女视频免费永久观看网站| 欧美激情极品国产一区二区三区| 国产成人av教育| 欧美人与性动交α欧美软件| 老熟女久久久| 人人妻人人澡人人爽人人夜夜| 最新在线观看一区二区三区 | 亚洲自偷自拍图片 自拍| 中文字幕制服av| 久久精品aⅴ一区二区三区四区| 久久99一区二区三区| bbb黄色大片| 亚洲色图综合在线观看| 晚上一个人看的免费电影| 少妇裸体淫交视频免费看高清 | 嫩草影视91久久| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 王馨瑶露胸无遮挡在线观看| 一本一本久久a久久精品综合妖精| 日韩伦理黄色片| 天天躁日日躁夜夜躁夜夜| 狠狠精品人妻久久久久久综合| 日本黄色日本黄色录像| 中文精品一卡2卡3卡4更新| 免费在线观看视频国产中文字幕亚洲 | 黄色怎么调成土黄色| 好男人电影高清在线观看| 大香蕉久久网| 日韩伦理黄色片| 久久久久网色| av线在线观看网站| 国产91精品成人一区二区三区 | 亚洲精品av麻豆狂野| 国产一区二区 视频在线| 中文字幕av电影在线播放| 国产成人a∨麻豆精品| 精品高清国产在线一区| 一本大道久久a久久精品| 免费日韩欧美在线观看| 黄色视频在线播放观看不卡| 亚洲精品第二区| 亚洲精品国产av蜜桃| 青春草视频在线免费观看| 欧美精品亚洲一区二区| 婷婷成人精品国产| 韩国高清视频一区二区三区| 男女之事视频高清在线观看 | xxxhd国产人妻xxx| 啦啦啦中文免费视频观看日本| 久久精品国产综合久久久| 老司机深夜福利视频在线观看 | 9色porny在线观看| 啦啦啦在线免费观看视频4| 欧美日韩国产mv在线观看视频| 亚洲欧美一区二区三区久久| 中文乱码字字幕精品一区二区三区| 成年av动漫网址| 亚洲精品一二三| 在线观看免费日韩欧美大片| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 国产av一区二区精品久久| 视频在线观看一区二区三区| 欧美另类一区| 久久精品久久久久久久性| 99国产精品免费福利视频| 在线av久久热| 2021少妇久久久久久久久久久| 美女午夜性视频免费| 一区在线观看完整版| 亚洲中文日韩欧美视频| 国产亚洲精品第一综合不卡| 又大又黄又爽视频免费| 99国产精品99久久久久| 国产精品麻豆人妻色哟哟久久| 免费人妻精品一区二区三区视频| 久久国产亚洲av麻豆专区| 成年av动漫网址| 免费看十八禁软件| av在线播放精品| 日韩中文字幕欧美一区二区 | av一本久久久久| 精品卡一卡二卡四卡免费| 少妇精品久久久久久久| 亚洲欧美精品自产自拍| 黑人猛操日本美女一级片| www.自偷自拍.com| 纵有疾风起免费观看全集完整版| 久久久精品94久久精品| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 久久99热这里只频精品6学生| 少妇粗大呻吟视频| 成人免费观看视频高清| 黑人猛操日本美女一级片| 亚洲成人免费电影在线观看 | 男女高潮啪啪啪动态图| 妹子高潮喷水视频| 久久精品国产亚洲av涩爱| 国产麻豆69| 亚洲欧美一区二区三区久久| 精品人妻熟女毛片av久久网站| 亚洲伊人久久精品综合| 亚洲专区国产一区二区| 狠狠婷婷综合久久久久久88av| 中文字幕最新亚洲高清| 国产日韩欧美亚洲二区| 色精品久久人妻99蜜桃| 中文字幕精品免费在线观看视频| 午夜福利乱码中文字幕| 亚洲伊人久久精品综合| 一二三四在线观看免费中文在| 狠狠精品人妻久久久久久综合| 亚洲精品久久久久久婷婷小说| 日韩av免费高清视频| 欧美 亚洲 国产 日韩一| 国产一区二区三区综合在线观看| 精品久久蜜臀av无| 精品卡一卡二卡四卡免费| 青春草视频在线免费观看| 亚洲av电影在线进入| 777米奇影视久久| 青春草亚洲视频在线观看| 丁香六月天网| 视频在线观看一区二区三区| 国产av国产精品国产| 91麻豆精品激情在线观看国产 | 欧美成狂野欧美在线观看| 我的亚洲天堂| 一区二区日韩欧美中文字幕| 国产精品一国产av| 精品高清国产在线一区| 久久鲁丝午夜福利片| 日韩制服骚丝袜av| 中文字幕精品免费在线观看视频| 日韩熟女老妇一区二区性免费视频| 国产成人免费观看mmmm| 一本大道久久a久久精品| 少妇被粗大的猛进出69影院| av国产久精品久网站免费入址| 亚洲,欧美精品.| 成人三级做爰电影| videos熟女内射| 日韩人妻精品一区2区三区| 中文字幕色久视频| 免费黄频网站在线观看国产| 超色免费av| a级毛片黄视频| 人人妻人人爽人人添夜夜欢视频| 丝袜人妻中文字幕| 亚洲欧美一区二区三区国产| 国产精品国产av在线观看| 男人操女人黄网站| 欧美人与性动交α欧美软件| 亚洲久久久国产精品| 久久这里只有精品19| 天天躁日日躁夜夜躁夜夜| 国产又色又爽无遮挡免| 国产麻豆69| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品成人av观看孕妇| 国产在线免费精品| 亚洲精品自拍成人| 久久99精品国语久久久| 午夜老司机福利片| 脱女人内裤的视频| 丝袜美足系列| 一边摸一边做爽爽视频免费| 国产在线观看jvid| 国产福利在线免费观看视频| 成人国语在线视频| 免费在线观看日本一区| 少妇 在线观看| 97在线人人人人妻| 男女床上黄色一级片免费看| 最近最新中文字幕大全免费视频 | 国产亚洲午夜精品一区二区久久| 99久久综合免费| 国产日韩欧美视频二区| 国产在线视频一区二区| 中文字幕色久视频| 国产日韩欧美在线精品| 女警被强在线播放| 乱人伦中国视频| 国产精品av久久久久免费| 少妇猛男粗大的猛烈进出视频| 青春草视频在线免费观看| 免费高清在线观看视频在线观看| 脱女人内裤的视频| av欧美777| 十分钟在线观看高清视频www| 亚洲欧美精品综合一区二区三区| 成在线人永久免费视频| 777米奇影视久久| 亚洲欧美一区二区三区国产| 久久人妻熟女aⅴ| 人人澡人人妻人| 国产日韩欧美亚洲二区| 男人操女人黄网站| 欧美 日韩 精品 国产| 电影成人av| 中文字幕高清在线视频| 精品国产一区二区三区四区第35| 欧美日韩一级在线毛片| 久久人人爽人人片av| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生| 日韩 欧美 亚洲 中文字幕| 亚洲精品日韩在线中文字幕| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| 五月天丁香电影| 国产在线视频一区二区| 亚洲,一卡二卡三卡| 久久亚洲精品不卡| 久久这里只有精品19| 成人国产一区最新在线观看 | 热re99久久国产66热| 在线观看免费视频网站a站| 国产成人一区二区三区免费视频网站 | 秋霞在线观看毛片| 国产免费福利视频在线观看| 中文字幕av电影在线播放| 91精品三级在线观看| 国产激情久久老熟女| 999久久久国产精品视频| 亚洲人成电影免费在线| 亚洲熟女毛片儿| 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲国产一区二区在线观看 | 国产片特级美女逼逼视频| 狠狠婷婷综合久久久久久88av| 免费在线观看视频国产中文字幕亚洲 | 高清黄色对白视频在线免费看| 久久久久网色| 国产黄色免费在线视频| 国产一级毛片在线| 在线观看免费午夜福利视频| 一级黄色大片毛片| 亚洲av电影在线进入| 两个人看的免费小视频| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 欧美av亚洲av综合av国产av| 久久性视频一级片| 男人舔女人的私密视频| 成年人黄色毛片网站| 国产高清视频在线播放一区 | 大型av网站在线播放| 久久午夜综合久久蜜桃| 国产欧美日韩精品亚洲av| 脱女人内裤的视频| 色94色欧美一区二区| 亚洲av片天天在线观看| 大话2 男鬼变身卡| 亚洲国产精品成人久久小说| 亚洲精品一二三| 两个人看的免费小视频| 男女无遮挡免费网站观看| 欧美少妇被猛烈插入视频| 天天躁狠狠躁夜夜躁狠狠躁| av网站免费在线观看视频| 国产精品 欧美亚洲| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 美女视频免费永久观看网站| 亚洲av国产av综合av卡| 欧美亚洲日本最大视频资源| 欧美亚洲 丝袜 人妻 在线| 丝袜美足系列| www.av在线官网国产| 激情五月婷婷亚洲| 亚洲成av片中文字幕在线观看| 男女边吃奶边做爰视频| 满18在线观看网站| 一区二区三区乱码不卡18| 人妻一区二区av| 久久精品久久精品一区二区三区| 亚洲成人国产一区在线观看 | 欧美国产精品va在线观看不卡| 欧美精品一区二区大全| netflix在线观看网站| 国产欧美日韩精品亚洲av| 亚洲av电影在线观看一区二区三区| 热re99久久精品国产66热6| 国产极品粉嫩免费观看在线| 亚洲精品在线美女| 一区二区av电影网| 美女国产高潮福利片在线看| 欧美日韩亚洲综合一区二区三区_| 久久久久网色| 欧美日韩亚洲综合一区二区三区_| 91精品伊人久久大香线蕉| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站| 精品亚洲乱码少妇综合久久| 黄色片一级片一级黄色片| videos熟女内射| 国产成人av教育| 免费黄频网站在线观看国产| 成人国语在线视频| 亚洲精品一卡2卡三卡4卡5卡 | 日本五十路高清| 99精品久久久久人妻精品| 色播在线永久视频| 日本五十路高清| 亚洲欧美一区二区三区黑人| 国产亚洲午夜精品一区二区久久| 国产精品国产三级国产专区5o| 午夜福利,免费看| av网站免费在线观看视频| 咕卡用的链子| 欧美97在线视频| 日韩av在线免费看完整版不卡| 亚洲国产精品一区三区| 热re99久久精品国产66热6| 天天操日日干夜夜撸| 亚洲精品国产色婷婷电影| 久久久久久久大尺度免费视频| 视频区欧美日本亚洲| 成人18禁高潮啪啪吃奶动态图| 婷婷成人精品国产| 久久人妻福利社区极品人妻图片 | 免费看av在线观看网站| 亚洲精品久久成人aⅴ小说| 久久久亚洲精品成人影院| 国产黄色免费在线视频| 久久久久久久久免费视频了| 精品亚洲成a人片在线观看| 久久久亚洲精品成人影院| 欧美精品亚洲一区二区| 成年美女黄网站色视频大全免费| 在线观看免费高清a一片| www.自偷自拍.com| 后天国语完整版免费观看| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 91精品国产国语对白视频| 99re6热这里在线精品视频| 亚洲精品国产色婷婷电影| 久久精品熟女亚洲av麻豆精品| 首页视频小说图片口味搜索 | 久久中文字幕一级| 真人做人爱边吃奶动态| 国产欧美日韩精品亚洲av| av天堂久久9| 国产xxxxx性猛交| 两人在一起打扑克的视频| 午夜精品国产一区二区电影| 两性夫妻黄色片| 丰满饥渴人妻一区二区三| 色精品久久人妻99蜜桃| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| 国产成人a∨麻豆精品| 黄频高清免费视频| 亚洲熟女毛片儿| 精品久久久久久电影网| 欧美变态另类bdsm刘玥| 中文字幕高清在线视频| 久久久亚洲精品成人影院| 午夜精品国产一区二区电影| 少妇的丰满在线观看| 国产日韩欧美在线精品| 青青草视频在线视频观看| 国产成人一区二区在线| 国语对白做爰xxxⅹ性视频网站| 久久久国产精品麻豆| 五月天丁香电影| 亚洲av日韩在线播放| 免费在线观看日本一区| 91精品三级在线观看| 精品国产国语对白av| 欧美激情高清一区二区三区| 高清不卡的av网站| 电影成人av| 国产黄色视频一区二区在线观看| 亚洲人成电影免费在线| 好男人电影高清在线观看| 久久久国产一区二区| 男女无遮挡免费网站观看| 视频区图区小说| 亚洲伊人色综图| 一区在线观看完整版| 免费少妇av软件| 一本—道久久a久久精品蜜桃钙片| 久久天堂一区二区三区四区| 丰满迷人的少妇在线观看| 叶爱在线成人免费视频播放| 国产精品久久久久成人av| 国产亚洲精品久久久久5区| 侵犯人妻中文字幕一二三四区| 久久中文字幕一级| 国产精品一二三区在线看| 18在线观看网站| 伊人久久大香线蕉亚洲五| 老司机午夜十八禁免费视频| a 毛片基地| av网站免费在线观看视频| 水蜜桃什么品种好| bbb黄色大片| 国产深夜福利视频在线观看| 成人影院久久| 一二三四社区在线视频社区8| 一级毛片我不卡| 精品免费久久久久久久清纯 | 亚洲国产日韩一区二区| 日日爽夜夜爽网站| 赤兔流量卡办理| av国产久精品久网站免费入址| 免费高清在线观看日韩| 亚洲专区国产一区二区| 两性夫妻黄色片| 美女午夜性视频免费| 欧美少妇被猛烈插入视频| 色网站视频免费| 最近手机中文字幕大全| 69精品国产乱码久久久| 丝瓜视频免费看黄片| 国产片内射在线| 中文字幕人妻丝袜一区二区| 男女边吃奶边做爰视频| 国产av精品麻豆| 最近手机中文字幕大全| 亚洲熟女精品中文字幕| 欧美+亚洲+日韩+国产| 午夜老司机福利片| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 日本五十路高清| 欧美黄色淫秽网站| 无遮挡黄片免费观看| 亚洲色图综合在线观看| 老司机影院毛片| 婷婷色av中文字幕| 亚洲av成人精品一二三区| 午夜免费男女啪啪视频观看| 国产日韩一区二区三区精品不卡| 在线看a的网站| 亚洲三区欧美一区| 午夜激情久久久久久久| 自线自在国产av| 欧美大码av| 欧美人与善性xxx| 成年人午夜在线观看视频| 满18在线观看网站| 人人妻,人人澡人人爽秒播 | 日本av免费视频播放| 日本五十路高清| 亚洲五月婷婷丁香| 日韩精品免费视频一区二区三区| 久久ye,这里只有精品| 成在线人永久免费视频| 热re99久久精品国产66热6| 亚洲国产av影院在线观看| 日韩电影二区| 国产高清国产精品国产三级| 亚洲中文av在线| 老司机影院成人| 免费高清在线观看视频在线观看| 91字幕亚洲|