基于Tsallis 熵的复杂网络节点重要性评估方法*

杨松青 蒋沅† 童天驰 严玉为 淦各升

1) (南昌航空大学信息工程学院,南昌 330063)

2) (南京理工大学自动化学院,南京 210094)

复杂网络中节点重要性的评估是网络特性研究中的一项重要课题,相关研究具有广泛的应用.目前提出了许多方法来评估网络中节点的重要性,然而大多数方法都存在评估角度片面或者时间复杂度过高的不足.为了突破现有方法的局限性,本文提出了一种基于Tsallis 熵的复杂网络节点重要性评估方法.该方法兼顾节点的局部和全局拓扑信息,综合考察节点的结构洞特征和K 壳中心性,并充分考虑节点及其邻域节点的影响.为了验证该方法的有效性,本文采用单调性指标、SIR 模型和Kendall 相关系数作为评价标准,在8 个来自不同领域的真实网络上与其他方法进行比较.实验结果表明,此方法能更有效和准确地评估网络节点的重要性,可以显著区分不同节点的重要性.此外,该方法的时间复杂度仅为 O(n2),适用于大型复杂网络.

1 引言

随着网络信息技术的发展与进步,对复杂网络的研究已经在许多领域受到关注.关键节点作为复杂网络的核心要素之一,对网络的结构和功能有着重要影响,被广泛的研究与应用[1-3].例如,通过对电力网络中某关键的调度中心和变电所进行维保,能保证整个电力网络正常运转的稳定性与高效性[4];通过对病毒传播网络中超级传播源实施隔离,能够大幅降低病毒的传播速度和减小扩散范围[5];通过对社交网络中核心的言论加以控制,能够有效促进或抑制信息的传播[6].因此,评估网络中节点重要性具有重大意义.

在不同研究背景下,研究人员提出了众多节点重要性评估方法.目前,度中心性[7]是评估网络节点重要性的常用方法,该方法简单但不够精准.介数中心性[8]和接近中心性[9]分别考虑的是节点间的最短路径数量和平均最短距离,此类方法从全局信息的角度出发,效果明显但计算时间复杂度较高,不适用于大型复杂网络.Kitsak 等[10]依据网络中节点处于网络的核心位置往往有较高影响力的思想,提出用K 壳分解法确定网络中节点的位置,此方法时间复杂度低,适用于大型复杂网络,然而其区分度过于粗粒化.为了提高节点区分度,Zeng 等[11]提出了混合度分解(mixed degree decomposition,MDD)法,以网络中剩余的邻居节点以及删除的邻居节点的混合度进行K 壳分解.Bae和Kim[12]在改进K 壳算法时提出了邻域核心的概念,利用节点自身及其邻居的K 壳值之和来评估网络中节点的重要性.Hou 等[13]考虑多指标评估,基于欧拉距离公式计算节点的度值、介数、K 壳等3 个不同指标的组合度量,更为综合地评估了节点的重要性.王凯莉等[14]提出了一种基于多阶邻居壳数的向量中心性方法,利用向量的形式来表示节点在网络中的相对重要性.Burt 等[15]提出经典社会学中的“结构洞”理论,设计了网格约束系数来量化节点形成结构洞时所受到的约束,该方法利用了节点的局部信息,具有良好的精度和较低的计算时间复杂度.苏晓萍和宋玉蓉[16]提出一种基于节点及其邻域结构洞的评估节点重要性的方法,该方法综合考虑了节点及其邻域的度值和拓扑结构,能有效评估网络中节点的重要性.韩忠民等[17]提出了一种融合结构洞等7 个度量指标的方法,该方法利用基于ListNet 的排序学习方法,能够较为全面地评估网络中节点的重要性.

当使用熵来评估复杂网络时,复杂网络的结构越有序,熵值就越小,反之亦然.同时,熵也可用来描述整个网络结构的复杂性和复杂网络的统计特性.例如,Chen 等[18]提出了结构熵来量化复杂网络的结构特征.Zhang 等[19]提出了一种基于非广延统计力学的结构熵来度量网络的复杂性.黄丽亚等[20]提出了一种基于网络节点在K步内可达的节点总数的K-阶结构熵,用于度量复杂网络的异构性.Wang 等[21]提出了一种基于K 壳和节点信息熵的改进算法,用于区分上层外壳和下层外壳中节点的重要性.

