時序超網(wǎng)絡(luò)上重要節(jié)點挖掘方法研究

詹秀秀　余小燕　劉闖　張子柯

摘要：在如何識別時序超網(wǎng)絡(luò)上的重要節(jié)點方面取得了一定的進(jìn)展。定義了該類網(wǎng)絡(luò)上度量節(jié)點重要性程度的 8 個中心性方法及隨機(jī)移除節(jié)點的基線方法，分別側(cè)重于網(wǎng)絡(luò)不同的拓?fù)浣Y(jié)構(gòu)性質(zhì)和時間特征，從多個角度綜合考慮了該類網(wǎng)絡(luò)上節(jié)點的重要性。同時，構(gòu)建了時序超網(wǎng)絡(luò)上的SI 傳播模型，基于該模型提出了新的評估方法來衡量所提出的中心性方法的有效性。研究表明，在時序超網(wǎng)絡(luò)上，基于最快到達(dá)路徑的介數(shù)中心性方法是評價該類網(wǎng)絡(luò)上節(jié)點重要性的良好指標(biāo)。此外，基于時間分辨率的度和超度中心性方法通過尋找網(wǎng)絡(luò)的最佳時間分辨率，可以進(jìn)一步優(yōu)化普通的度和超度中心性方法，彌補(bǔ)了普通方法不能有效考慮網(wǎng)絡(luò)時間信息的缺點，且在多個真實網(wǎng)絡(luò)上表現(xiàn)出與介數(shù)中心性方法相當(dāng)?shù)男阅堋?/p>

關(guān)鍵詞：時序超網(wǎng)絡(luò)；SI 傳播模型；重要節(jié)點；中心性方法

中圖分類號：N 94??????????? 文獻(xiàn)標(biāo)志碼：A

Mining methods of important nodes on temporal hyper-networks

ZHAN Xiuxiu1， YU Xiaoyan1， LIU Chuang1， ZHANG Zike2

（1. Alibaba Research Center for Complexity Sciences， Hangzhou Normal University， Hangzhou 311121， China; 2. College ofMedia and International Culture， Zhejiang University， Hangzhou 310058， China）

Abstract： Some? progress? has? been? made? on? how? to? identify? important? nodes? on? the? temporal? hypernetwork. Eight centrality methods for measuring the importance of nodes on this type of network and abaseline? method? for? randomly? removing? nodes? were? defined，? focusing? on? the? different? topologicalproperties and time characteristics of the network， and comprehensively considering the importance ofnodes on this type of network from multiple perspectives. At the same time， a SI spreading model ontemporal hyper-network was constructed. Based on this model， a new evaluation metric was put forwardto? measure? the? effectiveness? of? the? proposed? centrality? metrics.? The? main? results? are? as? follows.Betweenness metric which is based on the fastest arrival path performs better than the other baselines inidentifying? important? nodes? in? temporal? hyper-networks.? In? addition，? the? degree? and? hyper-degreecentrality which consider time resolution are superior to the degree and hyper-degree based on the staticnetwork topology， which makes up for the shortcomings that the static methods that can not effectivelyconsider? the? network? time? information.? In? addition，? they? show? comparable? performance? with? thebetweenness in multiple real networks.

Keywords： temporal hyper-network; SI spreading model; important nodes; centrality methods

網(wǎng)絡(luò)結(jié)構(gòu)的異質(zhì)性和復(fù)雜性使得在特定觀察角度下，某些節(jié)點相比于其他節(jié)點更具重要性意義，這種重要性可通過中心性來量化。如何衡量網(wǎng)絡(luò)中節(jié)點的中心性是網(wǎng)絡(luò)分析的基礎(chǔ)問題之一，該問題研究目標(biāo)的實質(zhì)是確定網(wǎng)絡(luò)中節(jié)點的重要性程度，并對網(wǎng)絡(luò)中的節(jié)點進(jìn)行重要性排序[1-2]。具有較高中心性值的節(jié)點往往是高影響力節(jié)點，擁有較強(qiáng)的信息擴(kuò)散或者疾病傳播的能力。因此，找到網(wǎng)絡(luò)中的高影響力節(jié)點具有重要的實際意義。以往的中心性研究成果在現(xiàn)實生活中也獲得了廣泛應(yīng)用，如疾病傳播控制[3-4]、廣告精準(zhǔn)投放[5]、交通堵塞疏導(dǎo)[6-8]等。

