Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2018: 147.55

Open Access
See all formats and pricing
More options …
Volume 15, Issue 1


Volume 13 (2015)

Research on Critical Nodes Algorithm in Social Complex Networks

Xue-Guang Wang
  • Corresponding author
  • Department of Computer Science, East China University of Political Science and Law, Shanghai 201620, China
  • Department of Computer Science, University College London (UCL), LondonWC1E 6BT, United Kingdom of Great Britain and Northern, Shanghai, Ireland
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2017-03-16 | DOI: https://doi.org/10.1515/phys-2017-0008


Discovering critical nodes in social networks has many important applications and has attracted more and more institutions and scholars. How to determine the K critical nodes with the most influence in a social network is a NP (define) problem. Considering the widespread community structure, this paper presents an algorithm for discovering critical nodes based on two information diffusion models and obtains each node’s marginal contribution by using a Monte-Carlo method in social networks. The solution of the critical nodes problem is the K nodes with the highest marginal contributions. The feasibility and effectiveness of our method have been verified on two synthetic datasets and four real datasets.

Keywords: Critical Node Problem; Complex Networks; Community Structure; Influence Maximization

PACS: 01.20.+x; 07.05.Mh

1 Introduction

The analysis of social networks has many important applications and has attracted more and more institutions and scholars. How to identify influential individuals from a network is called a critical node problem (CNP), and this is one of the basic problems in the field of social network analysis. An effective solution for this problem will have an important practical value [1]. For example, we can better target and isolate the source of the disease and block its spread and diffusion in the disease network; or companies may identify individuals with influence and let them recommend products to their friends so that a cascade spreads by the greatest extent.

Domingos and Richardson first studied the CNP as an algorithmic problem and presented a greedy approximation algorithm. Currently, the research on the CNP of social networks mainly concentrates on improving the greedy algorithm or developing new methods by graph theory and complex network theory. Considering the widespread community structure in social networks, this paper gives a solution for the CNP, which assigns a marginal contribution for every node in a community of a social network using the solution concept and union concept of cooperative games. Then we sort all nodes by their contribution and obtain critical nodes according to some rules.

2 Background

The models for information propagation in networks have been widely studied [2]. We consider two models in this paper: independent cascade model (ICM) and linear threshold model (LTM) [2].

1) ICM Model. In this model, a propagation probability pu,v is given for each edge (u, v) ∈ E, that is, vertex v is activated with probability pu,v by u. When an initial set A0 of active vertices is given, the diffusion process spreads up according to the following randomized rule. When a vertex u is activated at time-step i, it has a single chance of activating its neighbor v with pu,v. If u succeeds, v will become active at time-step i + 1. Here, if v has multiple parent vertices that become active at time-step i for the first time, then their activation attempts are sequenced in an arbitrary order. Whether or not u succeeds, it cannot make any further attempts to activate v in subsequent steps. The process runs until no more activation is possible.

2) LTM Model. In this model, vertex v is influenced by each neighbor u according to a weight wv,u such that ∑v wv,u ≤ 1. Each vertex has a predefined threshold θv ∈ [0, 1], which is chosen uniformly at random. When an initial set A0 with active vertices is given, the diffusion process unfolds according to the following randomized rule. All activated vertices at time-step i still keep active at time-step i + 1. Whether or not any inactivated vertex is activated is determined by its neighbors’ weights such that Σv wv,uθv. The process runs until no more activation is possible.

The difference between ICM and LTM is that (for ICM,?) each attempt of activation is independent of the attempts by all the other active individuals, while in the latter model each inactive individual is influenced by the aggregated weight of all its active neighbors.

A network is modeled as a graph G = (V, E) with vertices in V modeling the individuals in the network and edges in E modeling the relationship between individuals. A vertex has two states: active and inactive, which in this context means whether a product or idea is accepted by individuals or not. We assume that a vertex can only change from inactive to active and not vice versa; an inactive vertex can be activated by its active neighbor vertices and an active vertex can activate its inactive neighbors; the increment of activated vertices represents the dissemination of information.

