arxiv: v3 [cs.si] 27 Oct 2016

Transcription

1 Booean Gossiping Networs Bo Li, Junfeng Wu, Hongsheng Qi, Aexandre Proutiere, and Guodong Shi Abstract arxiv: v3 [cs.si] 27 Oct 2016 This paper proposes and investigates a Booean gossiping mode as a simpified but non-trivia probabiistic Booean networ. With positive node interactions, in view of standard theories from Marov chains, we prove that the node states asymptoticay converge to an agreement at a binary random variabe, whose distribution is characterized for arge-scae networs by mean-fied approximation. Using combinatoria anaysis, we aso successfuy count the number of communication casses of the positive Booean networ expicity in terms of the topoogy of the underying interaction graph, where remaraby minor variation in oca structures can drasticay change the number of networ communication casses. With genera Booean interaction rues, emergence of absorbing networ Booean dynamics is shown to be determined by the networ structure with necessary and sufficient conditions estabished regarding when the Booean gossiping process defines absorbing Marov chains. These resuts iustrate possibiities of reating detaied dynamica properties of Booean networs to graphica properties of the underying interactions. Keywords: Booean networs, Gossiping process, Marov chains, Communication casses 1 Introduction 1.1 Bacground A variety of random networ dynamics with nodes taing ogica vaues arises from bioogica, socia, engineering, and artificia inteigence systems [1 4]. In the 1960s, Kauffman introduced random Booean iteration rues over a networ [1] to describe proto-organisms as randomy aggregated nets of chemica reactions where the underying genes serve as a binary (on-off) device. Inspired by neuron systems, the so-caed Hopfied networs [2] provided a way of reaizing coective computation inteigence, where nodes B. Li and H. Qi are with the Key Laboratory of Mathematics Mechanization and with the Key Laboratory of Systems and Contro respectivey, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing , China. Emai: ibo@amss.ac.cn, qihongsh@amss.ac.cn. J. Wu and A. Proutiere are with ACCESS Linnaeus Centre, KTH Roya Institute of Technoogy, Stochom 10044, Sweden. Emai: junfengw@th.se, aepro@th.se. G. Shi is with the Research Schoo of Engineering, The Austraian Nationa University, Canberra 0200, Austraia. E-mai: guodong.shi@anu.edu.au. 1

2 having binary vaues behave as artificia neurons by a weighted majority voting via random or deterministic updating. Rumors spreading over a socia networ [3] and virus scattering over a computer networ [4] can be modeed as epidemic processes with binary nodes states indicating whether a peer has received a rumor, or whether a computer has been infected by a type of virus. Booean dynamica networs, consisting of a finite set of nodes and a set of deterministic or stochastic Booean interaction rues among the nodes, are natura and primary toos for the modeing of the above node dynamics with ogica vaues. The study of Booean networs received considerabe attention for aspects ranging from steady-state behaviors and input-output reations to imit cyce attractors and mode reduction, e.g., [6, 7, 10 12, 18, 21]. It has been we understood that deterministic Booean rues are essentiay inear in the state space [6, 18], whie probabiistic Booean networs are merey standard Marov chains [7, 19 21]. There however exist fundamenta chaenges in estabishing expicit and precise theoretica resuts due to computation compexity barriers as we as the ac of anaytica toos. In this paper, we propose and study a randomized Booean gossiping process, where Booean nodes pairwisemeet overanunderyinggraphinarandommannerateach timestep,andthenthetwointeracting nodes update their states by random ogica rues in a prescribed set of Booean operations. 1.2 The Mode We consider n nodes indexed in the set V = {1,...,n}. Time is sotted at t = 0,1,... Each node i hods a binary vaue from the set {0,1} at each time t, denoted x i (t). Note that, there are a tota of 16 Booean functions with two arguments mapping from {0,1} 2 to {0,1}. We index these functions in the set H := { 0,..., 15 }, where each specifies a binary Booean function in the way that a b is the vaue of the function with arguments (a, b). The underying interaction structure of the networ is modeed by an undirected graph G = (V, E), where E is the edge set with each entry being an unordered pair of two distinct nodes in V. Throughout our paper we assume that the graph G = (V,E) is connected. We use the foowing gossip mode [13] to describe the node interaction process. Each node meets other nodes at independent time instants defined by a rate-one Poisson process. This is to say, the inter-meeting times at each node foows a rate-one exponentia distribution. Without oss of generaity, we can assume that at most one node is active at any given instant. Node interactions are characterized by an n n stochastic matrix M = [m ij ], where m ij 0 and n j=1 m ij = 1 for a i,j = 1,...,n. The matrix M compies with the structure of the graph (V,E) in the sense that m ij > 0 if and ony if {i,j} E. Then independenty at each time t 0, a node i V is first drawn with probabiity 1/n, and then node i pics the pair {i,j} with probabiity m ij. Denote N i = {j : {i,j} E} as the neighbourhood of node i. 2

3 Let C be a subset of H specifying potentia node interaction rues aong the edges. Let q := C be the cardinaity of the set C. We index the eements in C by C1,..., Cq. Suppose the node pair {i,j} is seected at time t. Introducing q constants p > 0 satisfying q =1 p = 1, the evoution of the x m (t) is determined by x i (t+1) = x i (t) C x j (t), with probabiity p ; x j (t+1) = x i (t) C x j (t), with probabiity p ; (1) x m (t+1) = x m (t), m / {i,j}, where the updates of nodes i and j are independent with each other. 1.3 Reated Wor The proposed randomized Booean gossiping mode apparenty cover the cassica gossiping process [13 16] as a specia case. Compared to standard gossiping process [3] in which rumors are passed by from informed nodes to uninformed ones during pairs of interactions, the random Booean gossiping process (1) is a generaization in the sense that it taes into consideration node Booean decisions where nodes, informed or un-informed, decide whether a rumor hods true or not by random Booean rues. On the other hand, the random Booean gossiping process (1) is aso a specia case of the probabiistic Booean networ mode [10, 11], where random Booean interactions are posed pairwise. This means that, conceptuay the mode (1) under consideration can certainy be paced into the studies of genera probabiistic Booean networs, e.g., [12,19,20]. The graphica nature of the mode (1) however maes it possibe to go beyond these existing wor [12,19,20] for more direct and expicit resuts. Finay, majority Booean dynamics [24] and asynchronous broadcast gossiping [25] are reated to the mode (1) in the way that they describe Booean interactions between one node and a its neighbors at a given time instant, in contrast to the gossiping interaction rue which happens between one node and one of its seected neighbors. 1.4 Contributions The proposed random Booean gossiping mode is fuy determined by the underying graph G and the Booean interaction set C. Cassica (deterministic or probabiistic) Booean networs aso have graphica characterization [10] whereain appears if thestate of theend nodesdependon each other in the Booean updating rues. To the best of our nowedge, few resuts have been obtained regarding how the structure of the interaction graph infuences detaied networ state evoution in the study of Booean networs. First of a, we study a specia networ where the Booean interaction rues in the set C do not invove the negation, which is termed positive Booean networs. Using standard theories from Marov chains, we 3

