The domination number of on-line social networks and random geometric graphs

Anthony Bonato, Marc Lozier, Dieter Mitsche, Xavier Pérez-Giménez, Paweł Prałat

Introduction

On-line social networks (or OSNs) such as Facebook have emerged as a hot topic within the network science community. Several studies suggest OSNs satisfy many properties in common with other complex networks, such as: power-law degree distributions , high local clustering , constant or even shrinking diameter with network size , densification , and localized information flow bottlenecks . Several models were designed to simulate these properties , and one model that rigorously asymptotically captures all these properties is the geometric protean model (GEO-P) (see for models where various ranking schemes were first used, and which inspired the GEO-P model). For a survey of OSN models see , and for more general complex networks . A fundamental difference with GEO-P versus other models is that it posits an underlying feature or metric space. This metric space mirrors a construction in the social sciences called Blau space . In Blau space, agents in the social network correspond to points in a metric space, and the relative position of nodes follows the principle of homophily : nodes with similar socio-demographics are closer together in the space. We give the precise definition of the GEO-P model (actually, one of its variants, the so-called MGEO-P model) below. We focus on the MGEO-P model, since it simpler than GEO-P and generates graphs with similar properties.

The study of domination and dominating sets plays a prominent role in graph theory with a number of application to real-world networks. A dominating set in a graph GG is a set of nodes SS in GG such that every node not in SS is adjacent to at least one node in SS. The domination number of GG, written γ(G)\gamma(G), is the minimum cardinality of a dominating set in GG. Computing γ(G)\gamma(G) is a well-known NP-complete problem, so typically heuristic algorithms are used to compute it for large-scale networks. Dominating sets appear in numerous applications such as: network controllability , as a centrality measure for efficient data routing , and detecting biologically significant proteins in protein-protein interaction network . For more additional background on domination in graph theory, see .

In social networks, we consider the hypothesis that minimum order dominating sets contain agents with strong influence over the rest of the network. Our goal in the present paper is to consider the problem of finding bounds on dominating sets in stochastic models of OSNs, and also in real-world data derived from OSNs. We consider bounds on the domination number of a stochastic model (see next paragraph), and upper bounds for that model are well-correlated with real-world OSN data. We note that the domination number has been studied previously in complex network models, including preferential attachment , and recently in .

The OSN model we consider is called the memoryless geometric protean model (MGEO-P), first introduced in . The MGEO-P model depends on five parameters which consists of: the number of nodes nn, the dimension of the metric space mm, the attachment parameter 0<α<10<\alpha<1, the density parameter 0<β<1−α0<\beta<1-\alpha, and the connection probability 0<p≤1.0<p\leq 1.

With probability pp, the node vv is adjacent to an existing node uu satisfying D(v,u)≤I(ru)\mathcal{D}(v,u)\leq I(r_{u}), where the distances are computed with respect to the following metric:

and where ∥⋅∥∞\left\|\cdot\right\|_{\infty} is the infinity-norm. We note that this implies that the geometric space is symmetric in any point as the metric “wraps” around like on a torus. The volume of the space influenced by the node is rv−αn−βr_{v}^{-\alpha}n^{-\beta}. Then the next node arrives and repeats the process until all nn nodes have been placed. We refer to this model by MGEO-P(n,m,α,β,p)(n,m,\alpha,\beta,p).

We give rigorous bounds on the domination number of a typical graph generated by the MGEO-P. An event AnA_{n} holds asymptotically almost surely (a.a.s.) if it holds with probability tending to 11 as nn tends to infinity. Our main result on MGEO-P is the following.

If m=o(log⁡n)m=o(\log n), then a.a.s. the domination number of a graph GG sampled from the MGEO-P(n,m,α,β,p)(n,m,\alpha,\beta,p) model satisfies

where CC is any constant greater than 66. In particular, a.a.s. γ(G)=nα+β+o(1)\gamma(G)=n^{\alpha+\beta+o(1)}.

We defer the proof of Theorem 1.1 to Section 2. It is noteworthy that the domination number of the preferential attachment model is linear in the order of the graphs sampled; see ; this fact further demonstrates the differences between MGEO-P and other complex graph models.