一般情況下，每個網(wǎng)絡(luò)都有特定的節(jié)點重要性排名，不同的識別方法考慮網(wǎng)絡(luò)的不同結(jié)構(gòu)屬性，會給出不同的排名列表[9]。截至目前，研究者們已經(jīng)提出了許多用于識別網(wǎng)絡(luò)中重要節(jié)點的中心性方法。例如：Zareie等[10]使用熵對有影響力的節(jié)點進(jìn)行排序； Pu 等[11]提出了局部維度來識別重要節(jié)點。 Teng 等[12]追蹤了社交網(wǎng)絡(luò)中的真實信息流，找到了有影響力的傳播者；張喜平等[13]引入 m 階鄰居節(jié)點的概念，提出了一種基于 m 階鄰居節(jié)點重要度貢獻(xiàn)的復(fù)雜網(wǎng)絡(luò)節(jié)點重要度方法；顧亦然等[14]引入影響度概念，用節(jié)點自身的 K- Shell 值與對其鄰居節(jié)點的影響度來表征節(jié)點對鄰居節(jié)點的貢獻(xiàn)重要度。

然而，這些對網(wǎng)絡(luò)節(jié)點中心性的研究僅停留在表示成對節(jié)點交互的普通靜態(tài)網(wǎng)絡(luò)上，當(dāng)網(wǎng)絡(luò)的結(jié)構(gòu)發(fā)生動態(tài)變化時，靜態(tài)網(wǎng)絡(luò)并不能完全準(zhǔn)確地描述各式各樣的現(xiàn)實世界系統(tǒng)，即現(xiàn)實中的交互均是非靜態(tài)的、隨時間不斷發(fā)展和變化的，而時序網(wǎng)絡(luò)可以很好地表達(dá)這一概念[15-18]。另外，相比于普通網(wǎng)絡(luò)所表達(dá)的二元交互，現(xiàn)實世界更普遍存在的是表達(dá)多個個體間聯(lián)系的高階交互，如電子郵件的發(fā)送者和接收者、出版物的共同作者、購物網(wǎng)站共同購買一樣物品的用戶等。放松節(jié)點數(shù)量限制的超網(wǎng)絡(luò)是這一類型交互的理想載體[19-20]。因此，為更準(zhǔn)確地表達(dá)現(xiàn)實世界的復(fù)雜關(guān)系，本研究將這兩類網(wǎng)絡(luò)組合，構(gòu)造了一類全新的網(wǎng)絡(luò)——時序超網(wǎng)絡(luò)。該類網(wǎng)絡(luò)隨著時間的推移，網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)會隨時建立和瓦解，且可形成節(jié)點數(shù)大于2的超邊，其為研究真實復(fù)雜系統(tǒng)的動力學(xué)提供了良好的表示工具[21]。該框架尤其適用于描述社交系統(tǒng)，典型的如物理接觸網(wǎng)絡(luò)和在線社交網(wǎng)絡(luò)。在這兩種網(wǎng)絡(luò)所對應(yīng)的實際場景中，人群間的真實接觸和在線會話等交互會隨著時間的推移發(fā)生變化，此類交互往往是高階的，并在時間上具有突變性和持續(xù)性。已有研究表明，在研究人類時間動力學(xué)時，考慮社會互動的高階結(jié)構(gòu)是非常重要的[21]。時序超網(wǎng)絡(luò)能以超高的時間分辨率描述現(xiàn)實世界動態(tài)演化過程，并從中挖掘到具有高影響力的重要節(jié)點，該過程對于上述問題的解決具有非常重要的參考意義。然而，目前關(guān)于時序超網(wǎng)絡(luò)上節(jié)點中心性評估方法的研究還寥寥無幾，且還未形成系統(tǒng)化的定義。

本文的研究目標(biāo)是在時序超網(wǎng)絡(luò)上尋找重要的節(jié)點，通過構(gòu)建基于時序超網(wǎng)絡(luò)的 SI 傳播模型來模擬網(wǎng)絡(luò)的動態(tài)演化過程，定義時序超網(wǎng)絡(luò)上8個節(jié)點中心性評估方法。此外，通過攻擊不同中心性方法所尋找到的重要性排名較高的節(jié)點來評估所提出的中心性指標(biāo)的有效性，為挖掘時序超網(wǎng)絡(luò)上的重要節(jié)點提供了理論支撐。