4 show that the networ nodes asymptoticay converge to a consensus represented by a binary random variabe, whose distribution is studied for arge-scae networs in ight of mean-fied approximation methods. Moreover, by combinatoria anaysis the number of communication casses of positive Booean networs is fuy characterized with respect to the structure of the underying interaction graph G, where surprisingy oca cycic structures can drasticay change the number of communication casses of the entire networ. Next, we move to genera Booean interaction rues and study the reation between emergence of absorbing networ Booean dynamics and the networ structure. Particuary, necessary and sufficient conditions are provided for the proposed Booean gossiping process to be absorbing, where whether the underying graph G contains an odd cyce pays a crucia roe. 1.5 Paper Organization The remainder of this paper is organized foows. Section 2 investigates positive Booean dynamics in terms of steady-state distribution and communication casses. Section 3 further studies genera Booean dynamics with a focus on how the interaction graph determines absorbing Marov chains aong the random Booean dynamics. Finay Section 4 concudes the paper with a few remars. 2 Positive Booean Gossiping In this section, we consider a specia case where the Booean interaction rues in the set C do not invove the negation. Note that conventionay represents Booean AND operation, whie represents Booean OR operation. We term such types of Booean interaction as positive Booean dynamics, indicated in the foowing assumption. Assumption 1. There hods C = {, }. We denote C1 = and C2 =, and denote p = p 1 as the probabiity for operation in the dynamics (1). 2.1 State Convergence Let X t = (x 1 (t),...,x n (t)),t = 0,1,... be the random process driven by (1). It is nown that this random process ( ) X t defines a 2 n -state Marov chain M(G,C) = (S n,p), where S n = { (s 1,...,s n ) : s i {0,1},i V } is the state space, and P is the state transition matrix. An eement (s 1,...,s n ) S n is denoted in short as [s 1...s n ]. We denote P = [ ] P [s1...s n][q 1...q n] R 2 n 2 n 4

5 with its rows and coumns indexed by the eements in S n, i.e., ( ) P [s1...s n][q 1...q n] := P x(t+1) = [q 1...q n ] x(t) = [s 1...s n ]. Reca that a state in a Marov chain is caed absorbing if it is impossibe to eave this state [23]. A Marov chain is caed absorbing if it contains at east one absorbing state and it is possibe to go from any state to at east one absorbing state in a finite number of steps. In an absorbing Marov chain, the non-absorbing states are caed transient. It is not hard to find that the Marov chain M(G,C) is an absorbing chain with [0...0] and [1...1] being the two absorbing states. The state transition matrix P therefore wi have the form P = I 2 0, R Q where I 2 corresponds to the two absorbing states [0...0] and [1...1], R is a (2 n 2) 2 matrix describing transition from the 2 n 2 transient states to the two absorbing states, and Q is a (2 n 2) (2 n 2) matrix describing the transition between the transient states. We can concude the foowing resut from standard theories for absorbing Marov chains. Proposition 1 Let Assumption 1 hod. Given X 0 = X(0) S n. (i) There hods im t P ( X t = Y ) = [ (I 2 n 2 Q) 1 R ] X 0 Y for any Y S n, for which the convergence is exponentia at a rate determined by the spectrum radius of Q. (ii) There hods [ (I 2 n 2 Q) 1 R ] X 0 = 0 for a transient Y. In other words, the stationary distribution Y of the Marov chain entirey concentrates on the two absorbing states [0...0] and [1...1], with probabiities [ (I 2 n 2 Q) 1 R ] X 0 [0...0] and [ (I 2 n 2 Q) 1 R ] X 0 [1...1], respectivey. (iii) There exists a Bernoui random variabe x with mean vaue [ (I 2 n 2 Q) 1 R ] X 0 [1...1] such that im t x i(t) = x amost surey for a i V. 2.2 Continuous-Time Approximation Athough it has been cear from Proposition 1 that the imit of the node states is fuy characterized by [ (I 2 n 2 Q) 1 R ] [0...0]X 0 and [ (I 2 n 2 Q) 1 R ] [1...1]X 0, computing the exact vaue or even obtaining an approximation for the matrix (I 2 n 2 Q) 1 R is difficut for arge networs due to the exponentiay increasing dimension of the matrix. In this subsection, using mean-fied method [4, 22], we construct a continuous-time differentia equation to approximate the behavior of X(t) for arge scae networs. To this end, we assume that the x i (0) are i.i.d Bernoui random variabes. 5

6 2.2.1 Compete Graph Define n D(t) = x i (t)/n i=1 as the proportion of nodes that tae vaue 1 at time t. Assume the underying networ forms a compete graph and each node chooses the node that it interacts with from its neighbours with uniform probabiity at each step. Denote D(t) as the expected vaue of D(t), i.e., D(t) = E{D(t)}. The density D(t) evoves by the foowing rues: Let the two nodes in the seected pair {i,j} hod different vaues. When n is arge, and the graph is compete, this happens with an approximate probabiity 2D(t)(1 D(t)). The vaue D(t) wi increase by 1/n if the two seected nodes both use operations to update their vaues, an event with probabiity p 2. The vaue D(t) wi decrease by 1/n if the two seected nodes both appy operations, an event with probabiity (1 p ) 2. For a other cases, D(t) is unchanged. As a resut, we concude that E{D(t+1) D(t) D(t)} 1 n p 2 2D(t)(1 D(t)) 1 n (1 p ) 2 2D(t)(1 D(t)). (2) For a networ of n nodes, V{D(t)} = E{D 2 (t)} E 2 {D(t)} can be considered very sma for arge n. We further have D(t+1) D(t) 1 n p 2 2D(t)(1 D(t)) 1 n (1 p ) 2 2D(t)(1 D(t)). (3) Define s = t/n and D(s) = D(ns) = D(t). Then, (3) can be written as D(s+ 1 ) D(s) n 1 n p 2 2 D(s)(1 D(s)) 1 n (1 p ) 2 2 D(s)(1 D(s)) (4) We can therefore approximate (4) for arge n by the foowing differentia equation whose soution reads anayticay as d ds D(s) = p 2 2 D(s)(1 D(s)) (1 p ) 2 2 D(s)(1 D(s)), (5) D(s) = D(0) (1 D(0))e 2(1 2p )s + D(0). (6) 6