Theorem 1.1 suggests a sublinear bound on the domination number for OSNs, and we evidence for this in real-world data. In Section 3, we find bounds for graphs in the Facebook 100 data set, and compare these results to those for the stochastic models. We chose to work with the so-called Facebook 100 (or FB100) data set, as it provides representative samples from the network of increasing orders. Hence, we may consider trends for the domination number in the data. While the data presented is our first and initial study, the bounds we find for the domination number of FB100 are of sublinear order, and these bounds are well-correlated with those from MGEO-P. Sublinear domination results for other complex networks were also reported in ; our approach is distinct as we consider social networks of increasing orders.

In addition to the results above, we find rigorous bounds on the domination number for classical random geometric graphs. Given a positive integer nn, and a non-negative real rr, we consider a random geometric graph G=(V,E)∈G(n,r)G=(V,E)\in\mathscr{G}(n,r) defined as follows. The node set VV of GG is obtained by choosing nn points independently and uniformly at random in the square S=2\mathcal{S}=^{2}. (Note that, with probability 11, no point in S\mathcal{S} is chosen more than once, and hence, we may assume that ∣V∣=n|V|=n.) For notational purposes, we identify each node v∈Vv\in V with its corresponding geometric position v=(vx,vy)∈Sv=(v_{x},v_{y})\in\mathcal{S}, where vxv_{x} and vyv_{y} denote the usual xx- and yy-coordinates in S\mathcal{S}, respectively. Finally, the edge set EE is constructed by connecting each pair of nodes uu and vv by an edge if and only if dE(u,v)≤rd_{E}(u,v)\leq r, where dEd_{E} denotes the Euclidean distance in S\mathcal{S}.

Random geometric graphs were first introduced in a slightly different setting by Gilbert to model the communications between radio stations. Since then several closely related variants on these graphs have been widely used as a model for wireless communication, and have also been extensively studied from a mathematical point of view. The basic reference on random geometric graphs is the monograph by Penrose .

We note that our study is the first to explicitly provide provable bounds on the domination number of random geometric graphs. In particular, we derive the following result.

Let G∈G(n,r)G\in\mathscr{G}(n,r) and let ω=ω(n)\omega=\omega(n) be any function tending to infinity as n→∞n\to\infty. Then a.a.s. the following holds:

Denote by N(x)N(x) the minimal number of balls of radius xx needed to cover S\mathcal{S}. If r=Θ(1)r=\Theta(1), then

Define C=2π3/9≈1.209.C=2\pi\sqrt{3}/9\approx 1.209. If ωlog⁡n/n≤r=o(1)\omega\sqrt{\log n/n}\leq r=o(1), then

If 1/n≤r<ωlog⁡n/n1/\sqrt{n}\leq r<\omega\sqrt{\log n/n}, then

The proof of Theorem 1.2 is deferred to Section 4. The final section summarizes our results and presents open problems.

Proof of Theorem 1.1

For each node v∈[n]v\in[n], we consider the ball Bv={x∈[0,1)m:D(x,v)≤I(rv)}B_{v}=\{x\in[0,1)^{m}:\mathcal{D}(x,v)\leq I(r_{v})\}, which has volume bv=rv−αn−βb_{v}=r_{v}^{-\alpha}n^{-\beta}. The next lemma will be useful to estimate the sum of volumes of the balls corresponding to a set of nodes.

Let TT be a set of tt nodes (fixed before ranks are chosen) with ωnαlog⁡n≤t≤n\omega n^{\alpha}\log n\leq t\leq n, for a function ω\omega going to infinity with nn arbitrarily slowly.

Furthermore, given any integer ss such that 1≤s≤t1\leq s\leq t and s1−α≥ω(n/t)αlog⁡ns^{1-\alpha}\geq\omega(n/t)^{\alpha}\log n, a.a.s. all subsets S⊆TS\subseteq T of ss nodes satisfy

Observe that the sum in (1) is asymptotic to what one would expect. Indeed, if the ranks of the nodes in TT are distributed evenly, then one would obtain ∑i∈Ttri−α=∑i=1t(in/t)−α=tn−α/(1−α)+O(1)\sum_{i\in T}^{t}{r_{i}}^{-\alpha}=\sum_{i=1}^{t}(in/t)^{-\alpha}=tn^{-\alpha}/(1-\alpha)+O(1).