1方法

1.1時序超網(wǎng)絡(luò)的定義與表示

為方便討論，假設(shè)記錄網(wǎng)絡(luò)的時間有限，且開始記錄的時間tstart =0，結(jié)束記錄的時間tend = T。在時間窗口[0， T]上，根據(jù)不同時間戳的節(jié)點與超邊的隸屬關(guān)系，構(gòu)建時序超網(wǎng)絡(luò)H0（T）={H0 ， H1 ， ··· ， HT}，即時序超網(wǎng)絡(luò)可以由每個時刻的靜態(tài)超網(wǎng)絡(luò)Ht =（Vt ， Et）來表示，且t e [0， T]。其中，Ht =（Vt ， Et）表示屬于t 時刻的一組節(jié)點 Vt ={v1 ， v2 ， ··· ， vN }和一組超邊Et ={e1 ， e2 ， ··· ， em}構(gòu)成的超網(wǎng)絡(luò)， N表示時序超網(wǎng)絡(luò)中的總節(jié)點數(shù)， m表示Ht中的超邊數(shù)。值得注意的是，任意一條超邊ej可以包含多個節(jié)點，任意一個節(jié)點vi也可以隸屬于多條超邊。在時序超網(wǎng)絡(luò)中，假設(shè)每個時間片中的節(jié)點個數(shù)相同，而超邊的數(shù)量會隨著時間的推移而發(fā)生改變。時序超網(wǎng)絡(luò)若忽略不同交互所發(fā)生的時間，即為H0（T）所對應(yīng)的靜態(tài)超網(wǎng)絡(luò) H =（V， E），其中 V ={v1 ， v2 ， ··· ， vN }，超邊eje E，當(dāng)且僅當(dāng)ej在Et（t e [0， T ]）中至少出現(xiàn)過一次。

圖1（a）為一個由3個時間戳和6個節(jié)點組成的時序超網(wǎng)絡(luò)示例，根據(jù)不同時間戳發(fā)生的交互所形成的不同連接進(jìn)行著色。中間的縱軸線代表時間流動的路徑，曲線連接著一個交互所涉及到的多個節(jié)點，不同顏色代表該時刻的高階交互所形成的不同超邊：如在t1時刻，節(jié)點v1，v2，v3之間的兩兩交互形成一個擁有3個節(jié)點的超邊（3–超邊）；節(jié)點v2和v4之間的交互形成擁有兩個節(jié)點的超邊（2–超邊）。若忽略圖1（a）中時序超網(wǎng)絡(luò)每條超邊出現(xiàn)的時間信息，就可以得到其對應(yīng)的靜態(tài)超網(wǎng)絡(luò)，如圖1（b）所示。

1.2時序超網(wǎng)絡(luò)上節(jié)點中心性方法

為了挖掘時序超網(wǎng)絡(luò)上的重要節(jié)點，定義度量該類網(wǎng)絡(luò)上節(jié)點重要性的一系列中心性指標(biāo)，包括基于節(jié)點的度、超度、接觸數(shù)、節(jié)點之間的距離以及網(wǎng)絡(luò)的時序信息等。

1.3評估方法

針對上述8種時序超網(wǎng)絡(luò)上重要節(jié)點挖掘的中心性指標(biāo)，提出合理可靠的方法來評估這些度量指標(biāo)的有效性。

1.3.1 SI 傳播模型

構(gòu)建時序超網(wǎng)絡(luò)的易感–感染（ susceptible- infected ， SI ）傳播模型。在該模型中，個體被分成兩組：易感者（susceptible）和感染者（infected）。為簡單起見，讓所有傳播過程都從t =0時刻的一個初始被感染的種子節(jié)點開始，其余所有節(jié)點的初始狀態(tài)均被設(shè)置為易感狀態(tài)。按照時間t 推移，依次遍歷每一個時刻，若感染節(jié)點出現(xiàn)在t 時刻的超邊中，就會以β的概率隨機(jī)去感染其鄰居易感節(jié)點，該時刻的節(jié)點感染完成之后，推移到下一時刻進(jìn)行感染，傳播過程在t = T時刻結(jié)束。最后，依次迭代序列中的每一個節(jié)點作為初始被感染的種子節(jié)點，進(jìn)行上述傳播過程[24]。用φ（t）表示t 時刻以每個節(jié)點作為傳播源的平均被感染節(jié)點總數(shù)。