The definition of the CNP is as follow:

Input: G = (V, E) and an integer k;

Output: A = arg maxSV,|S|k σ(S)

For k - CNP, we hope to find out a set A with k elements and maximize σ(A) which is a NP-hard (define). But, Nemhauser et al. have approved that there is a greedy algorithm (see Algorithm 1) which approximates the optimum within a factor of (1 - 1/e)4.

Algorithm 1 Greedy Algorithm

1: A = Φ

2: for i = 1 to k

3:   vi = arg maxuV\A(σ(A ⋃ {u}) - σ(A))

4:   A = A ⋃{vi}

5: end for

There is a key problem in how to compute the value of σ(A) in Algorithm 1. Currently, we do not have any efficient method to obtain its exact solution. However, we can use a Monte-Carlo method to simulate the process of influence spread for obtaining approximate results by high probability.

3 The Algorithm

3.1 Introduction to the algorithm

Given a finite set of players N, a cooperative game with transferable utility as a pair (N, v), characteristic function v : 2N →R and v(∅) = 0. For ∀iN, if payoff vector satisfies xiv({i}) and i=1Nxi=v(N), then it is so-called an allocation of (N, v). The solution of a cooperative game is a kind of allocation rule and the allocated payoff for every player denotes a method to measure the negotiation strength of the players in the game. Shapely (cite) presented a solution concept which determines the only allocation distribution scheme from the solutions with different property, i.e., it assigns the player’s payoff according to the importance of every player for the game [5]. The Shapely value of the player i in the game (N, v) is Shi(v)={SN|iS}(ns)!(s1)!n!(v(S{i})v(S)),iN

where n = |N| and s = |S|.

However, the Shapely value does not consider the impact of coalition structure and Owen extends it [6]. Each union obtains its payoff from the game between the unions, and then the payoff is allocated by the internal game among the members of the union. All the payoffs are computed by the Shapely value.

Assuming that N = {1, 2,..., n} and M = {1,..., m}, a partition P = {N1, N2,..., Nm} is a coalition structure on N. Let Nk a union and ⋃1≤k≤m Nk = N. When lk, NiNk = ∅. For iN, k(i) denotes the index of the union containing player i, so k(i) is defined by the relation iNk(i). For kM and SNk, the game S is defined by v^S(Q)=v(hQNh),kQv(hQ{k}NhS),kQ

where QM.

The game v¯k is defined by v¯k(S)=Shk(v^S), then the Owen value of the player i ∈ N in the game v¯k(i) is Owi(v,P)=Shi(v¯k(i)).

3.2 A new algorithm

Because the community structure is prevalent in social networks [7], we respectively consider the community’s influence on the information diffusion and every node’s influence in the community. We take the nodes in the social network as the players in the cooperative game and information diffusion as coalition formulation. We can map the information diffusion into a cooperative game with transferable utility which is formalized in the literature [8]. So, we can use the Owen value to obtain the marginal contribution of every node. Because the Owen value can be seen as a two-step procedure in which the Shapely value applies twice, we firstly compute a node’s Shapely value.

Given a node i ∈ N and a subset SN such that iS,the marginal contribution of node i is v(S ⋃ {i}) - v(S), ∀SN\{i}. Consider the set of all possible permutations Ψ on N, let ψΨ and define Si(ψ) to be the set of all nodes appearing before node i in the permutation ψ. So, the average marginal contribution of node i to the given coalition game is 1n!ψΨ[v(Si(ψ){i})v(Si(ψ))].