受到上述研究工作的启发,本文提出了一种基于Tsallis 熵的复杂网络节点重要性评估方法(TSM).该方法不仅考虑了节点的结构洞特征和K 壳中心性,还充分考虑了节点自身及其邻域节点的影响,在兼顾节点的局部和全局拓扑信息的同时,采用Tsallis 熵来度量网络结构的复杂性.时间复杂度分析表明,该方法的时间复杂度仅为O(n2),适用于大型复杂网络.为了说明该方法的有效性和适用性,本文选用了8 个来自不同领域的真实网络,并采用了7 个现有的节点重要性评估方法作为比较方法.在此基础上,利用单调性指标、SIR 模型和Kendall 相关系数τ说明了该方法的优越性以及不同方法之间的关系.

2 理论基础

2.1 网络定义

设图G=(V,E) 是1 个无向无权网络,其中V={v1,v2,···,vn}表示G中节点的集合,|V|=n,E={(vi,vj)|vi,vj ∈V}表示G中边的集合,|E|=m.记An×n=(aij)n×n为G对应的邻接矩阵,如果节点vi和节点vj之间存在连边关系,则aij=1,否则aij=0 .

2.2 度中心性

度中心性(degree centrality,DC)[7]是衡量节点影响力的最简单指标.节点拥有的度数越大,该节点的影响力就越大.节点vi的DC 定义为

其中,n为网络中节点的∑数量.于是,节点vi的邻接度可以定义为,其中,Γ(i)为节点vi的邻居节点的集合.

2.3 介数中心性

介数中心性(betweenness centrality,BC)[8]是基于节点的最短路径数定义的全局指标.通常,节点的BC 越大,该节点的影响力就越大.节点vi的介数中心性定义为

其中,σst为网络中任意两个节点vs和vt之间的最短路径数目,为经过节点vi的最短路径数目.

2.4 混合度分解方法

混合度分解方法(mixed degree decomposition method,MDD)由Zeng 等[11]提出,该方法有效提高了K 壳分解方法的区分度.它通过同时考虑节点的剩余度值和节点的已移除度值对K 壳分解过程进行了改进,网络按照新的混合度值继续分层,节点vi的混合度定义为

其中,λ是0 到1 之间的可调参数.当λ=0 时,MDD方法与K 壳分解方法一致,但当λ=1 时,它等同于DC.在本文中,设定了参考文献[11]中使用的参数λ=0.7 .

2.5 全方位距离

Hou 等[13]提出了全方位距离(all-around distance,AAD)指标,其综合考虑了节点的度值、BC、K 壳等3 个不同指标的影响,采用欧拉距离公式计算了这3 个指标的组合度量,计算公式如下所示:

其中,KS(i) 为节点vi的K 壳值.

2.6 改进的K 壳方法

改进的K 壳方法(improved K-shell method,IKS)由Wang 等[21]提出,该方法结合了K 壳分解方法和节点信息熵.它按照K 壳分层从大到小的顺序进行迭代,在每层未被选中节点中选择节点信息熵最高的节点,直至网络中所有节点都被选中为止,节点vi的节点信息熵定义为

2.7 结构洞理论

结构洞理论(structural holes theory)是Burt等[15]在1992 年研究社交网络中的竞争关系时提出的.在社会网络中,结构洞是指个体之间为非冗余关系,且存在互补信息的差距.为了量化结构洞节点对这些关系的控制,Burt 采用约束系数来测量结构洞中节点形成的约束.当系数越小,节点越容易形成结构洞,节点的影响力越大.节点vi的约束系数为

其中,节点vq是节点vi和节点vj的共同邻居节点;Lij表示在节点vi的所有邻居节点中节点vj所占的权重比例.

在计算节点的约束系数时,考虑到节点及其邻域的度值和拓扑结构对节点重要性的影响将Lij定义为

2.8 Tsallis 非广延统计力学

对于1 个有限的离散概率集,Boltzmann-Gibbs熵[22]被定义为

其中,k是玻尔兹曼常数,pi满足标准化条件.

1988 年,基于Boltzmann-Gibbs 熵[22]和Shan non 熵[23],Tsallis[24]提出了一种新的熵,它是一种更一般的形式,能够描述广延的和非广延的系统,被定义为

(9)式中的q-对数函数为

根据(11)式可以将(9)式改写为