1.3.2節(jié)點中心性評估方法

采用不同的中心性方法作為節(jié)點攻擊策略來評估節(jié)點重要性方法的有效性。具體步驟如下： a.對于給定的節(jié)點中心性方法，首先用該方法對網(wǎng)絡(luò)中的節(jié)點進(jìn)行重要性排序； b.刪除重要性排名前k的節(jié)點序列，構(gòu)造攻擊節(jié)點序列后的時序超網(wǎng)絡(luò)； c.將原始時序超網(wǎng)絡(luò)和攻擊節(jié)點序列后的時序超網(wǎng)絡(luò)分別執(zhí)行傳播概率β下的 SI 傳播模擬，分別得到兩條傳播曲線φ0（t）和φk（t）；d.比較兩條傳播曲線的差異，即計算兩條傳播曲線間的面積差。如圖2所示，用φ0和φk之間面積差（藍(lán)色標(biāo)出）評估所提出的中心性指標(biāo)的有效性，評價指標(biāo)的數(shù)學(xué)定義如下：

E （k）越大，說明攻擊該中心性方法挖掘到的重要節(jié)點能有效抑制傳播，同時也驗證了該方法的有效性。除了上述8種中心性指標(biāo)，還將隨機(jī)刪除包含 k個節(jié)點的序列作為對比方法（用 Random 表示），隨機(jī)實驗的平均次數(shù)為50次。

2實驗

2.1數(shù)據(jù)集

選擇9個真實的時序超網(wǎng)絡(luò)數(shù)據(jù)集，個體的交互覆蓋醫(yī)院、會議、藝術(shù)館、工作場所、研究機(jī)構(gòu)及學(xué)校等。通過這些網(wǎng)絡(luò)來探究不同的中心性方法在具有不同社會背景和拓?fù)浣Y(jié)構(gòu)性質(zhì)的真實場景中的適用性。 EEU1和 EEU2網(wǎng)絡(luò)描述了歐洲大型研究機(jī)構(gòu)兩個不同部門所有成員之間電子郵件的時間網(wǎng)絡(luò)[25]。HT2009網(wǎng)絡(luò)是 ACM Hypertext 2009會議期間的時間網(wǎng)絡(luò)[26]。Gallery 網(wǎng)絡(luò)記錄的是一家藝術(shù)館的參觀人員的時間網(wǎng)絡(luò)[27]。Hospital 網(wǎng)絡(luò)描述的是法國里昂醫(yī)院病房中患者、患者和衛(wèi)生保健工作者以及醫(yī)護(hù)人員之間接觸的時間網(wǎng)絡(luò)[21]。Workplace 網(wǎng)絡(luò)描述的是2015年在法國辦公樓測量的個人之間接觸的時間網(wǎng)絡(luò)[28]。HS2011和 HS2012網(wǎng)絡(luò)描述的是2011年和2012年法國馬賽一所高中學(xué)生之間的聯(lián)系時間網(wǎng)絡(luò)[29]。PS 網(wǎng)絡(luò)是兒童和教師之間的聯(lián)系時間網(wǎng)絡(luò)[29]。這些網(wǎng)絡(luò)的基本拓?fù)湫再|(zhì)如表1所示。其中： N 為節(jié)點數(shù)；C為接觸數(shù)； T 為網(wǎng)絡(luò)的最大時間戳；< f >為節(jié)點平均度；< g >為節(jié)點平均超度；< d >為網(wǎng)絡(luò)的平均最快到達(dá)距離。表2展示的是時序超網(wǎng)絡(luò)中擁有不同度數(shù)的超邊規(guī)模。表中，超邊的度表示一條超邊中包含的節(jié)點數(shù)，這里只統(tǒng)計至網(wǎng)絡(luò)中超邊的度為10的超邊數(shù)量，“—”表示網(wǎng)絡(luò)中不存在該度數(shù)的超邊。從表2可以觀察出，隨著超邊的度的增加，其超邊數(shù)量是逐漸減少的，大多數(shù)的真實時序超網(wǎng)絡(luò)中的交互局限于度為2，3，4的超邊之間，而郵件網(wǎng)絡(luò)的超邊規(guī)模偏大。