Let ss and tt be integers satisfying all the conditions of the statement in part (b). Set ω^=ω1/4→∞\hat{\omega}=\omega^{1/4}\to\infty, so we have t≥ω^4(t/s)1−αnαlog⁡nt\geq\hat{\omega}^{4}(t/s)^{1-\alpha}n^{\alpha}\log n. This also implies ω^=o((sn/t)1−α)\hat{\omega}=o((sn/t)^{1-\alpha}). Let YjY_{j} be the number of elements in TT with rank at most jj. Observe that YjY_{j} has expectation jt/njt/n, and follows a hypergeometric distribution. For (sn/t)1−α/ω^≤j≤n(sn/t)^{1-\alpha}/\hat{\omega}\leq j\leq n, a Chernoff bound (see e.g. ) gives that

We apply a union bound over all jj, and conclude that a.a.s., for every (sn/t)1−α/ω^≤j≤n(sn/t)^{1-\alpha}/\hat{\omega}\leq j\leq n,

In order to estimate the sums in the statement, we assume w.l.o.g. that T=[t]T=[t] and r1<r2<⋯<rtr_{1}<r_{2}<\cdots<r_{t} (otherwise we permute the indices of the vertices in TT). It follows that a.a.s., for every (sn/t)1−α/ω^≤j≤n(sn/t)^{1-\alpha}/\hat{\omega}\leq j\leq n,

For the lower bound on rir_{i} below, we need to use the fact that ⌊11+1/ω^in/t⌋≥(sn/t)1−α/ω^\left\lfloor\frac{1}{1+1/\hat{\omega}}in/t\right\rfloor\geq(sn/t)^{1-\alpha}/\hat{\omega}, which is easily verified to be true since ω^=o((sn/t)1−α)\hat{\omega}=o\left((sn/t)^{1-\alpha}\right). Finally, we infer that a.a.s., for any choice of SS,

This proves statement (b). For statement (a), take s=ts=t and note that for this choice of ss, for any ωnαlog⁡n≤t≤n\omega n^{\alpha}\log n\leq t\leq n, the condition s1−α≥ω(n/t)αlog⁡ns^{1-\alpha}\geq\omega(n/t)^{\alpha}\log n is satisfied. Observe that then S=T=[t]S=T=[t], so the first inequality in the above equation is an equality.∎

Fix a constant K>1−αpK>\frac{1-\alpha}{p}, and let DD be the set containing the first t=⌊Knα+βlog⁡n⌋t=\lfloor Kn^{\alpha+\beta}\log n\rfloor nodes added in the process. We will show that a.a.s. DD is a dominating set. By Lemma 1, we may condition on the event that (1) holds for t=∣D∣=⌊Knα+βlog⁡n⌋t=|D|=\lfloor Kn^{\alpha+\beta}\log n\rfloor. Note that this assumption on the ranks does not affect the distribution of the location of the nodes in [0,1)m[0,1)^{m}. Therefore, given a node u>tu>t (appearing in the process later than nodes in DD), the probability that uu is not dominated by DD is

Taking a union bound over all nodes not in DD, we can guarantee that a.a.s. all nodes are dominated.

As an alternative and relatively simple approach, one may prove the same upper bound on the domination number as follows. First, show that a.a.s. the minimum degree δ\delta is at least (1+o(1))pn1−α−β(1+o(1))pn^{1-\alpha-\beta}. Then we may use Theorem 1.2.2 in , which states that for every graph GG with minimum degree δ\delta,

0.2 Lower bound:

Given a constant 0<ε<10<\varepsilon<1, we define T′T^{\prime} to be the set of nodes in TT with rank greater than (1−ε)n(1-\varepsilon)n. Note that ∣T′∣|T^{\prime}| has a hypergeometric distribution, so it follows easily from Chernoff’s bound (see ) that a.a.s. ∣T′∣≥(ε/2)nα+β|T^{\prime}|\geq(\varepsilon/2)n^{\alpha+\beta}. For convenience, we choose ε=ε(α)\varepsilon=\varepsilon(\alpha) to be the only real in (0,1)(0,1) satisfying ε=2(1−ε)α\varepsilon=2(1-\varepsilon)^{\alpha}. For every node i∈T′i\in T^{\prime}, the corresponding ball BiB_{i} has length at most