Note that, this method must work with n! permutations and its computational complexity is O((n/e)n)9. Therefore, we give the approximate method for computing the Shapely value. Randomly generated t-sets Ψt with t permutations, see Algorithm 2, let ψ ∈ Ψt and ψ(i) denotes the ith node in the permutation. The number of activated nodes after running the diffusion model when the node ψ(1) is activated is the contribution of ψ(1). Next, we consider the node ψ(2). If ψ(2) becomes active after ψ(1) is activated, then the contribution of ψ(2) is 0. Otherwise, the contribution of ψ(2) is the number of activated nodes by ψ(2). Therefore, we can obtain the contributions of ψ(3),..., ψ(n). For ψΨt, we repeat the above process R times. Then the average contribution of each node in the diffusion process can be calculated. We can obtain the top-k nodes by sorting by the greatest influence and ensure that they are not adjacent to each other. See Algorithm 3 and Algorithm 4.

Algorithm 2 CREAT_Ψt(n)

1: t=random(n)

2: m=random(n)

3: Select_t (n)

4: for i = 0 to m do

5:  arr[i] = random() % (n-i)

6:  for j=0 to i do

7:   if arr[j] <= arr[i] then

8:    arr[i]=arr[i]+1

9:   end if

10:  end for

11:  v = arr[i]

12:  for k = i-1 to j do

13:   arr[k+1] = arr[k]

14:   k=k-1

15:  end for

16:  arr[j] = v

17: end for

18: for h=2 to t do

19:  |Ψh -1| =arr

20:  n=n-m

21:  m=random(n)

22:  Select_t(n)

23: end for

24: | Ψh | =arr

Algorithm 3 Shapely Value(v, R)

1: n=|ψ|

2: t=|CREAT_Ψt|

3: tmp[1...n]=0

4: ShV[1...n]=0

5: for i=1 to t do

6:  for r=1 to R do

7:   for j=1 to n do

8:    tmp_b=v(Sj(ψi)⋃ {j}) - v(Sj(ψi))

9:    tmp[j]=α tmp[j] + (1-α) tmp_b

10:   end for

11:  end for

12: end for

13: for i=1 to n do

14:  ShV[i]=tmp[i]/t

15: end for

Algorithm 4 TOP_K(V)

1: TopK[1...k]=0

2: AsceSort(V);

3: TopK[1]= V[1]

4: i=1

5: j=2;

6: while i < k do

7:  if V[j] is not adjacent to TopK[1... i-1] then

8:   TopK[i]= V[j]

9:   i=i+1

10:  end if

11:  j=j+1

12: end while

According to the description in Section 3.1, we provide the calculation method of the Owen value. Roger Guimerà, et al studied node roles in a community according to within-module degree z and participation coefficient P, which are divided into seven categories Role = {R1, R2, R3, R4, R5, R6, R7}10. We consider two types of the roles: Non-hub connector node (z < 2.5 and 0.62 < P ≤ 0.80, R3) and Connector hub (z ≥2.5 and 0.30 < P ≤0.75, R6). Most of the nodes which belong to these two roles connect to other communities.

We use the CNM algorithm [11] to divide the network G = (V, E) into l communities C = {C1, C2,..., Cl|Ci = (Vi, Ei), i = 1,...,l} and assign the role rvihRole for every node vihVi in the community. Let Role’ = {R3, R6}, i = 1, ..., l, j = 1, ..., l, V={vih,vjg|(vih,vjg),vihVi,rvihRole,ij}, E={eih,(vih,vjg)|eih=(viha,vihb)Ei,vihaVi,vihbVi,rvihRole,ij},

define G′ = (V′, E′) as the community game network. So, we can obtain the Shapely value of every node in the network G′ and take the sum of the Shapely values of all nodes in the same community as the community’s payoff. Then we treat a community as a separate network and calculate the Shapely value of every node in the community. The node’s Owen value is assigned according to the normalized Shapely value of the node in the community and the community’s payoff. So, we can obtain the k critical nodes by Algorithm 4.

4 Experiments

We validate our method using two synthetic datasets and four real network datasets. All experiments are executed using a PC with 3.2GHz CPU, 4G Memory and Windows 7. The development tools are Matlab 2009 and Microsoft Visual Studio 2010.

4.1 Datasets