其中,N是子系统的数量.将非广延参数q用于描述子系统之间相互作用,当q=1 时,Tsallis 熵将退化为Boltzmann-Gibbs 熵.

3 基于Tsallis 熵的节点重要性评估方法

3.1 方法构造

基于Tsallis 熵的节点重要性评估方法构造的基本思想是在综合考虑节点的结构洞特征和K 壳中心性的同时,利用Tsallis 熵度量网络结构的复杂性,评估网络中节点的重要性.由于结构洞特征能体现节点的局部拓扑信息,而K 壳中心性表征了节点的全局拓扑信息,结合这两种拓扑信息,利用Tsallis 熵度量网络结构复杂性,并通过对每个节点的自身及其邻域节点的综合考察,可以给出1 个合理的评估网络中节点重要性的方法,即本文提出的TSM 方法.

基于Tsallis 非广延统计力学,可以得到Tsallis熵的如下形式:

因为结构洞约束系数是基于网络局部拓扑信息的评估指标,考虑到相邻节点的结构洞特征对该节点重要性的影响,可以得到Pij的定义:

由于K 壳中心性是描述节点在网络中所处位置的全局指标,因此它可以用来表示每个节点与整个网络之间的关系以及不同节点之间的关系.本文运用MDD 方法得到节点vi的K 壳中心性 Ksi及其与网络中所有节点的K 壳中心性数值之和的比值.进而,可得到节点vi的Tsallis 参数qi:

(17)式的目的是使参数qi＞1,以反映子系统对网络的次可加性影响.在信息论中,不同的q值表示不同子系统对Tsallis 熵的可加性不同,即q <1,q=1 ,q ＞1 分别代表超可加性、可加性和次可加性.扩展到复杂网络中,整个网络可以看作1 个系统,每个节点和相应的连边可以看作网络的子系统.

利用Tsallis 熵度量网络结构复杂性得到T(vi)后,再结合节点自身的结构洞特征,可得

Bae 和Kim[12]在改进K 壳算法时提出了邻域核心的概念,通过对所有邻居的K 壳值进行求和来衡量网络中节点的重要性.本文借鉴了这一概念,充分考虑邻域节点的影响.对于节点vi,将其所有相邻节点的IC 值进行相加,得到节点vi的邻域核心(核心邻域中心性) Cnc(vi):

在(19)式的基础上,对其所有相邻节点的邻域核心值Cnc 进行求和,可以得到节点vi的扩展邻域核心Cnc+(vi),最终得到了节点vi的重要性评估值TSM(vi):

3.2 时间复杂度分析

本节将对本文提出的TSM 算法进行计算时间复杂度分析,其具体计算步骤如表1 所列.由表1可知,TSM 要遍历网络中每个节点,在计算节点的结构洞约束系数值RCi时,考虑了节点的度及其与邻居拓扑结构间的关系,故时间复杂度应为O(n2);在计算Tsallis 参数qi时,采用了MDD 方法得到节点的K 壳重要性指标,其时间复杂度为O(n);在计算节点的Ti值时,需要考虑节点的邻居节点信息,则时间复杂度为O(n*k) ;而在对IC 进行累加时,时间复杂度为O(n) .由此可以看出,TSM 算法的时间复杂度取决于节点的结构洞约束系数值RCi的计算,因此本文提出的TSM 算法的时间复杂度为O(n2) .其中,n为网络中节点的总数,k为网络的平均度.

表1 计算步骤Table 1.Step of the calculation.

表2 列出了几种不同方法的时间复杂度,并按照评估方法所包含的网络信息将这些方法分成了局部、全局、混合等3 个类型.其中,m为网络中连边的总数.当n

图3 不同传播率下SIR 和各评估方法之间的Kendall 相关系数,图中黑色虚线为传播阈值 βth (a) Karate 网络;(b) Dolphins网络;(c) Polbooks 网络;(d) Jazz 网络;(e) USAir 网络;(f) C.Elegans 网络;(g) Email 网络;(h) PowerGrid 网络Fig.3.Kendall correlation coefficient between SIR and the evaluation methods under different transmission rates,the black dashed line in the figure is the propagation threshold βth :(a) Karate network;(b) Dolphins network;(c) Polbooks network;(d) Jazz network;(e) USAir network;(f) C.Elegans network;(g) Email network;(h) PowerGrid network.