2.2節(jié)點重要性指標(biāo)的性能比較

上文定義了不同時間分辨率下的度和超度，為了探索不同時間分辨率的中心性指標(biāo)的性能，將基于時間分辨率的中心性方法的性能定義為從 k =5到k =30范圍內(nèi)有效性E（k）曲線下的面積。以EEU2網(wǎng)絡(luò)為例，在圖3（a）中顯示了D 和W方法的有效性結(jié)果，并將這兩種方法的性能作為對照組（圖3（b）中的 static 組），改變網(wǎng)絡(luò)的時間分辨率?t ，從100增至1000，計算基于不同時間分辨率的度（ D_t）和超度中心性（ W_t）方法的性能。如在圖3（b）中，D_t方法的最佳時間分辨率為?t =100，W_t方法的最佳時間分辨率為?t =300。按照這種想法，分別在9個真實網(wǎng)絡(luò)中，尋找不同方法的最佳時間分辨率，結(jié)果如表3所示。如果在某一?t情況下D_t和W_t方法的性能值大于 static 組中的 D 和W方法的性能，則該網(wǎng)絡(luò)的最佳時間分辨率在表3中用粗體進(jìn)行標(biāo)注，可以發(fā)現(xiàn)所有網(wǎng)絡(luò)均可以通過該方法尋找到最佳的時間分辨率，同時在絕大多數(shù)網(wǎng)絡(luò)上D_t和W_t方法的性能比 D 和 W方法更好。因為這兩種方法側(cè)重于考慮不同時期網(wǎng)絡(luò)中節(jié)點的重要性，同時結(jié)合了網(wǎng)絡(luò)的時間性質(zhì)和高階性質(zhì)。

在傳播中選取β=1來評估不同中心性方法的效果，結(jié)果如圖4和圖5（a ）所示。值得說明的是，D_t和W_t方法選取的均是表3中所列出的最佳時間分辨率。在圖5（ a ）中，以 Hospital 網(wǎng)絡(luò)為例，將9個中心性方法在該網(wǎng)絡(luò)中所得到的性能值組成一個序列X，并將其按照最大最小值方法進(jìn)行歸一化，即Xi（*）=（xi-xmin）/（xmax-xmin）。式中： Xi表示序列 X 的第i個數(shù)據(jù)； Xi（*）表示序列 X 的第i個數(shù)據(jù)歸一化后得到的結(jié)果；Xmax和Xmin分別表示序列 X 中的最大值和最小值。歸一化后，得到對應(yīng)的序列X*，即為9個中心性方法的平均性能序列，在圖5（ a ）中以不同顏色的柱子進(jìn)行顯示。圖4和圖5（ a ）的結(jié)果表明，在時序超網(wǎng)絡(luò)上定義的中心性方法中， B 方法的有效性和平均性能均明顯優(yōu)于其他中心性方法?？赡茉蛴?個： a. B 方法的節(jié)點介數(shù)中心性值越高，說明節(jié)點更大范圍地連接著不同時刻的超邊，在網(wǎng)絡(luò)中充當(dāng)橋梁的作用越明顯； b. B 方法在計算距離時是基于節(jié)點的最快到達(dá)路徑的思想，因此在考慮網(wǎng)絡(luò)全局特征的同時也考慮了時間的全域性，更加適用于時序超網(wǎng)絡(luò)。另外，圖4出現(xiàn)了兩個有趣的現(xiàn)象。首先，在 Gallery 網(wǎng)絡(luò)中， Q方法表現(xiàn)優(yōu)異，這可能是因為 Gallery網(wǎng)絡(luò)的節(jié)點平均度較小，平均最快到達(dá)距離較長，如表1所示。但在其余真實的時序超網(wǎng)絡(luò)中該方法的有效性并不是很高，尤其在Workplace 網(wǎng)絡(luò)中， Q 方法的有效性甚至和對比方法的接近，這說明 Q 方法在度量時序超網(wǎng)絡(luò)上節(jié)點的重要性時并不具備普適性。其次，擬合的曲線會出現(xiàn)下降的趨勢，如 Gallery網(wǎng)絡(luò)中的V方法和 EEU2網(wǎng)絡(luò)中的Q方法，此現(xiàn)象意味著這些方法挖掘到的節(jié)點在網(wǎng)絡(luò)中的重要性程度并不高。而在圖5（ a ）中， Q ，U ，V等方法的效果不如其他方法。最后，可以清晰地觀察到，在大多數(shù)網(wǎng)絡(luò)中D_t和W_t方法的平均性能比D和W方法的更好。這說明通過尋找最佳時間分辨率，可以進(jìn)一步有效優(yōu)化普通的 D 和W方法。同時，驗證了在β=0.5的情況下，不同中心性方法及對比方法之間的平均性能與上述β=1情況下的實驗結(jié)果類似，如圖5（b）所示。