We consider a tessellation of [0,1)2[0,1)^{2} into large cells. At the centre of each large cell we consider a smaller cell. Small cells have side length ((2/ε)n−α−β)1/m((2/\varepsilon)n^{-\alpha-\beta})^{1/m} and large ones have side length 2((2/ε)n−α−β)1/m2((2/\varepsilon)n^{-\alpha-\beta})^{1/m}. There are

large cells fully contained in [0,1)m[0,1)^{m} (we discard the rest), and thus NN small cells inside of those. By construction, if a node in T′T^{\prime} falls into a small cell, then its ball is contained into in the corresponding large cell. Let X\mathcal{X} be the set of small cells that contain at least one node in T′T^{\prime}, and let T′′⊆T′T^{\prime\prime}\subseteq T^{\prime} be a set of X=∣X∣X=|\mathcal{X}| nodes such that each cell in X\mathcal{X} contains precisely one node in T′′T^{\prime\prime} (if a given small cell contains at least two nodes in T′T^{\prime}, then a node is selected arbitrarily to be placed in T′′T^{\prime\prime}). Vertices in T′′T^{\prime\prime} are potentially dangerous for the adversary since, one node in D2D_{2} can “out-dominate” at most one single node in T′′T^{\prime\prime}. However, she may in theory get lucky and in-dominate many of these nodes (in the next section we will show that this will not happen a.a.s.).

We want to show that a.a.s. X≥N/4X\geq N/4. The probability that there are at least 3N/43N/4 small cells containing no nodes in T′T^{\prime} is at most

In-domination:

Let ξ=ξ(α)\xi=\xi(\alpha) be a sufficiently small positive constant, and define

The adversary chooses a set D1⊆TD_{1}\subseteq T of s=⌊ξμ−mnα+β⌋s=\lfloor\xi\mu^{-m}n^{\alpha+\beta}\rfloor nodes in her attempt to in-dominate T′T^{\prime}. By Lemma 1(b), a.a.s. regardless of her choice,

We tessellate the space into cells of volume (2/ε)n−α−β(2/\varepsilon)n^{-\alpha-\beta} (same size as the small cells in the out-domination part, but now we have the whole space partitioned into cells of that size). Recall that, for each node i∈D1i\in D_{1}, the ball BiB_{i} has length bi1/m≥n−(α+β)/m{b_{i}}^{1/m}\geq n^{-(\alpha+\beta)/m}. Therefore, the volume of the set of cells intersected by BiB_{i} is at most

Combining this and (5), a.a.s. and regardless of the adversary’s choice, the total volume of the cells intersected by the balls of the nodes in D1D_{1} is at most

Let Y\mathcal{Y} be the set of cells intersected by the balls of the nodes in D1D_{1}, and put Y=∣Y∣Y=|\mathcal{Y}|. By (6), a.a.s.

Thus, in view of (4) we just need to make ξ\xi small enough so that ∣X∖Y∣|\mathcal{X}\setminus\mathcal{Y}| is larger than ∣D2∣=⌊ξμ−mnα+β⌋|D_{2}|=\lfloor\xi\mu^{-m}n^{\alpha+\beta}\rfloor. That is because dangerous cells in X∖Y\mathcal{X}\setminus\mathcal{Y} contain nodes in T′T^{\prime} that are not in-dominated by D1D_{1}, and each one of these cells requires one different node in D2D_{2} to out-dominate its nodes. Recall that our choice of ε∈(0,1)\varepsilon\in(0,1) depends only on α\alpha. Then picking ξ\xi sufficiently small so that εξ1−α2(1−α)+ξ<ε8\frac{\varepsilon\xi^{1-\alpha}}{2(1-\alpha)}+\xi<\frac{\varepsilon}{8}, we get

where we also used that μ>3λ>λ\mu>3\lambda>\lambda. Finally, distinguishing the cases m=O(1)m=O(1) and m→∞m\to\infty, we observe that

Domination in Facebook 100 graphs

Several algorithms were used to bound the domination number of the FB100 graphs, but one providing the smallest dominating sets is an adaptation of the DS-DC algorithm . In the algorithm, initially all nodes VV are in the dominating set SS. It then selects a node uu of minimum degree in SS, and deletes it only if the set S∖{u}S\setminus\{u\} remains dominating. The algorithm then repeats these steps for all nodes in SS in order of their increasing degrees. We considered other algorithms, such as greedy algorithms where high degree nodes are added to an empty dominating set sequentially, or by choosing a random dominating set, but DS-DC outperformed these algorithms. We omit a detailed discussion of the performance of other algorithms owing to space.