表2 不同方法的时间复杂度Table 2.Time complexity of different methods.

4 实验数据与结果分析

4.1 实验数据

本文选用了8 个来自不同领域的真实网络进行实验,分别是:空手道俱乐部社交网络Karate network[25]、海豚社交网络Dolphins network[26]、美国政治书籍网络Polbooks network[27]、爵士音乐家网络Jazzmusic network[28]、美国航空运输网络USAirline network[29]、秀丽隐杆线虫网络C.Elegans network[30]、大学电子邮件网络Email network[31]、美国电力网络PowerGrid network[1].各网络的统计特征如表3 所列.其中n为节点总数,m为连边总数,〈k〉为平均度,c为聚类系数,〈d〉为平均最短路径长度,βth为SIR 模型传播率阈值,β为实际传播率取值.

表3 网络的统计特征Table 3.Statistical characteristics of the networks.

4.2 有效性分析

为了验证本文提出的TSM 算法的有效性和准确性,首先以Karate 网络为例进行仿真分析,其拓扑图如图1 所示.选用了DC[7]、BC[8]、MDD[11]、邻域结构洞方法[16](N-Burt)、核心邻域中心性[12](Cnc+)、AAD[13]以及IKS[21]等7 种方法与本文提出的TSM 方法进行对比,得到网络前15 个关键节点的重要性评估结果如表4所列.

图1 Karate 网络的拓扑结构Fig.1.Topological structure of the Karate network.

通过结合Karate 网络拓扑结构(图1)与重要性评估结果(表4)分析可知:节点v4和节点v32拥有的邻居数量均为6 且其所处的网络信息传输位置相同,而在其他方法中节点v4和节点v32的重要性是不一样的;节点v9,v14,v24拥有的邻居数量相同均为5,而它们存在的邻域结构洞差值较大,因此,仅仅考虑节点拥有的邻居数量和其所处的网络信息传输位置无法进一步区分节点的重要性.采用N-Burt 方法评估节点重要性时,未考虑到相邻节点的结构洞特征对节点重要性的影响,如节点v32比节点v4存在更多的结构洞,但节点v4的相邻节点较节点v32的相邻节点在网络中是更为重要的“桥接”节点,信息能通过它们更快地扩散至网络.BC 是一种能够很好地体现节点在网络信息传播过程中的重要性的方法,但是对于一些不在任何一条最短路径上的节点以及信息负载能力相似的节点,如节点v6,v7,v8,v28等仍需要进一步考虑节点的局部信息来评估,而且此方法的计算时间复杂度较高.

表4 Karate 网络的节点重要度评估结果Table 4.Node importance evaluation results in the Karate network.

表5 不同评估方法的单调性指标MTable 5.Monotonicity index M of different evaluation methods.

本文所提出的TSM 方法克服了以上方法的不足,它从节点的局部和全局拓扑信息两个角度出发,综合考虑节点自身的结构洞特征和网络位置信息以及邻域节点的影响,并基于Tsallis 熵度量网络结构的复杂性,可提高节点重要性的区分度,准确地评估节点重要性.如表4 所列,通过TSM 方法评估出的前15 个重要节点与其他方法评估出的重要节点基本相同,说明TSM 方法能够有效地评估节点的重要性.在TSM 方法中,节点v4比节点v32更为重要,这很好地体现了相邻节点的结构洞特征对节点的影响,同时也说明本方法克服了BC的不足;节点v9,v14,v24和节点v8,v28,v30的重要性得到了进一步的区分,其中,节点v9的重要性高于节点v14,这说明了本方法既有效地识别了网络中的“桥接”节点,亦克服了DC 的不足.通过对比,本文提出的TSM 方法能有效地评估节点重要性,可进一步提高节点重要性的区分度.

4.3 评价标准

接下来介绍验证所提出的TSM 方法准确性的评价标准.它们分别是单调性指标[12]、SIR 模型[32]和Kendall 相关系数[33].

本文采用单调性指标[12]M(R) 来量化不同重要性评估方法的分辨率,如下所示:

其中,R为评估方法得到的每个节点的重要性排序向量,n为网络中节点的个数,nr为重要性相同的节点的数量.M(R) 的取值介于0 到1 之间.如果M(R)=1,则排序向量R是完全单调的,说明该评估方法能将网络中所有节点的重要性完全区分;如果M(R)=0,则表示网络中所有节点的重要性相同.