7 Here D(0) = D(0) = D 0 is the mean of the i.i.d Bernoui random variabes x i (0). Consequenty, we estabish the foowing approximate equation for D(t): From (7), the foowing hods. D(t) = D 0 (1 D 0 )e 2(1 2p )t/n +D 0. (7) Concusion. Assume G is a compete graph. For arge n, D(t) approaches zero when p < 1/2, and D(t) approaches one when p > 1/2, as time tends to infinity. To verify this concusion, we show some numerica resuts. Exampe 1. Consider a compete graph with n = 1000 nodes. Fix D 0 = 0.5, and we randomy distribute the vaues of nodes according to D 0 = 0.5. Consider T = time steps. For p = 0.49 and 0.51, et the nodes update their vaues according to (1) for 1 t T for 2000 rounds, respectivey. Then, we cacuate the average density of nodes with vaue one for every t. We compare the numerica simuation with the approximate soution given by (7). Figure 1 shows that (7) approximates the rea process (1) remaraby we Node Proportion with State Simuated Reaization p * =0.49 ODE Approximation p * =0.49 Simuated Reaization p =0.51 * ODE Approximation p * = Time t x 10 4 Figure 1: A compete graph with 1000 nodes is considered. The soid ines are the approximate soution given by (7); the dashed ines are drawn according to the simuated reaization of the agorithm (1). The continuous-time approximations match the numerica reaizations rather precisey Reguar Graph A reguar graph is a graph where nodes have equa degrees. Suppose node i is seected to initiaize a gossiping interaction at time t. Because i is uniformy seected from V, the probabiity that the seected 7

8 node i is at state 1 is D(t). If G is a reguar graph with a high degree, the distribution of the random variabe wi tend to have a simiar distribution with j N i x j (t) N i j i x j(t), n 1 which is approximatey a Bernoui random variabe with mean D(t). Therefore, D(t) evoves foowing simiar rue as compete graphs, and the differentia equation (7) wi continue to be a good approximation for high-degree reguar graphs. Exampe 2. Consider a reguar graph of degree 500 with n = 1000 nodes. Fix D 0 = 0.5, and consider T = time steps. For p = 0.49, we randomy distribute the vaues of nodes and et the nodes update their vaues according to (1) for 1 t T for 2000 times. We cacuate the average density of nodes with vaue 1 for a t. Figure 2 shows that (7) continues to provide an acceptabe approximation of the rea process (1) Node Proportion with State Simuated Reaization ODE Approximation Time t x 10 4 Figure 2: A reguar graph with 1000 nodes is considered where node degree is 500 and p = The soid ine is the approximate soution given by (7) and the dashed ine is drawn according to numerica simuation. We see that (7) continues to be a good approximation of (1). On the other hand, when the degree of a reguar graph is sma, the degree wi significanty affect the vaidity of (7), as confirmed in the foowing exampe. Exampe 3. Consider reguar graphs of degree 2,6,12,36 with n = 1000 nodes, respectivey. Fix D 0 = 0.5, p = 0.4, and consider T = time steps. Aso, we cacuate the average density of nodes with vaue 8

9 1 for a t. Figure 3 shows that the deviation between (7) and the rea process (1) increases as the node degree decreases Ratio of Nodes with Vaue Simuation: Degree 2 Simuation: Degree 6 Simuation: Degree 12 Simuation: Degree 36 ODE Approximation Time t x 10 4 Figure 3: For p = 0.4, we run numerica simuation of (1) on reguar graphs of degree 2,6,12 and 36 with 1000 nodes, respectivey. The soid ine is approximate soution given by (7) and the dashed ines are drawn according to numerica simuation. The resuts show that the deviation between (7) and the rea process (1) becomes greater as degree becomes smaer. 2.3 Communication Casses We continue to investigate the communication casses of M(G,C). Reca that a state [s 1...s n ] is said to be accessibe from state [q 1...q n ] if there is a nonnegative integer t such that P ( X t = [s 1...s n ] X0 = [q 1...q n ] ) > 0. It is termed that [s 1...s n ] communicates with state [q 1...q n ] if [s 1...s n ] and [q 1...q n ] are accessibe from each other [23]. This communication reationship forms an equivaence reation among the states in S n. The equivaence casses of this reation are caed communication casses of the chain M(G,C). The number of communication casses of M(G,C) is denoted as χ (G,C). The foowing theorem provides a fu characterization to χ (G,C). Theorem 1 Let Assumptions 1 hod. Then (i) χ (G,C) = 2n, if G is a ine graph; (ii) χ (G,C) = m + 3, if G is a cyce graph with n = 2m; χ (G,C) = m + 2, if G is a cyce graph with n = 2m+1; 9

10 (iii) χ (G,C) = 5, if G is neither a ine nor a cyce, and contains no odd cyce; (iv) χ (G,C) = 3, if G is not a cyce graph but contains an odd cyce. Estabished by constructive proofs that can overcome the fundamenta computationa obstace in anayzing arge-scae Booean networs, Theorem 1 iustrates how oca structures can drasticay change the number of communication casses as a goba property of networs. The detaied proof of Theorem 1 has been put in Appendix A. 3 Genera Booean Dynamics In this section, we discuss the evoution of (1) under genera Booean interaction set C 2 H, where 2 H denotes the set containing a subsets of H. To this end, we need to identify precisey the eements in C. A fu ist of the 16 Booean functions mapping {0,1} 2 to {0,1} is provided in Appendix B, indexed from 0 to Main Resuts { } Introduce B := C 2 H : 10 C { 2, 3, 10, 11 } or { 2, 11 } C { 2, 3, 10, 11 }. The foowing resut proves that, when interaction rue set C does not beong to B, absorbing property of the induced Marov chain becomes a property determined fuy by the connectivity of G. Theorem 2 Suppose C / B. Then, for any two connected graphs G 1 and G 2, M(G 1,C) is an absorbing Marov chain if and ony if M(G 2,C) is an absorbing Marov chain. Next, we estabish a necessary and sufficient condition regarding M(G 1,C) being an absorbing chain when interaction rue set C does not beong to B. Theorem 3 Suppose C / B. Then M(G 1,C) is an absorbing Marov chain if and ony if one of the foowing conditions hods (i) C { 2, 3 }; (ii) C { 3, 11 }; (iii) C { 0, 1, 2, 3, 4, 5, 6, 7 } and C { 0, 1, 4, 5 } ; (iv) C { 1, 3, 5, 7, 9, 11, 13, 15 } and C { 5, 7, 13, 15 }. 10