Figure 1 presents the DS-DC predicted upper bounds on γ(G)\gamma(G), where GG is a graph in the FB100 data set.

We plotted the upper bound predicted by the MGEO-P model in Theorem 1.1, and we note the close similarity between that bound and ones for FB100. Note that we ignore constants in the big Oh term in the upper bound from the model, and simply plot the bound generated by nα+βlog⁡nn^{\alpha+\beta}\log n. The values for α\alpha, β\beta, and the dimension parameter mm for each of the FB100 graphs are taken from tables provided in . (For example, in order to determine the power-law exponent, the Clauset-Shalizi-Newman power law exponent estimator was used; see for more details.) The MGEO-P bound seems well-correlated with the bounds provided in the kk-core, especially where k=3,4,5.k=3,4,5. See Table 1, which fits the domination number of the FB100 graphs to the curve y=nxlog⁡ny=n^{x}\log n.

To contrast the bounds provided in Figure 1 with the bound in (3), we plot them in Figure 2. We plotted the theoretical bound using δ=5\delta=5 (that is, the minimum degree of the 55-core).

The figure shows a significant over-estimate of the domination number of the bound in (3), further corroborating the claim that the domination numbers of the FB100 graphs are sublinear with respect to the order of the graph.

Proof of Theorem 1.2

Given a bounded subset of the plane MM, for ε>0\varepsilon>0 let N(ε)N(\varepsilon) be the minimum number of balls of radius ε\varepsilon that can cover MM. Then we have that

where M‾\overline{M} denotes the closure of MM.

Observe that (C−1)(C-1) can therefore, be seen as measuring the proportion of unavoidable overlapping. Moreover, shows that an optimal covering of the square S\mathcal{S} using balls of radius ε\varepsilon corresponds to arranging the balls in such a way that their centers are the centers of the cells of a hexagonal tiling of length ε\varepsilon. More precisely, consider the lattice

Then the set of balls of radius ε\varepsilon and centre in Lε\mathcal{L}_{\varepsilon} that intersect S\mathcal{S} form a covering of S\mathcal{S} that gives the limit in Theorem 4.1.

Note that for all GG with maximum degree Δ\Delta, we trivially have γ(G)≥n/(1+Δ(G))\gamma(G)\geq n/(1+\Delta(G)) (for further relations between γ(G)\gamma(G) and other graph parameters, see, for example, ). Given any constant c>0c>0, for G∈G(n,r)G\in\mathscr{G}(n,r) with r≥clog⁡n/nr\geq c\sqrt{\log n/n}, it is easy to show, by Chernoff bounds together with union bounds, that a.a.s. Δ(G)=O(r2n)\Delta(G)=O(r^{2}n). Therefore, a.a.s. we have γ(G)=Ω(r−2)\gamma(G)=\Omega(r^{-2}). On the other hand, we can trivially construct a dominating set of G(n,r)\mathscr{G}(n,r) by tessellating S\mathcal{S} into square cells of side length r/2r/\sqrt{2} and picking one node from each cell (if the cell is not empty). This holds deterministically for any geometric graph (not necessarily random), with no restriction on rr, and gives γ(G)=O(r−2)\gamma(G)=O(r^{-2}). It follows that, for G∈G(n,r)G\in\mathscr{G}(n,r) with r≥clog⁡n/nr\geq c\sqrt{\log n/n}, a.a.s. γ(G)=Θ(r−2)\gamma(G)=\Theta(r^{-2}).