目前,大多数学者采用标准的SIR 模型[32]来检测信息和病毒的传播能力.在SIR 模型中,网络中的每个节点可以处于以下3 种状态中的1 种:易感状态S (susceptible)、感染状态I (infected)以及恢复状态R (removed).在传播初始阶段,网络中只有节点vi处于感染状态,其余节点均处于易感状态.在每个时间步中,处于状态I 的节点以传播率β尝试感染处于S 状态的邻居节点,同时以概率μ恢复成状态R 后,不再被感染.传播过程结束的标志是网络中不再出现I 态节点.在整个SIR 传播过程结束时,将网络中处于状态R 的节点的数量S(i)看作节点vi的传播能力.为不失一般性,本文设定恢复率μ=1 .在设定传播率时,如果传播率β过小或过大,都会导致传播过程无法正常进行,从而无法准确衡量每个节点的传播能力.因此,本文设定传播率阈值为βth≈〈k〉/〈k2〉.其中〈k〉为网络平均度,〈k2〉为网络二阶邻居平均度.为保证传播过程正常进行,本文的传播率β在传播率阈值βth附近取值.

为了评估节点重要性评估方法的精准度,本文使用Kendall 相关系数[33]τ来度量每个评估方法得到的节点重要性排序列表R与从SIR 模型获得的节点传播能力排序列表σ之间的相关性.假设存在两个包含n个节点的序列X和Y,其中X=(x1,x2,···,xn)和Y=(y1,y2,···,yn) ,将序列X和序列Y中的元素一一对应构造1 个新的序列XY=((x1,y1),···,(xi,yi),···,(xn,yn)).而对于序列XY中的任意一对元素XYi=(xi,yi)和XYj=(xj,yj),若xi ＞xj且xi ＞yj,或xi

其中,C和D分别表示同序对和异序对的数量.τ的取值介于—1 和1 之间,τ=1 时表示两个序列完全正相关,即相关系数τ越接近1,精确度越高;反之,则表示完全负相关.

4.4 实验分析

本节记录并分析了不同节点重要性评估方法在不同真实网络中的实验结果.

首先选用了DC,BC,MDD,N-Burt,Cnc+,AAD 以及IKS 等7 种评估方法作为对比方法,利用单调性指标[12]M考察了7 种对比方法与本文所提出的TSM 方法的分辨率,即对节点重要性的区分度.表5 记录了不同评估方法在8 个真实网络中的分辨率,并在图2 中更为直观地呈现了它们之间的差异.从图2 可以看出,TSM方法在8 个真实网络中表现十分出色.此外,Cnc+和IKS 方法的区分度也很好.与其他7 种方法相比,TSM 方法大多数情况下可以给出更高M值,而且在所有网络中M(TSM)都非常接近甚至等于1.因此,TSM 方法能够更好地区分节点的重要性.

图2 不同评估方法的单调性指标M (a) Karate 网络、Dolphins 网络、Polbooks 网络和Jazz 网络;(b) USAir 网络、C.Elegans 网络、Email 网络和PowerGrid 网络Fig.2.Monotonicity index M of different evaluation methods:(a) Karate network,Dolphins network,Polbooks network and Jazz network;(b) USAir network,C.Elegans network,Email network and PowerGrid network.

其次,使用SIR 模型得到节点的传播能力排序σ,并计算了不同评估方法的排序与σ之间的Kendall 相关系数τ.相关性越高,对节点重要性的评估就越准确.由于SIR 模型迭代过程的随机性,本文对每个节点重复模拟该过程,然后取其平均值.模拟将遵循以下规则:对于网络|V| ＜ 100,模拟迭代104次;对于 100<|V|<104,模拟迭代103次.表6 展示了8 个真实网络在选定传播率下各评估方法的Kendall 相关系数,SIR 模型的实际传播率取值如表3 所列.从表6 可以看出,TSM 方法的相关系数值在0.8 到1 之间.同时,与其他方法相比,TSM 方法的相关性在8 个真实网络中都最为显著.因此,TSM 方法的准确性得到了初步的验证.

表6 选定传播率下SIR 和各评估方法之间的Kendall 相关系数Table 6.Kendall correlation coefficient between SIR and evaluation methods under a certain transmission rate.