11 For interaction rues C B, the foowing theorem further gives a tight condition regarding M(G 1,C) being an absorbing Marov chain. Theorem 4 Suppose C B. Then M(G,C) is an absorbing Marov chain if and ony if G does not contain an odd cyce. Note that, Theorem 2 can actuay be inferred from Theorem 3. Theorem 3 and Theorem 4 together present a comprehensive understanding of the absorbing property of the networ Booean evoution. 3.2 Key Lemma To simpify the discussion, we introduce some new notations. For any S = [s 1...s n ] S n, we denote [S0] as [s 1...s n 0] S n+1 and [S1] as [s 1...s n 1] S n+1. For any a {0,1}, denote a = 1 a. We categorize the states into the foowing five casses: C 1 (G) = {[s 1...s n ] : s i = 0, 1 i n}, C 2 (G) = {[s 1...s n ] : s i = 1, 1 i n}, C 3 (G) = {[s 1...s n ] : s i s j for any edge {i,j} of G}, C 4 (G) = {[s 1...s n ]: i, j,, s.t. {i,j} is an edge of G and 0 = s i = s j s }, and C 5 (G) = {[s 1...s n ]: i, j,, s.t. {i,j} is an edge of G and 1 = s i = s j s }. We may simpy write C i instead of C i (G) whenever this simpification causes no confusion. In this subsection, we estabish a ey technica emma regarding whether a state in the C i can be an absorbing state in terms of the seection of C. Lemma 1 (i) The state in C 1 is an absorbing state if and ony if C { 0, 1, 2, 3, 4, 5, 6, 7 }. (ii) The state in C 2 is an absorbing state if and ony if C { 1, 3, 5, 7, 9, 11, 13, 15 }. (iii) Any state in C 3 is an absorbing state if and ony if C { 2, 3, 10, 11 }. (iv) Any state in C 4 is an absorbing state if and ony if C { 2, 3 }. (v) Any state in C 5 is an absorbing state if and ony if C { 3, 11 }. Proof. (i) We denote the ony state in C 1 as S [1]. Note that S [1] is the state that any two nodes associated with a common edge hod the same vaue 0. According to the agorithm (1), S [1] is an absorbing state, if and ony if for any i C there hods 0 i 0 = 0. Thus, S [1] is an absorbing state, if and ony if C { 0, 1, 2, 3, 4, 5, 6, 7 }. (ii) The proof is simiar to that in (i), whose detais are omitted. (iii) For any state S C 3, any two nodes associated with a common edge share different vaues. According to (1), S is an absorbing state if and ony if for any i C, 0 i 1 = 0 and 1 i 0 = 1. That is, S is an absorbing state, if and ony if C { 2, 3, 10, 11 }. 11

12 (iv) By (1), S C 4 is an absorbing state if and ony if for any i C, there hod 0 i 0 = 0, 0 i 1 = 0 and 1 i 0 = 1. In other words, S is an absorbing state, if and ony if C { 2, 3 }. (v) The proof is simiar to that of (iv), which is, again omitted. 3.3 Proof of Theorem 3 This subsection focuses on the proof of Theorem 3. (Necessity.) Assume M(G, C) is an absorbing Marov chain. If any state in C 4 is absorbing, according to Lemma 1(iv), we have C { 2, 3 }. If any state in C 5 is absorbing, according to Lemma 1(v), we have C { 3, 11 }. If no state in C 4 and C 5 is absorbing, then no state in C 3 is absorbing, by Lemma 1(iii)(iv)(v) and the condition C / B. Thus, if no state in C 4 and C 5 is absorbing, there must be states in C 1 or C 2 absorbing. We discuss the two cases, respectivey. If the state S [1] in C 1 is absorbing and there exists at east one state that can transit to S [1] in finite steps with positive probabiity, it must be that there exists i C such that 1 i 0 = 0 according to (1). That is C { 0, 1, 4, 5 }. By Lemma 1, C { 0, 1, 2, 3, 4, 5, 6, 7 }. If the state S [2] in C 2 is absorbing and there exists at east one state that can transit to S [2] in finite steps with positive probabiity, it must be that there exists i C such that 0 i 1 = 1 according to (1). That is C { 5, 7, 13, 15 }. Again, by Lemma 1, C { 1, 3, 5, 7, 9, 11, 13, 15 }. We have now concuded the necessity statement. (Sufficiency.) We divide the proof into a few cases: Let C { 2, 3 }. If 2 C, there is a positive probabiity that 2 is chosen. Because = 0, according to (1), any state in C 2 or C 5 can transit to state in C 4 in finite steps with positive probabiity. Moreover, states in C 1, C 3, C 4 are absorbing states. Thus, M(G,C) is an absorbing Marov chain. On the other hand, if 2 / C, C = { 3 }. In this case a the states are absorbing. Of course, M(G,C) is an absorbing Marov chain. Let C { 3, 11 }. The proof foows from the same argument as the proof under C { 2, 3 }. Let C { 0, 1, 2, 3, 4, 5, 6, 7 } and C { 0, 1, 4, 5 }. The state S [1] in C 1 is absorbing by Lemma 1(i). Note that there is a positive probabiity that i { 0, 1, 4, 5 } is chosen. Because 1 i 0 = 0, any states other than S [2] C 2 can transit to state S [1] in finite steps with positive probabiity. We discuss the absorbing status of S [2]. 12