2.3肯德爾相關(guān)性

肯德爾相關(guān)性系數(shù)τ，可以衡量兩個不同變量之間的相關(guān)性[30]。假設(shè)對于兩個隨機(jī)變量 Y和Z ，它們的第i個組合用（Yi ， Zi）表示。當(dāng)Yi>Yj和Zi>Zj或Yi

式中： n+表示具有一致性的序列個數(shù)； n一表示不一致的序列個數(shù)。τ=1表示不同方法得到的節(jié)點重要性排序列表相同；τ=0則表示相反的情況。τ的值越大，說明兩種中心性方法的關(guān)系越相關(guān)。采用該方法來檢測定義的不同中心性方法之間的相關(guān)性。

圖6顯示了9個真實時序超網(wǎng)絡(luò)不同中心性方法之間的肯德爾相關(guān)性結(jié)果。首先，可以觀察到 D ，D_t，W ，W_t這4個方法之間的顏色最深，即它們之間的相關(guān)性最高。雖然這4個方法都是側(cè)重于度量網(wǎng)絡(luò)中節(jié)點關(guān)于度的信息來判斷節(jié)點的重要性，但是從圖5（ a ）可以看出，D_t，W_t方法會明顯比D ，W方法的平均性能更高，這也進(jìn)一步說明度量節(jié)點中心性時考慮網(wǎng)絡(luò)時間信息的重要性。其次， U和V方法都是從節(jié)點的接觸次數(shù)出發(fā)，判斷節(jié)點的重要性時考慮的方向是相似的，所以這兩種方法之間的相關(guān)性程度很高，且在圖4中，這兩種方法的曲線重合度也很高，再次說明了這兩種方法具有很高的相似性。但是，節(jié)點的接觸數(shù)和接觸鄰居數(shù)同節(jié)點的度和超度是有差別的，因為它們考慮了節(jié)點出現(xiàn)的時間信息，統(tǒng)計了節(jié)點在不同時刻出現(xiàn)的次數(shù)，所以U和V方法同D ，D_t，W ，W_t方法之間的相關(guān)性程度減弱。最后， B 和Q方法分別從節(jié)點位于網(wǎng)絡(luò)中的位置和節(jié)點到其他節(jié)點之間的最快到達(dá)距離等方面來衡量節(jié)點的重要性，因此，它們之間的相關(guān)性程度略低，同其他方法之間的相關(guān)性也類似。

3總結(jié)與展望

構(gòu)建了一種全新類型的網(wǎng)絡(luò)——時序超網(wǎng)絡(luò)，其同時兼具網(wǎng)絡(luò)的高階性質(zhì)和時間性質(zhì)，能更加準(zhǔn)確、有效地描繪真實世界的社會系統(tǒng)。此外，基于節(jié)點的最快到達(dá)路徑計算了時序超網(wǎng)絡(luò)上節(jié)點之間的最快到達(dá)距離。接下來，針對挖掘該類網(wǎng)絡(luò)上重要節(jié)點的中心性方法匱乏且不成體系等問題，進(jìn)一步提出了8個中心性方法和評估指標(biāo)，并在9個不同領(lǐng)域的真實網(wǎng)絡(luò)上進(jìn)行了實驗研究。結(jié)果表明：提出的所有中心性方法中，介數(shù)中心性方法總是具有優(yōu)越性的，其能有效地挖掘到時序超網(wǎng)絡(luò)上具有高影響力的重要節(jié)點；在基于時間分辨率的中心性方法中，通過尋找網(wǎng)絡(luò)的最佳時間分辨率，可以進(jìn)一步優(yōu)化普通的度和超度中心性方法，解決了在時序超網(wǎng)絡(luò)上計算度或者超度時不能有效地考慮網(wǎng)絡(luò)時間信息的問題，想法更加合理，實驗結(jié)果也更加有效。同時，基于時間分辨率的度和超度中心性方法也表現(xiàn)得極具競爭力。綜上所述，在真實世界中為了快速且有效地控制疾病和信息的傳播，一方面可以通過尋找具有“樞紐”作用的節(jié)點并中斷其傳播能力，另一方面節(jié)點出現(xiàn)的時間信息以及其在不同時間窗口的重要性也是不可忽略的關(guān)鍵因素。

