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

    Influence fast or later: Two types of influencers in social networks

    2022-06-29 09:22:10FangZhou周方ChangSu蘇暢ShuqiXu徐舒琪andLinyuan呂琳媛
    Chinese Physics B 2022年6期

    Fang Zhou(周方) Chang Su(蘇暢) Shuqi Xu(徐舒琪) and Linyuan L¨u(呂琳媛)

    1Yangtze Delta Region Institute(Huzhou)&Institute of Fundamental and Frontier Sciences,University of Electronic Science and Technology of China,Huzhou 313001,China

    2Beijing Computational Science Research Center,Beijing 100193,China

    Keywords: social networks,fast influencers,final influencers,spreading dynamics,degree assortativity

    1. Introduction

    In network science, vital nodes identification[1–5]is an intriguing and important problem, which has attracted longstanding scientific attention. Vital nodes are considered as nodes that play an important role in the structure or function of networks. In the scenario of information spreading,[3,6]vital nodes are also known as influencers (or “influentials”),namely, those who can influence a large fraction of nodes in the spreading dynamics. Identifying influencers have crucial meanings on controlling or speeding the spreading process in networks and can support various practical applications. For example, in advertisement marketing campaigns,[1]having a celebrity promote a product on social media can largely enhance the visibility of the product, thus driving the product sales and obtaining a good marketing effect.

    Researchers have made great efforts in the identification of influencers and have gained a series of achievements.[1–4,7–10]To determine whether a node in a network is an influencer, most research relies on quantifying the spreading capacity of the node in spreading dynamics.The spreading process is usually specified as follows: for a given network, we choose a node (typically referred to as the “seed node”) and employ a spreading model to simulate the spreading process; when the spreading stops, we measure the seed node’s spreading capacity, and the nodes with outstanding capacity are defined as influencers. Researchers have evaluated many methods’ ability to reproduce the ranking of nodes’ spreading capacity. One of the most effective and simplest methods is the centrality metric that considers the network structure of the focal node. Commonly-used centrality metrics that are introduced to rank nodes include local metrics such as degree,[11]H-index,[12]LocalRank,[3]collective influence,[4]and global metrics such as eigenvector centrality,[13,14]betweenness centrality,[15]closeness centrality,[16]dynamics-sensitive centrality,[17]etc. Besides,other iterative-based methods proposed by researchers include PageRank[18]and its parameter-free variant LeaderRank,[19]HITS,[20]and SALSA.[21]

    Based on the above process of influencer identification,researchers mainly focus on the final spreading capacity of nodes (i.e., the result when the spreading is finished), but ignore nodes with great spreading capacity at the early stage.Zhouet al.[7]provided the first step towards uncovering this type of nodes, which they called “fast influencers”. The authors revealed that fast influencers are different from traditional influencers that pursue the final big impact, and they proposed a metric called social capital to identify them. Fast influencers are of both theoretical and practical significance.Potential applications include, for example, maximizing a newly-released movie’s promotion effect in a limited time and cutting off the spreading of the latest rumor as soon as possible.

    Despite the fact that fast influencers and final influencers have a clear distinction on definitions, there lacks an understanding of the underlying differences between fast and final influencers,which is what we aim to achieve in this work. To this end, we focus on comparing the two types of influencers based on two network-based features: local structure[22]and degree assortativity,[23–25]which have been revealed to influence the spreading process in networks. Specifically, we analyze the differences between the fast and final influencers in terms of their spreading capacity and four local centrality metrics, degree,h-index,k-core and LocalRank, as well as study how the network degree assortativity influences the fraction of the two types of vital nodes.The structure of the paper is organized as follows,firstly in Section 2,we introduce the related approaches,definitions of fast-only and final-only influencers,and the details of eight studied datasets. Then in Section 3,we analyze the distinct spreading performance and local structures of fast-only influencers and final-only influencers, and discuss the relationship between the network degree assortativity and the fraction of fast-only/final-only influencers. In Section 4, we conclude the main results of this work and reveal insights for practical applications.

    2. Method

    2.1. Susceptible–infected–recovered model

    We employ the susceptible–infected–recovered (SIR)model to simulate the spreading dynamics in networks. In the SIR model,there are three types of individuals:susceptible individuals,infected individuals,and recovered individuals,and there are two model parameters, infect probabilityβand recover probabilityμ. The SIR model supposes that infected nodes infect their neighboring susceptible nodes with probabilityβ, infected individuals recover to the recovered state with probabilityμ, and recovered individuals cannot be infected again. When a seed node is given,we assume that it is infected, and all other nodes are in the susceptible state, then the spreading proceeds according to the above rules. When no individual is found in the infected state,the spreading process ends, and the total number of nodes in the recovered state is measured as the spreading capacity(also referred to as spreading size or spreading effect)of the seed node.

    2.2. Definition of the fast-only and final-only influencers

    Fast influencers are considered as individuals that can influence the top fraction of individuals at the early stage of the spreading process.[7]In comparison,final influencers are those individuals whose spreading capacity is ranked among the top(according to a top ratiox%) at the final stage. In this paper, we focus on the non-overlapping parts between fast influencers and final influencers, which we refer to as fast-only and final-only influencers,respectively. Specifically,fast-only influencers are the fast influencers that do not belong to final influencers. Similarly,the final-only influencers are those final influencers that do not overlap with fast influencers. According to the definitions,fast-only influencers are individuals that achieve high spreading capacity at the early stage but lose the superiority at the final stage, and on the contrary, finalonly influencers are individuals that fail to exhibit a prominent spreading capacity at the early stage but influence a large fraction of nodes at the final stage. According to the same top ratio, the numbers of fast-only and final-only influencers are the same.

    2.3. Assortativity coefficient

    Assortativity coefficient(r)is the Pearson correlation coefficient of degrees between pairs of linked nodes.[26,27]A positiverdenotes that the high(low)-degree nodes tend to connect with high(low)-degree nodes,while a negativerindicates that the high(low)-degree nodes tend to connect with low(high)-degree nodes. In formula, the definition of the assortativity coefficient of a network is as follows:

    wherewkis the distribution of the remaining degree(i.e., the probability that the neighboring node of a randomly-chosen node in the network has degreek) andejkis the joint probability distribution of the remaining degrees of the two nodes at either end of a randomly-chosen edge,andσ2wis the variance of the distributionwk,-1≤r ≤1. Whenr=1, the network is totally assortative(all nodes are only connected with nodes with the same degree)andr=-1 if the network is totally disassortative.

    2.4. Degree configuration model

    To reveal how degree assortativity influences the fraction of fast-only and final-only influencers in the network, we introduce the degree configuration model[28]into the analysis to adjust the degree assortativity of networks. Specifically,for a given network,the model keeps the degree sequence fixed and changes the assortativity of the network.[28–30]According to the degree configuration model,given a network,we randomly choose two edges,(v1,v2)and(v3,v4),which are guaranteed to be different. The four nodes are ranked with respect to their degrees.

    To generate an assortative network from the original one,with probabilityp,we rewire the two edges to let one of them connects the two nodes with smaller degrees and the other link connects the two nodes with larger degrees; with probability 1-p, the two edges are randomly rewired. To avoid multiple edges connecting two same nodes, we will select a new pair of nodes if the rewired edge has already existed in the network. Such a process is repeated until the maximum iterations reaches. In this model,pcontrols the extent of degree assortativity. Ifpis close to 1,after the rewiring process,the generated network tends to be assortative. Similarly, through rewiring network edges with probabilitypto let nodes with larger degrees are connected with nodes with smaller degrees,we can generate networks with tunable disassortativity.

    2.5. Centrality metrics

    Centrality metrics measure the relative importance of a node in the network,which capture a node’s position features in the network structure.In the analysis,we focus on the structure of final-only influencers and fast-only influencers in terms of four metrics: degree,h-index,k-core, and LocalRank. We use aN×Nadjacency matrixAto represent the studied undirected un-weighted network,whereNis the number of nodes in the network,whereAi j=1 if there is a edge between nodeiandj,andAij=0 otherwise.

    Degree Degree is the simplest centrality metric. In a network, nodei’s degreekiis defined as the number of its neighbors,[11]in formulaki=∑j Aij.

    Theh-index Theh-index was firstly introduced to measure the academic impact of scholars[12]based on their publications and the corresponding citations. Then, it was extended to quantify nodes’ significance based on their neighbors’ impact.[31]For a given nodei, itsh-index is defined as the maximum valuehisuch thatihas at leasthineighbors whose degrees are all no less thanhi.

    Thek-core The calculation of nodes’k-core is according to a network pruning process.[2]Specifically,the process starts by removing nodes with degrees equal to one, the pruning is continued until all nodes’degrees are larger than one,then thek-core of the removed nodes is one. Next,the pruning process is repeated for the remaining nodes with degreek=2,3,...,and all nodes will be assigned theirk-core scores.[2]

    LocalRank For a given nodei, its LocalRank[3]scoreLRiis defined as

    whereNkis the number of the nearest and the next nearest neighbors of nodek. LocalRank takes nodes’four-order local structure into account.

    2.6. Empirical datasets

    In the analysis, we use eight social network datasets to support our research. Email is a communication network at the University Rovira i Virgili in Tarragona in Spain,[22]nodes are users and each edge represents that at least one email between the corresponding two users was sent. Cond-mat is a collaboration network of scientists,[33]in which nodes represent scientists and each edge means the two linked scientists have co-authored at least one paper. Lastfm-asia is a social network of LastFM users,[34]nodes are LastFM users from Asian countries and edges are mutual following relationships between them. Musae-facebook is a page-page network of Facebook sites.[35]Nodes represent official Facebook pages (such as pages of politicians, governmental organizations, television shows and companies) and edges are mutual likes between sites. Deezer-HR, deezer-HU and deezer-RO are three friendship networks collected from the music streaming service Deezer from Croatia, Hungary, and Romania, respectively.[36]Nodes represent the users and edges are the mutual friendships. Facebook is a friendship network[37]in which nodes represent Facebook users,and edges represent the mutual friendships. We summarize the characteristics of the eight datasets in Table 1.

    Table 1. The characteristics of the eight studied datasets. Here N is the size of network(i.e.,number of nodes),E is the number of edges,〈k〉is the average degree,λc is the threshold in the SIR model,[38] c is the clustering coefficient,and ρ is the degree assortativity.

    3. Results

    3.1. Overlap between fast influencers and final influencers

    To identify the fast and final influencers and analyze the differences between them,we employ the SIR model[2,39](see Subsection 2.1 for details) to simulate the spreading dynamics. In the simulation,the recover probability holds for a fixed value ofμ=1, and the infect probabilityβchanges fromλcto 8λc. Given nodei, we set it as the seed infected node to start the spreading process and recordi’s spreading capacity(i.e.,number of infected individuals)at timetasqi(t). When the spreading process ends, we denotei’s final spreading capacity asqi. We implement 100 times simulation to obtain average results. After obtaining all nodes’spreading capacity at any given time,we rank all individuals’q(t)and define the topx%(x=1,2,5)as the fast influencers at timet=5. Similarly,ranking all individuals’q,and the topx%(x=1,2,5)of individuals are defined as final influencers.

    Fig. 1. The overlapping fraction between fast influencers and final influencers on four social networks.Different curves in each panel denote the results based on different λ/λc, where λc is the threshold value of SIR model.[2,39] The result shows that in the four networks, with the increase of t,the curves decrease at first,then increase.

    Under different time steps (t), we calculate the overlapping fraction between the fast influencers(top 5%byq(t))and final influencers(top 5%byq),which is denoted as precision(top 5%). Precision represents the fraction of fast influencers that become the final influencers at the final stage. The results on four social networks (see Subsection 2.6 for details) are presented in Fig.1. We find that as the time(t)increases,the overlap between fast influencers and final influencers exhibits an overall decreased trend firstly and then increased trend.The results are consistent under the condition of other top ratios, see the results of top 1% and 2% in Figs. A1 and A2 in Appendix A.And the results on the other four datasets are shown in Figs. A3–A5, which indicates a monotonically increasing trend. The discrepancy could be due to the difference in network size, namely, compared with the four networks in Fig.1,and Deezer and Facebook are larger,in which the same top ratio involves more nodes, which could distort the initial trend of the curve. In all studied datasets,the precision value reaches 1 within 60 time steps,which indicates the spreading ends beforet=60. Consistent with the prior study about fast influencers,[7]we define nodes with the top 5%q(5) as fast influencers in the following analysis, which also corresponds to a low overlap between the two types of influencer.

    Further,we analyze how the spreading parameterλ(λ=β/μ) influences the overlap between the fast influencers and the final influencers. As mentioned above, in the SIR model,we fix the recover probabilityμ=1, and change the infect probabilityβto obtain differentλ. Whenλis large thanλc(the threshold of the SIR model,[5,38]see Table 1), i.e.,λ/λc>1, the spreading outbreaks; whenλis smaller thanλc(λ/λc<1),the spreading will vanish soon. In our simulation,under different infect probabilities,for a given seed nodei,the spreading process is implemented 100 times to obtain average results. We plot the results under four differentλ/λcin Fig.1.The results demonstrate that whent=5,the overlap between the two types of influencers becomes smaller with a largerλ,revealing that in a scenario with high infected probability,one needs to carefully consider the specific type of influencers he wants to identify, because in this case, few final influencers would be fast influencers at the same time. The results are robust in terms of different top ratios and other four datasets(see Figs.A1–A5).

    3.2. Fast-only influencers versus final-only influencers

    Both fast-only influencers and final-only influencers are of importance to practical applications. For example, in advertisement marketing, the marketer sometimes is more concerned about maximizing the advertisement effect at the early stage, thus spending on the fast-only influencers is a better choice to balance the efficiency and cost. On the contrary, if the goal of the marketer is the long-term effect, focusing on final-only influencers is a better choice that will bring a longlasting marketing effect.Therefore,having a deep understanding of the difference between fast-only influencers and finalonly influencers is crucial for making rational use of them.

    3.2.1. Comparison of spreading capacity

    To measure the gap of spreading capacity between the fast-only influencers and the final-only influencers,we calculate the average spreading capacity of the two types of individuals at different time. The results on four datasets are shown in Fig.2(see Fig.A6 for the results of the other four datasets),in which each column denotes the results based on differentλ/λc. In the four datasets, we observed consistent patterns.Firstly, by their definitions, final-only influencers underperform fast-only influencers at the early stage, but they can exhibit high spreading capacity in the final stage. Secondly,fastonly influencers perform well at the early stage but lose their superiority at the late time. Thirdly,the average spreading capacity of fast-only influencers and final-only influencers is always better than that of the random-chosen(5%)individuals.

    Fig.2. The spreading capacity of fast-only influencers,final-only influencers and random-chosen individuals. The x-axis is the time t,and the y-axis is the average spreading capacity of the three types of nodes. Different columns denote results based on different λ/λc.

    Fig.3.The average local metric scores of fast-only influencers and final-only influencers on the four social networks.The local metrics include degree,h-index,k-core and LocalRank. The x-axis denotes four different λ/λc, and the y-axis denotes the average centrality scores of two types of influencers. In general,the results indicate that compared with final-only influencers,fast-only influencers have more strong local structures,as λ/λc increases,the score gap between fast-only influencers and final-only influencers becomes wider.

    Finally, the last intersection between the curves of finalonly influencers and fast-only influencers (the time that the spreading capacity of the former exceeds that of the latter)gradually moves backward as the increase ofλ, which indicates that under a condition of high spreading probability,final-only influencers need more time to exceed fast-only influencers.

    3.2.2. Comparison of local structure via centrality metrics

    Next, we compare the local structure of fast-only influencers and final-only influencers by computing their average centrality metric scores(see Subsection 2.5). Figure 3 shows the results of fast-only and final-only influencers. The metrics include degree,h-index,k-core and LocalRank(the results on the other four datasets are shown in Appendix Fig.A7). From Fig. 3, we find that compared with the final-only influencers,the fast-only influencers have higher scores with respect to all considered four metrics. Similar results are observed on the other four datasets, which reveal that nodes with strong local structures tend to have a powerful ability to infect more nodes at the early stage. The results indicate that to maximize the spreading effect in a limited time,we need to emphasize more on individuals’ local structure than global structure, which is consistent with the finding in previous work.[7]On the contrary, to achieve a long-time spreading effect, we should not only consider the local position of individuals, but also consider the global position of individuals.[5]

    Figure 3 also demonstrates that asλincreases,the score gap between fast-only influencers and final-only influencers becomes larger,which is consistent with the results in Fig.1,namely, in the case of large infected probability, the discrepancy between fast-only influencers and final-only influencers becomes larger, non-top-spreading individuals at the early stage have more chance to become final influencers in the final stage.

    3.3. The impact of degree assortativity on the fraction of fast-only/final-only influencers

    Previous research has claimed that the degree assortativity of networks influences the spreading dynamics.[23–25]In this section,we investigate how the variation of degree assortativity influences the fast-only and final-only influencers. To this end, we adjust the degree assortativity of given networks by the configuration model,[28]which keeps the degree sequence fixed and changes the degree assortativity of networks by rewiring the edges(see Subsection 2.4 for details). To obtain robust results, for each dataset, we run the configuration model 10 times and analyze the average results.

    Recall that in the configuration model, parameterpcontrols the extent of the network assortativity or disassortativity.In Figs. 4(a) and 4(b), we show the relationship betweenpand the assortativityron four studied networks. Whenr >0,the network is assortative, which indicates that nodes tend to be connected with nodes with similar degree values,otherwiser <0. Taking the result on email network as an example, in panel a, we can see that whenpincreases from 0 to 1, the assortativity of the network rises from-0.03 to 0.51,indicating that the network tends to be more assortative. Similarly,in panel b, withpincreases from 0 to 1, the assortativity of email network declines from-0.03 to-0.37, which means the network tends to be more disassortative.

    Fig.4. The influence of degree assortativity on the fraction of the fast-only and final-only influencers. In panels (a) and (b), the x-axis represents the parameter p in the configuration model,and the y-axis is the degree assortativity (r) of the network. In panel (a), with the increase of p, the curves go up, indicating that the generated networks become more assortative. Similarly, in panel (b), when p rises from 0 to 1, the curves go down, meaning that the generated networks become more disassortative. In panels (c) and(d),the y-axis denotes the ratio of the number of fast-only influencers to the number of fast influencers. In panel (c), when p increases, the fraction of fast-only influencers decreases on the four studied networks. In panel (d),when p increases, we do not observe a consistent pattern. Specifically, in cond-mat, lastfm-asia and muase-facebook, the increase of p is associated with the decrease of fast-only influencers,but in email network,the opposite is true.

    We further study the relation between degree assortativity and the fraction of the fast-only and final-only influencers.For each generated network by the configuration model,[28]we run the SIR model and record all individuals’spreading capacity att=5 and at the end.Then we calculate the fraction of fast-only influencers. Figure 4(c)demonstrates that with the increase ofp,the fraction of fast-only influencers decreases. Such a result indicates that when a network becomes more assortative, fast influencers are more likely become final influencers, namely,ordinary individuals that did not perform well at the early stage can hardly achieve a top spreading effect at the final stage,and individuals in the network can hardly break the hierarchy in information spreading,which is unfair for most individuals.The reason for such a result is that if most nodes tend to connect with nodes with similar degrees,individuals with low degrees can hardly transmit information to high-degree nodes. As for Fig.4(d),it shows that when the network becomes more disassortative, the variation of the fraction of fast-only influencers is distinct on different datasets. In lastfm-asia, cond-mat and musae-facebook datasets,the fraction of fast-only influencers decreases as the networks become more disassortative,but for the email dataset, with the increase ofp, the fraction of fastonly influencers increases.

    Following we discuss the possible reasons for this discrepancy. Based on the analysis of Fig. 4(b), we notice that compared with the other three networks, the email network is more disassortative withpincreases, indicating it involves a larger fraction of high-degree nodes that connect with lowdegree nodes for a givenp. And in our paper, the defined fast-only influencers are individuals with relatively high degree(see Fig.2 and Table 1 for details),which leads to the case that with the increase of degree disassortativity, comparing with other three networks,in email network,there have more high-degree nodes tend to connect with low-degree nodes,and the fractions of fast-only and final-only influencers increase as well. While for the other three networks, there are relatively fewer high-degree nodes that connect with low-degree nodes,the fraction of fast-only and final-only influencers would not increase. The results with respect to the other four datasets are shown in Appendix Fig. A8, which demonstrates that the deezer-HU and deezer-RO networks show similar results(the more disassortative the network is,the more chance that ordinary individuals have to become influencers finally).

    In summary, our results indicate that for assortative networks,it is hard for those individuals that perform non-top at the early stage to achieve a breakthrough at the final stage,and individuals who perform well at the early stage have more chance to keep their superiority at the final stage.

    4. Conclusion

    Vital nodes play an important role in the structure and function of networks. In this paper, we distinguish two types of vital nodes: fast influencers and final influencers. We study the overlap between the two types of influencers at different time and spreading parameterλ. Then we focus on the non-overlapping parts of fast influencers and final influencers, which we call the fast-only influencers and final-only influencers, respectively, and we analyze the differences between them. We study howλinfluences the spreading capacity of the fast-only influencers and final-only influencers,and the result shows that increasingλwill make those individuals with non-top spreading performance have more chance to have top performance at the final stage. Then we compare the local structures of fast-only influencers and final-only influencers through centrality metrics, and find that on the studied eight social networks,compared with final-only influencers, the fast-only influencers have a stronger local structure, represented by higher centrality scores. Finally, we discuss how the degree assortativity influences the fraction of fast-only influencers. The result indicates that when the network becomes more assortative,individuals that perform well at the early stage will keep their superiority at the final stage.And mediocre individuals at the early time in the network can hardly achieve a high spreading effect at the final stage. Our research shows that if people focus on the effect at the early stage, they need pay more attention to fast-only influencers to reach a balance between effect and cost. If they are more concerned about the long-time effect, they need to focus on final-only influencers to achieve a long-lasting spreading effect. Our research can also shed the light on the sociology field, to keep the social system fairer, we need to make the society more disassortative rather than assortative.

    Acknowledgements

    L.L. acknowledges the Science Strength Promotion Programme of UESTC,Chengdu.

    Project supported by the National Natural Science Foundation of China (Grant Nos. 61673150 and 11622538) and Special Project for the Central Guidance on Local Science and Technology Development of Sichuan Province,China(Project No.2021ZYD0029).

    Appendix A: Overlap between fast influencers and final influencers for different top ratios and datasets

    Fig.A1. The overlap between fast influencers and final influencers on four social networks(top 1%).

    Fig.A2. The overlap between fast influencers and final influencers on four social networks(top 2%).

    Fig.A3. The overlap between fast influencers and final influencers on deezer-HR,deezer-HU,deezer-RO,Facebook networks(top 1%).

    Fig.A4. The overlap between fast influencers and final influencers on deezer-HR,deezer-HU,deezer-RO,Facebook networks(top 2%).

    Fig.A5. The overlap between fast influencers and final influencers on deezer-HR,deezer-HU,deezer-RO,Facebook networks(top 5%).

    Appendix B:The comparison between fast-only and final-only influencers in terms of spreading capacity at different time

    Fig.B1. The spreading capacity of the fast-only influencers,final-only influencers and random-chosen(5%)individuals on deezer-HR,deezer-HR, deezer-HR, Facebook networks. The x-axis is the time t, and the y-axis is the average spreading capacity of the three types of nodes.Different columns denote results based on different spreading parameters λ. Similar to Fig.3 in the main text,but for the other four datasets.

    Fig.B1. The spreading capacity of the fast-only influencers,final-only influencers and random-chosen(5%)individuals on deezer-HR,deezer-HR,deezer-HR, Facebook networks. The x-axis is the time t, and the y-axis is the average spreading capacity of the three types of nodes. Different columns denote results based on different spreading parameters λ. Similar to Fig.3 in the main text,but for the other four datasets(continue).

    Appendix C:Centrality metrics scores of fast-only and final-only influencers on four social networks

    Fig.C1. The average local metric scores of the fast-only influencers and final-only influencers on the other four social networks. Similar to Fig.4 in the main text,but for the other four datasets.

    Appendix D:The impact of degree assortativity on the fraction of fast/final-only influencers

    Fig.D1. The influence of degree assortativity on the fraction of the fast-only and final-only influencers. Similar to Fig.4 in the main text,but for the other four datasets.

    夜夜骑夜夜射夜夜干| 亚洲综合色网址| 精品久久久久久电影网| 亚洲av欧美aⅴ国产| 亚洲成av片中文字幕在线观看 | 成年动漫av网址| 国产欧美日韩一区二区三区在线| 高清欧美精品videossex| 有码 亚洲区| 亚洲精品日韩在线中文字幕| 国产精品二区激情视频| 丝袜喷水一区| 丝袜在线中文字幕| 日韩 亚洲 欧美在线| 亚洲精品在线美女| www.自偷自拍.com| 日韩 亚洲 欧美在线| 亚洲国产最新在线播放| 男女国产视频网站| 视频区图区小说| 波野结衣二区三区在线| 国产精品香港三级国产av潘金莲 | 成人亚洲精品一区在线观看| 国产又色又爽无遮挡免| 午夜福利一区二区在线看| 久久av网站| 黄色 视频免费看| 我要看黄色一级片免费的| 中国国产av一级| 亚洲国产av影院在线观看| 国产 一区精品| 香蕉丝袜av| 国产精品不卡视频一区二区| 精品亚洲乱码少妇综合久久| 97精品久久久久久久久久精品| 国产片特级美女逼逼视频| 校园人妻丝袜中文字幕| 中文字幕色久视频| 2022亚洲国产成人精品| 波野结衣二区三区在线| 久久久久国产网址| 成年女人在线观看亚洲视频| 中文字幕亚洲精品专区| 亚洲一区中文字幕在线| 亚洲,一卡二卡三卡| 日韩av免费高清视频| 日本午夜av视频| 边亲边吃奶的免费视频| 99国产精品免费福利视频| www.自偷自拍.com| 少妇人妻精品综合一区二区| 丝袜美足系列| 制服人妻中文乱码| 中国三级夫妇交换| 亚洲精品,欧美精品| 叶爱在线成人免费视频播放| 免费在线观看完整版高清| av在线老鸭窝| 999精品在线视频| 狂野欧美激情性bbbbbb| 亚洲国产欧美网| 久久影院123| 久久久久国产一级毛片高清牌| 国产精品一国产av| 亚洲精品久久久久久婷婷小说| 熟女电影av网| 欧美+日韩+精品| 少妇精品久久久久久久| 999精品在线视频| 亚洲国产成人一精品久久久| 久久久久精品性色| 丝袜美足系列| 亚洲av电影在线观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 亚洲国产av新网站| 只有这里有精品99| 欧美xxⅹ黑人| 欧美激情极品国产一区二区三区| h视频一区二区三区| tube8黄色片| 18禁动态无遮挡网站| 少妇人妻精品综合一区二区| 中文字幕精品免费在线观看视频| 黄色怎么调成土黄色| 天美传媒精品一区二区| 多毛熟女@视频| 美女国产视频在线观看| 日本午夜av视频| 2022亚洲国产成人精品| 亚洲国产精品一区二区三区在线| 人人妻人人澡人人爽人人夜夜| 亚洲,欧美,日韩| 日韩不卡一区二区三区视频在线| 亚洲久久久国产精品| 一区福利在线观看| 国产乱人偷精品视频| 18禁动态无遮挡网站| 韩国av在线不卡| 久久久久精品性色| 欧美 亚洲 国产 日韩一| 久久韩国三级中文字幕| 久久影院123| 国产亚洲午夜精品一区二区久久| 中国三级夫妇交换| 日本wwww免费看| 精品一区二区免费观看| 波多野结衣一区麻豆| 日产精品乱码卡一卡2卡三| 男的添女的下面高潮视频| 久久人妻熟女aⅴ| 男女下面插进去视频免费观看| 欧美最新免费一区二区三区| 精品少妇黑人巨大在线播放| 久久久欧美国产精品| 久久精品aⅴ一区二区三区四区 | 看免费成人av毛片| 欧美日韩综合久久久久久| 一级,二级,三级黄色视频| 免费观看在线日韩| 国产精品一区二区在线观看99| 99re6热这里在线精品视频| 亚洲国产精品999| 国产野战对白在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产精品蜜桃在线观看| 18禁观看日本| 妹子高潮喷水视频| 久久久久国产精品人妻一区二区| 久久韩国三级中文字幕| 久久久久久久精品精品| 丰满乱子伦码专区| av福利片在线| 两性夫妻黄色片| 国产成人精品久久二区二区91 | 亚洲精品成人av观看孕妇| av国产久精品久网站免费入址| 日日爽夜夜爽网站| 看免费av毛片| 欧美日韩视频高清一区二区三区二| 欧美国产精品一级二级三级| 亚洲av.av天堂| 国产精品亚洲av一区麻豆 | 亚洲伊人久久精品综合| 成人影院久久| 国产精品久久久久久av不卡| 日本欧美国产在线视频| 免费看av在线观看网站| 一二三四在线观看免费中文在| 狂野欧美激情性bbbbbb| 一二三四在线观看免费中文在| 亚洲av在线观看美女高潮| 日韩一区二区视频免费看| 欧美+日韩+精品| av网站免费在线观看视频| 日韩精品有码人妻一区| 国产亚洲午夜精品一区二区久久| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区| 亚洲欧美色中文字幕在线| 老汉色av国产亚洲站长工具| 寂寞人妻少妇视频99o| 亚洲天堂av无毛| 自拍欧美九色日韩亚洲蝌蚪91| 丰满迷人的少妇在线观看| 亚洲精品国产av成人精品| 人妻系列 视频| 巨乳人妻的诱惑在线观看| 亚洲男人天堂网一区| 久久久久久久久久人人人人人人| 久久久久久久久久久久大奶| 亚洲美女视频黄频| 成人18禁高潮啪啪吃奶动态图| 中文字幕精品免费在线观看视频| 精品卡一卡二卡四卡免费| 国产精品99久久99久久久不卡 | 久久99热这里只频精品6学生| 亚洲国产最新在线播放| 一级黄片播放器| 亚洲国产成人一精品久久久| 国产探花极品一区二区| 国产熟女欧美一区二区| 26uuu在线亚洲综合色| 一级毛片电影观看| 午夜激情av网站| 国产一区有黄有色的免费视频| 日本黄色日本黄色录像| 久久久久久久久久人人人人人人| 91午夜精品亚洲一区二区三区| 亚洲精品久久久久久婷婷小说| 免费播放大片免费观看视频在线观看| 亚洲av福利一区| 我的亚洲天堂| 老熟女久久久| 婷婷色综合www| 国产亚洲精品第一综合不卡| 久热这里只有精品99| 黑丝袜美女国产一区| 午夜av观看不卡| 日韩三级伦理在线观看| 高清不卡的av网站| 欧美国产精品va在线观看不卡| 激情视频va一区二区三区| 18禁观看日本| 国产精品香港三级国产av潘金莲 | 99热网站在线观看| 亚洲欧洲日产国产| 精品卡一卡二卡四卡免费| 精品少妇一区二区三区视频日本电影 | 大陆偷拍与自拍| 日韩av不卡免费在线播放| 日本色播在线视频| 女人精品久久久久毛片| 久久综合国产亚洲精品| 人妻人人澡人人爽人人| 亚洲三区欧美一区| 国产又爽黄色视频| 久久久久久久久久人人人人人人| 久久精品亚洲av国产电影网| 亚洲国产av新网站| 成人漫画全彩无遮挡| videossex国产| 国产 精品1| 亚洲少妇的诱惑av| 亚洲精品一二三| 寂寞人妻少妇视频99o| 少妇人妻 视频| 观看美女的网站| www日本在线高清视频| 三上悠亚av全集在线观看| 午夜免费观看性视频| 国产成人午夜福利电影在线观看| 少妇的逼水好多| a级毛片在线看网站| 免费播放大片免费观看视频在线观看| 久久影院123| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品古装| 在线观看免费视频网站a站| 亚洲男人天堂网一区| 亚洲成人av在线免费| 成人国产麻豆网| 午夜激情久久久久久久| 色视频在线一区二区三区| 最近的中文字幕免费完整| 日本91视频免费播放| 精品国产露脸久久av麻豆| 搡女人真爽免费视频火全软件| 亚洲情色 制服丝袜| 亚洲国产精品成人久久小说| 亚洲伊人色综图| 精品酒店卫生间| 制服诱惑二区| 91在线精品国自产拍蜜月| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 丝袜美腿诱惑在线| 久久免费观看电影| 七月丁香在线播放| 9191精品国产免费久久| 欧美亚洲日本最大视频资源| 久久毛片免费看一区二区三区| 亚洲一码二码三码区别大吗| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品av麻豆狂野| 丝袜在线中文字幕| 国产精品国产av在线观看| 满18在线观看网站| 女人精品久久久久毛片| 精品国产国语对白av| 日韩一区二区三区影片| 成年av动漫网址| 午夜福利一区二区在线看| 永久免费av网站大全| 欧美日韩精品成人综合77777| www.精华液| 日韩欧美精品免费久久| 性色av一级| 精品人妻在线不人妻| 色网站视频免费| 男人操女人黄网站| √禁漫天堂资源中文www| 好男人视频免费观看在线| 美女视频免费永久观看网站| 精品一品国产午夜福利视频| 90打野战视频偷拍视频| 亚洲精品,欧美精品| 国产男女内射视频| 午夜福利影视在线免费观看| 久久久久久久大尺度免费视频| 国产精品不卡视频一区二区| 高清黄色对白视频在线免费看| av片东京热男人的天堂| 色视频在线一区二区三区| 侵犯人妻中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 亚洲成人av在线免费| 在线观看免费日韩欧美大片| 午夜激情久久久久久久| 亚洲欧美一区二区三区久久| 纯流量卡能插随身wifi吗| 伦精品一区二区三区| 久久久亚洲精品成人影院| 精品午夜福利在线看| 色婷婷久久久亚洲欧美| 三上悠亚av全集在线观看| 日韩在线高清观看一区二区三区| 日本av免费视频播放| 日韩制服骚丝袜av| 黄色怎么调成土黄色| 一级黄片播放器| 欧美变态另类bdsm刘玥| 国产成人精品无人区| 日韩三级伦理在线观看| 国产免费一区二区三区四区乱码| 成人二区视频| 1024香蕉在线观看| 亚洲欧美精品自产自拍| 日韩制服骚丝袜av| av卡一久久| 久久久久久人妻| 看免费成人av毛片| 免费高清在线观看视频在线观看| www.精华液| 久久久久精品性色| 天天躁狠狠躁夜夜躁狠狠躁| 大陆偷拍与自拍| 国产成人91sexporn| 啦啦啦视频在线资源免费观看| 亚洲精品成人av观看孕妇| 国产日韩欧美在线精品| 久久久久久久久久人人人人人人| 久久精品国产自在天天线| 黄片播放在线免费| 国产亚洲一区二区精品| 欧美变态另类bdsm刘玥| 欧美精品人与动牲交sv欧美| 午夜激情久久久久久久| 黄网站色视频无遮挡免费观看| 亚洲国产看品久久| 色网站视频免费| 国产一区二区在线观看av| 免费观看无遮挡的男女| 欧美 日韩 精品 国产| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 男女国产视频网站| 亚洲欧美一区二区三区国产| 日韩欧美一区视频在线观看| 午夜福利网站1000一区二区三区| 一级片'在线观看视频| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 最近最新中文字幕大全免费视频 | 免费黄色在线免费观看| 考比视频在线观看| 久久久亚洲精品成人影院| 色哟哟·www| 不卡av一区二区三区| 国产淫语在线视频| 纵有疾风起免费观看全集完整版| 中文精品一卡2卡3卡4更新| 一二三四中文在线观看免费高清| 国产精品久久久久久精品古装| 十八禁网站网址无遮挡| 亚洲天堂av无毛| 日日摸夜夜添夜夜爱| 九九爱精品视频在线观看| 午夜福利在线免费观看网站| av免费观看日本| 免费观看无遮挡的男女| 国产深夜福利视频在线观看| 国产精品.久久久| 成人国产av品久久久| 天天影视国产精品| 久久人妻熟女aⅴ| 一级片'在线观看视频| 亚洲欧美精品综合一区二区三区 | 亚洲欧美精品自产自拍| 午夜激情久久久久久久| 少妇 在线观看| 国产成人一区二区在线| 婷婷色综合大香蕉| 大片免费播放器 马上看| 新久久久久国产一级毛片| 一级毛片黄色毛片免费观看视频| av又黄又爽大尺度在线免费看| 99久国产av精品国产电影| 丝袜喷水一区| 纵有疾风起免费观看全集完整版| 午夜激情久久久久久久| 欧美国产精品va在线观看不卡| 一区二区日韩欧美中文字幕| 街头女战士在线观看网站| 免费观看无遮挡的男女| 97在线人人人人妻| 亚洲av日韩在线播放| 成人漫画全彩无遮挡| 中文字幕制服av| 精品国产露脸久久av麻豆| av网站在线播放免费| 少妇熟女欧美另类| 国产一区二区激情短视频 | 国产熟女午夜一区二区三区| 成人精品一区二区免费| 免费看a级黄色片| 国产精品久久久久久人妻精品电影| 中文字幕最新亚洲高清| 亚洲欧美激情在线| 成人免费观看视频高清| 国产亚洲欧美98| 水蜜桃什么品种好| 亚洲av美国av| 18禁黄网站禁片午夜丰满| 女性被躁到高潮视频| 少妇被粗大的猛进出69影院| 欧美精品一区二区免费开放| 一级毛片高清免费大全| 亚洲自拍偷在线| 久久人妻熟女aⅴ| 国产高清videossex| 亚洲中文字幕日韩| 久久久久久久午夜电影 | 三上悠亚av全集在线观看| 性色av乱码一区二区三区2| 久久香蕉激情| www日本在线高清视频| 久久精品aⅴ一区二区三区四区| 真人一进一出gif抽搐免费| 黑人操中国人逼视频| av视频免费观看在线观看| 天天躁夜夜躁狠狠躁躁| 日韩精品青青久久久久久| 久久久久久久久久久久大奶| 丁香六月欧美| 国产又爽黄色视频| 亚洲五月婷婷丁香| 亚洲午夜理论影院| 国产精品香港三级国产av潘金莲| 欧美成狂野欧美在线观看| 成人永久免费在线观看视频| 国产一卡二卡三卡精品| 美女 人体艺术 gogo| 亚洲欧美精品综合一区二区三区| 亚洲熟妇中文字幕五十中出 | 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 夫妻午夜视频| www.自偷自拍.com| 日韩人妻精品一区2区三区| 中文字幕av电影在线播放| 99re在线观看精品视频| 国产精华一区二区三区| 一级毛片精品| 在线观看免费午夜福利视频| 亚洲欧美激情在线| 少妇裸体淫交视频免费看高清 | 国产黄色免费在线视频| 搡老熟女国产l中国老女人| 国产免费av片在线观看野外av| bbb黄色大片| 黑人巨大精品欧美一区二区mp4| 午夜福利在线免费观看网站| 丰满人妻熟妇乱又伦精品不卡| 天天躁夜夜躁狠狠躁躁| 电影成人av| 亚洲午夜理论影院| 欧美激情 高清一区二区三区| 亚洲一区高清亚洲精品| 超碰97精品在线观看| av网站免费在线观看视频| 国产成人免费无遮挡视频| 国产蜜桃级精品一区二区三区| 成人黄色视频免费在线看| 欧美不卡视频在线免费观看 | 免费av毛片视频| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 成人国产一区最新在线观看| 日韩欧美免费精品| 欧美最黄视频在线播放免费 | 大陆偷拍与自拍| 久久国产精品人妻蜜桃| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 女性被躁到高潮视频| 国产一区二区在线av高清观看| 熟女少妇亚洲综合色aaa.| 久久欧美精品欧美久久欧美| 高清毛片免费观看视频网站 | 亚洲专区中文字幕在线| a级毛片在线看网站| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 国产单亲对白刺激| 女人被躁到高潮嗷嗷叫费观| 免费搜索国产男女视频| av电影中文网址| 免费高清视频大片| 国产激情久久老熟女| 免费看a级黄色片| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 90打野战视频偷拍视频| 最近最新中文字幕大全电影3 | 国产av精品麻豆| 精品国产乱子伦一区二区三区| 久热这里只有精品99| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲一级av第二区| 久久精品国产亚洲av高清一级| 亚洲av成人一区二区三| 老汉色av国产亚洲站长工具| 黄色丝袜av网址大全| 波多野结衣av一区二区av| 国产欧美日韩综合在线一区二区| 久久国产精品影院| cao死你这个sao货| 91九色精品人成在线观看| 19禁男女啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| 国产xxxxx性猛交| 国产麻豆69| 成人av一区二区三区在线看| 国产精品久久视频播放| 欧美乱色亚洲激情| 成人亚洲精品一区在线观看| 日韩大尺度精品在线看网址 | 日韩 欧美 亚洲 中文字幕| 热99re8久久精品国产| 热re99久久精品国产66热6| 一级,二级,三级黄色视频| 日韩国内少妇激情av| 亚洲av日韩精品久久久久久密| 久久影院123| 精品人妻在线不人妻| 狂野欧美激情性xxxx| 国产av在哪里看| 国产精品秋霞免费鲁丝片| 丝袜在线中文字幕| 国产精品成人在线| 女人被狂操c到高潮| 亚洲第一av免费看| 午夜久久久在线观看| 另类亚洲欧美激情| 亚洲情色 制服丝袜| 亚洲精品美女久久久久99蜜臀| 欧美久久黑人一区二区| 国产av一区二区精品久久| 一级片'在线观看视频| 国产野战对白在线观看| 99久久国产精品久久久| 999精品在线视频| 国产精华一区二区三区| 在线观看66精品国产| 老鸭窝网址在线观看| 亚洲国产精品合色在线| 久久人妻熟女aⅴ| 午夜福利影视在线免费观看| 视频区欧美日本亚洲| 免费观看精品视频网站| 99精品欧美一区二区三区四区| 国产精品一区二区在线不卡| 久久精品aⅴ一区二区三区四区| 中出人妻视频一区二区| av天堂在线播放| 精品一区二区三区四区五区乱码| 91字幕亚洲| 国产精品秋霞免费鲁丝片| av中文乱码字幕在线| 老司机午夜福利在线观看视频| 国产又色又爽无遮挡免费看| 男女高潮啪啪啪动态图| 1024香蕉在线观看| 高潮久久久久久久久久久不卡| 国产一区二区三区综合在线观看| 久久 成人 亚洲| 99国产综合亚洲精品| 欧美国产精品va在线观看不卡| 夜夜爽天天搞| 天堂俺去俺来也www色官网| 久久久国产成人免费| 韩国精品一区二区三区| 成人国语在线视频| 美女午夜性视频免费| 伦理电影免费视频| 99在线人妻在线中文字幕| 亚洲午夜理论影院| 在线看a的网站| 黑人巨大精品欧美一区二区mp4| 欧美精品一区二区免费开放| 男人舔女人的私密视频| 欧美中文日本在线观看视频| 久久精品国产亚洲av香蕉五月| svipshipincom国产片| 99热只有精品国产| 日日夜夜操网爽| av在线天堂中文字幕 | 嫩草影视91久久| 97超级碰碰碰精品色视频在线观看| 新久久久久国产一级毛片| 亚洲成国产人片在线观看| 亚洲第一av免费看| 久久久国产一区二区| 国产av一区二区精品久久| 久久人妻av系列| 国产一区二区在线av高清观看| 一区二区三区国产精品乱码| 成在线人永久免费视频| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 国产乱人伦免费视频| 国产色视频综合| 亚洲av成人不卡在线观看播放网|