We use the BA model [12] and the Forest Fire model [13] to generate two synthetic datasets with 5000 nodes. The BA model takes the power-law distribution and has two important features: growth and preferential attachment. The model can generate a scale-free network whose power exponent is 3 and community structure is not obvious. The Forest Fire model can generate a network with degree power law, densification law, shrinking diameter, and obvious community structure.

The real datasets include DBLP, Facebook, Enron and Youtube. The DBLP dataset [14] constructs a co-author network with 143276 nodes and 359812 edges according to the papers published in the conferences and journals of the computer field from 1997 to 2006. The Facebook dataset [15] provides the friendship network of the New Orleans area with 60567 users and 583766 connections obtained from Jan. 1, 2007 to Dec. 31, 2008. The Enron dataset [16-19] is an E-mail network with 36692 nodes and 367662 edges. The Youtube dataset [20, 21] consists of 35468 nodes and 261191 edges obtained from Jan. 1, 2007 and Jan. 15, 2007.

4.2 The results

We compare our method CNP-OV with the greedy algorithm GREEDY (Algorithm 1), the degree-heuristic algorithm DEGREE and the random algorithm RANDOM. DEGREE selects k nodes with the greatest degree from the network as the initial set; RANDOM randomly selects k nodes as the initial set. In order to obtain the accurate influence of every algorithm, we use the average number of the activated nodes after running ICM and LTM 10000 times for every initial set. In ICM, the propagation probability is set to 0.05; in LTM, the edge-weight of a node is the reciprocal of the node’s degree. The size of the initial set is respectively set from 1 to 20. In the experiments, we only discuss the case of using the ICM because we obtain the same conclusions for ICM and LTM.

First, we consider the influence of community structure on the CNP-OV. We respectively compute the Shapely value and Owen value of every node in the BA and FF datasets and obtain the initial set and the number of activated nodes after running the ICM. The process is repeated 100 times and the average number of the nodes activated by the initial set with different size is drawn in Figure 1. The dashed Sv denotes the method based on the Shapely value and the solid line Ov is the method based on the Owen value. The Sv is almost the same as the Ov for the BA dataset with an unobvious community structure (Figure 1a). In contrast, the Ov is significantly better than the Sv for the FF dataset with an obvious community structure (Figure 1b).

The influence of community structure on our method.
Figure 1

The influence of community structure on our method.

Secondly, we analyze our method’s accuracy. We respectively use GREEDY, DEGREE, RANDOM and CNP-OV to determine the initial sets from above for four real datasets. Figure 2 describes the number of nodes activated by the initial set with different size based on ICM. The accuracy of the CNP-OV is similar to GREEDY, sometimes even better than it(Figure 2b). Compared with DEGREE and RANDOM, the CNP-OV provides a large advantage.

Performance Analysis
Figure 2

Performance Analysis

5 Conclusions

Considering the widespread community structure in social networks, for solving the CNP, in this paper we present a new algorithm for discovering critical nodes based on cooperative games. It assigns a marginal contribution for every node in a community of social networks using the solution concept and union concept of cooperative games. We sort all nodes by their contribution and obtain critical nodes according to selected rules. We have validated its feasibility and effectiveness. The next step will be to study how to improve the time efficiency of the CNP-OV.