We first prove the lower bound in part (b). Fix an arbitrarily small constant δ>0\delta>0. Tessellate S\mathcal{S} into cells of side length α=ωlog⁡n/n=o(r)\alpha=\sqrt{\omega\log n/n}=o(r). By Chernoff bounds together with a union bound over all cells, we get that a.a.s. each cell contains at least one node. We may condition on this event, and proceed deterministically. For a contradiction, suppose that there exists a dominating set of size s=⌊(C/π−δ)r−2⌋.s=\lfloor(C/\pi-\delta)r^{-2}\rfloor. Consider then ss balls whose centers are at the nodes of the dominating set. When using radius rr for all balls, each cell is at least touched by some ball, since each cell is non-empty and each node is covered. Hence, by using radius r′=r+α2r^{\prime}=r+\alpha\sqrt{2}, each square is totally covered by some ball. Therefore, S\mathcal{S} can be covered by ⌊(C/π−δ)r−2⌋\lfloor(C/\pi-\delta)r^{-2}\rfloor balls of radius r′r^{\prime}. On the other hand,

since α=o(r)\alpha=o(r). This contradicts Theorem 4.1 and so γ(G)>s\gamma(G)>s. Since the argument holds for any δ>0\delta>0, we get the desired lower bound.

For the upper bound, we will show that we can find a covering of S\mathcal{S} with (C/π+o(1))r−2(C/\pi+o(1))r^{-2} balls of radius rr that are centered at some nodes of GG. Again, fix some arbitrarily small constant δ>0\delta>0. Let r′=(1−δ)rr^{\prime}=(1-\delta)r, and consider the lattice Lr′\mathcal{L}_{r^{\prime}}, as defined in (7). Let Lr′′\mathcal{L}^{\prime}_{r^{\prime}} be the set of all points x∈Lr′x\in\mathcal{L}_{r^{\prime}} such that the ball with centre xx and radius r′r^{\prime} intersects S\mathcal{S}. Recall that Lr′′\mathcal{L}^{\prime}_{r^{\prime}} gives the optimal covering of S\mathcal{S} with balls of radius r′r^{\prime}, and therefore, attains the bound given by Theorem 4.1

It might happen that some point x∈Lr′′x\in\mathcal{L}^{\prime}_{r^{\prime}} does not belong to S\mathcal{S}. In this case, we replace xx by the closest point x^\hat{x} on the boundary of S\mathcal{S} (this can be uniquely done, since S\mathcal{S} is closed and convex). Note that B(x,r′)∩S⊆B(x^,r′)∩S\mathcal{B}(x,r^{\prime})\cap\mathcal{S}\subseteq\mathcal{B}(\hat{x},r^{\prime})\cap\mathcal{S}. We denote ^Lr′\hat{}\mathcal{L}_{r^{\prime}} the modified set of points that we obtained. By construction, ^Lr′⊆S\hat{}\mathcal{L}_{r^{\prime}}\subseteq\mathcal{S}, and we can cover S\mathcal{S} using balls with centre in ^Lr′\hat{}\mathcal{L}_{r^{\prime}} and radius r′r^{\prime} or larger. Moreover, ∣^Lr′∣=s|\hat{}\mathcal{L}_{r^{\prime}}|=s. Clearly, if we can guarantee that for each x∈^Lr′x\in\hat{}\mathcal{L}_{r^{\prime}} there exists a node of GG inside B(x,δr)∩S\mathcal{B}(x,\delta r)\cap\mathcal{S}, then GG is dominated by these nodes, and hence, ss yields an upper bound for γ(G)\gamma(G).

Observe that for any point x∈^Lr′x\in\hat{}\mathcal{L}_{r^{\prime}} (and therefore, in S\mathcal{S}), the area of B(x,δr)∩S\mathcal{B}(x,\delta r)\cap\mathcal{S} is at least (δr)2π/4(\delta r)^{2}\pi/4, since at least a quarter of a ball must be inside S\mathcal{S}. The probability that there is no node of GG in B(x,δr)∩S\mathcal{B}(x,\delta r)\cap\mathcal{S} is at most

Since there are ss events that we need to investigate and clearly s≤ns\leq n, by a union bound, a.a.s., for every x∈^Lr′x\in\hat{}\mathcal{L}_{r^{\prime}}, the region B(x,δr)∩S\mathcal{B}(x,\delta r)\cap\mathcal{S} contains at least one node of GG. It follows that a.a.s. γ(G)≤s\gamma(G)\leq s and since the argument holds for any δ>0\delta>0, we derive the desired upper bound.