本文基于時序超網(wǎng)絡(luò)定義的一系列中心性方法，一方面重點考慮了節(jié)點的度、超度和接觸數(shù)、節(jié)點之間的路徑和距離等信息，另一方面結(jié)合了網(wǎng)絡(luò)的時間性質(zhì)和高階性質(zhì)，力求尋求一種能囊括這兩方面且效果優(yōu)異的方法。因此，該工作彌補(bǔ)了目前在時序超網(wǎng)絡(luò)上挖掘重要節(jié)點的中心性方法稀少且不成體系等問題，為后續(xù)深入研究時序超網(wǎng)絡(luò)奠定了較為堅實的基礎(chǔ)。

然而，在時序超網(wǎng)絡(luò)上挖掘重要節(jié)點的中心性方法仍存在許多潛在問題：首先，本文所提出的中心性方法還未結(jié)合節(jié)點的更多方面特征來綜合度量節(jié)點的重要性；其次，在尋找網(wǎng)絡(luò)的最佳時間分辨率時，本文粗略地將網(wǎng)絡(luò)按照某一數(shù)值進(jìn)行劃分，這可能會遺漏一些重要信息，故可以根據(jù)網(wǎng)絡(luò)時間訊息的不同，按小時、天數(shù)、月份、年份等信息來劃分網(wǎng)絡(luò)，當(dāng)然，這不可避免地會增加網(wǎng)絡(luò)信息收集工作的復(fù)雜度和任務(wù)量。未來研究中，可以根據(jù)此類網(wǎng)絡(luò)的特點改變或補(bǔ)充一些因素，使方法更具普適性和有效性。

參考文獻(xiàn)：

[1]郭曉成， 馬潤年， 王剛.復(fù)雜網(wǎng)絡(luò)中節(jié)點重要性綜合評價方法研究[J].計算機(jī)仿真， 2017， 34（7）：264–268.

[2] WEN T， DENG Y. Identification of influencers in complex networksbylocalinformationdimensionality[J]. Information Sciences， 2020， 512：549–562.

[3] BELL D C， ATKINSON J S， CARLSON J W. Centrality measuresfordiseasetransmissionnetworks[J]. Social Networks， 1999， 21（1）：1–21.

[4] ZENGQ，LIUY，TANGM，etal. Identifyingsuper- spreaders in information-epidemic coevolving dynamics on multiplex networks[J]. Knowledge-BasedSystems， 2021，229：107365.

[5]陳崢， 高翔， 梁恒，等.基于復(fù)雜網(wǎng)絡(luò)的在線廣告投放方法[P].中國， 201610741623.3.2017-01-04.

[6]李樹彬， 吳建軍， 高自友， 等.基于復(fù)雜網(wǎng)絡(luò)的交通擁堵與傳播動力學(xué)分析[J].物理學(xué)報， 2011， 60（5）：050701.

[7]周艷， 李妍羲， 江榮貴， 等.交通擁堵與預(yù)警信息交互傳播動力學(xué)分析[J].地球信息科學(xué)學(xué)報 ， 2017， 19（10）：1279–1286.

[8]邢國正， 杜萍萍， 方晗琦， 等.基于復(fù)合復(fù)雜網(wǎng)絡(luò)交通運(yùn)輸系統(tǒng)研究[J].南陽理工學(xué)院學(xué)報， 2017， 9（4）：54–57.

[9] LI Z， HUANG X Y. Identifying influential spreaders by gravity model considering multi-characteristics of nodes[J]. Scientific Reports， 2022， 12（1）：9879.

[10] ZAREIE A， SHEIKHAHMADI A， FATEM A. Influential nodesrankingincomplexnetworks： Anentropy-based approach[J]. Chaos，Solitons & Fractals， 2017， 104：485–494.