13 1. If S [2] is absorbing, nothing needs to prove that M(G,C) is absorbing. 2. If S [2] is not absorbing, it can transit to other state with positive probabiity, which in turn transits to S [1]. Thus, M(G,C) is absorbing. Assume C { 1, 3, 5, 7, 9, 11, 13, 15 } and C { 5, 7, 13, 15 }. The proof is simiar as the case above. The proof of Theorem 3 is now competed. 3.4 Proof of Theorem 4 In this subsection, we prove Theorem 4. If G contains an odd cyce, C 3 is empty. By Lemma 1, no state in C 1 C 2 C 4 C 5 is absorbing. As S n = C 1 C 2 C 3 C 4 C 5, no state is absorbing. Thus, M(G,C) is not absorbing. On the other hand, if G does not contain an odd cyce, there are two eements in C 3, and we proceed to prove by induction on the number n of nodes that M(G,C) is an absorbing Marov chain. For n = 1, it is obvious. AssumeM(G,C) is absorbingfor n =. Theremust beaspanningtree, denoted G T1, of G. We further find a subtree G T2 of G T1 with G T2 containing nodes of G T1. Without oss of generaity, et G T2 contain nodes 1,..., of G. By our induction assumption, M(G T2,C) is absorbing. Now any state in S +1 can be represented as [Su], where S S and u {0,1}. As M(G T2,C) is absorbing, there is a positive probabiity that in finite steps S transits to a state S in C 3 (G T2 ). Because G T2 is a subgraph of G T1, [Su] can transit to [S u] in finite steps in M(G T1,C). There wi be two cases. If [S u] C 3 (G T1 ), for G contains no odd cyce, [S u] C 3 (G). The proof is done. If [S u] / C 3 (G T1 ), there must be some node j associated with node +1 over graph G T1. Because C B, there is a positive probabiity that 10 or 11 is chosen. Note that = 1, and = 1. Thus, by (1), [S u] transits to [S u] with positive probabiity in M(G T1,C). Moreover, [S u] C 3 (G T1 ). For G contains no odd cyce, [S u] C 3 (G) eads to the desired resut. The proof of Theorem 4 is competed. 4 Concusions We proposed and investigated a Booean gossiping mode, which may be usefu in describing socia opinion evoution as we as serves as a simpified probabiistic Booean networ. With positive node interactions, using standard theories from Marov chains, it was shown that the node states asymptoticay converge to a consensus represented by a binary random variabe, whose distribution is studied for arge-scae compete networs in ight of mean-fied approximation methods. By combinatoria anaysis the number of 13

14 communication casses of the positive Booean networ was counted against the topoogy of the underying interaction graph. With genera Booean interaction rues, the emergence of absorbing networ Booean dynamics was expicity characterized by the networ structure. It turned out that oca structures in terms of existence of cyces can drasticay change fundamenta properties of the Booean networ. In future, it wi be interesting to oo into the possibiity of extending the graphica anaysis estabished in the current wor to more genera Booean networs such as majority voting. Appendix A. Proof of Theorem 1 For each n, we use Mod n (i) to denote the unique integer j satisfying 1 j n and i j mod (n). Reca that for any a {0,1}, we denote a = 1 a. We prove the statements of Theorem 1 in a few steps starting with a few fundamenta graphs. A.1 Line graph In this subsection we prove Theorem 1.(i) stating that χ (G,C) = 2n when G is a ine graph. Without oss of generaity we assume the edges of G are {i,i+1} for i = 1,...,n 1. The proof is outined as foows. We first introduce the notion of L-reduced state for each state in S n. Then, we prove that any two states communicate with each other if and ony if their L-reduced states are identica. Finay, we count the number of L-reduced states in the state space and therefore obtain the number of communication casses. Definition 1 (L-reduced states) Let [s 1...s n ] S n. There exists a unique partition of s 1,...,s n into s 1 = = s i1 = r 1, i 1 1; s i1 +1 = = s i2 = r 2, i 2 > i 1 ;... s id 2 +1 = = s id 1 = r d 1, i d 1 > i d 2 ; s id 1 +1 = = s n = r d such that r i r i+1 for a i = 1,...,d 1. Then [r 1...r d ] := L([s 1...s n ]) is termed the L-reduced states of [s 1...s n ]. Note that the vaues of any two consecutive eements in an L-reduced state are different. The foowing two emmas hod. Lemma 2 Suppose G is a ine graph. Then L([s 1...s n ]) is a subsequence of L([q 1...q n ]) if [s 1...s n ] is accessibe from [q 1...q n ]. More precisey, denoting L([s 1...s n ]) = [r 1...r d ], L([q 1...q n ]) = [h 1...h d ] 14

15 we have d d and there exists 1 τ 1 < τ 2 < < τ d d such that r i = h τi for a i = 1,...,d. Proof. By the definition of accessibiity, there is a nonnegative integer t such that P ( X t = [s 1...s n ] X 0 = [q 1...q n ] ) > 0. First we assume t = 1. According to the structure of (1), either [s 1...s n ] = [q 1...q n ], or there is u {1,...,n} such that s u q u and s i = q i for a i u. The desired concusion obviousy hods if [s 1...s n ] = [q 1...q n ]. For the atter case, there is q v with v = u + 1 or v = u 1 such that q u q v. Consequenty, the two states [s 1...s n ] and [q 1...q n ] differ with each other ony at s u and q u and satisfy s u q u, s u = s v, q u q v. Then it is easy to verify that L([s 1...s n ]) is a subsequenceof L([q 1...q n ]) from the definition of L-reduced states. Now we proceed to et t = 2. There wi be a state [w 1...w n ] such that [s 1...s n ] is one step accessibe from [w 1...w n ], and [w 1...w n ] is one step accessibe from [q 1...q n ]. Utiizing the above understanding for the case with t = 1 we now L([s 1...s n ]) is a subsequence of L([w 1...w n ]) and L([w 1...w n ]) is a subsequence of L([q 1...q n ]), which in turn impy L([s 1...s n ]) is a subsequence of L([q 1...q n ]). Therefore the desired concusion hods for t = 2. Apparenty the argument can be recursivey carried out and the resut hods for arbitrary integer t. We have now competed the proof of the emma. Lemma 3 Let G be a ine graph and consider S = [s 1...s n ],Q = [q 1...q n ] S n. S and Q communicate with each other if and ony if they have identica L-reduced states. Proof. The necessity part of this emma foows directy from Lemma 2. In the foowing we focus ony on the sufficiency part. Let the identica L-reduced state of [s 1...s n ] and [q 1...q n ] be [r 1...r ]. We carry out an induction argument on for any n. Let = 1. Then [0...0] n and [1...1] n are the two possibe states for [s 1...s n ] and [q 1...q n ]. The desired concusion hods straightforwardy. Now assume: Induction Hypothesis: The statement of the emma hods true for a and a n. We proceed to prove the statement for = +1 and n. Denote i 1 = max{h : r 1 = s i, 1 i h} and j 1 = max{h : r 1 = q i, 1 i h}. By symmetry we may assume i 1 j 1 and we use the foowing two observations: a) The state [q 1...q n ] communicates with the state [q 1...q i1 q i q j1 q j q n ] by the definition of j 1. 15