This work is supported by the National Social Science Foundation of China (No.11BFX125), PuJiang Talent Project, Peak of law subject construction Project and Public Security Discipline Construction Foundation.


  • [1]

    He N., Li D.Y., Gan W.Y., Zhu X., Mining vital nodes in complex networks, Computer Science, 2007, 34(12): 1-5. Google Scholar

  • [2]

    Even-Dar E., Shapira A., A note on maximizing the spread of influence in social networks, WINE 2007, LNCS 4858, 2007: 281-286. Google Scholar

  • [3]

    Shetty J., Adibi J., Discovering Important Nodes through Graph Entropy The Case of Enron Email Database, The 3rd international workshop on Link discovery, Chicago, Illinois, 2005:74-81.Google Scholar

  • [4]

    Nemhauser G., Wolsey L., Fisher M., An analysis of the approximations for maximizing submodular set functions, Mathematical Programming, 1978, 14(1): 265-294.CrossrefGoogle Scholar

  • [5]

    Shapley L.S., A value for n-person games, In: Kuhn H W, Tucker A W (Eds.), Contributions to the Theory of Games II, Princeton University Press, 1953, 307-317. Google Scholar

  • [6]

    Owen G., Values of games with a priori unions, In: Henn R, Moeschlin O (Eds.), Essays in mathematical economics and game theory, Springer-Verlag, Berlin, 1977, 76-88. Google Scholar

  • [7]

    Scott J., Social Network Analysis: A Handbook (2nd ed), Sage, London, 2000. Google Scholar

  • [8]

    Jean-Francois C., Network Games as TU Cooperative Games: The Core, the Shapley Value and Simple Network Games, 2009, http://centres.fusl.ac.be/CEREC/document/seminars/caulier_cerec_feb2009.pdf 

  • [9]

    Cormen T.H., Leiserson C.E., Rivest R.L., Stein C., Introduction to Algorithms, 2nd ed, Cambridge, MA: MIT Press, 2001.Google Scholar

  • [10]

    Guimerà R., Amaral L.A.N., Functional cartography of complex metabolic networks, Nature, 2005, 433:895-900.CrossrefGoogle Scholar

  • [11]

    Clauset A., Newman M.E.J., Moore C., Finding community structure in very large networks, PHYSICAL REVIEW E, 2004, 70(6): 066111 (6). CrossrefGoogle Scholar

  • [12]

    Barabási A.L., Albert R., Emergence of scaling in random networks, Science, 1999, 286(5439): 509-512. CrossrefGoogle Scholar

  • [13]

    Leskovec J., Kleinberg J., Faloutsos C., Graphs over time: Densification laws, shrinking diameters and possible explanations, The 11th ACM SIGKDD international conference on Knowledge discovery in data mining (KDD). Chicago, Illinois, USA, 2005: 177-187. Google Scholar

  • [14]


  • [15]

    Viswanath B., Mislove A., Cha M., Gummadi K.P., On the evolution of user interaction in Facebook, The 2nd ACM SIGCOMM Workshop on Social Networks(WOSN), Barcelona, Spain,2009: 37-42. Google Scholar

  • [16]


  • [17]

    Awadalla N.S., Hanna M.A., Ismail M.N., Period Variation Study and Light Curve Analysis of the Eclipsing Binary GSC 02013-00288, Applied Mathematics and Nonlinear Sciences, 2016, 1(2): 321–334.CrossrefGoogle Scholar

  • [18]

    Chen Q., Chang H., Govindan R., Jamin S., Shenker S.J., Willinger W., The origin of power laws in Internet topologies revisited, The 21st Annual Joint Conference of the IEEE Computer and Communications Societies, Los Alamitos, CA, USA, 2002: 608-617.Google Scholar

  • [19]

    Watts D.J., Strogatz S.H., Collective dynamics of “small-world” networks, Nature, 1998, 393: 440-442. CrossrefGoogle Scholar

  • [20]

    Rosa M., Gandarias M.L., Multiplier method and exact solutions for a density dependent reaction-diffusion equation, Applied Mathematics and Nonlinear Sciences, 2016, 1(2):311–320. CrossrefGoogle Scholar

  • [21]

    Mislove A., Marcon M., Gummadi K.P., Druschel P., Bhattacharjee B., Measurement and analysis of online social networks, The 7th ACM SIGCOMM conference on Internet measurement conference (IMC), San Diego, California, USA, 2007: 29-42. Google Scholar

About the article

Received: 2016-08-10

Accepted: 2016-09-14

Published Online: 2017-03-16

Citation Information: Open Physics, Volume 15, Issue 1, Pages 68–73, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2017-0008.

Export Citation

© 2017 X.-G. Wang. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. BY-NC-ND 3.0

Comments (0)

Please log in or register to comment.
Log in