[11] PU J， CHEN X W， WEI D J， et al. Identifying influential nodes basedonlocaldimension[J]. EurophysicsLetters，2014， 107（1）：10010.

[12] TENG X， PEI S， MORONE F， et al. Collective influence of multiplespreadersevaluatedbytracingrealinformation flow in large-scale social networks[J]. Scientific Reports，2016， 6（1）：36043.

[13]張喜平， 李永樹， 劉剛， 等.節(jié)點重要度貢獻(xiàn)的復(fù)雜網(wǎng)絡(luò)節(jié)點重要度評估方法[J].復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)， 2014， 11（3）：26–32，49.

[14]顧亦然， 王兵， 孟繁榮.一種基于 K-Shell 的復(fù)雜網(wǎng)絡(luò)重要節(jié)點發(fā)現(xiàn)算法[J].計算機(jī)技術(shù)與發(fā)展 ， 2015， 25（9）：70–74.

[15]李佳佳.動態(tài)網(wǎng)絡(luò)中心性方法分析[D].成都：西安電子科技大學(xué)， 2012：12–15.

[16]董永強(qiáng)， 祖倩倩， 陶樺.基于真實數(shù)據(jù)集的動態(tài)網(wǎng)絡(luò)時間中心性分析[J].華中科技大學(xué)學(xué)報（自然科學(xué)版）， 2015， 43（5）：109–113.

[17]付凱， 夏靖波， 趙小歡.動態(tài)融合復(fù)雜網(wǎng)絡(luò)節(jié)點重要度評估方法[J].哈爾濱工業(yè)大學(xué)學(xué)報， 2017， 49（10）：112–119.

[18]田亮.復(fù)雜網(wǎng)絡(luò)的統(tǒng)計性質(zhì)及其上動力學(xué)分析[D].南京：南京航空航天大學(xué)， 2012.

[19]胡楓， 趙海興， 馬秀娟.一種超網(wǎng)絡(luò)演化模型構(gòu)建及特性分析[J].中國科學(xué)：物理學(xué)力學(xué)天文學(xué)， 2013， 43（1）：16–22.

[20]王眾托.關(guān)于超網(wǎng)絡(luò)的一點思考[J].上海理工大學(xué)學(xué)報，2011， 33（3）：229–237.

[21] CENCETTI G， BATTISTON F， LEPRI B， et al. Temporal propertiesofhigher-orderinteractionsinsocial networks[J]. Scientific Reports， 2021， 11（1）：7028.

[22]汪小帆， 李翔， 陳關(guān)榮.復(fù)雜網(wǎng)絡(luò)理論及其應(yīng)用[M].北京：清華大學(xué)出版社， 2006.

[23] AKSOYSG，JOSLYNC，MARREROCO，etal. Hypernetwork science via high-order hypergraph walks[J]. EPJ Data Science， 2020， 9（1）：16.

[24] ZHANXX，HANJALICA，WANGHJ. Suppressing informationdiffusionvialinkblockingintemporal networks[C]//ProceedingsoftheEighthInternational Conference on Complex Networks and Their Applications. Cham： Springer， 2019：448–458.

[25] PARANJAPE A， BENSON A R， LESKOVEC J. Motifs in temporalnetworks[C]//ProceedingsoftheTenthACM International Conference on Web Search and Data Mining. Cambridge： ACM， 2017：601–610.

[26] ISELLA L，STEHL? J， BARRAT A， et al. What's in a crowd？ Analysisof face-to-facebehavioralnetworks[J]. Journal of Theoretical Biology， 2011， 271（1）：166–180.

[27] VAN DEN BROECK W， QUAGGIOTTO M， ISELLA L， et al. The making of sixty-nine days of close encounters at the science gallery[J]. Leonardo， 2012， 45（3）：285–285.

[28] ZHANGSL，ZHAOXY，WANGHJ. MitigateSIR epidemicspreadingviacontactblockingintemporal networks[J]. Applied Network Science， 2022， 7（1）：2.

[29] ZHAN X X， LIU.C，WANGZP，etal Measuring and utilizingtemporalnetworkdissimilarity[DB/OL].（2021–11–02）. https：//arxiv.org/abs/2111.01334.

[30] LVLY，CHENDB，RENXL，et.alVitalnodesidentificationincomplexnetworks[J]. PhysicsReports， 2016， 650：1–63.

（編輯：丁紅藝）