为了进一步验证所提方法的准确性和可靠性,对比了不同传播率下各评估方法与SIR 模型排序结果之间的Kendall 相关系数,如图3 所示.从图3可以看出,当传播率β在传播阈值βth附近取值时,TSM 方法的相关系数除了在Karate 网络中略低于Cnc+方法以外,在大多数情况下都高于其他方法,与SIR 模型得到的节点传播能力有显著的相关性.同时,可以看出当传播率β远大于或者远小于传播阈值βth时,在除Karate,USAir 和PowerGrid网络以外的5 个真实网络中,DC,MDD 以及AAD等方法表现出了比TSM 方法更好的相关性.而BC的相关性在8 个真实网络中表现最差.因此,TSM方法能够更为准确地评估节点的重要性.

然后,对Kendall 相关系数的考察节点范围进行了调整,将考察的节点从整体节点变为部分节点,进一步研究了各个评估方法得到的不同比例排序靠前的节点与节点传播能力排序的相关性结果.在8 个真实网络中进行实验,设置节点比例L的变化范围为0.05—1.如图4 所示,给出不同比例节点下各评估方法的相关系数值.从图4 可以看出,本文提出的TSM 算法在不同比例节点下都能取得不错的节点重要性评估结果,且在更大范围的L值下可以取得更好的评估结果.

图4 不同比例节点下各评估方法的Kendall 相关系数 (a) Karate 网络;(b) Dolphins 网络;(c) Polbooks 网络;(d) Jazz 网络;(e) USAir 网络;(f) C.Elegans 网络;(g) Email 网络;(h) PowerGrid 网络Fig.4.Kendall correlation coefficient of the evaluation methods under different proportion nodes:(a) Karate network;(b) Dolphins network;(c) Polbooks network;(d) Jazz network;(e) USAir network;(f) C.Elegans network;(g) Email network;(h) PowerGrid network.

最后,分析了TSM 方法与其他评估方法之间的关系,如图5 所示.图中每个点代表复杂网络中的1 个节点,点的颜色表示由SIR 模型得到的节点传播能力,横轴和纵轴分别为TSM 方法和其他方法得到的网络中各节点重要性数值.当图中的节点呈线性分布时,说明TSM 方法与该方法有强的相关性.从图5 中的节点分布情况可以发现,TSM与DC,MDD,N-Burt 之间存在强的相关性,这说明TSM 方法既充分考虑了节点的DC 和位置信息,而且有能力识别到网络中的“桥接”节点.此外,TSM 与Cnc+,AAD,IKS 之间节点分布呈较明显的线性关系,这说明节点重要性排名总体相似,亦证明TSM 方法能有效评估节点重要性.从图5 还可以发现,有一些节点尽管其BC 数值大,但其传播能力并不是很强,因此仅靠BC 评估节点的重要性并不准确.所以在评估网络中节点的重要性上,TSM 方法较其他方法更有优势.

图5 TSM 与其他评估方法之间的关系 (a) Karate 网络,β=0.25;(b) Dolphins 网络,β=0.15;(c) Polbooks 网络,β=0.1;(d) Jazz 网络,β=0.04Fig.5.Relationship between TSM and other evaluation methods (a) Karate network,β=0.25;(b) Dolphins network,β=0.15;(c) Polbooksnetwork,β=0.1;(d) Jazz network,β=0.04.

5 结论

本文提出了一种TSM 方法,有效地对复杂网络中节点重要性进行评估和排序.TSM 方法兼顾局部拓扑信息和全局拓扑信息,结合了节点的结构洞特征和K 壳中心性,利用Tsallis 熵度量网络结构的复杂性,弥补了现存方法评估角度片面的不足,优化了可用资源的利用,使信息得以有效传播.该方法在考虑节点的邻域结构和位置信息的同时,有能力识别到网络中的“桥接”节点,从而对复杂网络中节点的重要性进行有效的评估.在8 个真实网络上的实验结果表明,与其他评估方法如DC,BC,MDD,N-Burt,Cnc+,AAD 和IKS 相比,TSM 方法能更好地区分节点重要性的差异,在不同比例节点下都能准确地评估节点的重要性.此外,TSM 方法的时间复杂度仅为O(n2),适用于大型复杂网络.寻找一种结合网络结构和传播动力学的更有效的方法来评估网络中节点的重要性仍然是一个长期的挑战.未来的工作将集中在如何应用所提出的TSM 方法来评估多层网络中的节点重要性.