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

    Percolation transitions in edge-coupled interdependent networks with directed dependency links

    2023-10-11 07:56:42YanLiGao高彥麗HaiBoYu于海波JieZhou周杰YinZuoZhou周銀座andShiMingChen陳世明
    Chinese Physics B 2023年9期
    關(guān)鍵詞:周杰銀座

    Yan-Li Gao(高彥麗), Hai-Bo Yu(于海波), Jie Zhou(周杰), Yin-Zuo Zhou(周銀座), and Shi-Ming Chen(陳世明),?

    1School of Electrical and Automation Engineering,East China JiaoTong University,Nanchang 330013,China

    2School of Physics and Electronic Science,East China Normal University,Shanghai 200241,China

    3Alibaba Research Center for Complexity Sciences,Hangzhou Normal University,Hangzhou 311121,China

    Keywords: edge-coupled interdependent networks with directed dependency links, percolation transitions,cascading failures,robustness analysis

    1.Introduction

    Many facilities and systems in real life can be abstracted into networks, such as transportation networks, biological networks, financial networks, electric networks, and social networks.[1–9]These networks have brought convenience to human production and life,while the propagation phenomena occurring on the networks have also brought negative effects,and complex networks science provides theories and analytical methods for studying complex systems in the real world.

    Early research on complex networks was limited to the topological properties of the single-layer network,[10–13]until it was discovered that the interdependence within networks in the real world posed extremely serious systemic risks,such as large-scale power outages and traffic congestion.[14,15]This is due to the fact that the failures in a single-layer network could spread not only in itself but also in the other networks.Therefore,the study of interdependent networks has attracted much attention of many scholars.

    Buldyrevet al.[16]first developed the analytical framework of interdependent networks,where there are one-to-one dependent nodes between two networks, and found that interdependence greatly reduces the robustness of the network.Following this research, many results have emerged for the study of interdependent networks.According to the difference of interdependent patterns, Parshaniet al.[17]analyzed the node-coupled interdependent networks in which a fraction of the nodes are interdependent in both networks, and found that reducing the coupling between the networks leads to a change from a first-order percolation phase transition to a second-order percolation transition at a critical point, i.e., its robustness is better than the fully node-coupled interdependent networks.[16]Shaoet al.[18]proposed the coupled network systems with multiple support-dependence relations and investigated them analytically and numerically.Buldyrevet al.[19]and Parshaniet al.[20]studied the coupled networks with positively correlated degrees of mutually dependent nodes and found that as the networks become more inter-similar the system would become significantly more robust to random failure.Huet al.[21]found that the phase transition includes the hybrid transition of mixed first and second orders in study of a coupled networks system where both interdependent and interconnected links exist.This hybrid transition phenomenon was also found in the studies of Zhouet al.[22]and Zhanget al.[23]afterward.Donget al.[24]developed a interacting networks model with feedback-dependency links and found that the feedback-dependency links have no effect on robustness of the system for a weak coupling strength corresponding to the second-order transitions, whereas the system is more vulnerable for strong coupling strength corresponding to first-order or hybrid phase transition.Liuet al.[25]found that reducing the coupling strength between interdependent nodes can improve the system’s robustness by studying the interdependent networks with weak node-coupling strength.

    In addition to studying interdependent patterns, scholars have proposed numerous network models.Gaoet al.[26–28]extended interdependent networks to systems with more than two independent networks and analyzed the effect of topology and the number of networks on network robustness.Cellaiet al.[29]showed that the overlap among layers can improve the robustness of interdependent multiplex systems and change the critical behavior of the percolation phase transition in a complex way.Huanget al.[30]studied the effect of clustering on the system robustness of interdependent networks and found that clustering makes the system more vulnerable.Valdet al.[31]set the nodes with high degrees as autonomous nodes and found them to be effective in improving the robustness of the network.Liuet al.[32]developed a theoretical framework of interdependent directed networks and found that the directionality within the network increases its vulnerability and exhibits hybrid phase transitions.

    The above works address the case of node-coupled interdependent networks (NINs), i.e., the failures of nodes in one network lead to the failures of nodes in the other network.Considering that there are also many cases of edge-coupled interdependent networks in the real world, Gaoet al.[33]developed a new model of edge-coupled interdependent networks(EINs) and analyzed it using the self-consistent probabilities method[34]and found that it is stronger than the NIN.

    The existing studies of EINs have only considered the case where both connectivity links and dependency links are undirected.However, there are also some EINs in which the dependency links are directed in the real world.One example is power-transportation interdependent networks.Generally speaking, power system and rail transportation system interdepend on each other.However,by analyzing the interactions between these two networks we found that there are one-way dependencies instead of mutual dependencies for the specific lines.For example,the failures of part power lines can directly affect the functions of some part of rail transportation lines,whereas these transportation lines may not just right affect these power lines.In other words, there are directed dependency relations between transportation lines and power lines.

    Therefore, a model of edge-coupled interdependent networks with directed dependency links (EINDDLs) is proposed,and the theoretical framework for the analysis of EINDDLs under random attacks is developed by settingβto change the dependence relationship between two networks based on the self-consistent probabilities method.Both theoretical analysis and computer simulations are conducted on RR-RR and ER-ER EINDDLs(RR:random regular,ER:Erd¨os–Rényi)to explore the relationship betweenβand the phase transition behaviors of networks A and B under random attacks.Some characteristics in the phase transition behaviors of networks A and B are found and the corresponding prevention strategy is proposed.

    2.Model description

    The model of edge-coupled interdependent networks with directed dependency links (EINDDLs) in this paper refers to two-layer edge-coupled interdependent networks(EINs)consisting of networks A and B,where the dependency links are one direction instead of mutual direction.We assume that the number of nodes and the average degree of both networks are the same, so the total numbers of edges in both networks are equal.A connectivity link in network A is randomly selected to establish a directed dependency link with a connectivity link in network B, as shown by the red and blue dashed lines in Fig.1(a).Figure 1 gives the process of cascading failures for the EINDDL.Dependency links make the system follow such cascading rules,e.g.,failure ofa1?a8causes failure ofb2?b3while the failure ofb2?b3cannot lead to failure ofa1?a8.Here,βis defined to control the ratio of two directed types of dependency links, i.e., network A has a 1-βratio of connectivity links depending on network B,and network B has aβratio of connectivity links depending on network A.The dependencies of the two networks will also be changed by changing the value ofβ.

    Fig.1.Cascading failures of the EINDDL.

    According to the bond percolation theory, when a network is under attacks, it may be broken into many connected clusters.The biggest one referred as giant component(GC)is usually considered as functional part of a network.We randomly remove the 1-pfraction of connectivity links in network A in the initial state,this will result in a portion of nodes failing due to being out of GC, while the edges connected to these failed nodes will also fail.Since there are connectivity links in network B that depends on network A,it causes a part of the connectivity links in network B to fail and makes a part of the nodes fail by leaving the GC,which subsequently leads to the failure of connectivity links in network A that depends on network B and causes further damage to network A.Cascading failures repeatedly propagate between two networks as the above process until no more new connectivity links or nodes fail in both networks,which is the steady state.In terms of interdependent networks,the remainder of the network after the cascading failure composes an extensive component, referred as mutually connected giant component(MCGC).The size of the MCGC at a steady state reflects the robustness of the network.

    The solid black lines in Fig.1 represent the connectivity links within the network,the dashed black lines represent the failed connectivity links,and the dashed red and blue lines are the dependency links between two networks, whereβis 0.6.Due to the initial statea3?a4anda4?a5in network A failed,as shown in Fig.1(a).According to the cascading failure rules,b5?b6in network B will not fail andb6?b7will fail.Finally,a steady state is reached,as shown in Fig.1(b).

    3.Theoretical analysis framework

    3.1.Single-layer network

    The degree distributionp(k)is defined as the probability that a node of degreekis randomly selected in the network.The condition that a node of degreekdoes not belong to the GC is that allkconnectivity links of the node do not lead to the GC.We introduce an auxiliary parameterR,defined as the probability that a randomly selected link leads to the GC along one direction, then the probability that the node does not belong to the GC can be expressed as(1-R)k.For any node with degreeksatisfying the degree distributionp(k),the probability that a node does not belong to the GC is obtained by multiplying the probability (1-R)kwith the degree distributionp(k)and summing up,[34–38]whereSdenotes the probability that a randomly selected node in the network belongs to the GC,

    From the generating function theory,[36,37]it is known that the network degree distribution generating function isΦ0(x) = ∑k p(k)xkand the excess-degree distribution generating function isΦ1(x) = ∑k q(k)xk-1, whereq(k) denotes the excess-degree distribution of the nodes, i.e.,q(k)=p(k)·k/〈k〉,and〈k〉is the average degree of the network,thus Eq.(1)can be rewritten as

    Apart from Eq.(2),in order to solveSorR,it is also necessary to build a consistent equation of the probabilityRthat a connectivity link leads to the GC,i.e.,at least one of the remaining links of the node reached by this link leads to the GC,i.e.,

    Thus Eq.(3)can reduce to

    Using the generating functionΦ1(x), Eq.(4) can be rewritten as

    If a fraction 1-pof connectivity links is randomly removed from the network (bond-percolation), the probability of randomly selecting a node in the network that belongs to the GC remains unchanged since no node is removed from the network,and the expression of the order parameterSremains as Eq.(2).Since only a fractionpof connectivity links remains in the network,the probability of randomly selecting a link to reach a node with degreekbecomesptimes the original one, i.e.,p·q(k).Thus, the consistent equation for the auxiliary parameterRis

    The order parameterSand the auxiliary parameterRcan be solved by combining Eqs.(2)and(6).

    3.2.Edge-coupled interdependent networks with directed dependency links

    Following the single-layer network method,the EINDDL is a two-layer interdependent networks, where we introduce two auxiliary parametersRAandRB,which denote the probabilities that randomly selected link in network A and network B leads to the MCGC along one direction, respectively.Network A has 1-βfraction of connectivity links depending on network B,and network B hasβfraction of connectivity links depending on network A.For simplicity, assuming that both networks have the same degree distribution and the same number of nodes,their generating function expressions are also the same.In this case, when a connectivity link is randomly selected in network A, whether the node reached by this link belongs to MCGC depends on the following two factors.

    (i) If a randomly selected connectivity link does not depend on network B,i.e., the corresponding link in network B will not cause this link to fail.The fraction of this part of connectivity links isβ.It is only necessary to consider whether the end nodes of these links in network A belong to the GC.

    (ii) If a randomly selected connectivity link depends on network B, i.e., the failure of the corresponding connectivity link in network B will cause this link to fail.The ratio of this part of connectivity links is 1-β.It is necessary to consider the probability that the end nodes of this part of links in network A belong to GC, but also the probability that the end nodes of the corresponding links in network B belong to GC.

    Therefore, based on Eq.(5), the probabilityRAAthat a fraction ofβconnectivity links leads to the GC is

    Since a 1-βfraction of connectivity links depends on network B,at least one of the nodes to which the corresponding connectivity links the network B lead must be in the GC,as shown in Fig.2.

    Fig.2.The survival of connectivity link dependent on network B.

    The probability that one end ofg′not in the GC isΦ1(1-RB) and the probability that both ends are not in the GC is [Φ1(1-RB)]2.Thus the probabilityRLBthat at least one end of the connectivity linkg′in the GC is

    Similarly, the node to which the 1-βratio of connectivity links leads has to be in the GC of network A, and its probabilityRABis known according to Eq.(5),

    Accordingly, the probability that the nodes to which the connectivity links of the 1-βfraction lead are in the MCGC has to satisfy the conditions in both networks A and B, i.e.,RAB·RLB.

    Combining the two aspects considered above, the auxiliary parameterRAcan be expressed as

    Similarly, whether the nodes reached by a randomly selected connectivity link in network B belong to the MCGC requires consideration of the above two aspects.Using the same analysis,we can obtain that theβfraction of connectivity links in network B leading to the GC needs to satisfy both constraints in networks A and B,and the probability is

    The 1-βratio of connectivity links in network B leading to the GC only needs to meet the constraint in network B leading to the GC,and the probability is

    Thus, the expression for the auxiliary parameterRBcan be expressed as

    Following the above formulas,if the 1-pfraction of the connectivity links in network A are randomly removed in the initial state,the fraction of remaining connectivity links in network A isp,and the connectivity links ofβand 1-βfractions becomep·βandp·(1-β),respectively.From Eq.(10),the expression of the auxiliary parameterRAat this time is

    Since there are no attacks in network B,what affects network B is the failed links in theβscaled connectivity links.Therefore, the surviving connectivity links in network B are only related to the connectivity links of the 1-βscale and the connectivity links of thep·βscale.From Eq.(13), the expression of the auxiliary parameterRBis

    In addition,after randomly removing the 1-pfraction of connectivity links in network A,the order parameterSwill not be affected by the deletion of connectivity links because the nodes between the two layers are not directly related.From Eq.(2),the expressions for the order parametersSAandSBof networks A and B are,respectively,

    Given the values ofβandp,combing Eqs.(14)–(17),we can solve the values of the auxiliary parametersRA,RBand the order parametersSAandSB.

    In the percolation theory,[16,35]a fraction 1-pof nodes or links are randomly removed which may trigger cascading failures due to the dependency between networks.An abrupt collapse may occur where the size of GC changes to zero discontinuously upon a removal of a critical fraction 1-pcof nodes/links.This discontinuous percolation transition is named as the fist-order phase transition andpcis known as the percolation threshold.If the size of GC changes to zero continuously that is typical on single-layer networks, we call it the second-order phase transition.

    The type of phase transition and the phase transition thresholdpcof the EINDDL can be investigated by comparingSandp.The percolation behaviors of the networks vary with different values of the parameterβ.If the network undergoes a first-order phase transition,the curve of Eq.(14)tangentially will meet the curve represented by Eq.(15).If the network undergoes a second-order phase transition,the value of the auxiliary parameterRAorRBwill be 0,and the second-order phase transition threshold can be found by the above equations.

    3.3.Theoretical framework verification

    Figure 3 gives the plots ofSAandSBversuspfor the RRRR EINDDL and the ER-ER EINDDL with differentβfor the number of nodesNA=NB=100000 and the average degree〈kA〉=〈kB〉=4,where the solid lines are the theoretical results and the dots are the simulated values.

    Fig.3.Plots of SA and SB versus p with different β.

    It can be seen from Fig.3 that the theoretical values can match well with the computer simulation values.The theoretical framework can be verified.We apply this framework on the RR and ER networks to analyze their phase transition behaviors.

    4.Robustness analysis of the EINDDL

    4.1.The RR-RR EINDDL

    In the RR network, all node degrees are〈k〉.Therefore,the degree distribution is

    Assuming that the numbers of nodes of networks A and B are the same and the average degree is〈k〉, the generating function of the degree distribution of the RR network isΦ0(x)=x〈k〉and the generating function of the excess degree distribution isΦ1(x)=x〈k〉-1from the theory of generating functions.According to Eqs.(14)and(15),the expressions for the auxiliary parametersRAandRBof the RR-RR EINDDL are

    From Eqs.(16)and(17),the expressions of the order parametersSAandSBare

    From Eqs.(19)–(22), the values of the order parametersSAandSBare determined by the auxiliary parametersRAandRB, and only non-zero values ofRAandRBcan be obtained for non-zero values ofSAandSB.Obviously, Eqs.(19) and(20)possess a trivial solutionRA=RB=0,corresponding to the case in which the GC is 0 in the network.

    The first-order phase transition thresholdof the system can be obtained when the curves represented by Eqs.(19)and(20)are tangent.The GC of network A in the critical case will be 0 if the system undergoes a second-order phase transition.At this time, the auxiliary parameterRA→0.The nontrivial solutionRB0of the auxiliary parameterRBcan be found by substituting it into Eq.(20), and then the second-order phase transition thresholdpIIcis obtained as follows:

    Assuming that the number of nodes in both network A and network B in the RR-RR EINDDL isN=100000 and the average degree is the same,Fig.4 gives a plot of the order parametersSAandSBversuspat differentβfor〈k〉=3,4,and 5.

    4.2.The ER-ER EINDDL

    The degree distribution of a Erd¨os–Rényi (ER) network is Poisson distribution, i.e.,Similarly, assuming that the number of nodes of networks A and B are the same and the average degree is〈k〉,the generation function of the degree distribution of the ER network and the generation function of the excess degree distribution are the same from the theory of generating functions, that is,Φ0(x)=Φ1(x)=e〈k〉(x-1).According to Eqs.(14)and(15),the expressions for the auxiliary parametersRAandRBare,respectively,

    From Eqs.(16)and(17),the expressions of the order parametersSAandSBare

    Given the values ofβand the average degree〈k〉,Eqs.(24)and(25)can be rewritten as

    Eliminating the auxiliary parameterRBfrom the system of Eq.(28),we obtain

    When Eq.(29)satisfies the following conditions

    Equations (24) and (25) appear as a nontrivial solution,corresponding to the critical case when the system undergoes a first-order phase transition, wherepis the first-order phase transition threshold.

    If the system undergoes a second-order phase transition, the same analytical method of substitutingRA→0 into Eq.(25) can give the solutionRB0of the auxiliary parameterRB,which leads to the second-order phase transition threshold,

    Figure 5 gives a plot of the order parametersSAandSBversuspat differentβfor〈k〉=3,4 and 5.

    Fig.4.Plots of SA and SB versus p for the RR-RR EINDDL at different β and〈k〉.

    4.3.Robustness analysis

    4.3.1.Robustness analysis of network A

    It can be seen that, no matter whether it is the RR-RR EINDDL or the ER-ER EINDDL, network A gradually converges to the first-order phase transition from the second-order phase transition asβincreases,as seen in Figs.4 and 5.Network A will change from the first-order to the hybrid phase transition and then from the hybrid to the second-order phase transition whenβkeeps increasing,as shown in Fig.6.

    This process can be explained by the graphical method.SAis a monotonically increasing function aboutRAwhenkandβare inevitable, so the value ofRAcan directly reflect the phase transition.Taking the RR-RR EINDDLk=4 as an example, make the plot of the curve represented by Eqs.(19)and (20) and label the solution of Eq.(19) (i.e., the value ofRA),as shown in Fig.7.It can be seen thatRAchanges from zero to non-zero value continuously aspincreases in Fig.7(a).Therefore,the second-order phase transition occurs in network A.RAjumps from zero to a non-zero value aspincreases in Fig.7(b), and the first-order phase transition occurs in network A.In Fig.7(c),RAchanges continuously from zero to non-zero values and jumps when the two curves are tangent so that network A undergoes the hybrid phase transition.In Fig.7(d),RAstill changes from zero to non-zero continuously aspincreases,andRAdoes not jump because the tangent(intersection)point of the two curves is precisely on theRAaxis.Thus,network A undergoes the second-order phase transition.

    Fig.5.Plots of SA and SB versus p for the ER-ER EINDDL at different β and〈k〉.

    Fig.6.Plots of SA versus p for the EINDDL when β >0.8.

    Fig.7.Plots of curves represented by Eqs.(19)and(20).

    4.3.2.Robustness analysis of network B

    From Figs.4 and 5, it can be found thatSBin network B gradually changes from a more significant non-zero value to zero whenp <asβincreases,because network A fails completely whenp <and network B is equivalent to the failure of a single network with the initial removalβfraction of connectivity links.

    Next,we investigate the conditions,on which network B completely collapses.WhenSBis non-zero, Eq.(20) intersects theRBaxis in the plane (RA,RB) as shown in Fig.9(a),where the red line indicates Eq.(19) and the blue line indicates Eq.(20).Asβgradually increases to a certain value,Eq.(20)gradually goes from intersecting theRBaxis to passing the(0,0)point, when the order parameterSBis exactly 0 at a certainp, as shown in Fig.9(b), we define theβandpat this case as the critical valuesandwhen network B collapses completely.Subsequently, atβ ≥andp ≤,the value of the order parameterSBis always 0.Equation(20)no longer intersects with theRBaxis,as shown in Fig.9(c).

    Fig.8.Plots of phase transition threshold pAc versus β for network A.

    Fig.9.Plots of Eqs.(19)and(20)with identical p and different β.

    Fig.10.Plots of SB as a function of p-pc for the EINDDL when β >0.8.

    Then we study what type of phase transition occurs when network B is able to completely collapses.The plots ofSBas a function ofp-pcare made whenβ ≥, as shown in Fig.10.It can be seen thatSBturns from discontinuous to continuous with the increase ofβfrom Fig.10.Hence, asβincreases,network B turns from the first-order to second-order phase transition whenβ ≥.

    The process is similarly explained graphically by labeling the solution of Eq.(20) (i.e., the value ofRB) in Figs.7(b)–7(d).Asβincreases,the value ofRBchanges from discontinuously to a non-zero value to continuously to a non-zero value.Hence,network B changes from the first-order to the secondorder phase transition.Also,from Fig.7(d),it can be seen that when the tangent(intersection)point of the two curves is exactly on theRAaxis,the values ofRAandRBno longer jump,and the second-order phase transition occurs in both network A and network B at this time.

    It is clear that the first-order phase transition occurs in network B whenand the second-order phase transition occurs in network B whenfrom the above study.

    4.3.3.The prevention strategy

    From Subsections 4.3.1 and 4.3.2, it is clear that the robustness of both network A and network B become weaker and then stronger asβincreases from 0 to 1 with the phase transition point of view.The mechanism of this phenomenon is the way of spreading the damages influenced byβ.To illustrate this,we defineFas the ratio of failed links of networkA(B) spreading to networkB(A) to the total number of links of cascading failure in the second and third steps,and make a plot ofFversusβfor EINDDLk=4,p=0.3, as shown in Fig.11.

    Fig.11.Ratio of failure links of A(B)spreading to B(A)with different β.

    The number of failed links of network A (B) cascade to network B(A)increases and decreases asβincreases from 0 to 1.The more failed links there are in each cascading failure step,the more fragile the network is,that is the reason why the robustness of networks A and B decreases and then increases.Therefore,the value ofβcan directly affect the robustness of the EINDDL.Thus the prevention strategy is to find the value ofβwhen the EINDDL is most robust.

    We useSABsimilar to the order parameterSto denote the size of MCGC.Figure 12 shows the plots ofSABas a function ofp.The robustness of the EINDDL can be measured by the phase transition thresholdpcand the jump valuePJat the transition threshold.In order to find the prevention methods to optimize the network design,we study the relationship between robustness measured by both phase transition thresholdpcand the jump valuePJandβ.As we know,the smaller thepcandPJare,the more robust the system is.rising-falling-rising trend with the increase ofβ,and there exists peaks which reflect the strongest cases.Therefore,to prevent mass cascading failures,the value ofβshould be chosen around the high peaks and away from the low peak.

    Fig.12.Plots of SAB as a function of p for the EINDDL.

    5.Conclusions

    In summary, we have established a model for edgecoupled interdependent networks with directed dependency links (EINDDLs), and the directed dependency strength between two networks is determined byβ.The theoretical analysis framework of an EINDDL is developed based on the selfconsistent probability method,and the theoretical analysis and simulation of RR-RR and ER-ER networks are carried out to prove the theoretical analysis.The results show that under random attacks,network A exhibits different types of phase transitions(first-order,second-order and hybrid)with the increase ofβ, and the critical values ofβfor the whole phase transition type change process of network A are obtained.Network B is rarely affected by cascading failure whenβis small,and asβincreases, network B will gradually converge from the first-order to the second-order phase transition.The robustness of both network A and network B becomes weaker and then stronger asβincreases from 0 to 1 with the phase transition point of view.Finally,by measuring the robustness of the EINDDL with the phase transition thresholdpcand the jump valuePJ,we obtain the value ofβwhere the robustness of the EINDDL is best,and propose the prevention strategy.

    Fig.13.Plots of SR versus β for the EINDDL.

    For differentβ,we get two vectorsand.To evaluate the robustness of networks,we forward transit and normalize the vectorsndto getSpcandSPJ,so the comprehensive robustness evaluation metricSRof the EINDDL is expressed in the following,where the weightsωcandωJof the two metrics are received by the entropy method.

    Due to the forward transition and normalized vectorsp′candthe larger theSRis, the stronger the networks are.The plots ofSRversusβfor differentkis shown in Fig.13.It can be seen that the robustness of the EINDDL shows an

    Acknowledgments

    Project supported by the National Natural Science Foundation of China (Grant Nos.61973118, 51741902,11761033, 12075088, and 11835003), Project in JiangXi Province Department of Science and Technology (Grant Nos.20212BBE51010 and 20182BCB22009), and the Natural Science Foundation of Zhejiang Province (Grant No.Y22F035316).

    猜你喜歡
    周杰銀座
    Experimental realization of quantum controlled teleportation of arbitrary two-qubit state via a five-qubit entangled state
    銀座集團-山東商院商貿(mào)流通虛擬仿真實訓(xùn)基地建設(shè)研討會在我校召開
    “一根筋”周杰,特立獨行的娛樂“圈外人”
    May You Always
    山東銀座集團的財務(wù)戰(zhàn)略分析
    新營銷(2019年14期)2019-12-18 01:34:40
    銀座東急購物廣場日本東京
    中鐵·銀座尚筑 迭代公寓 國匠品質(zhì)
    齊魯周刊(2017年30期)2017-08-09 00:39:46
    中鐵·銀座尚筑新品發(fā)布會暨銀座簽約儀式
    齊魯周刊(2017年30期)2017-08-09 00:36:55
    Characteristics of Meteorological Factors over Different Landscape Types During Dust Storm Events in Cele,Xinjiang,China
    臨沭縣玉山鎮(zhèn):小小紫薯“入駐”銀座,變身市民“新寵”
    長江蔬菜(2013年23期)2013-03-11 15:09:39
    欧美精品一区二区免费开放| avwww免费| 母亲3免费完整高清在线观看| 97在线人人人人妻| 亚洲精品国产av蜜桃| 97精品久久久久久久久久精品| 一本大道久久a久久精品| 精品国产乱子伦一区二区三区 | 国产精品一区二区在线观看99| 国产成人系列免费观看| 免费在线观看视频国产中文字幕亚洲 | 黄色视频在线播放观看不卡| 一本色道久久久久久精品综合| 成年美女黄网站色视频大全免费| 91精品三级在线观看| avwww免费| 欧美+亚洲+日韩+国产| 少妇精品久久久久久久| av不卡在线播放| 一区二区av电影网| 宅男免费午夜| 男女高潮啪啪啪动态图| 51午夜福利影视在线观看| 中文字幕人妻丝袜一区二区| 亚洲精品一二三| 亚洲精品国产av成人精品| 亚洲美女黄色视频免费看| 永久免费av网站大全| 啪啪无遮挡十八禁网站| 亚洲熟女精品中文字幕| 欧美午夜高清在线| 日韩视频一区二区在线观看| 精品少妇一区二区三区视频日本电影| 首页视频小说图片口味搜索| 久久人妻熟女aⅴ| 久久久水蜜桃国产精品网| tube8黄色片| 视频区图区小说| 欧美日韩亚洲综合一区二区三区_| 久久狼人影院| 亚洲熟女毛片儿| 一本大道久久a久久精品| 久久青草综合色| 亚洲国产中文字幕在线视频| 日韩制服骚丝袜av| 一区二区av电影网| 亚洲人成电影免费在线| 亚洲人成电影观看| 俄罗斯特黄特色一大片| 99精品欧美一区二区三区四区| 精品久久久久久久毛片微露脸 | 香蕉丝袜av| 最近最新中文字幕大全免费视频| 九色亚洲精品在线播放| 色视频在线一区二区三区| av网站在线播放免费| 精品福利永久在线观看| 美女高潮到喷水免费观看| 国产一卡二卡三卡精品| 老熟妇仑乱视频hdxx| 成人av一区二区三区在线看 | 久久精品熟女亚洲av麻豆精品| 成年人午夜在线观看视频| 91精品国产国语对白视频| 亚洲精品日韩在线中文字幕| 永久免费av网站大全| 国产成人精品在线电影| 欧美激情极品国产一区二区三区| 国产欧美日韩一区二区三 | 性高湖久久久久久久久免费观看| 久久人人爽人人片av| 久久这里只有精品19| 嫁个100分男人电影在线观看| cao死你这个sao货| 99国产精品一区二区蜜桃av | 老司机影院毛片| 国产av一区二区精品久久| 午夜视频精品福利| 久久国产精品影院| 无限看片的www在线观看| 黄色怎么调成土黄色| 国产亚洲欧美精品永久| 国产精品久久久久久人妻精品电影 | 亚洲第一欧美日韩一区二区三区 | 成年人免费黄色播放视频| 国产亚洲欧美在线一区二区| 99精国产麻豆久久婷婷| 性色av乱码一区二区三区2| 捣出白浆h1v1| 搡老岳熟女国产| 99国产极品粉嫩在线观看| 老汉色∧v一级毛片| 亚洲欧美日韩另类电影网站| 亚洲三区欧美一区| 国产亚洲一区二区精品| 精品卡一卡二卡四卡免费| 久久人人爽av亚洲精品天堂| 黄色怎么调成土黄色| 一本综合久久免费| 又黄又粗又硬又大视频| 精品福利观看| 99热全是精品| 亚洲精品中文字幕在线视频| 一区二区三区激情视频| 精品少妇久久久久久888优播| 婷婷丁香在线五月| 欧美精品一区二区大全| 国产黄色免费在线视频| 中亚洲国语对白在线视频| 一级片免费观看大全| 交换朋友夫妻互换小说| 国产野战对白在线观看| 成人影院久久| 亚洲av日韩在线播放| cao死你这个sao货| 亚洲av国产av综合av卡| 国产黄频视频在线观看| 免费在线观看黄色视频的| 少妇人妻久久综合中文| 久久精品亚洲av国产电影网| 国产一区二区三区综合在线观看| 美女高潮喷水抽搐中文字幕| av视频免费观看在线观看| 乱人伦中国视频| 免费在线观看影片大全网站| 激情视频va一区二区三区| 色婷婷久久久亚洲欧美| 午夜日韩欧美国产| 亚洲国产欧美网| 一区二区三区乱码不卡18| 在线观看舔阴道视频| 亚洲专区字幕在线| 亚洲精品成人av观看孕妇| 伊人久久大香线蕉亚洲五| 最新在线观看一区二区三区| 一区二区三区精品91| www.999成人在线观看| 久久av网站| 亚洲少妇的诱惑av| 国产精品麻豆人妻色哟哟久久| av在线老鸭窝| 中文字幕av电影在线播放| 亚洲色图综合在线观看| 精品国产国语对白av| 一级黄色大片毛片| 一二三四在线观看免费中文在| 欧美精品亚洲一区二区| 亚洲 欧美一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 午夜福利影视在线免费观看| 亚洲第一av免费看| 欧美成狂野欧美在线观看| 大片免费播放器 马上看| 亚洲精品美女久久av网站| 亚洲第一青青草原| 亚洲专区中文字幕在线| 国产成人欧美| 婷婷成人精品国产| 国产亚洲一区二区精品| 免费在线观看影片大全网站| 妹子高潮喷水视频| 国产视频一区二区在线看| 99热全是精品| 在线观看www视频免费| 亚洲国产欧美日韩在线播放| 国产亚洲精品一区二区www | 国产精品久久久久久精品电影小说| 老司机亚洲免费影院| 91精品伊人久久大香线蕉| 欧美乱码精品一区二区三区| 亚洲专区字幕在线| 亚洲国产欧美日韩在线播放| 欧美变态另类bdsm刘玥| 成在线人永久免费视频| av超薄肉色丝袜交足视频| 国产欧美日韩一区二区精品| 精品久久久久久电影网| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费午夜福利视频| 满18在线观看网站| 一级,二级,三级黄色视频| 在线观看免费日韩欧美大片| 两人在一起打扑克的视频| av超薄肉色丝袜交足视频| 国产欧美日韩综合在线一区二区| 人人妻人人澡人人看| 伦理电影免费视频| 夜夜骑夜夜射夜夜干| 777久久人妻少妇嫩草av网站| 自线自在国产av| 国产不卡av网站在线观看| 亚洲精华国产精华精| 午夜激情久久久久久久| 一级毛片女人18水好多| 亚洲专区国产一区二区| 精品亚洲乱码少妇综合久久| 成人国产一区最新在线观看| 天天影视国产精品| 国产精品.久久久| 国产精品成人在线| 女警被强在线播放| 色视频在线一区二区三区| 国产精品久久久久久精品电影小说| 午夜视频精品福利| 国产有黄有色有爽视频| 亚洲国产av影院在线观看| 黑人欧美特级aaaaaa片| 狂野欧美激情性bbbbbb| 国产真人三级小视频在线观看| 国产欧美日韩精品亚洲av| 久久久久国内视频| 日本一区二区免费在线视频| 9191精品国产免费久久| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 精品国产一区二区久久| netflix在线观看网站| 婷婷色av中文字幕| 国产精品自产拍在线观看55亚洲 | 激情视频va一区二区三区| 深夜精品福利| 69精品国产乱码久久久| 午夜福利视频精品| 老熟妇仑乱视频hdxx| www.av在线官网国产| 免费观看a级毛片全部| 国产人伦9x9x在线观看| 蜜桃国产av成人99| 一区福利在线观看| 国产有黄有色有爽视频| 麻豆国产av国片精品| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜添小说| 极品少妇高潮喷水抽搐| 最近最新免费中文字幕在线| 免费久久久久久久精品成人欧美视频| 欧美另类一区| 啦啦啦视频在线资源免费观看| 国产在视频线精品| 一边摸一边做爽爽视频免费| 精品国产一区二区三区四区第35| videosex国产| 亚洲熟女精品中文字幕| 亚洲欧洲日产国产| 99久久综合免费| 亚洲中文字幕日韩| 久久久国产欧美日韩av| 19禁男女啪啪无遮挡网站| 久久性视频一级片| 亚洲国产av影院在线观看| 王馨瑶露胸无遮挡在线观看| 嫩草影视91久久| 老司机影院毛片| 亚洲伊人久久精品综合| 中国国产av一级| 国产精品熟女久久久久浪| 国产成人欧美| 欧美亚洲日本最大视频资源| 国产av一区二区精品久久| 老鸭窝网址在线观看| 国产视频内射| 亚洲中文日韩欧美视频| 人妻夜夜爽99麻豆av| 丰满的人妻完整版| 成在线人永久免费视频| 久久午夜亚洲精品久久| 亚洲第一电影网av| 老熟妇乱子伦视频在线观看| 美女黄网站色视频| 狠狠狠狠99中文字幕| 真人一进一出gif抽搐免费| 日本一区二区免费在线视频| 久9热在线精品视频| 亚洲成a人片在线一区二区| 国产成人av激情在线播放| 一本大道久久a久久精品| 午夜福利在线在线| 天堂动漫精品| 午夜福利视频1000在线观看| www.熟女人妻精品国产| 亚洲欧美日韩东京热| 国产三级在线视频| 黄色视频不卡| 久久香蕉激情| 午夜精品一区二区三区免费看| 女同久久另类99精品国产91| 草草在线视频免费看| 少妇裸体淫交视频免费看高清 | 中文字幕精品亚洲无线码一区| 男女那种视频在线观看| 国产乱人伦免费视频| 久久亚洲精品不卡| 国产成人精品无人区| 欧美中文日本在线观看视频| 中文资源天堂在线| 精品人妻1区二区| 亚洲无线在线观看| 久久香蕉激情| 欧美高清成人免费视频www| 欧美最黄视频在线播放免费| 欧美黄色淫秽网站| 成人av在线播放网站| 高清在线国产一区| 国产精品98久久久久久宅男小说| av在线播放免费不卡| 国产高清有码在线观看视频 | 最新美女视频免费是黄的| 欧美激情久久久久久爽电影| 亚洲第一电影网av| 国产探花在线观看一区二区| 国产爱豆传媒在线观看 | 国产成人av教育| 久久精品影院6| 99久久精品国产亚洲精品| 中文字幕av在线有码专区| 国产真人三级小视频在线观看| 久久久精品欧美日韩精品| 性欧美人与动物交配| 99国产精品一区二区三区| 亚洲国产精品sss在线观看| 天天添夜夜摸| 欧美成人免费av一区二区三区| 国内精品一区二区在线观看| 搡老岳熟女国产| 久久久久久久久久黄片| 国产精品久久视频播放| 一本精品99久久精品77| 日韩免费av在线播放| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 欧美成人午夜精品| videosex国产| 国产日本99.免费观看| 在线播放国产精品三级| 成人高潮视频无遮挡免费网站| 又粗又爽又猛毛片免费看| 久久人人精品亚洲av| 99国产精品一区二区蜜桃av| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 亚洲国产欧美人成| 又大又爽又粗| 99国产精品99久久久久| av天堂在线播放| 听说在线观看完整版免费高清| 操出白浆在线播放| 黄色毛片三级朝国网站| 国产三级在线视频| 国内毛片毛片毛片毛片毛片| 在线观看免费午夜福利视频| netflix在线观看网站| 免费看a级黄色片| 69av精品久久久久久| 国产97色在线日韩免费| 亚洲欧美日韩东京热| 他把我摸到了高潮在线观看| 国产精品久久视频播放| 欧美另类亚洲清纯唯美| 亚洲国产欧美网| 老汉色∧v一级毛片| av有码第一页| 国产成人av教育| 久久精品国产综合久久久| 亚洲成人精品中文字幕电影| 欧美性猛交╳xxx乱大交人| 国内精品久久久久久久电影| 韩国av一区二区三区四区| 少妇人妻一区二区三区视频| 久久亚洲真实| 亚洲欧洲精品一区二区精品久久久| 国产精品1区2区在线观看.| 人妻夜夜爽99麻豆av| 欧美午夜高清在线| 桃红色精品国产亚洲av| 国产精品电影一区二区三区| 国产一级毛片七仙女欲春2| 国内少妇人妻偷人精品xxx网站 | 老熟妇乱子伦视频在线观看| 精品国产乱子伦一区二区三区| 一边摸一边做爽爽视频免费| 久久香蕉国产精品| 午夜老司机福利片| 成年版毛片免费区| АⅤ资源中文在线天堂| 国产成人aa在线观看| 色老头精品视频在线观看| 亚洲18禁久久av| 欧美丝袜亚洲另类 | 久久九九热精品免费| 男女那种视频在线观看| 亚洲av成人不卡在线观看播放网| 久久中文字幕人妻熟女| 国产激情久久老熟女| 国产高清视频在线播放一区| 欧美性猛交╳xxx乱大交人| 国产激情久久老熟女| 日韩欧美 国产精品| 久久久久久久精品吃奶| 丁香六月欧美| 成人三级黄色视频| a在线观看视频网站| 脱女人内裤的视频| 村上凉子中文字幕在线| 日韩大码丰满熟妇| 在线观看舔阴道视频| avwww免费| 国产成人aa在线观看| 国内精品一区二区在线观看| 中文字幕人妻丝袜一区二区| 久久久久久亚洲精品国产蜜桃av| 国产精品99久久99久久久不卡| 亚洲国产欧洲综合997久久,| 夜夜夜夜夜久久久久| 成人高潮视频无遮挡免费网站| 男女做爰动态图高潮gif福利片| 精品欧美国产一区二区三| 国产伦在线观看视频一区| 久久99热这里只有精品18| 可以在线观看的亚洲视频| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机午夜十八禁免费视频| 手机成人av网站| 欧美黑人欧美精品刺激| 三级毛片av免费| 亚洲一区中文字幕在线| 九色成人免费人妻av| 精品久久蜜臀av无| 免费在线观看影片大全网站| 叶爱在线成人免费视频播放| 国产在线精品亚洲第一网站| 免费看a级黄色片| 免费看十八禁软件| 国产亚洲精品综合一区在线观看 | 国产精品精品国产色婷婷| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| 两个人免费观看高清视频| 午夜激情福利司机影院| 国产免费男女视频| 亚洲国产欧洲综合997久久,| 日韩精品青青久久久久久| 亚洲黑人精品在线| 麻豆国产av国片精品| 曰老女人黄片| 久久中文字幕一级| 欧美丝袜亚洲另类 | 免费在线观看完整版高清| 桃色一区二区三区在线观看| 精品少妇一区二区三区视频日本电影| 淫秽高清视频在线观看| 中文字幕最新亚洲高清| 精品人妻1区二区| 99riav亚洲国产免费| 99久久99久久久精品蜜桃| 亚洲成人国产一区在线观看| 久久久久性生活片| 小说图片视频综合网站| 久久中文字幕一级| 中文字幕久久专区| av欧美777| 久久久久性生活片| 亚洲一区二区三区色噜噜| 男女之事视频高清在线观看| 日韩大尺度精品在线看网址| 国内揄拍国产精品人妻在线| 午夜福利视频1000在线观看| 中文资源天堂在线| www.999成人在线观看| 免费看日本二区| 大型黄色视频在线免费观看| 色精品久久人妻99蜜桃| 好看av亚洲va欧美ⅴa在| АⅤ资源中文在线天堂| 动漫黄色视频在线观看| 美女大奶头视频| 久久久久九九精品影院| 在线观看日韩欧美| 色精品久久人妻99蜜桃| 中文字幕最新亚洲高清| 日韩欧美 国产精品| 亚洲av熟女| 久久久精品大字幕| 最新在线观看一区二区三区| 激情在线观看视频在线高清| 亚洲精品美女久久久久99蜜臀| 成人亚洲精品av一区二区| 男人舔女人的私密视频| 亚洲精品久久国产高清桃花| 国内精品一区二区在线观看| av在线天堂中文字幕| 亚洲成av人片在线播放无| 亚洲国产欧洲综合997久久,| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 久9热在线精品视频| 在线观看66精品国产| 十八禁网站免费在线| 高清在线国产一区| 亚洲成人精品中文字幕电影| 欧美日韩福利视频一区二区| av免费在线观看网站| 91成年电影在线观看| 亚洲在线自拍视频| 丁香欧美五月| 国产日本99.免费观看| 亚洲专区中文字幕在线| 在线观看免费日韩欧美大片| av免费在线观看网站| 免费搜索国产男女视频| 日韩国内少妇激情av| 天天添夜夜摸| 在线国产一区二区在线| 国产高清视频在线观看网站| 国产aⅴ精品一区二区三区波| 五月伊人婷婷丁香| 免费在线观看亚洲国产| 午夜福利高清视频| 午夜精品在线福利| 亚洲成人免费电影在线观看| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 成在线人永久免费视频| 丁香六月欧美| 中文字幕人妻丝袜一区二区| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩无卡精品| 国产成人av教育| 老司机靠b影院| 欧美成人一区二区免费高清观看 | 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 久久婷婷成人综合色麻豆| 久9热在线精品视频| 久热爱精品视频在线9| 成人午夜高清在线视频| 欧美乱码精品一区二区三区| av有码第一页| 欧美最黄视频在线播放免费| 亚洲国产欧美人成| 日韩欧美在线乱码| 欧美日韩乱码在线| 老司机在亚洲福利影院| 久久亚洲真实| 一级黄色大片毛片| 欧美乱妇无乱码| 久久婷婷成人综合色麻豆| 久9热在线精品视频| 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 久久人人精品亚洲av| 国产69精品久久久久777片 | 神马国产精品三级电影在线观看 | www.自偷自拍.com| 国产区一区二久久| 99re在线观看精品视频| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品第一综合不卡| 女人被狂操c到高潮| 18禁黄网站禁片免费观看直播| 久久久久久亚洲精品国产蜜桃av| 欧美av亚洲av综合av国产av| av片东京热男人的天堂| 国产成人系列免费观看| 首页视频小说图片口味搜索| 少妇被粗大的猛进出69影院| 亚洲在线自拍视频| 九色成人免费人妻av| 久久久久久久久中文| 亚洲欧美一区二区三区黑人| 亚洲午夜精品一区,二区,三区| 18禁裸乳无遮挡免费网站照片| 国内少妇人妻偷人精品xxx网站 | 99热6这里只有精品| 巨乳人妻的诱惑在线观看| 级片在线观看| 老汉色av国产亚洲站长工具| 一级毛片精品| 亚洲av电影在线进入| 又黄又爽又免费观看的视频| 久久性视频一级片| 99国产极品粉嫩在线观看| 校园春色视频在线观看| 婷婷六月久久综合丁香| 欧美乱妇无乱码| 久久欧美精品欧美久久欧美| 亚洲欧美精品综合久久99| 18禁国产床啪视频网站| 日日夜夜操网爽| 一夜夜www| 中文亚洲av片在线观看爽| 亚洲欧美一区二区三区黑人| 视频区欧美日本亚洲| 老司机深夜福利视频在线观看| 亚洲国产精品成人综合色| 精品国产超薄肉色丝袜足j| 99久久精品国产亚洲精品| 国产免费av片在线观看野外av| 欧美最黄视频在线播放免费| 99国产综合亚洲精品| 手机成人av网站| 久9热在线精品视频| 天天躁狠狠躁夜夜躁狠狠躁| 校园春色视频在线观看| 久久久水蜜桃国产精品网| 搞女人的毛片| 香蕉国产在线看| 亚洲国产精品999在线| 天堂影院成人在线观看| 色综合站精品国产| 欧美精品亚洲一区二区| 99久久国产精品久久久| 草草在线视频免费看| 看片在线看免费视频| 久久亚洲真实|