Irreversible -threshold processes: Graph-theoretical threshold models of the spread of disease and of opinion
Keywords
1. Introduction
Mathematical modeling of the spread of infectious disease has a long history, going back to the work of Bernoulli on smallpox in 1760. Much of the modern disease modeling literature studies the spread of disease using methods of dynamical systems. However, as diseases are spread through social networks, we can gain insight by using simple graph-theoretical models. The models we will describe have a variety of other applications, and in particular fall into an extensive literature on the spread of opinions through networks. Sample early work in this literature can be found in the work of Poljak and Sûra [23], DeGroot [3] and French [8].
We will consider models in which a network is represented by a graph , where is a set of vertices representing people and means that the people are connected in the network. We will assume that each vertex in a network is in a state at time . In our models, we will consider these states as changing only at discrete times, and we will assume that each state is either 0 or 1. We will consider different models for how these states change over time, depending on the states of neighboring vertices in the graph . In the disease application, we think of 1 as corresponding to a vertex having a given disease and 0 as corresponding to it not having the disease. In the case of opinion, we consider either agreement (1) or disagreement (0) with a given position.
Our major emphasis will be on a model we call an irreversible -threshold process, where a vertex never leaves state 1 (is never cured once infected) and enters state 1 from state 0 if at least of its neighbors are in state 0. A variation on the model will allow “vaccination” of a vertex, preventing it from ever entering state 1. We will be especially concerned with identifying sets of vertices we will call irreversible -conversion sets, sets of vertices with the property that if they are all in state 1 at the beginning, then eventually all vertices end up in state 1. Such a set of vertices represents a set of individuals a bioterrorist could infect so as to be sure to infect everyone eventually, or a set of individuals we could recruit to have a given opinion so as to eventually “infect” everyone with this opinion. Defense against infection would then amount to finding good “vaccination” or “quarantine” strategies or to designing the network topology so as to prevent the network from having small irreversible -conversion sets. In a -threshold process, we fix a number and assume that a vertex changes its state at time if at least of its neighbors are in the opposite state at time . The number is a threshold for state change. In an irreversible -threshold process, a vertex never changes from state 1 to state 0, but changes from state 0 to 1 if at least of its neighbors are in state 0 at the previous time. In terms of diseases, the idea is that an infection spreads to a vertex if sufficiently many of its neighbors are infected. In some models, will suffice. It does not seem likely that ordinary -threshold processes make sense for disease spread, since it is hard to understand the motivation for switching from 1 to 0. In this paper, we will concentrate on irreversible -threshold processes. Extensive analogous results on -threshold processes are contained in [4], as are results about the monotonic -threshold processes where the set of vertices in state 1 at time is a subset of the set at time . Mention should be made of majority processes, where at each time step, a vertex changes its state if and only if at least half of its neighbors have the opposite state. There are variations of the majority process which include the state of in the determination of majority, along with various rules for breaking ties (always choose 1 or always 0, always keep the current state, always switch states). These processes do not seem to be relevant to the spread of disease, though they do seem interesting for opinion change. We mention them here because they motivate our interest in a key notion of “dynamo” or “conversion” set that we shall define below. Majority processes (sometimes also called repetitive polling processes) have had applications in problems in distributed computing (see [5], [17]). A concrete example of an application in distributed computing, in particular, maintaining data consistency, was suggested by Peleg in [22]. Assume that all of the processors in a distributed network store an identical bit of data and some of these bits get changed for whatever reason (rounding error, for example). If a processor discovers that the values (states) of the bits of some of its neighbors are different, a majority process could be used to update the values of all of the processors until every processor has the same value for its bit. In a general context, though, there is no guarantee that every processor will return to the same value, much less the correct one. Peleg also referred to other applications of majority processes to distributed computing in system-level diagnosis and resource allocation. The problem of maintaining data consistency raises an interesting question: for what graphs does the majority process always lead to the situation of all of the vertices having the same state as the initial global majority state, regardless of the initial state configuration? In [19], Mustafa and Pekeč analyzed democratic consensus computers (d.c.c.’s), graphs that exhibit precisely this behavior. They produced a variety of results, including that a graph that is a d.c.c. has diameter at most four, has a trivial min-cut and a non-unique max-cut, and that for a graph on vertices with maximum degree at least , if is a d.c.c., then must have a vertex (called a master) that is adjacent to every other vertex. They conjectured that every d.c.c. contains a master, and the question of a full characterization of d.c.c.’s remains open. For a given graph, one might wonder what subsets of the vertex set can “force” a majority process to the situation with every vertex in the same state. More formally, a set of vertices is a dynamic monopoly, or dynamo for short, if when a majority process is started with all of the vertices in in state 1, then the system reaches the fixed point with all vertices being in state 1. Dynamos can be incredibly small; Berger [1] proved that for every , there exists a graph with at least vertices containing a dynamo of size 18. In addition, in [21], Peleg showed that for all graphs, the size of a dynamo that converts all vertices to state 1 in a single step (Peleg called such a dynamo a (static) monopoly) must be , and that there are graphs with static monopolies of size . We shall consider a similar concept for irreversible threshold processes. One can also consider special classes of dynamos. For example, if is a dynamo and denotes the set of vertices in state 1 at time (), then as in monotone -threshold processes, is a monotone dynamo if for all . In addition, one can consider a different update rule where the state of a vertex changes from 0 to 1 if and only if more than half of its neighbors are in state 1, but not vice versa. These irreversible majority processes (and the associated irreversible dynamos) appear in the context of fault-tolerant computing, where a vertex entering state 1 corresponds to permanent faults occurring in the associated processor. In [5], [6], [17], bounds are shown for both monotone and irreversible dynamos in toroidal graphs of various types and in butterflies, two graph structures which appear in parallel computation. In addition, Peleg [22] showed that for most variants of majority processes, every monotone dynamo is of size . There are, however, some majority process variants with constant size monotone dynamos; in the majority process where a vertex’s state is counted in its own majority determination and in case of ties, the vertex enters state 1, any pair of adjacent vertices in the cycle is a monotone dynamo. Majority processes are an interesting model for opinion change in social networks. However, in the case of disease, they do not seem to be applicable. Motivated by the disease application, we can think of modifying an irreversible threshold process by allowing us to “vaccinate” certain vertices that are in state 0. Once a vertex is vaccinated, it can never change into the infected state 1. There are several variants of the vaccination problem that one can consider. One is to only allow vaccination at the beginning (time ), say of vertices (where reflects a limitation on the amount of vaccine available). An alternative is to allow us to vaccinate vertices at time . The vaccination problem has been studied in the context of firefighting. We can think of the vertices as trees in a forest, and the state 1 corresponding to a tree being on fire. A fire spreads to neighboring trees under different rules, e.g., in an irreversible 1-threshold process. In such a process, if a tree is burning at time , all its neighbors burn at time except those that are protected by “vaccination” or, what is equivalent, by placement of firefighters on them at a time no later than time . There is a rather robust literature about this problem, particularly on finite and infinite planar grids and trees ([7], [10], [11], [18], [20], [24], to list a few references). In the rest of this paper, we will not allow vaccinations. However, we will allow epidemics to break out at more than one vertex. For the rest of this paper, we consider a -irreversible threshold process on a graph . Let be the set of all vertices in state 1 at time and let be similarly defined. Motivated by the notion of dynamo discussed in Section 1.2, we define a set of vertices to be an irreversible -conversion set if when an irreversible -threshold process is started with , then the process reaches the situation with all of the vertices in state 1. We can define -conversion sets for ordinary -threshold processes in an analogous manner, and in such processes, it is also of interest to study -conversion sets in monotone -threshold processes. These two concepts are studied in detail in [4]. For a graph , let denote the size of the smallest irreversible -conversion set for . We will be interested in determining and in finding irreversible -conversion sets of minimum size. Irreversible conversion sets have clear and interesting interpretations in various applications. In the case of opinion, they correspond to sets of individuals with the property that if they agree with an opinion at the beginning, then the entire group is guaranteed to agree with that opinion eventually. In the case of disease, we are looking for a set of vertices that an intentional attacker (bioterrorist) can infect to guarantee to infect all vertices eventually. In terms of bioterrorism defense, the goal is to design social networks so that irreversible conversion sets are either impossible to find or always large enough to offset the ability of an attacker to infect that large a set. Irreversible conversion sets also have an economic interpretation. We think of achieving saturation of a new product in a market by placing it with vertices of an irreversible conversion set in a process where one buys a product if sufficiently many of our neighbors own it. The rest of this paper is organized as follows: In the next section, we show that it is NP-complete to determine if a graph has an irreversible -conversion set of size at most , at least for . Section 3 considers special graphs: complete multipartite graphs, trees, and various kinds of grids. It gives exact values and bounds on the size of the smallest irreversible -conversion set for such graphs. Section 4 gives a general lower bound for the size of such a set in a general graph. We close with a discussion of open problems and future research.1.1. Threshold processes
1.2. Majority processes
1.3. Vaccination
1.4. -conversion sets
2. NP-completeness of the problem of finding the smallest irreversible -conversion set
The Irreversible -CONVERSION SET problem is as follows:
Irreversible-CONVERSION SET (-CS)
Instance: A graph and a positive integer .
Question: Does have an irreversible -conversion set where ? (Equivalently, is ?)
We will show that this problem is NP-complete for a fixed , but we must first prove a result about irreversible -conversion sets in -regular graphs that will also be useful later.
Lemma 1 Let be a connected -regular graph, and let . Then is an irreversible -conversion set if and only if is independent.
Proof In an irreversible -conversion process, a vertex will switch from state 0 to state 1 if and only if all of its neighbors are in state 1. If is independent, then every vertex in will switch to state 1 at time . If is not independent, then two vertices starting in state 0 share an edge, and those adjacent vertices will never enter state 1. Hence, cannot be an irreversible conversion set. □
Corollary 1 Let be a connected graph with no vertex of degree greater than , let be the set of vertices of degree less than , and let . Then is an irreversible -conversion set if and only if and is independent.
Although we do not use this result immediately, we will now note a similar result about irreversible -conversion sets in -regular graphs. In a graph , a feedback vertex set is a subset of the vertices such that the subgraph induced by the vertices in is cycle-free (equivalently, is a forest).
Proposition 1 Let be an -regular graph, and let . Then is an irreversible -conversion set if and only if is a feedback vertex set.
Proof If is not a feedback vertex set, then contains a cycle. Suppose . At every time step, every vertex in the cycle will have at least two neighbors in state 0. If is a feedback vertex set, then is a forest, so at each time step, all of the isolated vertices and leaves in the forest induced by the vertices in state 0 will enter state 1, leaving a smaller forest, and this process will eventually reach the fixed point with the entire vertex set in state 1. □
In a graph , a detour is a simple path (no repeated vertices) of maximum length, and the detour number, denoted , is the length of such a path. The detour number was first studied in [14]. The transient length of an irreversible threshold process is that value of time beyond which no vertex enters state 1.
Lemma 2 For an irreversible -threshold process on a connected graph of order , the transient length is at most ( if ).
Proof Assume we have an irreversible -threshold process on a graph , and assume that there exists a vertex which enters state 1 at a time . Let denote the set of vertices that enter state 1 at time (), and assume that for all . The are disjoint and every vertex in is adjacent to some vertex in . Thus, every vertex in is the endpoint of a simple path of length . If , then there is a simple path in of length , a contradiction. This bound can be improved if . If , every vertex in has at least two neighbors in . We will show that every vertex in is in a simple path of length that includes a vertex in each of . Let . It is the endpoint of a simple path of length , , , with . If any of ’s neighbors are not on , then clearly has a simple path of length containing . If all of ’s neighbors are on , we construct a simple path of length as follows. We will show by induction that every edge in is contained in a simple path of length consisting solely of vertices in . The claim is clearly true for ; has another neighbor in besides , so there is a path of length two containing the edge . Assume the claim is true for all . Consider the edge in , with . If has a neighbor not on , then is a simple path of length containing the edge . If all of the neighbors of in are on , then let be a neighbor of on , . By our inductive assumption, we know there is a simple path of length containing the edge consisting of vertices in . Let , where and are disjoint simple paths in . Consider the path . The vertices that were in are distinct, and the vertices are in for , so they are distinct from the vertices in and from each other, so the path is simple, of length , and consists solely of vertices in . Thus, if , , or . □
Remark This result gives a sharp bound on transient length. Consider an irreversible 1-conversion process in a path of vertices. If consists of one of the endpoints of the path, then the process completes entering all vertices in state 1 at time . This can be generalized for any irreversible -conversion process by starting with a path , adding leaves to and leaves to each of . If consists of all of the leaves added to the path, then vertex enters state 1 at time . Therefore, the process reaches the situation of all vertices in state 1 at time , and the longest path in this graph is of length , involving all of the vertices on the path and at most two leaves.
Theorem 1 For a fixed integer , IRREVERSIBLE -CONVERSION SET is NP-complete.
Proof By Lemma 2, we know that the transient length of an irreversible threshold process is polynomial in the number of vertices, so IR-CS is in NP. Our reduction will be from the INDEPENDENT SET problem, given below. INDEPENDENT SET Instance: A graph and an integer . Question: Is there an independent set in of size ? Fricke, Hedetniemi, and Jacobs [9], using a reduction from NOT-ALL-EQUAL-3SAT, show that for a fixed , INDEPENDENT SET is NP-complete for the class of -regular, non-bipartite graphs (the problem is in for the class of bipartite graphs). Given an instance of INDEPENDENT SET for a -regular, non-bipartite graph : “Does have an independent set of size ?”, we construct an instance of IR-CS: “Does have an irreversible -conversion set of size ?” By Lemma 1, a vertex set is an irreversible -conversion set in a -regular graph if and only if its complement is independent. Therefore, there is an irreversible -conversion set of size if and only if there is an independent set of size . So, there is a ‘yes’ answer to the instance of IR-CS if and only if there is a ‘yes’ answer to the instance of INDEPENDENT SET. □
The complexity of IR2-CS still remains open. We cannot use the above arguments, since it is trivial to find independent sets in 2-regular graphs (i.e., vertex-disjoint unions of cycles). In addition, Li and Liu [15] have developed a polynomial algorithm for finding a minimum feedback vertex set in 3-regular graphs, so we cannot use Proposition 1.
3. Irreversible -conversion sets for special graphs
In this section, we give constructions for minimum -conversion sets for some classes of graphs. For other classes, we present constructive upper bounds for the value of and give proofs of lower bounds.
Proposition 2 For the path and cycle on vertices, and .
Proof For simplicity, we will denote the vertices of () in order as . By Lemma 1, if is an irreversible 2-conversion set in , must be independent. Thus, and . Picking the complement of an independent set of size in will be an irreversible 2-conversion set of size . For the path , by Corollary 1, both and must be in any irreversible 2-conversion set . The largest independent set in not containing or has size , hence . Constructing irreversible 2-conversion sets of size is straightforward. □
In this section, we consider the complete multipartite graph consisting of classes of vertices, with vertices in the th class , and with edges between all vertices in different classes. We will assume that . Proposition 3 For the complete bipartite graph , : Proof Let and be the two partite sets of , with . Any subset of of size will be an irreversible -conversion set if . Clearly is the lower bound for the size of any irreversible -conversion set. If and , is an irreversible -conversion set. There are not enough vertices in to change the state of any vertex in , so all of must be in state 1 at time . If both partite sets are smaller than , then is the only -conversion set. □ Consider the complete -partite graph with . For any subset , let . A union of partite sets has size just above if and if is created by removing the partite set of highest index from (), then . Proposition 4 Let be a complete -partite graph with . Suppose that there is a union of partite sets such that: has size just above , and , where is the largest element of . Proof Let be a union of partite sets of satisfying conditions (1) and (2). If has size exactly , then clearly is an irreversible -conversion set for . If does not have size exactly , let , where is the largest element of , and let consist of and vertices in . At time , all vertices will be in state 1 except for the state 0 vertices in , which will enter state 1 at time by condition (2). Thus, . Since must be , equality is proven. Assume that no union of partite sets satisfies conditions (1) and (2). Let be a minimum irreversible -conversion set for , and let be the smallest union of partite sets that contains . At time , every vertex in will be in state 1. If no , contains elements of , then has size exactly and the result is straightforward. Thus, let us assume that some , , contains elements of . We will reach a contradiction. We may assume that at most one partite set has vertices in . If two partite sets contain vertices in , let and . Removing vertices in from and adding the same number of vertices from to gives a new irreversible -conversion set for . This is because nothing changes in the other partite sets during the threshold process with as compared to the one with , because will be in state 1 at time , and because a vertex in will have at least as many state 1 neighbors in at time as it did in the original process. We can also assume that the partite set in containing members of is the partite set in of highest index. To see this, let be the partite set in containing vertices in . If , consider the new set which removes of the vertices in from and includes all of . Clearly, is an irreversible -conversion set. If , then consider the new set that removes all of the vertices in from and includes more vertices from . We let and repeat the process of removing vertices in from the partite set of highest index in and adding vertices from to until either no partite set contains vertices in and , or the only partite set that contains vertices in and is the partite set of highest index in . Suppose that is the partite set of highest index in and is the only partite set in that contains vertices in . If , the 0 vertices in will never switch to 1 and this contradicts being an irreversible -conversion set. Suppose . Then does not have size just above . Clearly, since is an irreversible -conversion set. This implies that where . But now is an irreversible -conversion set. This is because if , at time all vertices in become 1 and now . This contradicts being a minimum irreversible -conversion set since is a proper subset of . □ Proposition 5 Let be the complete -partite graph with and . Then . Proof Let be any subset of of size , and let for . We claim that is an irreversible -conversion set. Since every vertex in has neighbors in state 1 in , . If , then at time , has at least neighbors in state 1 in . Since , has at least neighbors in , so . Hence, is an irreversible -conversion set. □ Remark Let be the complete -partite graph with . If and , then . For a tree , let denote the set of leaves of , i.e., vertices of degree 1. Clearly, must be contained in any irreversible -conversion set for , and in some circumstances, is a minimum irreversible -conversion set. We will use the term internal vertex in a tree to be any vertex that is not a leaf and use to denote the subtree induced by vertices in set . Proposition 6 If is a tree with every internal vertex of degree for , then . Proof Let be an internal vertex of , and let , where the eccentricity is the greatest distance between and any other vertex. The result follows by proving the following statement by induction on : “At time , in a -threshold process with , .” □ For a graph , a set is a vertex cover for if for every edge , or . It is easily shown that is a vertex cover if and only if is independent. Proposition 7 Let be a tree with every internal vertex of degree . Then is an irreversible -conversion set if and only if , where is a vertex cover of . Proof If is an irreversible -conversion set for , it must contain , and must be independent, since if two 0 vertices shared an edge, they would never become 1. Hence, for some vertex cover of . It remains to show that for any vertex cover of , is an irreversible -conversion set. Since is independent, all of the vertices in will become 1 at time . □ To extend these results to all trees, we make a few simplifying assumptions. We will restrict our analysis to trees with all internal vertices of degree . If there is an internal vertex with , it will have to be in any irreversible -conversion set. Thus, to every neighbor of , acts exactly like a leaf, so we can break into subtrees with as a leaf in each subtree. We define a -region of to be a maximal connected subset of vertices of degree . We can partition the vertex set of a tree into leaves, -regions, and non--regions (maximal connected subsets of vertices of degree ). For a -region , define the outer boundary of (denoted ) to be the set of vertices not in that have a neighbor in (or, ) and the inner boundary of (denoted ) to be the set of vertices in adjacent to vertices not in (or, ). Proposition 8 Let be a tree with every internal vertex of degree . If is a vertex cover for every subtree , where is a -region of , then is an irreversible -conversion set for . Proof If is a vertex cover for every subtree , where is a -region, then after one time step, every vertex in every -region will be 1, since the set of 0 vertices in a -region will be independent. At time , every non--region will be surrounded by vertices in state 1, and will proceed to turn to 1 as in the proof of Proposition 6. □ If is an irreversible -conversion set for a tree where every internal vertex has degree , then must be independent for every -region . This is guaranteed by Proposition 8, but there is no guarantee that the irreversible -conversion set described in Proposition 8 is minimum. If is an inner boundary vertex of some -region , and is its neighbor on the outer boundary of , then there is no need for to cover the edge , as long as other neighbors of eventually enter state 1. We now present an algorithm to compute that runs in time, where is the number of vertices in the tree. For simplicity, we will present the algorithm for and then discuss how to implement it for larger . For a tree rooted at a vertex , a set of vertices is an almost irreversible 2-conversion set (a-I2CS) if when an irreversible 2-conversion process is started with the vertices in in state 1, then all of the vertices in enter state 1. The root vertex may be in to convert other vertices in , but it is not necessary for to be in state 1 when the process reaches its fixed point. However, if the root vertex of has degree 2 or greater, then any a-I2CS for is also an irreversible 2-conversion set. If is not in , then will enter state 1 after two of its neighbors enter state 1. We define the composition of a tree with a tree and the composition of tree-subset pairs and with being an a-I2CS. We use the rule of composition that joins a tree rooted at and a tree rooted at (with ) by adding the edge and specifying as the root of the resulting larger tree . Let and be two trees rooted at vertices and , respectively. Also, let and be a-I2CS’s for and , respectively. If is the composition of with , then consider the set . might not be an a-I2CS for , since if is not converted to state 1 by in , then has at most one neighbor in . If is not an isolated vertex in and is converted to state 1 in by , then will eventually have two neighbors in in state 1, so will enter state 1 as well. If not, then either or must be added to for it to be an a-I2CS. Let be the resulting a-I2CS for . We will use the notation to denote this composition of two trees and their associated a-I2CS’s. To construct an algorithm to compute for any tree , we characterize the class of possible tree-subset pairs that can occur in a tree with root vertex , where is an a-I2CS of . For this problem, there are four classes of tree-subset pairs: If , we can determine the class of if is of class and is of class for each pair , . The classes of resulting from each composition are summarized in Fig. 1. The verification of the class of for each pair is straightforward. Note that in some cases it is necessary to add or to in order to convert . For example, if one merges any tree-subset pair with a tree-subset pair of class , since the degree of in the resulting tree is one, it must be added to since it cannot be converted otherwise. Fig. 1. The number in brackets in the th entry in this table gives the class of that results from composing of class with of class . Examples are given for each pair . The black vertices denote the vertices that appear in . A dotted edge represents the edge joining the roots of and . If it was necessary to add a copy of or to to convert , the added vertex is shown in grey, and the number of vertices that are added to is noted in parentheses. In some cases, the class of the resulting composition is different depending on whether or is added to . Using Fig. 1, we can write out a system of recurrence relations that capture how we can obtain a tree-subset pair of class after composition. It is impossible to compose two tree-subset pairs to get a tree-subset pair of class , because when two trees are composed, the root vertex has degree greater than zero. The superscripts and on some of the compositions denote that or has been added to to make an a-I2CS for . (1)(2)(3) For example, to see recurrence (2), we note that if is of class , then must be of class or and can be of any class. However, if is of class or if both and are of class , then must be added to to produce an a-I2CS. Given a tree of order , we root it at an arbitrary vertex that we label . We then label the remaining vertices using the well-known depth first search procedure so that if is ’s parent, then . For each vertex , let equal the subscript of the label of the parent of (assume that ). The algorithm will build the tree one edge at a time. At step , the algorithm builds a graph with vertex set and designates a root for each maximal connected subtree of . is obtained from by composing the subtree rooted at with the subtree rooted at . In graph there will be a 4-place class vector associated to each root vertex. If is the root of a subtree in , then the th entry of its class vector at step equals the size of the smallest a-I2CS for such that is of class . The size can be if there is no such a-I2CS. consists of isolated vertices, each the root of a one vertex tree. The class vector associated with each vertex in is , where ‘’ is the entry for class and , as the degree of an isolated vertex is zero, so it is impossible for a tree-subset pair to be of class or at the start of the algorithm. For classes and , we associate the a-I2CS’s and of sizes 1 and 0, respectively. At step of the algorithm, the only class vector that will change will be the one associated with the vertex . This is updated using the table in Fig. 1 and the recurrences (1)–(3). is the tree rooted at . Therefore, the th entry of the class vector associated to is the size of the smallest a-I2CS of such that is of class . If is of class or , then is also an irreversible 2-conversion set, but not if is of class or , since is not converted by . Thus, equals the minimum of the first two entries of the class vector associated to . We now have all the ingredients for an algorithm to compute for a tree . The input to the algorithm is an -place vector for , whose th entry is , as defined above. The output of the algorithm is the value of . Every vertex has an associated class vector whose th entry is denoted . Every class vector is initialized with the vector . At step , where runs from 0 to , the tree rooted at is composed with the tree rooted at producing a tree . The vector is updated by a procedure called combine that is derived directly from the recurrences (1)–(3) summarized in Fig. 1. equals the size of the smallest a-I2CS for such that the pair is of class . For example, let be the composition of a tree rooted at with a tree rooted at . Let be a minimum a-I2CS for such that is of class . By relation (3), if is of class , then is the union of the smallest a-I2CS for such that is of class and the smallest a-I2CS for such that is of class or , or the smallest a-I2CS for such that is of class or , plus the vertex . Since equals the size of the smallest a-I2CS for such that is of class , and is similarly defined for , then will be updated to equal the minimum of , , , and . The other values of are updated similarly. After step , all of the edges of the tree rooted at have been included. Therefore, equals the size of the smallest a-I2CS for such that is of class . We want to be an irreversible 2-conversion set for , so will be the minimum of and . An example of the algorithm is illustrated in Fig. 2. Fig. 2. An illustration of the algorithm to compute . At step , the root vertex of each maximal connected subtree in is denoted with a filled in circle, and the class vector associated to each root vertex at step is noted next to the vertex. The bold edge denotes the edge added by the composition of the subtree rooted at with the subtree rooted at . The filled in circles in the lower right tree give the minimum irreversible 2-conversion set for . Note that is of class , as the root vertex is not in . Also, the minimum of the first two entries of the class vector associated to the root vertex of is the entry associated with class . Pseudo-code for the algorithm is given below. Since the procedure combine takes time to compute, it is clear that the procedure runs in time. To summarize, we have the following theorem. Theorem 2 For a tree of order , the procedure (Parent) computes in time. To extend this algorithm to compute for , one simply defines additional classes of trees: one where the root is in , one where the root is not in but its degree is at least , and classes where the root is not in and the degree is , . For a fixed , running the combine procedure takes time, and the procedure is executed times, producing an algorithm. Let be the standard grid. Here, the vertices will be denoted , , , and there is an edge between and if and only if . We will frequently refer to the parity (even or odd) of as the parity of a vertex . For simplicity, we will refer to even vertices or odd vertices, instead of vertices of even (odd) parity. We also introduce the circulant grid , which is the standard grid with the edge added for every , and the toroidal grid which is the circulant grid with the edge added for every , . We will refer to the rows and columns of a grid by and . The toroidal grid is 4-regular, which allows us to take advantage of some of our earlier results on conversion sets in regular graphs. Proposition 9 For the toroidal grid graph : Proof By Lemma 1, is an irreversible 4-conversion set in a 4-regular graph if and only if is independent. must be a vertex cover for each , and must be a vertex cover for each , so each row must contain at least state 1 vertices and each column must contain at least state 1 vertices. Hence, (4) We can construct a 4-conversion set as follows. Given a row of a grid and a set of vertices in that row, is a left cyclic shift of in row if for and . Similarly, is a right cyclic shift of in row if for and . Assume that is odd, and let be a vertex cover for of size . If is even, then let be the right cyclic shift of in , and then let be the left cyclic shift of in , and so on. If is odd, then assume that . For the first rows (or rows if ), generate a vertex cover for via a right cyclic shift of the cover for (where is the vertex cover for ). If , then the cover for is the left cyclic shift of , and if , then the cover for is the right cyclic shift of . If , in , let be the left cyclic shift of the vertex cover in , and in , let be the right cyclic shift of , and continue to alternate cyclic shifts to the left and right so that the cover for row is the right cyclic shift of . Let be the union of all of these vertex covers. Since we have assumed that is odd and that if is odd, then , . We claim that forms an irreversible 4-conversion set of size , which together with the lower bound in Eq. (4) gives us the desired result for the case when is odd. No vertex in has a neighbor in in its row, since is a vertex cover for each row. Also, no vertex in has a neighbor in in its column, since both of its neighbors in the row are in , so when the vertex cover in its row is cyclically shifted, a vertex in will have a neighbor in in the row above and below it. Finally, it is easy to see that either the first or the last vertex in each row or columns is in . Hence, is independent. If both and are even, then simply start with a vertex cover of row by vertices and generate each subsequent row’s covering by a right cyclic shift. The proof that is independent is straightforward. □ Proposition 10 For the standard grid graph : Proof Because all of the vertices on the border of the grid have degree less than 4, they must be in any 4-conversion set . The vertices in the interior of the grid have degree 4, and no two vertices in can be neighbors. Hence, , corresponding to the set of even non-border vertices or the set of odd non-border vertices. Let be the larger of these two sets. It can be shown (as in Corollary 1) that every vertex within distance of any border vertex will be in state 1 for all , as long as remains independent for all , which is guaranteed by the independence of . Thus, the smallest irreversible 4-conversion set for has size . □ The following results on irreversible 3-conversion sets in grids are derived from the work on “irreversible dynamos” by Flocchini, et al. [5], [6], and from the work on minimum feedback sets in grids by Luccio [16]. Proposition 11 For the toroidal grid graph and the standard grid graph , assume that . Then: . The results in this proposition come from the work of Flocchini, et al. and Luccio, where they provided constructions of strong dynamos and feedback vertex sets in toroidal grid graphs, which correspond to irreversible 3-conversion sets. With some small modifications, their results can also be modified for standard grid graphs (which are no longer 4-regular, hence requiring some additional border vertices in ). Some sample constructions are given in Fig. 3. Fig. 3. Some constructions of conversion sets based on work from Flocchini et al., and Luccio. (a) A minimum irreversible 2-conversion set for . (b) A minimum irreversible 3-conversion set for . We close this section with some results on and grids. Theorem 3 For the toroidal grid graph , the standard grid graph , and the circulant grid graph : Proof For the toroidal grid graph , by Proposition 1, is an irreversible 3-conversion set if and only if is cycle-free. Every column must contain a vertex in , so . However, if every column has exactly one vertex in , then every column has two adjacent vertices in . Therefore, in every column , some vertex in has a neighbor in and some vertex in has a neighbor in . Therefore, the subgraph induced by contains a cycle which contains at least one vertex from each column in the grid, which is a contradiction. Hence, . Let be and all of the even vertices in rows . is cycle-free, so is an irreversible 3-conversion set of size . For the standard grid graph , if is an irreversible 3-conversion set, must contain the four corners, since they have degree two. In addition, must be a vertex cover of and , since all of the non-corner border vertices have degree three. Finally, there cannot be a path in connecting two border vertices in , as none of the vertices along such a path would ever enter state 1, since all of the non-border vertices would always have two state 0 neighbors, and the two border vertices would always have a state 0 neighbor. So, every path between two border vertices in must contain a vertex in in . If a column contains two border vertices in , then must be in . In addition, we know that , , , and must be in , since there cannot be two consecutive border vertices in . Let be a maximal collection of consecutive columns where each column has two border vertices in . We can assume that no vertex , , is in . To see this, let be a minimum irreversible 3-conversion set that contains for some , . is also a minimum irreversible 3-conversion set, because any shortest path between two border vertices one of which is in a column for and that contains also contains . Similarly, is a minimum irreversible 3-conversion set. If either or was already in , then was not minimal, since could be removed. Let be a maximal collection of consecutive columns where each column has one border vertex in . Since must be independent in the border vertices, the border vertices in will have to alternate rows in consecutive columns. In addition, columns and must have two border vertices in , since if is in , must be in , and since the collection of columns is maximal, must also be in . For every consecutive pair of columns and in this collection, at least one of and must be in , or else there would be a path of length four between two border vertices in . Therefore, every column that contains only one vertex in has two neighboring columns that each contain two vertices in . If is minimal, we can assume that no column contains three vertices in . Thus, if is the number of columns with vertices in , then , , and . If is even, then , and . If is odd, then , and . If is odd, letting be the even vertices of gives an irreversible 3-conversion set of size . If is even, letting be the even vertices of except , and including and gives an irreversible 3-conversion set of size . For the circulant grid graph , we apply the same argument as above, that every column with one vertex in has two neighboring columns each containing two vertices of . Again, , , and . If is even, then , so . If is odd, then , so . In both cases, letting consist of the even vertices in gives an irreversible 3-conversion set of the desired size. □ Some of the conversion sets constructed in the proofs above are shown in Fig. 4. Fig. 4. Minimum irreversible conversion sets (c.s.) for various and grids. (a) An irreversible 3-c.s. for . (b) An irreversible 3-c.s. for . All edges that extend beyond the border of the grid are assumed to wrap around to the opposite side.3.1. Complete bipartite and multipartite graphs
Then . Otherwise, .3.2. Trees
3.3. Grids
[5], [6], [16]
4. Lower bounds on
For the irreversible threshold process, we have the trivial lower bound , and there are several graphs that achieve this bound.
The following lower bound for irreversible conversion set sizes for -regular graphs is obtained using a method suggested by Eli Berger [2].
Proposition 12 If is an -regular graph with vertices, for .
Proof We define a slightly different method for updating an irreversible -threshold network. Instead of switching every 0 vertex with or more 1 neighbors simultaneously, at each time step, we pick any 0 vertex with or more 1 neighbors and switch it to 1, and repeat this process until no more vertices can be switched. If is an irreversible -conversion set under the normal update rule and is the set of vertices that enter state 1 at time , then under the new update rule, one can switch all of the vertices in one at a time, then the vertices in , and so on. Therefore, if is an irreversible -conversion set under the normal update rule, will also be an irreversible -conversion set under this new update rule. Define the energy of a network to be the number of edges between vertices of opposite state. Let be the energy of the network at time . For every vertex that we switch from 0 to 1, we lose an amount of energy of size at least and create an amount of energy of size at most . Hence, switching a vertex from 0 to 1 reduces the amount of energy in the network by at least , so . If is an irreversible -conversion set, then . Therefore, , or . However, is at most , so , and . □
If we consider irreversible 4-conversion on the toroidal grid , this result provides a sharp bound () for irreversible 4-conversion sets when and are even.
5. Discussion and open problems
The results described here have been formulated mostly in the language of spread of disease. As models of disease spread, they are somewhat oversimplified. However, we feel that models are tools for reasoning about problems and, as such, these models provide ways to formalize concepts such as effect of the topology of social networks on the spread of disease; alternative vaccination strategies; and other concepts. As we have noted, the concepts and ideas developed here also have potential interest in other applications. They were originally developed for understanding how opinions might spread through social networks. Clearly the work also has interest in the context of the “firefighter problem”. We hope that the results and ideas are also of interest as mathematical concepts and problems.
There are several questions that would provide avenues for further research.
- 1.
We would like a characterization of irreversible minimum conversion sets for all complete multipartite graphs and trees, as well as other classes of graphs (chordal graphs, interval graphs, planar graphs, other kinds of grids, etc.), or possibly algorithms for finding such sets.
- 2.
We would like to know the complexity of the IR2-CS problem.
- 3.
In our work, updates to the vertices occur in parallel at each time step. We could update individual vertices serially, which would certainly have an effect on the dynamics of the network. In addition, we have considered a deterministic update rule. In [12], [13], Hassin and Peleg replace the deterministic update rule with a random one, where the probability that a vertex switches signs is equal to the fraction of its neighbors with the opposite sign. We could also replace this with any probability function, perhaps one based on some threshold (for example, a vertex switches sign with probability where is the number of neighbors with the opposite sign: if the vertex switches sign automatically).
- 4.
As we have noted, analogous results on the -threshold processes where one can switch back from state 1 to state 0 would be of interest. There are extensive results about this process and also the process we have called monotone in [4]. However, exact results for minimum conversion sets are still missing for the most part. Other models of interest could allow switching back from state 1 to state 0 automatically after a given number of time steps (given number of days for a disease to run its course). We could also bring in “immunity” after a vertex has changed from state 1 to state 0. Probabilistic variants of all of these models would also be of interest.
- 5.
We have described vaccination strategies in Section 1.3. A number of questions and research directions remain open here. For instance, the bulk of the work described in Section 1.3 involves rectangular grids and trees. The same problems are of interest for chordal graphs, interval graphs, planar graphs, and other kinds of grids such as annular grids. There is very little known about vaccination strategies where we have a varying number of doses of vaccine each time period. Also, the models considered so far do not consider the possibility that vaccination will only be successful with a certain probability and they do not consider the possible side effects of vaccination, which could be brought into the models in various ways.
- 6.
Other public health interventions could be made precise in the language of the irreversible -threshold processes. For example, can we make precise the notion of “ring quarantine” — quarantining all individuals within a certain graph distance from newly-infected individuals? What are “optimal” values of ? Could we have quarantines only lasting a certain number of time steps? Could we bring in non-compliance with quarantines?
- 7.
We would like to find graph topologies that protect against having small irreversible -conversion sets. For example, one could ask, given and , what are graphs on vertices that have no irreversible -conversion sets of or fewer vertices?
In short, the results in this paper are just a beginning. They illustrate the types of analyses one can make using graph-theoretical models of the spread of disease or opinion. There is much more work to do.
Acknowledgements
The authors would like to thank the National Science Foundation for supporting this work through grants CCR 91-19999, BIR 94-12594, DBI 99-82983, SBR 97-09134, SES-92-11492 and EIA-02-05116 to Rutgers University.
References
- [1]
- E. Berger, Dynamic monopolies of constant size, Technical Report, Los Alamos National Laboratories, 1999. Available at http://xxx.lanl.gov/abs/math/9911125
- [2]
- E. Berger, Personal communications via e-mail, May 1999
- [3]
- M.H. DegrootReaching a consensusJournal of the American Statistical Association, 69 (1974), pp. 167-182
- [4]
- P.A. Dreyer, Applications and variations of domination in graphs. Ph.D. Thesis, Rutgers University, 2000
- [5]
- P. Flocchini, E. Lodi, F. Luccio, L. Pagli, N. Santoro, Irreversible dynamos in tori, in: Proceedings of EUROPAR, 1998, pp. 554–562
- [6]
- P. Flocchini, E. Lodi, F. Luccio, L. Pagli, N. Santoro, Monotone dynamos in tori, in: Proceedings of the 6th Colloquium on Structural Information and Communication Complexity, 1999
- [7]
- P. Fogarty, Catching the fire on grids, Masters Thesis, University of Vermont, 2003
- [8]
- J.R.P. FrenchA formal theory of social powerPsychology Review, 63 (1956), pp. 181-194
- [9]
- G.H. Fricke, S.T. Hedetniemi, D.P. JacobsIndependence and irredundance in -regular graphsArs Combinatorica, 49 (1998), pp. 271-279
- [10]
- S.G. Hartke, Graph-theoretic models of spread and competition, Ph.D. Thesis, Rutgers University, 2004
- [11]
- B. Hartnell, Q. LiFirefighting on trees: How bad is the greedy algorithm?Congressus Numerantium, 145 (2000), pp. 187-192
- [12]
- Y. Hassin, D. Peleg, Distributed probabilistic polling and applications to proportionate agreement, in: Proceedings of the 26th International Colloquium on Automata, Languages, and Programming, 1999, pp. 402–411
- [13]
- Y. Hassin, D. Peleg, Extremal bounds for probabilistic polling in graphs, in: Proceedings of the 7th International Colloquium on Structural Information and Communication Complexity, 2000
- [14]
- S.F. Kapoor, H.V. Kronk, D.R. LickOn detours in graphsCanadian Mathematics Bulletin, 11 (1968), pp. 195-201
- [15]
- D. Li, Y. LiuA polynomial algorithm for finding the minimum feedback vertex set of a 3-regular simple graphActa Mathematica Scientia, 19 (1999), pp. 375-381
- [16]
- F. LuccioAlmost exact minimum feedback vertex set in meshes and butterfliesInformation Processing Letters, 66 (1998), pp. 59-64
- [17]
- F. Luccio, L. Pagli, H. Sanossian, Irreversible dynamos in butterflies, in: Proceedings of the Sixth Colloquium on Structural Information and Communication Complexity, 1999
- [18]
- G. MacGillivray, P. WangOn the firefighter problemJournal of Combinatorial Mathematics and Combinatorial Computing, 47 (2003), pp. 83-96
- [19]
- N. Mustafa, A. Pekeč, Democratic consensus and the majority rule, Technical Report, BRICS, 2000. Available at http://www.brics.dk/RS/00/8/
- [20]
- K.L. Ng, P. RaffA generalization of the firefighter problem onDiscrete Applied Mathematics, 156 (2008), pp. 730-745
- [21]
- D. Peleg, Local majorities, coalitions, and monopolies in graphs: A review, in: Proceedings of the Third Colloquium on Structural Information and Communication Complexity, 1996
- [22]
- D. PelegSize bounds for dynamic monopoliesDiscrete Applied Mathematics, 86 (1998), pp. 263-273
- [23]
- S. Poljak, M. SûraOn periodical behavior in society with symmetric influencesCombinatorica, 3 (1983), pp. 119-121
- [24]
- P. Wang, S.A. MoellerFire control on graphsJournal of Combinatorial Mathematics and Combinatorial Computing, 41 (2002), pp. 19-34