16 b) Thetwostates [q i q j1 q j q n ]and[s i1 +1s i s n ]havethesamel-reducedstate[r 2...r ]. Therefore by our induction hypothesis, [q i q j1 q j q n ] and [s i s n ] communicate with each other, which in turn yieds that [s 1...s n ] communicates with [q 1...q i1 q i q j1 q j q n ]. Combining a) and b) we immediatey now that [s 1...s n ] communicates with [q 1...q n ]. By the principe of mathematica induction we have competed the proof of the emma. We are now ready to count the number of communication casses for the ine graph, which equas to the number of L-reduced states according to Lemma 3. For each m = 1,...,n, there are two different L-reduced states with ength m, i.e., [r 1...r m ] with r 1 = 0 or r 1 = 1. Consequenty, there are a tota of 2n different L-reduced states. This concudes the proof for Theorem 1.(i). A.2 Cyce graph In this subsection, we prove the case with G being a cyce graph. Without oss of generaity, et G be the cyce graph with edges { i,mod n (i+1) }, i = 1,...,n. We introduce some usefu notations that wi be used subsequenty. For any, we use σ to denote the permutation on set {1,...,} with σ (i) = Mod (i +1) for i = 1,...,. We further define P σ as a mapping over S by P σ ([s 1...s ]) = [s σ (1)...s σ ()] for a [s 1...s ] S. Intuitivey, if we pace these nodes uniformy on a cyce and denote the vaue of each nodeon them, then the resutof P σ on astate is obtained by rotating a the vaues countercocwise. We aso define a mapping h [1, 2 ] over S n by that for any [t 1...t n ], h [1, 2 ]( [t1...t n ] ) = [r 1...r n ] satisfies r i = t i for i 2 and r i = t 1 for i = 2. Definition 2 (C-reduced states) Let [s 1...s n ] S n with [r 1...r d ] = L([s 1...s n ]) being its L-reduced states. The C-reduced states of [s 1...s n ] S n, denoted C([s 1...s n ]), is defined as foows: If d = 1, then C([s 1...s n ]) := [r 1 ]. If d > 1 and r d r 1, then C([s 1...s n ]) := [r 1...r d ]. If d > 1 and r d = r 1, then C([s 1...s n ]) := [r 1...r d 1 ]. Let C(S) be the number of digits in C(S) for S S n. According to the definition, the vaues of any two consecutive eements of C-reduced states are different. Moreover, if there are at east two entries of C-reduced states, the first entry is different from the ast one. The foowing emma can be estabished using a simiar anaysis as we used in Lemma 2. 16

