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

    Identifying vital edges in Chinese air route network via memetic algorithm

    2017-11-21 12:54:14DuWenboLiangBoyuanYanGangOriolLoranCaoXianbin
    CHINESE JOURNAL OF AERONAUTICS 2017年1期
    關(guān)鍵詞:醋酸鈉緩沖溶液蒸餾水

    Du Wenbo,Liang Boyuan,Yan Gang,Oriol Loran,Cao Xianbin,*

    aSchool of Electronic and Information Engineering,Beihang University,Beijing 100191,China

    bBeijing Key Laboratory for Network-based Cooperative Air Traffic Management,Beijing 100191,China

    cSchool of Physics Science and Engineering,Tongji University,Shanghai 200092,China

    dUniversitat Polite`cnica de Catalunya-BarcelonaTech,C/Colom no.11,Terrassa 08222,Spain

    Identifying vital edges in Chinese air route network via memetic algorithm

    Du Wenboa,b,Liang Boyuana,b,Yan Gangc,Oriol Lordand,Cao Xianbina,b,*

    aSchool of Electronic and Information Engineering,Beihang University,Beijing 100191,China

    bBeijing Key Laboratory for Network-based Cooperative Air Traffic Management,Beijing 100191,China

    cSchool of Physics Science and Engineering,Tongji University,Shanghai 200092,China

    dUniversitat Polite`cnica de Catalunya-BarcelonaTech,C/Colom no.11,Terrassa 08222,Spain

    Air route network;Air transport network;Memetic algorithm;Robustness;Vital edges

    Due to rapid development in the past decade,air transportation system has attracted considerable research attention from diverse communities.While most of the previous studies focused on airline networks,here we systematically explore the robustness of the Chinese air route network,and identify the vital edges which form the backbone of Chinese air transportation system.Specifically,we employ a memetic algorithm to minimize the network robustness after removing certain edges,and hence the solution of this model is the set of vital edges.Counterintuitively,our results show that the most vital edges are not necessarily the edges of the highest topological importance,for which we provide an extensive explanation from the microscope view.Our findings also offer new insights to understanding and optimizing other real-world network systems.

    1.Introduction

    With the increasing people and goods transport demand during the accelerating globalization process,the air transportation system plays a more important role than ever before due to its high-speed and high-security advantages.For example,the air transport volume of China grows at an average annual speed of over 10%in the past decades,and now it possesses over one seventh of the total comprehensive transport volume(including roadways,railways,shipping and air transport),which was only 7.9%in 2000.Hence the air transportation system has been drawing much attention from different research communities.Oneofthe mostinteresting directions is to analyze the structure and function of air transportation systems within the framework of complex network theory.

    The air transportation system can be represented as a network,in which nodes denote airport and an edge will be created if there is a direct flight between two airports.In the vast majority of previous literature,the air transport network(ATN)was primarily classified into two scales:worldwide and national.

    For the worldwide scale,Amaral et al.firstly found that worldwide ATN is a small-world network with a power-lawdegree distribution,and the highest-degree airport is not necessarily the most central node,prompting them to propose a network model where both geographical and political factors are taken into account.1,2Barrat et al.investigated the worldwide ATN from a perspective of complex weighted networks and found the nonlinear positive correlation between flight flow and topology properties.3,4They proposed a weighted network model,enlightening the understanding of weighted feature of complex systems.Verma et al.decomposed the worldwide ATN into three distinct layers via k-core decomposition and found that this network is robust to the removal of long distance edges,but fragile to the disconnectivity of short and apparently insignificant edges.5,6

    For the national scale,ATNs of several major nations,such as US,Brazil,India and China,are extensively studied3,7–11,and the national ATNs usually exhibit different features from the worldwide ATN.Gautreau et al.studied US ATN during 1990–2000.3A remarkable result they presented is that although most statistical properties are stationary,an intense activity takes place at the local level.Fleurquin et al.proposed a delay propagation model via quantifying the network congestion for US ATN,revealing that even under normal operating condition the systemic instability risk is non-negligible.11Rocha investigated the Brazilian ATN during 1995–2006,and found that it shrank in topology but grew in traffic volume.7Bagler studied the Indian ATN,and found its signature of hierarchy feature.12As the most active economy,the Chinese aviation industry ranks second to US in the past decade and keeps a high increase rate.Consequently,Chinese ATN attracts continuous attention in different aspects from topology to dynamics and evolution,8–10,13,14one of which is to investigate the backbone of ATN,the air route network(ARN).

    ATN is actually a logic network with origin-destination(OD)relationships.In real air traffic operation,a flight does not straightly fly from departure airport to landing airport,but along some air route waypoints.ARN consists of air route waypoints and connections between them.In 2012,Cai et al.firstly investigated the Chinese ARN15and found that the degree distribution of Chinese ARN is homogeneous but the traf fic flow is rather heterogeneous.Vitali et al.then investigated the horizontal deviation and delays in Italian ARN.16The analysis of ARN is quite a novelty in the literature.However,the network robustness,which is an important issue for infrastructure systems17and has been extensively studied in ATN,18,19is still rare in ARN.In the typical network robustness model,edges are removed by different targeted attack strategies and the size of giant component estimates the robustness of the network.20When a small amount of edges are removed,the size of giant component is of a very small change.In this paper,we focus on identifying the vital edges in Chinese ARN by examining the robustness of the new network after removing an edge set via memetic optimization.Remarkably,we find that the most vital edges are not necessarily the edges of the highest topological importance.

    The rest of this paper is organized as follows.In the next section,we demonstrate Chinese air route network and its basic properties.Section 3 describes the optimization model and the memetic algorithm.Section 4 presents the simulation results and corresponding analysis.Finally,the paper is concluded in Section 5.

    2.Chinese air route network

    In Ref.15,the authors found that the topology structure of the Chinese ARN is homogeneous,yet its distribution of flight flow is quite heterogeneous.If we compare the Chinese ATN with the Chinese ARN,we found signi ficant differences.On one hand,the Chinese ATN is a typical small-world with low average shortestpath length and large clustering coefficient.On the other hand,the Chinese ARN is not a small-world network due to its low clustering coefficient,large average shortest path length and exponential spatial distance distribution.

    3.Model

    3.1.Optimization model

    The static robustness of complex networks has been extensively studied in the past decades.In Ref.21,it is quantified by the relative size of the largest connected componentG=N′/NwhereNis the total number of nodes in initial network andN′is the number of nodes in the largest component after attack.The larger value ofGrepresents a more robust network.Based on the largest connected component,Schneider et al.proposed a measureRto evaluate the robustness against targeted attack on nodes.17

    wheres(Q)is the fraction of nodes in the largest component after removingQnodes.For calculating the robustness of a network,we will follow a degree adaptive strategy:the highest degree nodes will be systematically removed one by one.It is a more comprehensive measure of network robustness.Obviously,a network with higherRhas a stronger resistance to targeted attacks.

    In the Chinese ARN,the closure of air route segment will decrease the connectivity of the whole network.If the vital edges can be recognized,we can prevent the cascading effect induced by the remove of the edges.It is of great significance to identify vital edges that lead to the vulnerability of Chinese ARN.In Ref.22,Freeman proposed a global metric edgebetweenness to measure the importance of an edge,which can identify influential edges effectively.It is defined as follows:

    where Γ is the set of nodes,njkthe number of the shortest paths fromjtok,andnjk(i)the number of the shortest paths fromjtokvia edgei.

    In this work,we formulate a combinational optimization problem to identify the vital edges within network robust model in the Chinese ARN.The objective is to minimize the network robustness after removing certain edges,i.e.closing certain air route segments.Therefore,these edges play an important role in maintaining network robustness.The optimization model is formulated as follows:

    A樣品空白對照——25 μl樣品+25 μl底物(37 ℃孵育10 min)+50 μl蒸餾水(37 ℃下孵育 60 min)+100 μl醋酸-醋酸鈉緩沖溶液;

    whereVis the total number of edges in the network and x is aVdimensional binary variable,e(k)∈ {0,1} (e(k)represents thekth edge in the network).e(k)=1 represents that the edgee(k)is removed,otherwise it remains in the network.The total number of removed edges is certain and denoted as cost(C).Thus,we can identify critical edges in the network and minimize its robustness.

    3.2.Memetic algorithm

    For solving this optimization model,we will use the memetic algorithm(MA),a useful tool for dealing with large-scale combinational problem.23–26Coming from the concept of meme,MA is defined as a part of local improvement in the process of cultural evolution.It is a hybrid metaheuristic of global search and heuristic local search with three operations:crossover,local search and tournament selection.

    (1)Initialization

    In MA,the population is composed ofPnindividuals.Each individualxrepresents a scheme of removing edges in the network and was generated randomly.

    (2)Crossover

    The crossover operator works on two parent individuals and can search in a large area.Suppose thatxp1andxp2are two parent individuals,andxc1andxc2are two child individu-als.First,we assignxp1toxc1andxp2toxc2and obtain the following sets of edges:

    Table 1 Pseudocode of MA.

    In summary,only the set of non-common edges that we want to remove will be swapped betweenxc1andxc2(Fig.3).

    (3)Local search

    The local search operator is an important part in MA that can accelerate the convergence speed.Based on previous edge importance evaluations,27,28we adopted a local search in the direction of removing more important edges.We first select an individual from parent and child population using the roulette wheel selection based on their fitness.Then for each edge of the selected individual,we conduct a local search with probabilityPl.For example,edgeeij(eijrepresents the edge between nodeiand nodej)will mutate into a randomly selected existing edgeelmbut not in its individual.This mutation will be accepted when the following formula is satisfied:

    Table 2 Parameters of memetic algorithm.

    where μ is a formula parameter in the range[0,1]andkl,km,kiandkjrepresent the degree of nodesl,m,iandjrespectively.

    (4)Tournament selection

    In this part,two individuals are respectively chosen from parent population and child population to run a tournament for the tournament selection.The population with the best fitness is selected for the next generation,and the total number of generation isPm.

    To conclude,the pseudocode of MA proposed is presented in Table 1.

    4.Results

    In some previous papers,under different kinds of malicious attacks on edges,the strategy based on edge-betweenness is actually a commonly adopted attack strategy.20,29,30Here,in order to identify the vital edges,we examine the robustness of the new networks after removing edges via MA,and compare it with the highest edge-betweenness adaptive strategy(BEAS).As many real-world networks are of scale-free properties,such as air transportation network,9World-Wide Web,31Internet32and social network,33the experiments are carried out not only on Chinese ARN but also on Baraba′si-Albert(BA)scale-free network to demonstrate the universality of our method.The Chinese ARN has 1499 nodes and 2242 edges.The BA scale-free network is generated withm0nodes and a new node is added withmedges at each time step,which connect the new node withmdifferent existing nodes.Here,it is set thatm0=2 andm=2 and the BA network is of 1000 nodes and 2000 edges.The cost denoting total number of removing edges is set from 0 to 300 and the network without edges removed represents the initial network.Table 2 shows the configurations of the memetic algorithm parameters used on the optimization model.

    Fig.4 shows the simulation results of the MA and BEAS of identifying the vital edges for the BA network and Chinese ARN.Looking at both networks via MA,we can see that the robustnessRdecreases and the costCincreases when removing edges(Fig.4(a)and(b)).It can also be noticed that the MA is significantly better than the BEAS.The difference between the two methods is especially high.It is obvious that the critical edges in the network are not extremely related with the edge-betweenness.Moreover,the memetic algorithm works better in the Chinese ARN and decreases its robustnessRfaster when a few edges are removed.The reason is that the Chinese ARN is not a small-world network.And the Chinese ARN is vulnerable because of its small clustering coefficient and large average shortest path length.In detail,Fig.4(c)and Fig.4(d)separately show 10 edges identified by the MA and BEAS methods in the Chinese ARN.It is found that all the top 10 highest edge-betweenness edges in the network are located at the middle China,which are almost completely different from the 10 edges identified by MA.

    We have seen that the MA works on both Chinese ARN and BA network.In order to reveal the underlying mechanism clearly,we examine a toy model with a network containing 10 nodes and 17 edges(Fig.5).In Fig.5,the blue lines are the existing edges and gray dotted lines are edges removed in that step.In the same way,the gray nodes mean that these nodes are removed and yellow nodes still exist.When no edge is removed from the network(Fig.5(a)),the robustness of the initial network is 0.35.Here,three edges are removed to measure the criticality of these edges using the MA and BEAS methods.Since edgese1,4,e2,7ande4,9have the highest edgebetweenness in the network,they are removed in the BEAS(Fig.5(b))and now we have a new network namedA.Similarly,edgese1,3,e4,9ande7,9are identified as the most important ones in the MA,and we now have a new network B(Fig.5(f)).For estimating which group of edges is critical,we compare the robustness of networkAandB(Fig.5(c)–(e))and Fig.5(g)–(i)).Table 3 illustrates the corresponding solutions and the robustness of the solution for both methods.

    In networkA(BEAS),we first remove node 2 with the highest-degree together with all edges connected with it:e2,3,e2,4,e2,5,e2,8ande2,10.s(Q=1)of the new network is 0.9,which is the fraction of nodes in the giant component after removing 1 node(Fig.5(c)).However in networkB(MA),the value ofs(Q=1)quickly decreases to 0.7(Fig.5(g)).Then,after the second nodes are removed,s(Q=2)in both networksAandBare 0.4 and 0.3 respectively(Fig.5(d)and Fig.5(h)),which is reduced to 0.4 and 0.2 after the third node is removed(Fig.5(e)and Fig.5(i)).At the end,all nodes are removed from the network and the robustness of networkAandBis 0.37 and 0.29 respectively.Thus,as previous results revealed,the critical edges in the network are not extremely related to the edge-betweenness,which apparently contradicts common intuitions.

    Table 3 illustrates the corresponding solutions and the robustness of the solution to both methods of the toy model in Fig.5.In the network after removingQnodes and the edges connected with them,Nis the number of nodes,Mis the number of edges,ands(Q)is the fraction of nodes in the largest component.The results demonstrate that the most vital edges are not necessarily the edge with the highest topological importance.Thus these edges identified by MA are important for the network robustness and should be protected to ensure the survivability of the network.

    Table 3 Illustration of vital edges identified by MA and BEAS.

    5.Conclusions

    It is of great importance to improve the robustness of real networks.In this paper,we identified the vital edges in Chinese air route network,which lead to fast breakdown after targeted attacks.Our results reveal that the edge-betweenness,an index to measure the importance of edges in short paths,is of little relevance to this problem.Furthermore,we demonstrate that the memetic algorithm is able to pinpoint the edges that have been proven more important than edges of high edgebetweenness.We also confirm these findings in scale-free model networks,hence offering novel insights of edge essentiality in various real networks.Thus,we think the vital edges identified by memetic algorithm should be especially protected to ensure a good performance of a network.In Chinese ARN,this means that air traffic managers should foresee complex solutions when considering the closure of one vital air route segment.

    Acknowledgements

    This paper is supported by the National Natural Science Foundation of China(Nos.91538204,61425014,61521091),National Key Research and Development Program of China(No.2016YFB1200100),andNationalKeyTechnology R&D Program of China(No.2015BAG15B01).

    1.Amaral LAN,Scala A,Barthe′le′my M,Stanley HE.Classes of small-world networks.ProcNatlAcadSciUSA2000;97(21):11149–52.

    2.Guimera`R,Mossa S,Turtschi A,Amaral LAN.The worldwide air transportation network:Anomalous centrality,community structure,and cities’global roles.Proc Natl Acad Sci USA2005;102(22):7794–9.

    3.Barrat A,Barthe′lemy M,Pastor-Satorras R,Vespignani A.The architecture of complex weighted networks.Proc Natl Acad Sci USA2004;101(11):3747–52.

    4.Gautreau A,Barrat A,Barthe′lemy M.Microdynamics in stationary complex networks.Proc Natl Acad Sci USA2009;106(22):8847–52.

    5.Verma T,Arau′jo NAM,Herrmann HJ.Revealing the structure of the world airline network.Sci Reports2014;4:5638.

    6.Verma T,Russmann F,Arau′jo NAM,Nagler J,Herrmann HJ.Emergence of core–peripheries in networks.Nat Commun2016;7(1):10441.

    7.da Rocha LEC.Structural evolution of the Brazilian airport network.J Stat Mech2009;2009(4):125–36.

    8.Wang JE,Mo HH,Wang FH.Evolution of air transport network of China 1930–2012.J Trans Geo2014;40:145–58.

    9.Zhang J,Cao XB,Du WB,Cai KQ.Evolution of Chinese airport network.Physica A2010;389:3922–31.

    10.Li W,Cai X.Statistical analysis of airport network of China.Phys Rev E2004;69(4 Pt 2):046106.

    11.Fleurquin P,Ramasco JJ,Eguiluz VM,Victor M.Characterization of delay propagation in the US air-transportation network.Transp J2014;53(3):330–44.

    12.Bagler G.Analysis of the airport network of India as a complex weighted network.Physica A2008;387(12):2972–80.

    13.Du WB,Zhou XL,Lordan O,Wang Z,Zhao C,Zhu YB.Analysis of the Chinese Airline Network as multi-layer networks.Trans Res Part E2016;89:108–16.

    14.Liu HK,Zhou T.Empirical study of Chinese city airline network.Acta Phys Sin2007;56(1):106[Chinese].

    15.Cai KQ,Zhang J,Du WB,Cao XB.Analysis of the Chinese air route network as a complex network.Chin Phys B2012;21(2):028903.

    16.Vitali S,Cipolla M,Gurtner G,Lillo F,Beato V,Pozzi S.Statistical regularities in ATM:Network properties,trajectory deviations and delays.Second SESAR innovation days2012.

    17.Schneider CM,Moreira AA,Andrade Jr JS,Havlin S,Herrmann HJ.Mitigation of malicious attacks on networks.Proc Natl Acad Sci2011;108(10):3838–41.

    18.Lordan O,Sallan JM,Simo P,Gonzalez-Prieto D.Robustness of airline alliance route networks.Commun Nonlin Sci Numer Simu2015;22(1–3):587–95.

    19.Lordan O,Sallan JM,Escorihuela N,Gonzalez-Prieto D.Robustness of airline route networks.Physica A2016;445:18–26.

    20.Zeng A,Liu WP.Enhancing network robustness against malicious attacks.Phys Rev E2012;85(6):066130.

    21.Motter AE,Lai YC.Cascade-based attacks on complex networks.Phys Rev E2002;66(6):065102.

    22.Freeman L.Set of measures of centrality based on betweenness.Sociometry1977;40(1):35–41.

    23.Menc?′a R,Sierra MR,Menc?′a C,Varela R.Memetic algorithms for the job shop scheduling problem with operators.Appl Soft Comput2015;34:94–105.

    24.Gong MG,Cai Q,Li YY,Ma JJ.An improved memetic algorithm for community detection in complex networks.2012 IEEE world congress on evolutionary computation;2012 June 10–15,Brisbane,Australia.Piscataway(NJ):IEEE Press;2012.

    25.Neri F,Cotta C.Memetic algorithms and memetic computing optimization: A literature review.SwarmEvolComput2012;2:1–14.

    26.Bhuvana J,Aravindan C.Memetic algorithm with preferential local search using adaptive weights for multi-objective optimization problems.Soft Comput2016;20(4):1365–88.

    27.Tan F,Xia YX,Zhang WP,Jin XY.Cascading failures of loads in interconnected networks under intentional attack.EPL2013;102(2):28009.

    28.Peng XZ,Yao H,Du J.Load-induced cascading failures in interconnected networks.Nonlin Dyna2015;82(1–2):97–105.

    29.Wang JW,Rong LL.Robustness of the western United States power grid under edge attack strategies due to cascading failures.Safety Sci2011;49(6):807–12.

    30.Mirzasoleiman B,Babaei M,Jalili M,Safari M.Cascaded failures in weighted networks.Phys Rev E2011;84(4):046114.

    31.Albert R,Jeong H,Baraba′si AL.Diameter of the World-Wide Web.Nature1999;401(6):130–1.

    32.Faloutsos M,Faloutsos P,Faloutsos C.On power-law relationships of the internet topology.Comput Commun Rev1999;29(4):251–63.

    33.Newman MEJ.The structure and function of complex networks.SIAM Rev2003;45(2):167–256.

    28 July 2016;revised 12 September 2016;accepted 8 October 2016

    Available online 21 December 2016

    ?2016 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is anopenaccessarticleundertheCCBY-NC-NDlicense(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    *Corresponding author.

    E-mail address:xbcao@buaa.edu.cn(X.Cao).

    Peer review under responsibility of Editorial Committee of CJA.

    猜你喜歡
    醋酸鈉緩沖溶液蒸餾水
    無水醋酸鈉結(jié)構(gòu)及熱穩(wěn)定性
    幾種緩沖溶液簡介及應(yīng)用*
    Analysis of Wastewater Membrane Pollutants in Joint Station and Research on Biological Control Technology
    基礎(chǔ)化學(xué)緩沖溶液教學(xué)難點(diǎn)總結(jié)
    科技視界(2017年25期)2017-12-11 20:30:32
    二氯醋酸鈉提高膠質(zhì)母細(xì)胞瘤U251細(xì)胞的放療敏感性
    剖宮產(chǎn)術(shù)前輸注醋酸鈉鉀鎂鈣葡萄糖注射液對新生兒血?dú)狻⒀?、乳酸及電解質(zhì)的影響
    用于蒸餾水機(jī)高溫測量的DPI系列智能測量儀表
    多效蒸餾水機(jī)冷凝水的熱能回收利用
    緩沖溶液法回收置換崗位中二氧化硫尾氣
    河南科技(2014年15期)2014-02-27 14:12:29
    電位滴定法測定NH3-NH4Cl緩沖溶液的緩沖容量
    国产欧美日韩精品一区二区| 精品午夜福利在线看| 免费大片18禁| 日韩 亚洲 欧美在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲无线观看免费| 国内少妇人妻偷人精品xxx网站| 美女国产视频在线观看| 亚洲五月天丁香| 欧美成人精品欧美一级黄| 免费av观看视频| av.在线天堂| www.色视频.com| 日日摸夜夜添夜夜爱| 成年免费大片在线观看| 国产亚洲精品av在线| 精品熟女少妇av免费看| 亚洲av.av天堂| 国产成人a∨麻豆精品| 欧美3d第一页| 日韩欧美在线乱码| 男女啪啪激烈高潮av片| 亚洲欧洲日产国产| 国产精品一区www在线观看| 男女那种视频在线观看| 在线观看免费视频日本深夜| 婷婷六月久久综合丁香| 久久6这里有精品| 久久国内精品自在自线图片| 国产在线精品亚洲第一网站| 99热6这里只有精品| 一级黄片播放器| 青春草国产在线视频 | 国产av一区在线观看免费| 观看免费一级毛片| 中文字幕人妻熟人妻熟丝袜美| 蜜桃久久精品国产亚洲av| 一级av片app| 国产亚洲欧美98| 久久精品国产亚洲av天美| 精品久久久久久久久久久久久| 嘟嘟电影网在线观看| 国产日本99.免费观看| 中文字幕制服av| 亚洲,欧美,日韩| 国产亚洲5aaaaa淫片| 亚洲精品国产成人久久av| 乱系列少妇在线播放| 最近最新中文字幕大全电影3| 在线观看免费视频日本深夜| 免费不卡的大黄色大毛片视频在线观看 | 在线观看午夜福利视频| 国产91av在线免费观看| 亚洲av二区三区四区| 少妇丰满av| 岛国在线免费视频观看| 蜜桃久久精品国产亚洲av| 久久久久久久久久久免费av| 国产精品av视频在线免费观看| 夜夜爽天天搞| 亚洲人成网站在线播放欧美日韩| 国内揄拍国产精品人妻在线| 嫩草影院入口| 午夜久久久久精精品| 波野结衣二区三区在线| 久久99热这里只有精品18| 性插视频无遮挡在线免费观看| 日本三级黄在线观看| 人人妻人人澡人人爽人人夜夜 | 成人亚洲欧美一区二区av| 九九久久精品国产亚洲av麻豆| 91久久精品国产一区二区成人| 精品99又大又爽又粗少妇毛片| 欧美不卡视频在线免费观看| 久久精品久久久久久噜噜老黄 | 在线播放国产精品三级| .国产精品久久| 久久99热这里只有精品18| 午夜老司机福利剧场| 级片在线观看| 日本-黄色视频高清免费观看| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 黄色欧美视频在线观看| 日本黄大片高清| 在线观看一区二区三区| 搡女人真爽免费视频火全软件| 国产精品不卡视频一区二区| 伊人久久精品亚洲午夜| 久久久久国产网址| 婷婷亚洲欧美| 夜夜看夜夜爽夜夜摸| 成人三级黄色视频| 长腿黑丝高跟| 免费看美女性在线毛片视频| 国模一区二区三区四区视频| 欧美色欧美亚洲另类二区| 日韩欧美国产在线观看| 国产亚洲av片在线观看秒播厂 | 欧美不卡视频在线免费观看| 亚洲精华国产精华液的使用体验 | 国内揄拍国产精品人妻在线| .国产精品久久| 在线观看66精品国产| 熟妇人妻久久中文字幕3abv| 麻豆精品久久久久久蜜桃| 国内精品美女久久久久久| 国产在线男女| 一级av片app| 性插视频无遮挡在线免费观看| 国产精品爽爽va在线观看网站| 自拍偷自拍亚洲精品老妇| 在线播放无遮挡| 亚洲不卡免费看| 天天躁日日操中文字幕| 日韩欧美一区二区三区在线观看| 老熟妇乱子伦视频在线观看| 亚洲三级黄色毛片| 午夜激情福利司机影院| videossex国产| 国产淫片久久久久久久久| 久久久久网色| 成人二区视频| 亚洲在线观看片| 欧美日本视频| 美女脱内裤让男人舔精品视频 | 日韩中字成人| 婷婷六月久久综合丁香| 日韩强制内射视频| 女的被弄到高潮叫床怎么办| 成人亚洲欧美一区二区av| 国产一区亚洲一区在线观看| 三级毛片av免费| 中文字幕av在线有码专区| 国语自产精品视频在线第100页| 两个人的视频大全免费| 自拍偷自拍亚洲精品老妇| 中文字幕av成人在线电影| 国产精品久久久久久精品电影| 亚洲一级一片aⅴ在线观看| 日韩成人av中文字幕在线观看| 一个人观看的视频www高清免费观看| 免费观看的影片在线观看| 99热6这里只有精品| 乱码一卡2卡4卡精品| 中国美白少妇内射xxxbb| 九九在线视频观看精品| 日日啪夜夜撸| 国产精品av视频在线免费观看| 麻豆国产97在线/欧美| 久久综合国产亚洲精品| 男人的好看免费观看在线视频| 久久精品国产99精品国产亚洲性色| 国内精品久久久久精免费| 夫妻性生交免费视频一级片| 国产精品人妻久久久影院| 国产美女午夜福利| 日韩,欧美,国产一区二区三区 | 欧美精品一区二区大全| 成人三级黄色视频| 久久人人精品亚洲av| 麻豆久久精品国产亚洲av| 大香蕉久久网| 内射极品少妇av片p| 成人欧美大片| 国产日韩欧美在线精品| 中文字幕免费在线视频6| 久久久久性生活片| 最近手机中文字幕大全| 美女内射精品一级片tv| 色播亚洲综合网| 一区二区三区免费毛片| 99久国产av精品| 国产激情偷乱视频一区二区| 日本-黄色视频高清免费观看| avwww免费| 欧美性感艳星| 在线免费观看不下载黄p国产| 美女高潮的动态| 中文字幕制服av| 欧美一级a爱片免费观看看| 观看免费一级毛片| 亚洲av成人精品一区久久| 国产伦理片在线播放av一区 | 精品免费久久久久久久清纯| 久久久久久久亚洲中文字幕| 日本与韩国留学比较| 天天躁日日操中文字幕| 午夜免费男女啪啪视频观看| 亚洲美女视频黄频| 国产美女午夜福利| 国产高清不卡午夜福利| 综合色av麻豆| 精品人妻一区二区三区麻豆| 日日啪夜夜撸| 美女cb高潮喷水在线观看| 国产精品人妻久久久久久| 日韩中字成人| 乱系列少妇在线播放| 国产精品.久久久| 丰满人妻一区二区三区视频av| 三级毛片av免费| 婷婷精品国产亚洲av| 久久久精品94久久精品| 99久久人妻综合| 色哟哟哟哟哟哟| 69人妻影院| 国产成人精品一,二区 | 日本一本二区三区精品| 亚洲国产精品sss在线观看| 亚洲最大成人手机在线| 成人高潮视频无遮挡免费网站| 欧美精品国产亚洲| 午夜精品国产一区二区电影 | 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 色吧在线观看| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产| 91久久精品国产一区二区成人| 可以在线观看毛片的网站| 最新中文字幕久久久久| 最近的中文字幕免费完整| 免费无遮挡裸体视频| 国产成人福利小说| 久久久久久久久久久丰满| 亚洲精品亚洲一区二区| 日韩高清综合在线| 免费看光身美女| 麻豆国产97在线/欧美| 亚洲精品乱码久久久v下载方式| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 麻豆成人av视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 白带黄色成豆腐渣| 变态另类丝袜制服| 内地一区二区视频在线| 国产亚洲5aaaaa淫片| 一本久久精品| 搞女人的毛片| 九色成人免费人妻av| 蜜桃久久精品国产亚洲av| ponron亚洲| 日韩成人av中文字幕在线观看| 欧美一级a爱片免费观看看| 国产不卡一卡二| 精品欧美国产一区二区三| 久久久久久久久中文| 美女大奶头视频| 婷婷精品国产亚洲av| 神马国产精品三级电影在线观看| av又黄又爽大尺度在线免费看 | 成人漫画全彩无遮挡| 国产色爽女视频免费观看| 国产一区二区三区av在线 | 最好的美女福利视频网| 亚洲最大成人手机在线| 亚洲av二区三区四区| av在线天堂中文字幕| 国产精品麻豆人妻色哟哟久久 | 2021天堂中文幕一二区在线观| 99久国产av精品国产电影| 国产午夜精品一二区理论片| 亚洲欧美精品专区久久| 欧美在线一区亚洲| 一级av片app| 99久久成人亚洲精品观看| 精品不卡国产一区二区三区| 熟女电影av网| 丰满的人妻完整版| 美女内射精品一级片tv| 国产乱人偷精品视频| 女同久久另类99精品国产91| 色综合亚洲欧美另类图片| 久久午夜福利片| 别揉我奶头 嗯啊视频| 久久久久久久久久久免费av| 99热网站在线观看| 亚洲高清免费不卡视频| 亚洲在线自拍视频| 神马国产精品三级电影在线观看| 深爱激情五月婷婷| 人妻久久中文字幕网| 小说图片视频综合网站| 国产精品人妻久久久影院| 欧美另类亚洲清纯唯美| 亚洲精品久久久久久婷婷小说 | 国产中年淑女户外野战色| 天堂中文最新版在线下载 | 熟妇人妻久久中文字幕3abv| 国产精品爽爽va在线观看网站| 久久久精品大字幕| 免费电影在线观看免费观看| 亚洲久久久久久中文字幕| 国产午夜精品久久久久久一区二区三区| 国产人妻一区二区三区在| 三级经典国产精品| 国产黄色小视频在线观看| 成人欧美大片| 国产精品人妻久久久影院| 国产精品永久免费网站| 久久精品国产自在天天线| 久久精品久久久久久久性| 日韩欧美 国产精品| 亚洲精品国产av成人精品| 日韩国内少妇激情av| 久久精品综合一区二区三区| 男女那种视频在线观看| 我的女老师完整版在线观看| 在线播放国产精品三级| 久久婷婷人人爽人人干人人爱| 91精品国产九色| 最好的美女福利视频网| a级毛色黄片| 精品人妻偷拍中文字幕| 精品人妻熟女av久视频| 三级国产精品欧美在线观看| 亚洲四区av| .国产精品久久| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩国产亚洲二区| 欧美性猛交╳xxx乱大交人| 国产在线男女| 一级毛片电影观看 | 国产精品99久久久久久久久| 免费大片18禁| 亚洲无线观看免费| 深爱激情五月婷婷| 成人永久免费在线观看视频| 黄色视频,在线免费观看| 看片在线看免费视频| 国产亚洲91精品色在线| 青春草视频在线免费观看| 高清日韩中文字幕在线| 国产 一区 欧美 日韩| 日韩欧美一区二区三区在线观看| 99久久成人亚洲精品观看| 亚洲欧美成人综合另类久久久 | 老熟妇乱子伦视频在线观看| 卡戴珊不雅视频在线播放| 日韩国内少妇激情av| 国产精品女同一区二区软件| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看| 麻豆成人av视频| 乱人视频在线观看| av在线亚洲专区| 91久久精品国产一区二区三区| 亚洲丝袜综合中文字幕| 亚洲内射少妇av| 日韩视频在线欧美| 天堂√8在线中文| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看| 天美传媒精品一区二区| 黄色欧美视频在线观看| 18禁黄网站禁片免费观看直播| 亚洲av免费高清在线观看| 国产成人一区二区在线| 色吧在线观看| 麻豆乱淫一区二区| 熟女电影av网| 18禁在线播放成人免费| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 给我免费播放毛片高清在线观看| 不卡视频在线观看欧美| 国产亚洲91精品色在线| 精品日产1卡2卡| 中文字幕av成人在线电影| 97热精品久久久久久| 日日摸夜夜添夜夜添av毛片| 波多野结衣高清作品| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 日韩av在线大香蕉| 亚洲丝袜综合中文字幕| 欧美日韩精品成人综合77777| 99久国产av精品国产电影| 大香蕉久久网| 亚洲av成人av| 午夜福利成人在线免费观看| 精品一区二区三区人妻视频| 丝袜美腿在线中文| 国产av麻豆久久久久久久| 国产不卡一卡二| 99热全是精品| 三级男女做爰猛烈吃奶摸视频| 激情 狠狠 欧美| 最近最新中文字幕大全电影3| 免费看a级黄色片| 久久午夜亚洲精品久久| 国产成人freesex在线| 变态另类丝袜制服| 国产在视频线在精品| 亚洲欧美精品专区久久| 亚洲国产精品成人久久小说 | 日本爱情动作片www.在线观看| 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 成人亚洲欧美一区二区av| 神马国产精品三级电影在线观看| 国产日本99.免费观看| 最好的美女福利视频网| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女视频黄频| 亚洲在线观看片| 国产一区二区激情短视频| a级一级毛片免费在线观看| or卡值多少钱| 国内精品一区二区在线观看| 日韩高清综合在线| 91久久精品电影网| 国产亚洲欧美98| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片| 热99在线观看视频| 99热这里只有是精品50| 变态另类丝袜制服| 尾随美女入室| 国产高清有码在线观看视频| 伦理电影大哥的女人| 18禁在线无遮挡免费观看视频| 你懂的网址亚洲精品在线观看 | 国产伦精品一区二区三区四那| 变态另类丝袜制服| 国产午夜精品论理片| av视频在线观看入口| 久久久久国产网址| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站 | 我的老师免费观看完整版| 国产成人影院久久av| 国产午夜精品论理片| 99在线人妻在线中文字幕| 成人二区视频| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 日韩视频在线欧美| 搡女人真爽免费视频火全软件| 美女 人体艺术 gogo| 亚洲第一区二区三区不卡| 校园人妻丝袜中文字幕| 成人欧美大片| 爱豆传媒免费全集在线观看| 国产久久久一区二区三区| 村上凉子中文字幕在线| 51国产日韩欧美| 久久久久久九九精品二区国产| 国产三级中文精品| 久久久色成人| 国产老妇女一区| 性插视频无遮挡在线免费观看| 久久精品国产亚洲av涩爱 | 成人美女网站在线观看视频| av免费在线看不卡| 国内久久婷婷六月综合欲色啪| 26uuu在线亚洲综合色| 亚洲一区高清亚洲精品| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 国产成人一区二区在线| 国产一区二区在线av高清观看| 久久人妻av系列| 国产不卡一卡二| 男人狂女人下面高潮的视频| 国产女主播在线喷水免费视频网站 | 精品99又大又爽又粗少妇毛片| 男女啪啪激烈高潮av片| 男人和女人高潮做爰伦理| 岛国毛片在线播放| 日本免费a在线| kizo精华| 国产一区二区激情短视频| 99久久成人亚洲精品观看| 亚洲真实伦在线观看| 99热6这里只有精品| 久久热精品热| 麻豆国产av国片精品| 美女脱内裤让男人舔精品视频 | 久久久精品大字幕| 在线免费观看的www视频| 精品人妻一区二区三区麻豆| 国内揄拍国产精品人妻在线| 成年av动漫网址| 亚洲美女视频黄频| 亚洲不卡免费看| 久久精品国产亚洲av香蕉五月| 亚洲不卡免费看| 国产蜜桃级精品一区二区三区| 中文字幕免费在线视频6| 哪个播放器可以免费观看大片| 99在线视频只有这里精品首页| 变态另类成人亚洲欧美熟女| 一个人免费在线观看电影| 在线播放无遮挡| 噜噜噜噜噜久久久久久91| 免费人成在线观看视频色| 深夜精品福利| 高清毛片免费看| 一区福利在线观看| 中文字幕免费在线视频6| 男人舔女人下体高潮全视频| 欧美另类亚洲清纯唯美| 午夜a级毛片| 欧洲精品卡2卡3卡4卡5卡区| www.色视频.com| 一本精品99久久精品77| 国产精品一及| 在线播放国产精品三级| 久久99热这里只有精品18| 国产精品久久电影中文字幕| 亚洲成人中文字幕在线播放| 成年免费大片在线观看| 午夜福利在线在线| 色视频www国产| 国产亚洲5aaaaa淫片| 午夜爱爱视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲高清免费不卡视频| 在线a可以看的网站| 亚洲精品国产av成人精品| 日本黄大片高清| 成人性生交大片免费视频hd| 精品人妻偷拍中文字幕| 久久久久性生活片| 亚洲国产精品sss在线观看| 99久国产av精品| 久久亚洲精品不卡| 国产av一区在线观看免费| av免费在线看不卡| 国产日韩欧美在线精品| 有码 亚洲区| 亚洲精华国产精华液的使用体验 | 国产av麻豆久久久久久久| 久久久久久大精品| 人人妻人人看人人澡| 乱人视频在线观看| 久久久久久久久久成人| 女人十人毛片免费观看3o分钟| 99久久九九国产精品国产免费| 国产精品久久视频播放| 精品不卡国产一区二区三区| 一级毛片我不卡| 国产一级毛片七仙女欲春2| 国产精品嫩草影院av在线观看| 日韩欧美 国产精品| 国产中年淑女户外野战色| 村上凉子中文字幕在线| 亚洲久久久久久中文字幕| 久久久久久久午夜电影| 日本成人三级电影网站| 日本撒尿小便嘘嘘汇集6| 99久久中文字幕三级久久日本| 成人永久免费在线观看视频| 老司机影院成人| 五月伊人婷婷丁香| 亚洲丝袜综合中文字幕| 18禁在线无遮挡免费观看视频| 国内少妇人妻偷人精品xxx网站| 99热这里只有精品一区| 狂野欧美激情性xxxx在线观看| 亚洲人成网站在线播放欧美日韩| 日本黄色片子视频| 国产高清视频在线观看网站| 伊人久久精品亚洲午夜| 1024手机看黄色片| 一区福利在线观看| 亚洲成av人片在线播放无| 日韩欧美精品免费久久| 中文欧美无线码| 免费观看精品视频网站| 亚洲高清免费不卡视频| 欧美日本亚洲视频在线播放| 真实男女啪啪啪动态图| 国产黄片视频在线免费观看| 99久久久亚洲精品蜜臀av| 亚洲欧美精品综合久久99| 精品久久久久久久久av| 亚洲成人av在线免费| 国产单亲对白刺激| 国产精品三级大全| 成人亚洲欧美一区二区av| 日本黄色片子视频| 国产av不卡久久| 亚洲国产精品合色在线| 成人性生交大片免费视频hd| 丰满的人妻完整版| 亚洲在久久综合| 成人av在线播放网站| 久久这里只有精品中国| 欧美日韩国产亚洲二区| 日本免费一区二区三区高清不卡| 99久久人妻综合| 国产精品女同一区二区软件| 悠悠久久av| 高清在线视频一区二区三区 | 最后的刺客免费高清国语| 亚洲精品日韩av片在线观看| 亚洲av中文字字幕乱码综合| 一本久久中文字幕| 久久精品久久久久久噜噜老黄 | 91aial.com中文字幕在线观看| 国产探花在线观看一区二区| 精品久久久久久久末码| 亚洲成人精品中文字幕电影| av专区在线播放| 午夜免费激情av| 男女下面进入的视频免费午夜| 欧美日本视频| 日韩亚洲欧美综合|