For the proof of part (a), note that N(x)N(x) is non-decreasing function of xx, and N(x)=1N(x)=1 for x≥1/2x\geq 1/\sqrt{2}. Fix r=Θ(1)r=\Theta(1). Tessellate S\mathcal{S} into cells of side length α=(ω/2)log⁡n/n\alpha=\sqrt{(\omega/2)\log n/n}. For the lower bound, suppose for contradiction that γ(G)≤N(r+α2)−1\gamma(G)\leq N(r+\alpha\sqrt{2})-1. By Chernoff bounds together with a union bound over all cells, a.a.s. there is at least one node in each such cell. Now place N(r+α2)−1N(r+\alpha\sqrt{2})-1 many balls with centers at the nodes of the dominating set. Since by using radius rr each cell is at least touched by some ball, by using radius r+α2r+\alpha\sqrt{2} each cell is totally covered by a ball, and therefore S\mathcal{S} is covered by N(r+α2)−1N(r+\alpha\sqrt{2})-1 balls of radius r+α2r+\alpha\sqrt{2}, contradicting the definition of N(x)N(x). Therefore, a.a.s. γ(G)≥N(r+α2)\gamma(G)\geq N(r+\alpha\sqrt{2}).

For the upper bound, consider an optimal arrangement of N(r−β)N(r-\beta) balls of radius r−βr-\beta, where β=ω/n\beta=\omega/\sqrt{n}. As before, if the centre pp of a ball is outside S\mathcal{S}, but B(p,r−β)∩S≠∅\mathcal{B}(p,r-\beta)\cap\mathcal{S}\neq\emptyset, we may shift the centre of the ball towards its closest point p′p^{\prime} on the boundary of S\mathcal{S}. Since B(p,r−β)∩S⊆B(p′,r−β)∩S\mathcal{B}(p,r-\beta)\cap\mathcal{S}\subseteq\mathcal{B}(p^{\prime},r-\beta)\cap\mathcal{S}, we still preserve the covering property, and therefore, we can obtain an optimal covering of S\mathcal{S} with balls of radius r−βr-\beta and centered at points inside of S\mathcal{S}. As in part (a), it suffices to show the existence of a node v∈Vv\in V inside B(c,β)∩S\mathcal{B}(c,\beta)\cap\mathcal{S} for any centre cc in this optimal arrangement of balls. Since N(r−β)=O(1)N(r-\beta)=O(1), the probability that there exists a centre cc such that (B(c,β)∩S)∩V=∅(\mathcal{B}(c,\beta)\cap\mathcal{S})\cap V=\emptyset is at most

and hence, a.a.s. for all centers cc, we have that B(c,β)∩S\mathcal{B}(c,\beta)\cap\mathcal{S} contains at least one node of GG. These nodes form a dominating set, and so a.a.s. γ(G)≤N(r−β)\gamma(G)\leq N(r-\beta).

Finally, note that the lower bound in part (b) can be easily adopted to show that a.a.s. γ(G)=Ω(r−2)\gamma(G)=\Omega(r^{-2}) as a.a.s. a positive fraction of cells contain at least one node for the range of rr considered in part (c). As already mentioned, the upper bound of O(r−2)O(r^{-2}) holds for any (deterministic) geometric graph and any rr. Hence, part (c) follows. For part (d), the upper bound is trivial. The lower bound comes from the fact that a.a.s. there will be Θ(n)\Theta(n) isolated nodes, and a dominating set has to contain all of them. The proof of the theorem is finished. ∎

Conclusions and open problems

We considered the domination number of a stochastic model for OSNs, the MGEO-P model. Theorem 1.2 shows a sublinear bound on the domination number of OSNs, which is well correlated with estimates for the domination number taken for the Facebook 100 data set. In addition, we provided bounds for the domination number of random geometric graphs.

In future work, we would like to broaden our analysis of the domination number to other data sets, and to test larger samples of OSNs. We will contrast the estimates provided by other heuristic algorithms for computing minimum order dominating sets, and provide a fitting of the data to bounds provided by the model.

So-called “elites”, those who exert strong influence on the ambient network, are studied extensively in the sociology literature (see for an overview of the literature on this topic). One approach to detecting elites is via their relatively high degree; hence, the use of kk-cores in . A different approach to detecting elites is to search for them within a minimum order dominating set, as these sets reach the entire network. Further, if minimum order dominating sets have much smaller order than the network (as we postulate), then that reduces the computational costs of finding elites. We plan on considering this approach to finding elites via dominating sets in future work.

References