17 Lemma 4 Suppose G is a cyce graph, (i) C(S) is either 1 or an even integer; (ii) If d is one or an even integer, then there is S S n with C(S) = d. (iii) If S is accessibe from T, then C(S) C(T). Lemma 5 Consider S,T S n. If 1 < C(S) = C(T) < n, then S and T communicate with each other. Proof. Denote S = [s 1...s n ] and T = [t 1...t n ]. We prove this emma in a few steps. Step 1. We first prove that S communicates with Pσ n (S) for any integer if C(S) < n. Note that if C(S) = 1, S must be [0...0] n or [1...1] n. The caim hods straightforwardy. Now we assume C(S) > 1. Since C(S) < n, the set I := {i : s i = s Modn(i+1), 1 i n} is nonempty. Moreover, because C(S) > 1, we can find j I such that s Modn(j+1) s Modn(j+2). By the structure of (1), the state h [Modn(j+2),Mod n(j+1)](s) is accessibe from S. By the definition of j, there hods h [j,modn(j+1)]h [Modn(j+2),Mod n(j+1)](s) = S. Thatistosay,S isaccessibefrom h [Modn(j+2),Mod n(j+1)](s).therefore,s communicateswith h [Modn(j+2),Mod n(j+1)](s Appyingthisargumentrecursivey,weobtainthatS communicateswith h [Modn(j+n),Mod n(j+n 1)]...h [Modn(j+2),Mod n which is equa to P σn (S). It is then convenient to concude that S communicates with P σ n (S) for any integer. Step 2. In this step, we prove that if S = [s 1...s n ] and T = [t 1...t n ] have identica C-reduced states, then S and T communicate with each other. Let [c 1...c d ] = C(S) = C(T). If d = 1 or n, it is easy to see S = T. Now assume 1 < d < n. Because d > 1, the sets {i : s i s 0 } and {i : t i t 0 } are not empty. Denote j 1 = max{i : s i s 0 }, and j 2 = max{i : t i t 0 }. Without oss of generaity we assume j 1 > j 2. By (1), T communicates with h [j2,j 2 +1](T). Further we now that T communicates with T = [t 1...t n] := h [j1 1,j 1 ]...h [j2 +1,j 2 +2]h [j2,j 2 +1](T). Moreover, we can concude that j 1 = max{i t i t 0 }, and T and T have the same C-reduced state. So T and S have the same C-reduced state. By the definition of C-reduced state and the fact that j 1 = max{i s i s 0 } = max{i t i t 0 }, we now that the L-reduced state of S is equa to the L-reduced state of T. Define a new ine graph G, whose nodes are the nodes of G with edges being {i,i + 1} for i = 1,...n 1. According to Lemma 3, S aso communicates with T in M( G,C). Therefore, S communicates with T in M(G,C), because G is a subgraph of G. Thus, S and T communicate with each other. Step 3. This step wi compete the proof. 17

18 Let d = C(S). If C(S) = C(T), we have nown that S and T communicate with each other. We ony need to consider the case C(S) C(T). Because C(S) = C(T), there must hod that C(T) = P σd (C(S)). For d > 1, the set {i : s i s 0 } is nonempty. Define j = min{i : s i s 0 }. According to Step 1, S communicates with P j 1 σ n (S). By the definition of C-reduced states, we now that the C-reduced state of P j 1 σ n (S) is P σd (C(S)), i.e., C(T). Therefore, P j 1 σ n (S) communicates with T, impyingthat S communicates with T. Now, we are ready to count the number of communication casses. According to Lemma 4, the digit number d of the C-reduced states of a the states in the same communication cass are identica. Moreover, d can be 1 or even numbers. If n = 2m, there are three cases: (i) For d = 1, there are two communication casses {[0...0]} and {[1...1]}. (ii) For each d = 2,4,...,2m 2, according to Lemma 4 and Lemma 5, there is a unique communication cass whose eements have C-reduced with d digits. (iii) For d = 2m, the two states S 0 := [s 1...s 2m ] and T 0 := [ s 1... s 2m ] with s 2i 1 = 1 and s 2i = 0 for i = 1,...,m, are the ony states whose C-reduced states are of ength 2m. Moreover, either S 0 or T 0 cannot be accessibe from any other state. That is to say, they form two communication casses. As a resut, there are a tota of m+3 communication casses. We have competed the proof for the case n = 2m. The case with n = 2m + 1 can be simiary anayzed, whose detaied proof is omitted. This concudes the proof of Theorem 1(ii). A.3 Star graph In this subsection, we prove that χ (G,C) = 5 if G is a star graph with n( 4) nodes. Note that a connected graph is caed a star graph if there is a node such that a the edges of the graph contain this node. This particuar node is caed the center node of the graph. The foowing proposition characterizes the communication casses for M(G, C) over a star graph G. Proposition 2 Let G be a star graph with n( 4) nodes. Then χ (G,C) = 5. Moreover, etting node 1 be the center node, the five casses are F 1 n = {[s 1...s n ] : s i = 0, 1 i n}, F 2 n = {[s 1...s n ] : s i = 1, 1 i n}, F 3 n = {[s 1...s n ] : s 1 = 0, s i = 1, 2 i n}, F 4 n = {[s 1...s n ] : s 1 = 1, s i = 0, 2 i n}, F 5 n = {[s 1...s n ] : i, j, 2 i, j n, s i = 0, s j = 1}. 18

19 Proof. Denote S [1] n,...,s [4] n as the singeton state in F 1 n,...,f 4 n, respectivey. Moreover, any other state cannotbeaccessibefroms [1] n ors [2] n,whies [3] n ors [4] n cannot be accessibe from any other state. Thus, they do form communication casses, respectivey. We ony need to prove a the eements in F 5 n communicate with each other. We prove this by induction. First, et n = 4. Thereare 12 eements of F 5 4, isted as [0100], [0010], [0001], [0011], [0101], [0110], [1100], [1010], [1001], [1011], [1101], [1110]. It is easy to verify that they are in the same communication cass. Assume that for n = 4, a the eements in F 5 n communicate with each other. Now we prove the case for n = + 1. Let G be a star graph with + 1 nodes with node 1 being its center node. Let G be the subgraph of G with nodes 1,..., and a edges containing them in G. In fact, G is a star graph with nodes. By our induction assumption, a eements in F 5 communicate with each other in M(G,C). Because G is a subgraph of G, a eements in A := {[S0] S +1 : S F 5 } communicate with each other, and a eements in B := {[S1] S +1 : S F 5 } communicate with each other. Note that F 5 +1 = A B {[S[1] 1],[S[2] 0],[S[3] 0],[S[4] 1]}. Introduce U [a] = [ ], U [b] = [ ], U[c] = [ ] and U[d] = [ ]. They are eements of F 5. It is easy to verify that [U[a] 0] A is accessibe from [U[a] from [U [b] 0] A. Therefore, a eements in A B communicate with each other. [S [3] It is straightforward to verify that [S [2] 0] and [S[2] 0] communicates with [S[3] 0] is accessibe from [U[b] 0] A. Thus, a eements in A {[S[2] with each other. Moreover, [S [1] [S [1] 1] communicates with [S[4] 1] B. Moreover, [U[b] 1] B is accessibe 1], [U[d] 1] is accessibe from [U[a] 1] B. Therefore, a eements in B {[S[1] each other. Summarizing a these reations we now a eements in F 5 +1 = A B {[S[1] 1],[S[2] 0],[S[3] 0],[S[4] 1]} 0]. Aso, [U[c] 0] A is accessibe from 0],[S[3] 0]} communicate 1] B is accessibe from [S[4] 1], and 1],[S[4] 1]} communicate with communicate with each other. This competes the proof of this proposition. A.4 Tree The foowing resut presents a characterization of the number of communication casses for tree graph that is not a ine. Proposition 3 Let G be a tree, having at east one node with degree greater than 2, i.e., G is not a ine graph. Then χ (G,C) = 5. The five communication casses can be described as foows: J 1 n = { [s 1 s 2...s n 1 s n ] s i = 0, 1 i n }, J 2 n = { [s 1 s 2...s n 1 s n ] s i = 1, 1 i n }, J 3 n = { [s 1 s 2...s n 1 s n ] s 1 = 0, s i s j for any edge {i,j} of G }, 19

20 J 4 n = { [s 1 s 2...s n 1 s n ] s 1 = 1, s i s j for any edge {i,j} of G }, and J 5 n = { [s 1 s 2...s n 1 s n ] i, j,, s.t. {i,j} is an edge of G and s i = s j s }. Proof. It is straightforward to verify that any of J 1 n, J 2 n, J 3 n, J 4 n contains a unique eement, and forms a communication cass. We now prove J 5 n is a communication cass using an induction argument. For n = 4, G is a star graph which is proved in Proposition 2. Now assume that this proposition hods for n = 4. For any tree G with +1 nodes that is not a ine graph, there is a subgraph G with nodes which is sti a tree. Without oss of generaity, we denote the node not in G as node v = +1 V. We use v 0 to denote the node with the highest degree in G (If there are more than one such nodes, we just choose one of them arbitrariy). There is a path (v 0, v 1,..., v h, v ) connecting node v 0 and node v in G, where h 0 is an integer. By the induction assumption, the communication casses of M(G,C) are J 1,...,J5 with each J defined by repacing n with in J n. Denote A = {[S0] S +1 : S J 5 } and B = {[S1] S +1 : S J 5 }. Note that J 5 +1 = A B {[S[1] 1],[S [2] 0],[S [3] 0],[S [4] 1]}. Because G is a subgraph of G, a eements in A communicate with each other, and a eements in B communicate with each other. Note that if G is a star graph with the v h being the center node, G wi be a star graph. This fas to the case discussed in Proposition 2. We assume G is not a star graph for the remainder of the proof. Introduce U [a] = [ ], U [b] = [ ]. v h v h We have U [a], U [b] J 5. Thus,[U[a] 1] Bisaccessibe from[u [a] 0] A. Moreover, [U [b] 0] Ais accessibe from [U [b] 1] B. Therefore, a eements in A B communicate with each other. We further denote S [4] = [γ 1...γ ] and S [3] = [β 1...β ], and then U [c] := [γ 1...γ vh 1γ vh γ vh +1...γ ], U [d] := [β 1...β vh 1β vh β vh +1...β ]. It is straightforward to verify that [S [1] 1] communicates with [U [a] 1] B, [S [2] 0] communicates with [U [b] 0] A, [S [3] 0] communicates with [U [d] 0] A, and [S [4] 1] communicates with [U [c] 1] B. Thus, any two eements in J 5 +1 communicate with each other. A.5 Competion of Proof of Theorem 1 The statements (i) and (ii) in Theorem 1 have been proved for the cases of ine and cyce graphs. We are now in a pace to prove (iii) and (iv) based on our resuts for tree graphs. According to Proposition 3, for tree graphs without being a ine graph, there are five communication casses C 1 n, C 2 n, C 3 n, C 4 n and 20

21 C 5 n. Since any connected graph contains a spanning tree, the communication casses of M(G,C) for any connected graph G that is not a ine or cyce, can ony be unions of the C j n, j = 1,..., = = = = = = = = 0 0 : 1 : 2 : 3 : = = = = = = = = = = = = = = = = 1 4 : 5 : 6 : 7 : = = = = = = = = = = = = = = = =0 8 : 9 : 10 : 11 : = = = = = = = = = = = = = = = = 1 12 : 13 : 14 : 15 : = = = = = = = = 1 Figure 4: The 16 Binary operators mapping from {0,1} 2 to {0,1}. Proofof Theorem1(iii): SupposeG isneither ainegraphnoracyce graphandit contains nooddcyce. There is a spanning tree of G, denoted G T. For M(G T,C), C 1 n, C2 n, C3 n, C4 n and C5 n are communication casses. C 1 n and C 2 n are absorbing states in G. Because there is no odd cyce, C 3 n and C 4 n are the states that any pair of nodes associated with a common edge G share different vaues. That is to say, C 3 n and C 4 n cannot be accessibe from any other states in M(G,C). Thus, C1 n, C2 n, C3 n, C4 n and C5 n are sti communication casses in M(G,C), i.e. χ (G,C) = 5. Proof of Theorem 1(iv): Now suppose G contains an odd cyce. Again, there is a spanning tree G T of G. C 1 n, C2 n, C3 n, C4 n and C5 n are communication casses in M(G T,C). Aso, C 1 n and C2 n are absorbing states in M(G,C). For states in C 3 n and C 4 n, there is an edge e beonging to the odd cyce such that the pair nodes of this edge tae different vaues. Now, by choosing another spanning tree G T containing the edge e, we can prove that eements in C 3 n, C4 n and C5 n communicate with each other in M(G,C). In turn, χ (G,C) = 3. 21

22 Appendix B. Booean Functions There are a tota of 16 Booean functions mapping from {0,1} 2 to {0,1}. We index them in Figure 4. References [1] S. A. Kauffman, Metaboic stabiity and epigenesis in randomy constructed genetic nets, Journa of Theoretica Bioogy, 22: , [2] J. J. Hopfied, Neura networs and physica systems with emergent coective computationa abiities, Proceedings of the Nationa Academy of Sciences, vo. 79, no. 8, pp , Apri [3] R. Karp, C. Schindehauer, S. Shener, and B. Vöcing, Randomized rumor spreading, Proc. Symp. Foundations of Computer Science, , [4] P. van Mieghem, J. Omic, and R. Kooij, Virus spread in networs, IEEE/ACM Transactions on Networing, vo. 17, no. 1, pp 1-14, [5] S. A. Kauffman. The origins of order: Sef-organization and seection in evoution. New Yor: Oxford University Press [6] D. Cheng and H. Qi, Controabiity and observabiity of Booean contro networs, Automatica, 45: , [7] L. Tournier and M. Chaves, Interconnection of asynchronous Booean networs, asymptotic and transient dynamics, Automatica, 49(4), pp , [8] D. G. Green, T. G. Leishman, and S. Sadedin, The emergence of socia consensus in Booean networs, Proceedings of the IEEE Symposium on Artificia Life, pp , [9] H. Kitano. Foundations of Systems Bioogy. Cambrige: MIT Press, [10] I. Shmuevich, E. R. Dougherty, S. Kim et a., Probabiistic Booean networs: a rue-based uncertainty mode for gene reguatory networs, Bioinformatics, 2: , [11] I. Shmuevich, R. F. Hashimoto, E. R. Dougherty, et a., Steaty-state anaysis of genetic reguatory networs modeed by probabiistic Booean networs, Comp Funct Genomics, 4: , [12] M. Brun, E. R. Dougherty, and I. Shmuevich, Steady-state probabiities for attractors in probabiistic Booean networs, Signa Process, 85: , [13] S. Boyd, A. Ghosh, B. Prabhaar and D. Shah, Randomized gossip agorithms, IEEE Trans. Information Theory, vo. 52, no. 6, pp ,

23 [14] D. Acemogu, A. Ozdagar and A. ParandehGheibi, Spread of (Mis)information in socia networs, Games and Economic Behavior, vo. 70, no. 2, pp , 2010 [15] B. Doerr, M. Fouz, and T. Friedrich, Why rumors spread so quicy in socia networs? Communications of ACM, 55(6): [16] D. Shah, Gossip Agorithms, Foundations and Trends in Networing, vo. 3, no. 1, pp , [17] L. Tournier and M. Chaves, Uncovering operationa interactions in genetic networs using asynchronous Booean dynamics, Journa of Theoretica Bioogy, 260(2), , [18] D. Cheng and H. Qi, A inear representation of dynamics of Booean networs, IEEE Transactions on Automatic Contro, 55: , [19] R. Pa, Context-sensitive probabiistic Booean networs: steady-state properties, reduction, and steady-state approximation, IEEE Transactions on Signa Processing, 58(2): , [20] R. Pa and S. Bhattacharya, Characterizing the effect of coarse-scae PBN modeing on dynamics and intervention performance of genetic reguatory networs represented by stochastic master equation modes, IEEE Transactions on Signa Processing, 58(6): , [21] M. Chaves and A. Carta, Attractor computation using interconnected Booean networs: testing growth rate modes in E. Coi, Theoretica Computer Science, [Onine] Ju [22] J. O. Kephart and S. R. White, Directed-graph epidemioogica modes of computer viruses, Proc. IEEE Comput. Soc. Symp. Research in Security and Privacy, pp , pp., 20-22, May [23] C. M. Grinstead and J. L. Sne. Introduction to Probabiity: Second Revised Edition. American Mathematica Society, [24] Y. Kanoria and A. Montanari, Majority dynamics on trees and the dynamic cavity method, The Annas of Appied Probabiity, vo. 21, no. 5, , [25] H. Amini, M. Draief, and M. Learge, Fooding in weighted sparse random graphs, SIAM Journa on Discrete Mathematics, 27(1): 1-26,