July 10, 2026
Real-world networks like the internet share patterns like a power law degree distribution and a high clustering coefficient. Many of these properties are captured by the generative model of hyperbolic random graphs (HRGs), which provides a theoretical framework for studying such networks. Motivated by the observation that several algorithms perform better on real-world networks than their worst-case guarantees suggest, we design and analyse distributed algorithms under the assumption that the input graph is an HRG. Indeed, prior work has shown that the classical symmetry-breaking problem of \(\Delta+1\) colouring, where \(\Delta\) is the maximum degree of the graph, can be solved in 2 rounds on HRGs [Maus and Ruff; SODA’26].
In stark contrast to this 2-round algorithm for \(\Delta+1\) colouring, we prove that the related symmetry-breaking problems of maximal independent set (MIS) and maximal matching (MM) are substantially harder: we establish a lower bound of \(\Omega\left(\frac{\log\log n}{\log\log\log n}\right)\) for MIS and MM on HRGs. Our lower bound techniques rely on new structural insights that may be of independent interest: we show that HRGs contain \(d\)-ary trees with large height and degree which enables us to adapt and lift prior impossibility results for distributed algorithms to the setting of HRGs.
We also show that these lower bounds are polynomial tight: we design algorithms tailored to HRGs that solve MIS and MM in \(\tilde{\mathcal{O}}(\log^{5/3}\log n)\) rounds with high probability in the \(\mathsf{LOCAL}\) model, improving over the general worst-case lower bound of \(\Omega\left(\min\left\{\log \Delta, \sqrt{\log n}\right\}\right)\) rounds [Khoury and Schild; FOCS’25].
Finally, we show that access to geometric information can significantly reduce the complexity of maximal matching: if vertices know their respective geometric embeddings, then MM can be solved in \(\mathcal{O}(\log\log\log n)\) rounds on HRGs. This reveals a separation between the standard \(\mathsf{LOCAL}\) model and an embedding-aware variant on hyperbolic random graphs.
Symmetry breaking problems play a key role in the theory of distributed computing. In this paper, we study the fundamental symmetry breaking problems maximal independent set (MIS) and maximal matching (MM). A maximal independent set is an independent set that cannot be extended by adding any vertex, while a maximal matching is a matching where no additional edge can be added. We study these problems through the lenses of the classic \(\mathsf{LOCAL}\) model and the \(\mathsf{CONGEST}\) model [1], [2]: A communication network is represented by a graph \(G=(V,E)\) with \(|V|=n\). In synchronous rounds, vertices exchange messages with their neighbours, and the round complexity of an algorithm is the number of rounds until all vertices have produced their output. We consider both the \(\mathsf{LOCAL}\) model, where messages have unbounded size, and the \(\mathsf{CONGEST}\) model, where each message is limited to \(\mathcal{O}(\log n)\) bits.
One of the classic results from the 1980’s is Luby’s Algorithm, giving a randomised algorithm that finds an MIS within \(\mathcal{O}(\log n)\) rounds with high probability [3], [4]. The same approach can then also be extended to an MM. Remarkably, despite much progress on distributed computing since then, \(\mathcal{O}(\log n)\) remains the best-known round complexity as a function of \(n\) for general graphs [5]–[19]. Meanwhile, Khoury and Schild recently showed a lower bound of \(\Omega(\sqrt{\log n})\) for MIS and MM on worst-case general graphs [9], improving upon an earlier bound of \(\Omega(\sqrt{\log n/\log\log n})\) rounds by Kuhn, Moscibroda, and Wattenhofer [20].
Worst-case graphs, however, may be a poor proxy for the communication networks that motivate many distributed problems. This raises the question of whether symmetry breaking becomes easier when the communication network is drawn from a model that captures structural features commonly observed in real-world networks. A recent step in this direction was taken by Maus and Ruff, who showed that a distributed colouring problem requiring \(\Omega(\log^* n)\) rounds on worst-case graphs [21], [22] can be solved in just two rounds on hyperbolic random graphs, a common abstraction of real-world networks [23]. Their simple algorithm significantly outperforms the state of the art for general graphs. Their work raises the broader question of whether this dramatic speed-up is specific to colouring, or whether it reflects a more general phenomenon. Thus, in this paper, we further explore the following research question:
How does the structure of real-world networks affect the complexity of symmetry breaking?
While different real-world networks have different structures, e.g., the internet graph has a different topology than a road network, many real-world networks have been observed to share several topological features. For example, many networks exhibit a heterogeneous degree distribution, close to a power-law.3 Examples include, among others, social networks [24], communication networks like the internet [25], and biological networks [26]. Another property observed in real-world networks is a high clustering coefficient [27]–[29]. We follow the lines drawn in [30]–[33] using hyperbolic geometry to encapsulate both of these properties. In particular, we use the model of hyperbolic random graph (HRG) as introduced by Krioukov, Papadopoulos, Kitsak, Vahdat, and Boguñá [31]. At this point, HRGs are well established as a model for complex real-world networks, and many different centralised algorithmic problems have been studied [34]–[40]. Moreover, the excellent empirical work by Bläsius and Fischbeck suggests that the theoretical analysis of algorithms on HRGs provides a good predictor for performance on real-world networks [41], and HRGs also serve as the underlying model in the work of Maus and Ruff [23]. A threshold hyperbolic random graph is obtained by sampling \(n\) vertices into a hyperbolic disk with radius \(R \approx 2\log n\), and a pair of vertices share an edge if the hyperbolic distance between the two vertices is at most \(R\). We provide a formal definition in 4 and refer the interested reader to [42], [43] for an in-depth introduction to HRGs.
We show that the complexity landscape of MIS and MM on HRGs differs dramatically from that of the distributed graph colouring problem studied in [23]; see also 1 for visualisiation of our results: On the algorithmic side, we show that MIS and MM can be solved exponentially faster on HRGs than on worst-case graphs by giving \(\log^{\mathcal{O}(1)}\log n\)-round algorithms ([thm:mainpolylog]). But, on the lower-bound side, we show that this speed-up has inherent limits. In contrast to graph colouring, MIS and MM do not collapse to constant time on HRGs; in fact, our lower bounds show that any further improvement to our upper bounds is limited to a small polynomial improvement ([thm:lowerbound]). This impossibility result relies on new structural insights into HRGs: we prove that their giant component contains substructures (trees) on which faster distributed algorithms are impossible ([pro:d-regular-trees]). We believe the existence of these substructures is of independent interest, as it contributes to the structural understanding of HRGs as models of real-world networks.
The lower bound above applies in the standard \(\mathsf{LOCAL}\) model, where nodes do not know their position in the underlying hyperbolic geometry. This naturally raises the question of whether such geometric information can be exploited algorithmically. In the classic \(\mathsf{LOCAL}\) model, the initial knowledge of nodes is limited to their own ID and potentially some coarse upper bounds on global parameters such as the maximum degree or the number of nodes. However, in real-world networks, nodes may have additional information about the network that may be exploited algorithmically. For maximal matching, we show that such extra information can indeed be useful: we design an algorithm that beats the aforementioned lower bound when nodes have access to their geometric HRG coordinates, giving a doubly-exponential improvement over the complexity on general graphs ([thm:geometric]).
We next present our results in detail and provide further background.
We obtain the following algorithms for MIS and MM on HRGs that are exponentially faster than the best algorithms for general graphs [3], [6], [7].
theoremmaintheorem There are randomised distributed algorithms that, for threshold hyperbolic random graphs, compute Maximal Independent Set and Maximal Matching
All probabilistic statements in our theorems are with respect to both the random generation of the HRG and the random choices made by the algorithm. See 1 for a visual representation of our results in comparison to runtimes of existing work. All algorithms of [thm:mainpolylog] use a two-step approach. For MIS, the first step consists of two iterations of a Luby-style algorithm in which certain nodes activate themselves, draw a random number, and nodes whose drawn numbers are local minima join the independent set. As a result, we obtain that the graph shatters into exponentially smaller remaining unsolved connected components, on each of which we can then run the best deterministic algorithm exploiting the component’s small size, namely the algorithm of Ghaffari and Grunau for MIS/MM in the \(\mathsf{LOCAL}\) model [13] and the fastest deterministic algorithms in the \(\mathsf{CONGEST}\) model [44]. Such a two-step approach is usually referred to as the shattering technique and has been used in essentially all state-of-the-art symmetry breaking algorithms, e.g., [6], [7], [45], and one also knows that the randomised complexity is lower bounded by the deterministic complexity on exponentially smaller instances [46]. A central difference in our work and the main contribution of [thm:mainpolylog] is that we use a fundamentally different argument than all these existing works to show that such a shattering phenomenon emerges. This is interesting as recently discovered flaws and missing details in existing shattering approaches ask for more different shattering approaches [47].
One might hope that ordinary Luby iterations already yield such a shattering. We show that this is not the case: for any constant number of Luby iterations, both large components and large degrees remain; see 8 for details.
[thm:mainpolylog] is not only faster than the known fastest algorithms for MIS and MM but also exponentially faster than the lower bound of \(\Omega\left(\min\{\log \Delta, \sqrt{\log n}\}\right)\) by Khoury and Schild [9] for general graphs, since the maximum degree \(\Delta\) of an HRG is polynomial in \(n\). Combined with the fact that the other central symmetry-breaking problem of \(\Delta+1\)-colouring, and even colourings with fewer colours, can be solved in constant time [23], suggests that the algorithms in [thm:mainpolylog] may be far from optimal. We show that this is not the case. Any improvements to [thm:mainpolylog] can only be polynomial in its runtime. More precisely, we prove the following theorem, providing a strong separation between the runtime of MIS/MM and colouring problems on HRGs.
theoremtheoremlowerbound Asymptotically almost surely6, the randomised complexity7 for computing an MIS or an MM for the giant component of a hyperbolic random graph is \(\Omega\left(\frac{\log\log n}{\log\log\log n}\right)\) rounds in the \(\mathsf{LOCAL}\) model.
To prove [thm:lowerbound] we show that with high probability any HRG contains relatively large induced \(d\)-ary trees8 and then leverage the fact that distributed algorithms for MIS/MM cannot be fast on tree-like graphs [9], [20], [48]–[52].
::: {#thm:trees .theorem} Theorem 1 (Informal version of [pro:d-regular-trees]). In an HRG, asymptotically almost surely, there exist polynomially many induced \(d\)-ary trees with \(n' \approx \sqrt{\log n}\) vertices, any degree \(d\), depth \(h \approx \log_d\log n\) and where the root of each tree has one additional edge to a vertex that connects it to the giant component.* :::
Then [thm:mainpolylog] implies that [pro:d-regular-trees] is in some sense almost-tight, i.e., significantly larger trees cannot exist. For example, \(d\)-ary trees with degree \(d \approx e^{\log^2\log n}\) and height \(h \approx \log^2\log n\) cannot exist, as the resulting lower bound for computing MIS on HRGs would contradict the existence of the algorithm in [thm:mainpolylog]. In fact, we consider [pro:d-regular-trees] as a main technical contribution of our work that is of independent interest. More generally, these structural insights contribute to existing literature of analysing the structure of HRGs, like the use of “dangling paths” to obtain bounds on the cover and hitting times of random walks [39], bounds on the degeneracy to obtain efficient colouring approximations [40], the construction of multi-commodity flows to bound the spectral gap [53] or more broadly the connectivity of an HRG [54]–[59].
The lower bounds of [thm:lowerbound] apply in the standard distributed setting, where nodes only interact through the graph and have no access to the underlying hyperbolic embedding. This raises the question of whether this hidden structure can be exploited when it is made available to the algorithm. We therefore also study a coordinate-aware variant of the model, in which each node initially knows its position in the hyperbolic disk and can communicate this information to its neighbours in one round. This is a strong assumption, but geometric information of this kind is available or can be inferred in several network settings [32]. Similar assumptions have also been considered for distributed MIS in Euclidean unit disk graphs, where they lead to deterministic \(\mathcal{O}(1)\)-round algorithms in the \(\mathsf{CONGEST}\) model [60] and to efficient maximum independent set approximations for graphs embedded in the hyperbolic plane [36], [61].
Under this additional geometric information, the complexity landscape changes again. For maximal matching, we break the lower-bound barrier of [thm:lowerbound] and obtain a randomised \(\mathcal{O}(\log\log\log n)\)-round algorithm in the \(\mathsf{CONGEST}\) model. This shows, similar to quantum algorithms which yield an advantage for some distributed problems [62], [63], that one can beat lower bounds if additional computational resources like that of a geometric embedding are available. For maximal independent set, we obtain a \(\mathcal{O}(\log\log n)\)-round \(\mathsf{CONGEST}\) algorithm on HRGs with known geometric coordinates.
theoremtheoremgeometric There are distributed algorithms that, for threshold hyperbolic random graphs given in their geometric representation, compute
Maximal Independent Set in \({\mathcal{O}}(\log\log n)\) rounds w.h.p.9 of the \(\mathsf{CONGEST}\) model;
Maximal Matching in \({\mathcal{O}}(\log\log\log n)\) rounds w.h.p. of the \(\mathsf{CONGEST}\) model.
While we do not break the lower bound barrier of [thm:lowerbound] for MIS, we obtain a runtime of a “flat” \(\log\log n\) for a deterministic algorithm. Moreover, the techniques do not rely on shattering as does [thm:mainpolylog] and, as such, provide a different approach for solving symmetry-breaking problems on real-world networks.
Finally, we remark that [thm:lowerbound] together with [thm:geometric] imply that a distributed computation of the embedding for an HRG requires \(\Omega\left(\frac{\log\log n}{\log\log\log n}\right)\) rounds.
A hyperbolic random graph combines the non-vanishing clustering coefficient of random geometric graphs (RGGs) [64] and the power-law degree distribution exhibited by models like Chung–Lu graphs [65], [66] or the Barabási–Albert model [67]. Other properties like that of a giant component10 [54], [55], [68], [69] and the small-world phenomena11 [57]–[59], [70], [71] follow. From the algorithmic side problems like Shortest Path [34], Maximum Clique [35], Maximum Independent Set and Maximum Matching [36], Vertex Cover [37], [38], Random Walks [39] or Colouring [23], [40] have been studied.
A closely related model which captures the same properties is Geometric Inhomogeneous Random Graphs (GIRGs) introduced by Bringmann, Keusch and Lengler [72]. While HRGs and GIRGs are often regarded as roughly equivalent [73], several differences between the two models have been identified [40], [74], [75]. GIRGs have been used to analyse algorithms like routing [76], shortest paths [77], rumour spreading [78], community testing [79], geometry detection [80] or first passage percolation [81].
Another related model is the class of hyperbolic uniform disk graphs (HUDGs) [61], [82], [83]; previously pioneered by Kisfaludi-Bak [84]. In this model, vertices are placed arbitrarily in a hyperbolic disk of radius \(R\), and two vertices are adjacent whenever their hyperbolic distance is at most \(R\). Unlike hyperbolic random graphs, no probability distribution is used to generate the vertex positions. When \(R=2\log n\), the model can be viewed as a deterministic or worst-case analogue of a hyperbolic random graph: the same geometric threshold is used, but the vertex positions are chosen adversarially rather than sampled from the hyperbolic random graph distribution.
While Luby’s algorithm with \(\mathcal{O}(\log n)\) rounds remains the fastest known algorithm on general graphs, substantial progress has been made on MIS and MM along other directions. From the algorithmic point of view, deterministic algorithms have been developed where the current state of the art for the \(\mathsf{LOCAL}\) model is \(\tilde{\mathcal{O}}(\log^{5/3} n)\) due to Ghaffari and Grunau [13], and in \(\mathsf{CONGEST}\) \(\mathcal{O}(\log^2 \Delta\cdot \log\log \Delta \cdot \log n)\) rounds for MIS [44] and \(\mathcal{O}(\log^2\Delta \cdot \log n)\) for MM [85]. Moreover, Barenboim, Elkin, Pettie and Schneider introduced a randomised algorithm that runs faster on graphs with small maximum degree \(\Delta\) [6]. In his seminal paper, Ghaffari, improved on their MIS algorithm and accomplished a round complexity of \(\mathcal{O}\left(\log \Delta + \operatorname{\text{{\rm poly}}}\log\log n \right)\) in the \(\mathsf{LOCAL}\) model [7], [13], [19] and \(\mathcal{O}\left(\log \Delta\log\log n + \operatorname{\text{{\rm poly}}}\log\log n \right)\) in the \(\mathsf{CONGEST}\) model [86]. Recently, it was shown by Khoury and Schild that these algorithms solve MIS on trees in \(o\left(\sqrt{\log n}\right)\) rounds [10] while MM on trees requires \(\Theta(\sqrt{\log n})\) rounds [6], [7], [9]. This shows that MM is harder than MIS on trees, while on general graphs, MIS is just as hard as MM.12
Within the broader context of lower bounds, due to a classic construction by Kuhn, Moscibroda and Wattenhofer there is an \(\Omega\left(\min\left\{\frac{\log \Delta}{\log\log \Delta}, \sqrt{\frac{\log n}{\log\log n}} \right\}\right)\) lower bound for randomised MIS and MM on certain irregular trees [20] (see also [48] for a simplified version). For Maximal matching, the breakthrough work by Balliu, Brandt, Hirvonen, Olivetti, Rabie and Suomela introduced an \(\Omega\left(\min\left\{\Delta, {\frac{\log n}{\log\log n}} \right\}\right)\) deterministic lower bound on regular trees13 [49]. Recently Khoury and Schild improved the randomised lower bound to \(\Omega\left(\min\left\{\log \Delta, \sqrt{\log n} \right\}\right)\) for randomised MM on regular trees [9]. For MIS on trees, Balliu, Ghaffari, Kuhn and Olivetti established a lower bound of \(\Omega\left(\min\left\{\frac{\log \Delta}{\log\log \Delta}, \sqrt{\frac{\log n}{\log\log n}} \right\}\right)\) for randomised algorithms [51]. As previously mentioned, MIS can be solved in \(o(\sqrt{\log n})\) on trees, and thus, it is not a coincidence that this bound is weaker than that for MM on trees.
For more specialised graph classes, further algorithms have been developed. If the input graph is a growth-bounded graph, Kuhn, Moscibroda, Nieberg and Wattenhofer showed that MIS can be solved in deterministic \(\mathcal{O}(\log \Delta \cdot \log^* n)\) rounds in \(\mathsf{LOCAL}\)[87]. Moreover, Barenboim and Elkin developed a \(o(\sqrt{\log n})\) rounds algorithm for MIS in the \(\mathsf{LOCAL}\) model if the input graph has bounded arboricity [11], [88]. MIS on bounded-independence graphs requires \(\mathcal{O}(\log^* n)\) rounds for the \(\mathsf{LOCAL}\) model [18] and for unit disk graphs, MIS can even be solved in constant rounds by a deterministic \(\mathsf{CONGEST}\) algorithm [60].
The distinguishing feature of hyperbolic random graphs as compared to general graphs is their underlying hyperbolic geometry that causes two nodes to be adjacent if and only if they are also geometrically close. Thus, our results are based on three geometric mechanisms. First, for the upper bounds, we use short randomised procedures to create geometric separators in the hyperbolic disk; these separators confine every remaining connected component to a narrow angular sector and hence to polylogarithmic size. Second, for the lower bounds, we show that HRGs contain many induced regular trees attached to the giant component by a cut edge, allowing lower bounds from tree-like instances to transfer to HRGs. Third, when geometric coordinates are available, we replace shattering by an explicit tiling of the disk, which lets nodes process well-separated tiles or annuli in parallel. We discuss these three ideas in turn.
The algorithms for MIS and MM in [thm:mainpolylog] follow the same high-level strategy. We first run a constant-round randomised procedure that shatters the graph into connected components of size \(\mathop{\mathrm{polylog}}n\). We then solve each remaining component independently using the best available deterministic algorithms. In the \(\mathsf{LOCAL}\) model, this gives the \(\tilde{\mathcal{O}}(\log^{5/3}\log n)\)-round bound via the algorithm of Ghaffari and Grunau [13], and in the \(\mathsf{CONGEST}\) model, it gives the \(\tilde{\mathcal{O}}(\log^3\log n)\)-round bound via the deterministic \(\mathsf{CONGEST}\) algorithms of [44].
Why don’t we use off-the-shelf algorithms and their analysis? Essentially all known sublogarithmic-time distributed algorithms for MIS and MM rely on such a shattering framework: an initial randomised step solves most of the graph, leaving only small connected components to be handled deterministically. Existing shattering-based algorithms are not directly suitable for our setting for two reasons. First, their randomised step inherently takes \(\Theta(\log \Delta)\) rounds, which is too slow on HRGs, where the maximum degree \(\Delta\) is polynomial in \(n\). Second, in reality, the standard shattering guarantees do not leave components of \(\mathop{\mathrm{polylog}}n\) size; rather, the residual components may have size \(\operatorname{\text{{\rm poly}}}(\Delta)\log n\), together with additional structure that can be exploited algorithmically. Relying on this additional structure would lead to weaker final round complexities in our setting. We therefore develop a different geometric argument that establishes the desired shattering phenomenon directly.
The key geometric target of the randomised step is the separator pattern shown in 2. We want to remove all vertices in an inner disk (hatched area in 2) and, in the remaining outer annulus, remove a collection of neighbourhoods (green areas in 2) that cut the annulus into narrow angular sectors. If these remaining unsolved sectors have angular width at most \(\phi\), for \(\phi \approx \mathop{\mathrm{polylog}}(n)/n\), then each sector contains only \(\mathop{\mathrm{polylog}}n\) vertices with high probability (blue regions in 2b). Moreover, by the choice of the removed regions, any two vertices lying in different sectors are at hyperbolic distance larger than \(R\), and hence are non-adjacent. Therefore, once this separator pattern is realised, every remaining connected component is contained in a single narrow sector and has size \(\mathop{\mathrm{polylog}}n\).
Figure 2: Shattering: (a) Vertices in the blue and red annuli are active in steps 1 and 2, respectively. (b) The hatched area and the green area together form the separator.. a — Overview of first two rounds., b — Result after the first two rounds.
The main challenge is to realise this separator pattern while producing a locally valid partial solution for MIS or MM. The way this is done differs substantially between the two problems.
The main challenge in designing a fast shattering algorithm, ideally a constant-time algorithm, is dealing with the huge amount of dependencies that occur in algorithms for HRGs. On the positive side for MIS, a selected vertex removes all its neighbours, which makes it well suited for creating geometric separators: if we can make vertices in suitable annuli join the independent set, then their neighbourhoods carve out the desired removed regions, see the green disks in 2a.
The first step creates separators in the outer disk. We activate vertices in an annulus at radius roughly \(\approx R-8\log\log n\) (see red annulus in 2a), equivalently, vertices whose degrees lie in an appropriate polylogarithmic range. We then perform one Luby-style step on these active vertices: each active vertex draws a random value, and active vertices that are local winners join the independent set. The neighbourhoods of these winners remove angular intervals in the outer disk. With high probability, every angular interval of width about \(\phi \approx \mathop{\mathrm{polylog}}n /n\) contains such a removed neighbourhood. Intuitively, this stems from the fact that active vertices have polylogarithmic degree and thus, given that within an interval of width \(\phi\) there \(\mathop{\mathrm{polylog}}n\) vertices located, in expectation there will be also a winner vertex. Consequently, the outer annulus is split into angular sectors of width at most \(\phi\) (see 2b).
The second step removes the inner disk (hatched region in 2). Here we activate remaining low-degree vertices close to the boundary (blue annulus in 2b). A vertex \(u\) in the inner disk has many such active neighbours close to the boundary (that also survived the first step), and for the sake of analysis we can determine \(\omega(\log )\) such neighbours for which the events of joining the independent set are independent. Each of these active neighbours has a constant probability of joining the independent set and thereby removing \(u\). Since \(u\) has \(\omega(\log n)\) such independent chances, \(u\) remains uncovered with probability \(n^{-\omega(1)}\). A union bound then shows that the entire inner disk is removed with high probability.
After these two steps, the remaining vertices lie only in the narrow sectors of the outer annulus. As argued above, different sectors are disconnected, and each sector contains at most \(\mathop{\mathrm{polylog}}n\) vertices with high probability. We then solve MIS independently on each remaining component using the deterministic algorithm for the corresponding models of \(\mathsf{LOCAL}\) and \(\mathsf{CONGEST}\).
For maximal matching, the same separator pattern is needed, but it is harder to realise. In MIS, one selected vertex removes its whole neighbourhood. In MM, a matched edge removes only its two endpoints, and a vertex can be matched only once. Thus, removing an entire region requires assigning distinct matching partners without creating conflicts.
We first remove the inner disk (2 hatched area). Vertices close to the boundary and of sufficiently small degree are activated (blue annulus in 2a). Each active vertex samples one neighbour uniformly at random and proposes to it. We then consider the vertices in the inner disk, which have a sufficiently large degree. With high probability, every such vertex receives at least one proposal from an active neighbour. Each vertex in the inner disk accepts one proposal, and these accepted proposals form matching edges. This matches all vertices in the inner disk and removes them from further consideration.
It remains to create separators in the outer annulus. As in the MIS case, neighbourhoods of vertices in a suitable annulus have the right geometry to separate the outer disk into narrow angular sectors. We activate such vertices (vertices in the red annulus in 2a) and let them draw random values and select local maxima as in one iteration of Luby. By choosing the annulus carefully, we ensure that the relevant outer-disk portions of the \(2\)-hop neighbourhoods of activated nodes are pairwise disjoint. Thus, in the \(\mathsf{LOCAL}\) model, each activated vertex can learn its \(2\)-hop neighbourhood and compute a maximal matching. Removing the matched nodes from the graph disconnects the unsolved parts into small connected components (blue areas in 2b). In the \(\mathsf{CONGEST}\) model, \(2\)-hop neighbourhoods cannot be learned efficiently; hence, we use additional structural information to compute the desired maximal matching efficiently.
The lower bounds are based on showing that HRGs contain hard substructures aka trees. The goal is not merely to find trees as abstract subgraphs, but to find induced regular trees that are attached to the giant component in a controlled way. More precisely, we show that, asymptotically almost surely, the giant component contains polynomially many induced \(d\)-ary trees (trees with degree \(d\)) of depth \(h\). In particular, we use \(d\)-ary trees with parameters \[\begin{align} d \approx {\log\log n} \qquad\text{and}\qquad h \approx \log \log n / \log\log\log n. \end{align}\] Each such tree is attached to the giant component by a cut edge incident to its root, while all other vertices of the tree have no additional edges leaving the tree. Once such trees exist, known lower bounds for MIS and MM on regular trees can be extended and lifted to the stated lower bounds for HRGs.
The trees are located close to the boundary of the hyperbolic disk; see 3 for an illustration. We consider a narrow angular sector of width \(\phi = o\left({\log n}/{n}\right)\). Inside such a sector, we prescribe small regions, called boxes (see yellow area in 3a for a sketch), where the vertices of the tree should appear. The root is placed at a distance about \(\log\log n\) from the boundary (see vertex \(u\) in 3a). Its children are placed slightly closer to the boundary, at a radial offset of about \(2\log d\). This construction is then repeated recursively \(h-1\) times: a vertex at level (i) has its \(d\) children in boxes at the next radial level.
The boxes are positioned so that three geometric conditions hold. First, each parent is adjacent to every vertex in its child boxes, i.e., within hyperbolic distance at most \(R\). Second, vertices in child boxes of the same parent are not adjacent to one another. Third, boxes belonging to different branches are sufficiently separated so that no unintended edges are created. These conditions follow from the hyperbolic distance formula and the choice of the radial and angular spacing between boxes.
Thus, if every prescribed box contains exactly one vertex and the remainder of the sector is empty, then the graph induced by the sector is precisely the desired \(d\)-ary tree of depth \(h\). We call such a sector nice; later, we also ensure that the tree in a nice sector is also connected to the giant component in the desired way.
We show that a single sector is nice with probability \(n^{-o(1)}\): to give intuition for why this is, we first remark that each of the \(n' \approx \sqrt{\log n}\) prescribed boxes in a sector has at least constant expected occupancy. So we categorise two different types of boxes, namely
boxes where the expected number of vertices is constant, i.e., \(\lambda = \Theta(1)\) and;
boxes where the expected number of vertices is larger than constant, i.e., \(\lambda = \omega(1)\).
For [item:constantbox], observe that the number of vertices \(X_i\) in the \(i\)-th box is known to be exactly distributed as a Poisson distribution, i.e., \(P(X_i=k)=\lambda^k/k! \cdot e^{-\lambda}\). Here, \(\lambda\) is the expected number of vertices in the box and scales with the size of the box. For boxes where \(\lambda =\Theta(1)\) we obtain \(P(X_i=1)=\lambda\cdot e^{-\lambda}=p=\Omega(1)\). Thus, since we have at most \(n' \approx \sqrt{\log n}\) boxes, the probability that all these boxes have exactly one vertex is given by \(p^{n'} = n^{-o(1)}\) as desired.
For [item:largerbox], recall that all boxes lie in a sector that has angle \(\phi = o(\log n / n)\). Focussing on boxes with expected number of vertices larger than 1, fix any such box \(\mathcal{B}_i\), and let \(X_i\) be the random variable with which we count the number of vertices in \(\mathcal{B}_i\). Using that we have a Poisson distribution, it follows that \[\mathbf{P}\left(X_i = 1\right) = \lambda \cdot \exp(-\lambda) \geq \exp(-\lambda) = \mathbf{P}\left(X_i = 0\right),\] since \(\lambda > 1\). Hence, the probability that a box has exactly one vertex is basically lower bounded by the probability that a box is empty. Now, let \(k\) be the number of boxes we consider, and we get from the equation above \(\prod_{i=1}^k \mathbf{P}\left(X_i = 1\right) \geq \prod_{i=1}^k\mathbf{P}\left(X_i = 0\right).\) On the other hand, since all boxes lie within a sector, the probability that all boxes are empty is lower bounded by the probability of the event that the entire sector is empty. Using that the sector has angle \(\phi = o(\log n / n)\), the expected number of vertices in a sector is \(\lambda_{\text{sector}} \approx n\cdot \phi = o(\log n)\). Thus, the probability of the event that all boxes are empty is \(\prod_{i=1}^k\mathbf{P}\left(X_i = 0\right) \geq \exp{(-\lambda_{\text{sector}})} = n^{-o(1)}\) by another application of a Poisson distribution. Putting everything together, it now follows that \(\prod_{i=1}^k \mathbf{P}\left(X_i = 1\right) \geq \prod_{i=1}^k\mathbf{P}\left(X_i = 0\right) = n^{-o(1)}\) as desired.
Of course, finding a nice sector is not enough: the resulting tree must be induced in the full HRG, and it must be attached to the giant component in a controlled way. To ensure this, we embed each candidate sector \(\Phi\) into a larger buffer sector \(\Psi\). The buffer is chosen large enough so that, conditioned on \(\Phi\) being nice and the relevant part of \(\Psi \setminus \Phi\) being empty, vertices of the tree have no additional neighbours outside the tree, except for a designated edge from the root to the giant component. We again show that this event holds with probability \(n^{-o(1)}\) and we use the standard fact that HRGs contain no vertices too close to the origin asymptotically almost surely to rule out long-range interference from the central part of the disk (see hatched area in 3b).
There are polynomially many disjoint buffer sectors \(\Psi\), and the events of the \(\Psi\)-sectors are independent of each other since vertices in two distinct “nice sectors” have disjoint neighbourhoods given that there is no vertex in the hatched area of 3b. Since each \(\Psi\)-sector succeeds with probability \(n^{-o(1)}\), a Chernoff bound implies that polynomially many of them contain the desired tree structure with high probability. Since we conditioned on the event that the hatched area of 3b is empty, an event that holds a.a.s., this proves the structural theorem and, via the lower-bound transfer from regular trees, yields [thm:lowerbound].
We finally outline the algorithms of [thm:geometric], where nodes are given access to geometric coordinates. This additional information changes the algorithmic approach. Instead of first producing random separators and then solving small residual components, the algorithms use the geometry directly.
For the MIS algorithm, we first partition the outer part of the hyperbolic disk into tiles arranged in \(\mathcal{O}(\log\log n)\) radial layers. The tiling is chosen to have two useful properties. First, every tile induces a clique (see red tiles 4b). Second, within each fixed layer, the tiles can be processed according to a constant-size schedule: in each step of the schedule, all currently active tiles are mutually non-adjacent and can therefore be handled in parallel; for a sketch, see red tiles in 4a.
The algorithm is then a simple layer-by-layer sweep. For each of the \(\mathcal{O}(\log\log n)\) layers, we execute the constant-size schedule for that layer. Whenever a tile is processed, we solve it locally by selecting one still-eligible vertex, if such a vertex exists, and adding it to the independent set. Since a tile is a clique, this is sufficient to resolve all the vertices of the tile. Since each layer requires only \(\mathcal{O}(1)\) rounds, the total runtime is \(\mathcal{O}(\log\log n)\).
Correctness for the outer disk follows directly from the sweep: every processed tile is solved when it is considered. It remains to argue that the inner disk (hatched area in 4b) is also dominated. Every inner-disk vertex has neighbours in many processed tiles; we show that by the tiling construction (with high probability over the draw of the HRG), any even adversarially chosen MIS (the algorithm is deterministic) within these tiles still contains at least one neighbour of the vertex. Hence, after the \(\mathcal{O}(\log\log n)\)-round sweep, every vertex is either selected or has a selected neighbour, and the resulting independent set is maximal.
For MM, a layer-by-layer sweep would only give an \(\mathcal{O}(\log\log n)\)-round algorithm, which is slower than our target runtime. Moreover, the MIS strategy does not transfer directly: matching one vertex in a tile only removes one neighbour, rather than dominating an entire neighbourhood. We therefore first remove the inner disk (hatched area in 4b) with the same constant-round randomised proposal step as for [thm:mainpolylog].
Afterwards, we use a bottom-up approach on the \(\mathcal{O}(\log\log n)\) remaining outer layers. Initially, we compute a maximal matching inside each single layer in \(\mathcal{O}(1)\) rounds using a tiling schedule that covers each intra-layer edge. We then repeatedly merge solved layers: when two neighbouring layers have already been solved separately, we may still need to add matching edges between them to obtain maximality for their union. Such merging steps can be performed in parallel for disjoint pairs of layers. The main technical work is to show that each individual merge can be implemented in \(\mathcal{O}(1)\) rounds as sketched in 11.
After \(t\) merging steps, the algorithm has solved blocks of \(2^t\) consecutive layers. Since there are only \(\mathcal{O}(\log\log n)\) layers, \(\mathcal{O}(\log\log\log n)\) steps suffice to solve the entire outer disk. Including the constant-round preprocessing of the inner disk, this gives the desired \(\mathcal{O}(\log\log\log n)\)-round algorithm for MM.
In this work, we investigate classical symmetry-breaking problems through the lens of hyperbolic random graphs. In particular, we resolve an open question posed by [23] by showing that maximal independent set and maximal matching can be computed in \(\operatorname{\text{{\rm poly}}}\log\log n\) rounds on hyperbolic random graphs. Thus, when the input graph is drawn from the hyperbolic random graph model, these problems admit substantially faster algorithms than the \(\Omega(\sqrt{\log n})\) round complexity known for worst-case instances [9]. Our result for the \(\mathsf{LOCAL}\) model relies on a constant-round algorithm that shatters an HRG into polylogarithmic-size components. The runtime then follows from the state-of-the-art deterministic algorithm [13]; any improvements for the deterministic algorithm in [13] would also imply a faster algorithm for hyperbolic random graphs.
Conversely, we show that MIS and MM are inherently harder than \((\Delta+1)\)-colouring on hyperbolic random graphs. While \((\Delta+1)\)-colouring can be solved in just 2 rounds [23], MIS and MM remain substantially more difficult, and we prove a lower bound of \(\Omega\left(\frac{\log\log n}{\log\log\log n}\right)\) rounds for both MIS and MM. Our lower bounds rely on new structural insights of hyperbolic random graphs that may be of independent interest. Specifically, we prove that the giant component of a hyperbolic random graph is likely to contain polynomially many logarithmic-size \(d\)-ary trees for a wide range of different degrees \(d\) and height \(h\).
Finally, for maximal matching, we demonstrate that access to the underlying geometric coordinates enables an exponentially faster algorithm than implied by this lower bound.
Our work opens up several research questions and the identified structural insights may also improve the design and analysis of distributed algorithms beyond the HRG model.
How does Luby’s algorithm perform on HRGs?
Is there a separation between MIS and MM on HRGs, as is the case on trees?
For embedding-aware graphs, can we overcome the lower bound barrier of [thm:lowerbound] for MIS?
How can we establish lower bounds in the setting of embedding-aware graphs?
What is the complexity of distributed approximation algorithms for Maximum Matching and Maximum Independent Set on HRGs?
Is the complexity landscape of other symmetry-breaking problems, such as edge colouring or vertex cover, similar to that of MIS/MM or \((\Delta+1)\)-colouring on HRGs?
The remainder of this article is structured as follows. In 4, we give a formal definition for the model of hyperbolic random graphs, and we collect important lemmas we make use of throughout this work. The analyses of our shattering algorithms can be found in 5. Our findings of \(d\)-ary trees and the lower bound implications are stated in 6.1. We conclude with our approach for algorithms where the embedding of an HRG is given (7).
In the following, we introduce the model of hyperbolic random graphs and follow the lines drawn by Papadopoulos et al. [30]. Points are represented by their native representation in the hyperbolic plane \(\mathbb{H}^2= [0, \infty) \times[0,2\pi)\). Thus, a point \(x \in \mathbb{H}^2\) is identified by a radial coordinate \(r(x)\) and an angular coordinate \(\varphi(x)\). We use exclusively the curvature -1 by which the hyperbolic distance \(\mathop{\mathrm{d_h}}(\cdot,\cdot)\) for two points \(x,y \in \mathbb{H}^2\) is given by \[\begin{align} \label{eq:hyperbolic95distance} \cosh(\mathop{\mathrm{d_h}}(x,y)) = \cosh(r(x))\cosh(r(y)) - \sinh(r(y))\sinh(r(x))\cos(\varphi(x) - \varphi(y)), \end{align}\tag{1}\] and consequently, we use for \(\mathbb{H}^2\) the topology induced by \(\mathop{\mathrm{d_h}}\).
For a hyperbolic random graph, we consider a subregion of the hyperbolic plane \(\mathbb{H}^2\): we operate on a disk with radius \(R\) denoted by \(\mathcal{D}_R= [0,R] \times [0,2\pi)\), centred at the origin \((0,0)\). Note that the set of points with distance \(r\) to \(x\in\mathcal{D}_R\) is given by the set \(\mathcal{B}_x(r) = \{y\in\mathcal{D}_R: \mathop{\mathrm{d_h}}(x,y)\leq r\}\), as we restrict the hyperbolic plane to \(\mathcal{D}_R\subset \mathbb{H}^2\).
In order to obtain a power-law degree distribution, we introduce the parameter \(\alpha \in (1/2, 1)\) (see also [30], [89]). Then, the probability measure \(\mu\) on \(\mathcal{D}_R\) for measurable \(\mathcal{S}\subseteq \mathcal{D}_R\) is defined by \[\begin{align} \mu(\mathcal{S}) = \int_S \rho(x) \,\mathrm dx, \qquad \rho(x) = \frac{\alpha\sinh(\alpha x)}{2\pi(\cosh(\alpha R) - 1)}, \end{align}\] where \(\rho\) is the density of \(\mu\) with respect to the Lebesgue measure on \(\mathcal{D}_R\).
Vertices \(V\) are identified by their point coordinates in \(\mathcal{D}_R\): for a subset of vertices that are identified by a set of points \(\mathcal{S} \subseteq \mathcal{D}_R\), we write \(V \cap \mathcal{S}\). Throughout this work, we exclusively work on the Poissonised version of (threshold) hyperbolic random graphs (see also [39], [53], [56]). This model allows us to analyse the number of vertices in disjoint areas of the hyperbolic disk independently. Let \(n \in \mathbb{Z}^+\) and \(N\) be a Poisson random variable with expectation \(\mathbf{E}\left[N\right] = n\). Then, let \(R:= 2\log(n) + C\) for some constant \(C \in \mathbb{R}\) and we use an inhomogeneous Poisson point process on \(\mathcal{D}_R\) where for any \(\mathcal{S} \subseteq \mathcal{D}_R\) it holds that the number of expected vertices is given by \(n\cdot\mu(\mathcal{S}) = \mathbf{E}\left[|V \cap \mathcal{S}|\right]\). That is, the intensity function at polar coordinates \((r, \varphi)\) for \(0\leq r < R\) is \(\lambda(r, \varphi):=n\rho(x)\) and the set of vertices is a random variable \(V = \{(r_1,\varphi_1), (r_2,\varphi_2),.., (r_N,\varphi_N)\}\). For a pair of vertices \(u,v \in V \cap \mathcal{D}_R\), there exists the edge \(\{u,v\} \in E\) in a threshold hyperbolic random graph, if and only if it holds \(\mathop{\mathrm{d_h}}(u,v) \leq R\).
In the following, we collect some results we make use of throughout this work. It will be convenient to characterise the connection of vertices in terms of their angular distance \(\delta_\varphi(\cdot,\cdot)\). Since vertices are connected if and only if their distance is at most \(R\), we define \[\begin{align} \label{eq:angle-func} \theta_R(r_1,r_2) = \arccos\left(\frac{\cosh(r_1)\cosh(r_2) - \cosh(R)}{\sinh(r_1)\sinh(r_2)}\right), \end{align}\tag{2}\] which denotes by 1 the angle distance for two points with radii \(r_1\) and \(r_2\) such that their hyperbolic distance is exactly \(R\). Thus, two vertices with an angle distance of at most \(\theta_R\) are connected. Throughout this work, we use for \(\theta_R(\cdot,\cdot)\) the following bounds, which are due to [83].
Lemma 1 (Angle distance). Let \(x,y \in \mathcal{D}_R\) and \(1 \leq r(x), r(y) \leq R\) and \(r(x)+r(y) \geq R\). Then it holds \[\sqrt{e^{R-r(x)-r(y)}} \leq \theta_R(r(x), r(y)) \leq \pi \sqrt{e^{R-r(x)-r(y)}}.\]
As remarked in [56], the function \(\theta_R(\cdot, \cdot)\) is decreasing in both arguments.
\(\theta_R(\cdot, \cdot)\) is strictly decreasing in both arguments.
Next, we look into the measure for different areas in \(\mathcal{D}_R\). First, we consider a ball of radius \(r\) with the origin as its centre point, as sketched by the red area in 5a ([89]).
Lemma 2 (Measure of inner disk). For any \(0 \leq r \leq R\) we have \(\mu\left(\mathcal{B}_0(r)\right) = (1 + o(1))e^{-\alpha(R- r)}~.\)
We write \(N(u) := \{v \in V : \{u,v\} \in E(G)\}\) for the neighbourhood of \(u\) and \(\deg(u) := |N(u)|\) for the degree of \(u\). Since a vertex \(u\) has an edge to every vertex within distance \(R\), it is adjacent to all vertices in \(\mathcal{B}_u(R) \cap \mathcal{D}_R\) (see blue area in 5a). For the expected degree of a vertex, we use the following bounds ([23]).
Lemma 3 (Vertex Degree). Let \(u \in V \cap \mathcal{D}_R\) be a vertex with radius \(r\geq 1\). Then, the expected degree of \(u\) in a threshold HRG is \[n e^{-r/2}\cdot\frac{(1-o(1))\alpha }{\pi(\alpha-1/2)} \leq \mathbf{E}\left[\deg(u)\right] \leq ne^{-r/2}\cdot \frac{ (1+o(1))\alpha}{\alpha-1/2}.\]
Another sub-area of the disk we use are layers as used by Friedrich and Krohmer [58]. For \(\ell \in [\lfloor R \rfloor]\), a layer is defined by \(\mathcal{L}_{\ell}:= \mathcal{B}_0(R-\ell)\setminus \mathcal{B}_0(R-\ell-1)\) (see 5b for a sketch). Using 2 and 3, we obtain the following.
Lemma 4 (Layer properties). Let \(\ell \in [\lfloor R \rfloor]\) and \(u \in V \cap \mathcal{L}_{\ell}\). Then \[\label{eq:layer95measure} \mu(\mathcal{L}_{\ell}) = (1-e^{-\alpha} +o(1)) e^{-\alpha \ell},\qquad{(1)}\] and \[\label{eq:layer95expected95degree} e^{\ell/2} \cdot \frac{ (1-o(1))\alpha e^{C/2}}{\pi(\alpha-1/2)} \leq \mathbf{E}\left[\deg(u)\right] \leq e^{\ell/2} \cdot \frac{ (1+o(1))\alpha e^{C/2}}{\alpha-1/2}.\qquad{(2)}\]
For notational convenience, we sometimes write \(V_\ell := V \cap \mathcal{L}_{\ell}\) for the set of vertices in layer \(\ell\).
Finally, the following statement will be useful: it states that if two vertices share a neighbour with larger radius, then they form a connected triangle with that neighbour.
Lemma 5. Let \(x,y,z \in \mathcal{D}_R\) such that \(1 \leq r(x) \leq r(y) \leq r(z) - 2\log(2\pi) \leq R - 2\log(2\pi)\). Then, if \(\mathop{\mathrm{d_h}}(x,z) \leq R\) and \(\mathop{\mathrm{d_h}}(y,z) \leq R\), it holds that \(\mathop{\mathrm{d_h}}(x,y) \leq R\).
Proof. By our hypothesis that \(\mathop{\mathrm{d_h}}(x,z) \leq R\) and \(\mathop{\mathrm{d_h}}(y,z) \leq R\), it follows for the angular distance between \(x\) and \(y\) by 1 that \[\begin{align} \delta_\varphi(x,y) \leq \theta_R(r(x), r(z)) + \theta_R(r(y), r(z)) \leq 2\pi\sqrt{e^{R-r(x)-r(z)}}, \end{align}\] since \(r(x) \leq r(y)\) and using [rmk:theta-monotonicity]. Consequently, using our hypothesis that \(r(z) \geq r(y) + 2\log(2\pi)\), it follows \[\begin{align} \delta_\varphi(x,y) \leq \sqrt{e^{R-r(x)-r(y)}} \leq \theta_R(r(x), r(y)), \end{align}\] and thus, by 1 we have \(\mathop{\mathrm{d_h}}(x,y) \leq R\) as desired. ◻
We write \([k]\) (where \(k \in \mathbb{Z}^+\)) for the set \(\{0,1,\dots, k-1\}\). For the randomness over the graph distribution of a hyperbolic random graph, we use \(\mathbb{P}\left(\cdot\right)\) and \(\mathbb{E}\left[\cdot\right]\). Conversely, we write \(\mathbf{P}\left(\cdot\right)\) and \(\mathbf{E}\left[\cdot\right]\) when we deal with the randomness of an algorithm on a sampled hyperbolic random graph \(G\). A sector \(\Phi \subseteq \mathcal{D}_R\) with angle \(\phi\) and bisector \(\varphi\), is defined by the set of points \(\Phi := \{x \in \mathcal{D}_R: \varphi - \phi / 2 \leq \varphi(x) \leq \varphi + \phi / 2 \}\). For a vertex \(u \in V\), we write \(E(u) := \{\{u,v\} \in E : v \in N(u)\}\) for the set of incident edges to \(u\).
The goal of this section is to prove the following theorem. All algorithms follow a shattering-based approach that exploits structural properties of hyperbolic random graphs to identify separators that shatter the graph into small components. These can then be solved independently and in parallel using any deterministic algorithm. The overall running time is determined by the algorithm used to solve the remaining components, as the shattering itself takes constant rounds. The proof for the MIS part can be found in 5.1. The MM part we show in 5.2.
We use Luby’s algorithm as a subroutine. In Luby’s algorithm, all (active) vertices draw a random ID. If a vertex \(u\) has the largest ID among its neighbours \(N(u)\), then \(u\) is added to the independent set \(I\). An iteration of Luby’s algorithm has \(2\) rounds: first, all active vertices send their random ID to all neighbours. Then any active vertex \(u\) informs all its neighbours based on the IDs if \(u\) is (a) in the independent set, (b) has a neighbour that is part of the independent set, or (c) neither of the two. For \(t \in \mathbb{Z}^+\), we denote by \(V_{(t)}\) the set of vertices that are neither part of the independent set nor have a neighbour in the independent set after the \(t\)-th step of HRG-Shattering-MIS, where by HRG-Shattering-MIS we refer to the following process.
Step 1: Activate all vertices of the set \(U_{(1)} : = \{v \in V : \left\lfloor\frac{\log^{4}(n)}{1000}\right\rfloor\;\leq \deg(v) \leq \lceil1000\log^{4}(n) \rceil\}\). Execute one iteration of Luby’s algorithm on active vertices \(U_{(1)}\).
Step 2: Activate all vertices of the set \(U_{(2)} := \{v \in V_{(1)} : \deg(v) \leq \lceil \log^{3/2}(n)\rceil\}\). Execute one iteration of Luby’s algorithm on active vertices \(U_{(2)}\).
We analyse step 1 in 6 and step 2 is addressed in 7. We then use the two lemmas to show that HRG-Shattering-MIS shatters a hyperbolic random graph into components that are all of size at most \(\operatorname{\text{{\rm poly}}}\log n\); see also [pro:shattering-mis].
The following lemma tells us that when we consider the independent set \(I\) generated after the first step of HRG-Shattering-MIS, that (1) all vertices of \(I\) are contained in an annulus \(\mathcal{A}\) with constant thickness and (2) for any sector with a large enough angle, we find a vertex that is part of the independent set \(I\); see also 6a for a sketch.
Lemma 6 (Step 1 of HRG-Shattering-MIS). Let \(G\) be a threshold hyperbolic random graph and consider the set of vertices \(I_{(1)}\) that is part of the independent set after step 1 of HRG-Shattering-MIS. Moreover, let \(c(C, \alpha)\) be a constant large enough and consider annulus \(\mathcal{A} \coloneq \mathcal{B}_0(R - {8\log\log n} + c(C, \alpha))\setminus \mathcal{B}_0(R - {8\log\log n} - c(C, \alpha))\). Then the following holds for \(I_{(1)}\) with probability \(1 - n^{-\omega(1)}\):
All vertices contained in the independent set after the first step are in \(\mathcal{A}\); \(I_{(1)} \subseteq U_{(1)} \subseteq V \cap \mathcal{A}\).
In any sector with angle \(\frac{\log^{30\alpha}n}{n}\) there exists a vertex \(v\) such that \(v \in I_{(1)}\).
Proof. We prove the first part of the lemma by showing that, w.e.h.p., no vertex outside the annulus has the same degree as an active vertex \(U_{(1)}\). Consequently, no such vertex joins the independent set in the first step.
To this end, let \(c(C, \alpha) \mathrel{\vcenter{:}}= 100 + |C| + 2\log(1/(\alpha - 1/2))\) chosen with hindsight, and consider any vertex \(u \in V \cap \mathcal{B}_0(R - {8\log\log n} - c(C, \alpha))\). Then, for \(u\) it holds via 3 that \[\begin{align} \mathbf{E}\left[\deg(u)\right] &\geq \frac{(1-o(1))\alpha}{\pi(\alpha - 1/2)} e^{(c(C, \alpha) - C)/2} \log^{4}n\\ &\geq 10000 \log^{4}n, \end{align}\] by our choice of \(c\). Thus, a Chernoff-bound (24) and union bound yield that for any \(u \in V \cap \mathcal{B}_0(R - {8\log\log n} - c(C, \alpha))\), that \(\deg(u) > \lceil1000\log^{4} n \rceil\) w.e.h.p.
Analogously, we get for any vertex \(u \in V \cap \mathcal{D}_R \setminus \mathcal{B}_0(R - {8\log\log n} + c(C, \alpha))\) that \[\begin{align} \mathbf{E}\left[\deg(u)\right] & \leq \frac{(1+o(1))\alpha}{(\alpha - 1/2)} e^{-(c(C, \alpha) + C)/2} \log^{4}n\\ & \leq \frac{\log^{4} n}{10000} . \end{align}\] Another combination of Chernoff bound with a union bound then reveals that for any \(u \in V \cap \mathcal{D}_R \setminus \mathcal{B}_0(R - {8\log\log n} + c(C, \alpha))\) that \(\deg(u) < \lfloor\frac{\log^{4} n}{1000}\rfloor\) w.e.h.p. This concludes the proof for [item:ring].
We now turn to the second point of the lemma. Consider a sector \(\Phi\) with angle \(\phi = \log ^{16\alpha }(n) /n\). W.l.o.g. let \(\Phi\) have bisector \(0\). Then, let \(\hat{\ell} = \lceil 8\log\log n + 2\log(\alpha - 1/2) - C \rceil\), \(\hat{\phi} = \log^{10\alpha}(n)/n\) and consider the area \(\mathcal{S} := \{x \in \mathcal{L}_{\hat{\ell}} : 2\pi - \hat{\phi} /2 \leq \varphi(x) \leq \hat{\phi}/2 \}\) (see also 6b for a sketch). We say \(\Phi\) is active if \(I_{(1)} \cap \mathcal{S} \neq \emptyset\), i.e., there is a vertex in the area \(\mathcal{S}\) that joins the independent set after the first step of HRG-Shattering-MIS.
A sector \(\Phi\) with angle \(\phi = \log ^{16\alpha }(n) /n\) is active with probability \(\Omega(1/\log^{4} n)\).
Proof of claim. We lower bound the probability that \(\Phi\) is active as follows: first, we show that any vertex in \(\mathcal{S}\) is active in the first step w.e.h.p. Then, we show that \(\mathcal{S}\) contains at least one vertex w.e.h.p. The bound then follows using a union bound over the complementary of these two events and since any active vertex \(u\) in step 1 of HRG-Shattering-MIS has degree \(\deg(u) \in \Theta(\log^{4} n)\).
By our choice of \(\hat{\ell}\) and \(\mathcal{S} \subset \mathcal{L}_{\hat{\ell}}\), we get for any \(u \in V \cap \mathcal{S}\) by using 4 that \[\frac{ (1-o(1))\alpha}{\pi} \cdot \log^4 n \leq \mathbf{E}\left[\deg(u)\right] \leq (1+o(1))\alpha \cdot \log^4 n.\] Since \(\alpha \in (1/2, 1)\), using Chernoff bounds then yields that for \(u \in V \cap \mathcal{S}\) it holds that \(\left\lfloor\frac{\log^{4}(n)}{1000}\right\rfloor\;\leq \deg(u) \leq \lceil1000\log^{4}(n) \rceil\) w.e.h.p. Hence, \(u\) is active for the first step of Luby in HRG-Shattering-MIS and a union bound shows that any vertex in \(\mathcal{S}\) is active w.e.h.p.
Next we use that \(\mathcal{S}\) spans an angle of \(\hat{\phi}=\log^{10\alpha}(n)/n\). This in conjunction with \(\mathcal{S} \subset \mathcal{L}_{\hat{\ell}}\) and \(\hat{\ell} \in 8\log\log n + \mathcal{O}(1)\) then yields via 4 that \[\mu(\mathcal{S}) = \hat{\phi}/2\pi \cdot \mu(\mathcal{L}_{\hat{\ell}}) =\Theta(1)\frac{\log^{10\alpha}n}{n} \cdot e^{-8\alpha \log\log n} \in \omega\left(\frac{\log n}{n}\right),\] since \(\alpha > 1/2\). Consequently, by another Chernoff bound, \(\mathcal{S}\) is non-empty w.e.h.p. The claim now follows since any vertex in by our outline of the proof for Claim [claim:active-sector] since any active vertex \(u\) joins \(I\) with probability at least \(1/(\deg(u) +1) \in \Omega(1/\log^4 n)\). ◻
Now, partition the disk \(\mathcal{D}_R\) into \(\lceil 2\pi n \cdot \log ^{-28\alpha } n\rceil =: k\) sectors such that each sector has angle \(\psi \in \Theta( \log ^{28\alpha } (n)/n)\). Note that if each such sector contains at least one vertex that joins the independent set \(I_{(1)}\), then it also holds for any sector with angle at least \(\frac{\log^{30\alpha}n}{n}\) that there exists a vertex that is in \(I_{(1)}\) as desired for [item:active-sector].
To show that this holds with the desired probabilistic guarantee, fix any sector \(\Psi\) with angle \(\psi\) and partition \(\Psi\) into \(k' \in \Theta(\log ^{12\alpha }n)\) sectors, such that each sector has angle \(\phi = \log ^{16\alpha }(n) /n\). For the \(i\)-th sector with angle \(\phi\), we write \(\Phi_i\) and let \(X_i\) be the indicator random variable that is \(1\) if the \(i\)-th sector is active. Moreover, let \(X = \sum_{i}^{k'}{X_i}\) so that by linearity of expectation the expected number of active sectors with angle \(\phi\) in \(\Psi\) is \[\mathbb{E}\left[X\right] = k' \cdot\mathbb{E}\left[X_i\right] \in \Omega(\log n^{12\alpha - 4}n) \in \omega(\log n),\] by applying Claim [claim:active-sector] and using that \(\alpha > 1/2\). Next, for any pair of sectors \(\Phi_i\) and \(\Phi_j\), we show that \(X_i\) and \(X_j\) are independent w.e.h.p.
To see this, recall that \(\Phi_i\) is active if the area \(\mathcal{S}_i \subset \Phi_i\) with angle \(\hat{\phi} = \log^{10\alpha}(n)/n\) contains a vertex that draws the largest ID among active neighbours. W.l.o.g. let \(\Phi_i\) have bisector \(0\). Moreover, using [item:ring], any active vertex \(v\) has radius \(r(v) \in R - 8\log\log n - \Theta(1)\) w.e.h.p. Hence, using 1 in conjunction with [rmk:theta-monotonicity], any active vertex \(v\) that is a neighbour of \(u \in V \cap \mathcal{S}_i\) has angle \[\varphi(v) \leq \hat{\phi} + \theta_R(r(u) , r(v)) \in \mathcal{O}\left(\frac{\log^8 n +\log^{10\alpha} n }{n} \right) \in o(\phi) \text{ w.e.h.p.},\] where the last step follows since \(\psi = \log ^{16\alpha }(n) /n\) and \(\alpha > 1/2\). Consequently the event only depends on the randomness in \(\Phi_i\) w.e.h.p. A union bound over \(k'\) sectors with angle \(\phi\) shows that this also holds for all \(k'\) \(\Phi\)-sectors in sector \(\Psi\) w.e.h.p.
Thus, using that for any pair of sectors \(\Phi_i, \Phi_j \subset \Psi\) the respective random variables \(X_i\) and \(X_j\) are independent w.e.h.p. and using \(\mathbb{E}\left[X\right] \in \omega(\log n)\), we obtain via a Chernoff bound that \(\Psi\) contains at least \(1\) vertex that is in \(I_{(1)}\) w.e.h.p. A union bound over the \(k \in o(n)\) \(\Psi\)-sectors with angle \(\psi\) wraps up the proof since every sector with angle \(\psi \in \Theta(\log^{28\alpha}(n)/n) \in o(\log^{30\alpha}(n)/n)\) contains a vertex of \(I_{(1)}\) w.e.h.p. ◻
Next, we show that after step two of HRG-Shattering-MIS, all vertices of degree larger than \(\approx \log^{7/2} n\) are removed by having a neighbour in the independent set \(I_{(2)}\) (except for the vertices we included in the independent set after step one). For the area of these vertices, see also the hatched area in 7a.
Lemma 7 (Step 2 of HRG-Shattering-MIS). Let \(r \mathrel{\vcenter{:}}= R - {7\log\log n}\). Then after the second step of HRG-Shattering-MIS, for any vertex \(u \in V_{(1)} \cap \mathcal{B}_0(r)\) it holds with probability \(1 - n^{-\omega(1)}\) that \(u \not\in V_{(2)}\).
Proof. Consider the annulus \(\mathcal{A}' := \mathcal{D}_R\setminus \mathcal{B}_0(R - 4 \log\log n - C)\). Then, by 3, any vertex \(v \in V \cap (\mathcal{D}_R\setminus\mathcal{A}' )\) has expected degree \(\mathbf{E}\left[\deg(v)\right] \in \Omega(\log^2 n)\). Using a Chernoff- and a union bound, all vertices outside the annulus \(\mathcal{A}'\) have a degree larger than \(\lceil \log^{3/2}(n)\rceil\) w.e.h.p. Hence, all vertices that are active in the second step of HRG-Shattering-MIS are within the annulus \(\mathcal{A}'\), i.e., \(U_{(2)} \subseteq V_{(1)} \cap \mathcal{A}'\) w.e.h.p.
Now, using \(r= R- \log\log n\) as defined in our lemma statement, fix a vertex \(u \in V_{(1)} \cap \mathcal{B}_0(r)\), i.e., a vertex in the inner disk that was not removed in the first round of HRG-Shattering-MIS. We prove our desired statement by the following case distinction: (1) there exists at least one vertex \(v\) in the neighbourhood of \(u\), such that \(v\) lies in the annulus \(\mathcal{A}'\) and \(v\) was removed after the first step of HRG-Shattering-MIS, i.e., \((N(u) \cap \mathcal{A}' )\setminus (V \setminus V_{(1)}) \neq \emptyset\) (see also 7a), and (2) non of the neighbours of \(u\), that lie in annulus \(\mathcal{A}'\) was removed after the first step of HRG-Shattering-MIS, i.e., \((N(u) \cap \mathcal{A}' )\setminus (V \setminus V_{(1)}) = \emptyset\).
Case 1 [\((N(u) \cap \mathcal{A}' )\setminus (V \setminus V_{(1)}) \neq \emptyset\)]: Consider a vertex \(v \in (N(u) \cap \mathcal{A}' )\setminus (V \setminus V_{(1)})\). Since \(v \notin V_{(1)}\), either (a) \(v\) is in the independent set after the first step, \(v \in I_{(1)}\), or (b) \(v\) has a neighbour \(w \in I_{(1)}\) that is in the independent set after the first step. Using 6, no vertex in the set \(I_{(1)}\) is in the annulus \(\mathcal{A}'\) w.e.h.p. and thus, case (a) does not occur w.e.h.p. Hence, only case (b) remains (see 7a for a sketch). Using again 6, it follows that vertex \(w \in N(v) \cap I_{(1)}\) has radius \(r(w) \leq R - 8\log\log n + c\) w.e.h.p. Then, since \(v \in V \cap \mathcal{A}'\), it holds \(r(v) \geq R - 4\log\log n\) - C. Thus, using 5, \(u \in N(w)\) since \(r(u) \leq R - 7\log\log n\). It follows that \(u \not\in V_{(1)}\) w.e.h.p. since \(w \in I_{(1)}\); \(u\) was removed after the first step of HRG-Shattering-MIS, and thus, we conclude by (a) and (b) that case 1 does not occur w.e.h.p.
Case 2 [\((N(u) \cap \mathcal{A})\setminus (V \setminus V_{(1)}) = \emptyset\)]: Consider the sector \(\Psi\) with angle \(\psi = e^{C/2} \cdot \log^{7/2} n /n\) and bisector \(\varphi(u)\). Note that, since \(r(u) \leq R - 7\log\log n\), we obtain via 1 in conjunction with [rmk:theta-monotonicity] that \(\Psi \subset \mathcal{B}_u(R)\), and thus, if a vertex \(v \in V \cap \Psi\) joins the independent set, \(u\) is removed since \(v \in N(u)\). We partition \(\Psi\) into \(k:= e^{C}\cdot\log^{3/2}n / 100\) sectors, such that for \(i \in [k]\) each sector \(\Phi_i \subset \Psi\) has angle \(\phi = 100\cdot e^{-C/2}\log^2 n /n\). Let \(\varphi_i\) be the bisector of \(\Phi_i\) and let \(\mathcal{S}_i := \{x \in \Phi_i \cap \mathcal{L}_{0} : \varphi_i - 1/n \leq \varphi(x) \leq \varphi_i + 1/n \}\), i.e., a sub-area of sector \(\Phi_i\) in layer \(\mathcal{L}_{0}\), with the same bisector as \(\Phi_i\) that spans an angle of \(2/n\) (see also green area in 7b). We say that \(\Phi_i\) is active if there exists a vertex \(v \in V \cap \mathcal{S}_i\) such that \(v\) joins the independent set in the second step of HRG-Shattering-MIS, i.e., \(I_{(2)} \cap \mathcal{S}_i \neq \emptyset\).
A sector \(\Phi_i\) is active with non-vanishing probability.
Proof of claim. Recall that due to case 2, we condition on the event that no neighbour \(v \in N(u) \cap \mathcal{A}'\) has a neighbour \(w \in N(w) \cap I_{(1)}\) in the independent set after the first step of HRG-Shattering-MIS.
Thus, using that \(|V \cap \mathcal{S}_i|\) follows a Poisson distribution and \[\mu(\mathcal{S}_i) = \frac{\mu(\mathcal{L}_{0})}{\pi n} \in \Omega(1/n),\] by 4 and \(\mathcal{S}_i\) spanning an angle of \(2/n\), it holds \(|v \in \mathcal{S}_i \cap V_{(1)}| = 1\) with non vanishing probability. We write for this event \(\mathcal{E}_i := \{|V \cap \mathcal{S}_i| =1\}\).
Moreover, since \(v \in V \cap \mathcal{L}_{0}\) it holds \(\mathbf{E}\left[\deg(v)\right] \in \mathcal{O}(1)\) by 4. Clearly, the number of active neighbours of \(v\) is upper bounded by \(\deg(v)\). Thus, again it follows from a Poisson distribution that the number of active neighbours of \(v\) is \(\mathcal{O}(1)\) with non-vanishing probability when conditioned on event \(\mathcal{E}_i\) (i.e., \(\mathcal{S}_i\) is empty except for \(v\)).
The claim then follows since, given that \(v\) has at most constant active neighbours, \(v\) joins the independent set with non-vanishing probability. ◻
To finish the proof, let \(X_i\) be the indicator random variable that is \(1\) if sector \(\Phi_i\) is active and let \(X = \sum_{i=0}^{k-1}X_i\) be the number of “active \(\Phi\)-sectors in sector \(\Psi\)”. By linearity of expectation and Claim [claim:active-sector2] in conjunction with \(k \in \Omega(\log n^{3/2})\) it holds \(\mathbb{E}\left[X\right] \in \omega(\log n)\). Thus, if we can show for any \(X_i\) and \(X_j\) that they are independent, our desired statement follows from a Chernoff bound in conjunction with a union bound over the set of vertices \(V \cap \mathcal{B}_0(R - 7\log\log n)\). This is since an active \(\Phi\)-sector implies that a neighbour \(v \in N(u)\) joins the independent set which removes \(u\).
We show that this independence holds for any two indicator random variables \(X_i\) and \(X_j\) w.e.h.p. Recall that \(X_i\) is \(1\) if and only if there is a vertex \(v \in V \cap \mathcal{S}_i\) that joins the independent set. Moreover, recall that w.e.h.p. only the area \(\mathcal{A}' = \mathcal{D}_R\setminus \mathcal{B}_0(R-4\log\log n -C)\) contain active vertices. W.l.o.g., let sector \(\Phi_i\) have bisector \(0\). Thus, considering \(v \in V \cap \mathcal{S}_i\) with radius \(r(v) \geq R - 1\), for any active neighbour \(w \in N(v)\) we have \(r(w) \geq R - 4\log\log n\) w.e.h.p., and it holds \[\varphi(w) \leq 1/n + \theta_R(r(v), r(w)) \leq \frac{\pi \cdot e^{(1-C)/2} \log^2 n}{n} < \phi \text{ w.e.h.p.,}\] using 1 and [rmk:theta-monotonicity] in conjunction with \(\phi = 100 \cdot e^{-C/2} \log ^2 n /n\) and \(\mathcal{S}_i\) spanning at most an angle of \(1/n\). Thus, all active neighbours of the vertex \(v \in V \cap \mathcal{S}_i\) are contained in sector \(\Phi_i\) w.e.h.p., yielding the independence w.e.h.p. as desired. A Chernoff bound yields \(X \in \omega(\log n)\) w.e.h.p. such that \(u\) is removed w.e.h.p. and a union bound over all the set of vertices \(|V \cap \mathcal{B}_0(R - 7\log\log n)| \in o(n)\) w.e.h.p. wraps up the proof. ◻
We now have all ingredients to show that HRG-Shattering-MIS shatters a hyperbolic random graph into components of size at most \(\operatorname{\text{{\rm poly}}}\log n\).
For a threshold hyperbolic random graph \(G\), it holds with probability \(1 - n^{-\omega(1)}\) that, after the second step of HRG-Shattering-MIS, the largest connected component in \(G[V_{(2)}]\) is of size at most \(\log^{\mathcal{O}(1)}(n)\).14 The round complexity of HRG-Shattering-MIS for \(\mathsf{CONGEST}\) is \(\mathcal{O}(1)\).
Proof. The runtime \(\mathcal{O}(1)\) directly follows from the fact that one iteration of Luby’s algorithm requires \(\mathcal{O}(1)\) rounds in \(\mathsf{CONGEST}\). Now, consider the set of vertices \(I_{(1)}\), i.e., the set of vertices in the independent set after the first step of HRG-Shattering-MIS. Let \(k := |I_{(1)}|\) and consider for \(I_{(1)}\) an ascending ordering by angular coordinates \(v_0, v_1,\dots, v_{k-1}\). For \(i \in [k]\), let \(\Phi_i\) be the sector defined by \(\{x \in \mathcal{D}_R: \varphi (v_{i}) \leq \varphi(x) \leq \varphi(v_{i+1})\}\) where \(v_k = v_0\). Note that by 6, for \(i \in [k]\) it holds for any consecutive pair \(v_i, v_{i+1}\) that \(\delta_\varphi(v_{i},v_{i+1}) \leq 2 \log^{30\alpha}(n)/n\) w.e.h.p. and thus, the angle for any sector \(\Phi_i\) is at most \(\mathcal{O}(\log^{30\alpha}(n)/n)\) w.e.h.p. Thus, it holds for any \(i \in [k]\) that \(\mu(\Phi_i) \in \mathcal{O}(\log^{30\alpha}(n)/n)\). Via Chernoff bound for \(|V \cap \Phi_i|\) and a union bound over \(k \in o(n)\) sectors, we then obtain for all \(i \in [k]\) that \(|V \cap \Phi_i| \in \mathcal{O}(\log^{30\alpha}n)\) w.e.h.p. Clearly, this implies that \(|V_{(2)} \cap \Phi_i| \in \mathcal{O}(\log^{30\alpha}n)\) w.e.h.p. Thus, if we can show for any pair \(i \neq j \in [k]\) that for all pairs \(u \in V_{(2)} \cap \Phi_i\) and \(v \in V_{(2)} \cap \Phi_j\) it holds \(\mathop{\mathrm{d_h}}(u,v) > R\) w.e.h.p., this finishes the proof.
To this end, fix any pair of sectors \(\Phi_i\) and \(\Phi_{j}\), and w.l.o.g. let \(j = i+1\) since non-adjacent sectors are separated by an even larger angular interval. Consider any pair \(u \in V_{(2)} \cap \Phi_i\) and \(v \in V_{(2)} \cap \Phi_{i+1}\). Using 7, it holds w.e.h.p. for both vertices \(u, v \in V_{(2)}\) that \(r(u), r(v) \geq R - 7 \log\log n\). Moreover, let \(s = v_{i+1}\) be the vertex that lies on the intersecting ray of \(\Phi_i\) and \(\Phi_j\). Then, by 1 in conjunction with [rmk:theta-monotonicity], it holds for the angular distance between \(u\) and \(v\) that \[\delta_\varphi(u,v) \geq \theta_R(r(v), r(s)) + \theta_R(r(u), r(s)) > \theta_R(r(u), r(v)) \text{ w.e.h.p.,}\] since \(u, v \not\in N(s)\) and \(r(v), r(u) \in r(s) + \Omega(\log \log n)\) w.e.h.p. by 6. Thus, for all pairs \(u \in V_{(2)} \cap \Phi_i\) and \(v \in V_{(2)} \cap \Phi_j\) it holds \(\mathop{\mathrm{d_h}}(u,v) > R\) w.e.h.p.; taking a union bound over all pairs of sectors, the desired result follows: every connected component of \(G[V_{(2)}]\) is contained in a single sector \(\Phi_i\). Since each sector contains at most \(\mathcal{O}(\log ^{30\alpha}n)\) vertices w.e.h.p., every connected component has poly-logarithmic size. This proves the proposition. ◻
Proof of [thm:mainpolylog] (MIS part). \(\mathsf{LOCAL}\): For \(\mathsf{LOCAL}\) MIS, we obtain a shattering where each connected component is of size at most \(\mathop{\mathrm{polylog}}n\) w.e.h.p. after \(\mathcal{O}(1)\) using [pro:shattering-mis]. Applying in parallel for each component [13], we obtain an MIS after \(\tilde{\mathcal{O}}(\log^{5/3}\log n)\) rounds w.e.h.p.
\(\mathsf{CONGEST}\): For \(\mathsf{CONGEST}\) MIS, [pro:shattering-mis] implies that after \(\mathcal{O}(1)\) rounds each connected component induced by undecided nodes is of at most poly-logarithmic size w.e.h.p. Consequently, applying in parallel for each connected component [44], we obtain an MIS after \(\tilde{\mathcal{O}}(\log^{3}\log n)\) rounds w.e.h.p. ◻
In this section, we turn to our shattering algorithm, which we use to obtain a maximal matching \(M \subseteq E\) in \(\mathop{\mathrm{polylog}}\log n\) rounds. For \(t \in \mathbb{Z}^+\), let \(M_t\) be the set of edges included in our matching after step \(t\) of our algorithm HRG-Shattering-MM(see description below). We denote by \(V_{(t)}\) the set of unmatched vertices after step \(t\): a vertex is unmatched if it has at least one incident edge that can be included without violating the matching condition after step \(t\) of HRG-Shattering-MM, i.e., \(V_{(t)} :=\{u \in V : \exists v \in N(u) \text{ such that } E(v) \cap M_t = \emptyset \text{ and } E(u) \cap M_t = \emptyset\}\). The following is a brief description of our algorithm, which we refer to as HRG-Shattering-MM:
Step 1 (Inner disk removal): Activate all vertices of the set \(U_{(1)} := \{v \in V_{(1)} : \deg(v) \leq \lceil \log^{3/2}(n)\rceil\}\). In parallel, each active vertex \(u \in U_{(1)}\) marks one incident edge uniformly at random. Then, in parallel, each vertex \(u \in V\) chooses a marked edge \(\{u,v\} \in E\) uniform at random from its incident marked edges, and we add \(\{u,v\}\) to the matching \(M\).
Step 2 (Separating the outer disk): Activate all vertices of the set \(U_{(2)} : = \{v \in V : \left\lfloor\frac{\log^{4}(n)}{1000}\right\rfloor\;\leq \deg(v) \leq \lceil1000\log^{4}(n) \rceil\}\). In parallel, each vertex \(u \in U_{(2)}\) draws a random ID, and we add \(u\) to the set \(S\) if \(u\) has the largest ID among its neighbours in \(U_{(2)}\).
If \(\mathsf{LOCAL}\): In parallel, each vertex \(s \in S\) collects the adjacency matrix of the induced subgraph \(G_s := G[N_2(s) \cap V_{(1)}]\).15 In parallel, each vertex \(s \in S\) computes a maximal matching \(M(s) \subseteq E(G_s)\) and we add \(M(s)\) to the matching \(M\).
If \(\mathsf{CONGEST}\): In parallel, each vertex \(s \in S\) activates the set of vertices \(U(s) := N_2(s) \cap V_{(1)}\). In parallel, for each \(s \in S\), compute a maximal matching \(M(s)\) for \(G[U(s)]\) via [85] and we add \(M(s)\) to the matching \(M\).
The rest of this section is dedicated to proving that HRG-Shattering-MM shatters an HRG into polylogarithmic connected components; see also [pro:shattering-mm]. We start by showing in 8 that after step 1, in constant rounds, all vertices of the inner disk \(\mathcal{B}_0(R - 6\log\log n)\) are matched w.e.h.p. We then continue by analysing step 2 and the properties it yields for set \(S\) in 9.
Lemma 8 (Step 1 of HRG-Shattering-MM). Let \(G\) be a threshold hyperbolic random graph. Then in \(\mathsf{CONGEST}\) after step 1 of HRG-Shattering-MM, the set of vertices \(V \cap \mathcal{B}_0(R- 6\log\log n)\) is matched in constant rounds w.e.h.p.
Proof. In the step 1 of HRG-Shattering-MM, we first activate all vertices with degree at most \(\lceil\log^{3/2} n\rceil\). In parallel, all active vertices mark one incident edge uniformly at random.
Fix a vertex \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) and let \(X_u\) be the random variable with which we count the number of marked edges incident to \(u\). Note that if we can show that \(X_u \geq 1\) for all \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) with probability \(1 - n^{-\omega(1)}\), then this proves the claim since then any vertex in \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) can simply select one of the marked edges so that \(u\) is matched.
To show this, consider the set of vertices which form an \(e\) edge with \(u\) such that \(e\) is potentially marked. Let us write \(N'(u) := N(u) \cap \{v \in V : \deg(v) \leq \lceil\log^{3/2} n\rceil\}\) for this set and note that for \(v \in N'(u)\), the edge \(\{u,v\}\) is marked with probability \(1/\deg(v)\), as \(v\) samples the edge uniform at random among its neighbours. Hence, it holds \[\begin{align} \label{eq:prob-bad-event} \mathbb{P}\left(X_u = 0\right) = \prod_{v \in N'(u)}(1-1/\deg(v)). \end{align}\tag{3}\] In the following, we lower bound \(|N'(u)|\). Since for any vertex \(v \in V \cap \mathcal{L}_{0}\), it holds via 3 that \(\mathbf{E}\left[\deg(v)\right] \in \mathcal{O}(1)\), such that a Chernoff-bound with a subsequent union bound over at most \(\mathcal{O}(n)\) vertices yields that \(\deg(v) \leq \log^{3/2}n\) and all vertices in layer \(\mathcal{L}_{0}\) are active w.e.h.p. and participate in our marking process. Using this bound and that \(r(u) \leq R - 6\log\log n\) in conjunction with the fact that \(\theta_R(\cdot, \cdot)\) is a monotonic decreasing function in both arguments ([rmk:theta-monotonicity]), it holds via law of total expectation that \[\begin{align} \mathbf{E}\left[|N'(u)|\right] \geq (1-o(1)) n\cdot \theta_R(R - 6\log\log n, R)\cdot \mu(\mathcal{L}_{0}) \in \Omega(\log^3(n)), \end{align}\] where in the last step we used ?? and 1. Another combination of Chernoff and union bounds reveals that for any \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) it holds \(|N'(u)| \in \Omega(\log^3 n)\) with probability \(1 - n^{-\omega(1)}\).
Let \(\mathcal{E}\) be the event that for all \(v \in V \cap \mathcal{L}_{0}\) it holds \(\deg(v) \in \mathcal{O}(\log^{3/2}n)\) and that for all \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) it holds \(|N'(u)| \in \Omega(\log^3(n))\). By previous discussion it follows via union bound that \(\mathbb{P}\left(\mathcal{E}\right) \in 1 - n^{-\omega(1)}\). Note that for our fixed \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) it holds \(\mathbb{P}\left(X_u = 0 |\mathcal{E}\right) \in e^{-\Omega(\log^{3/2}n)}\) by 3 . Hence, we obtain using a union bound \[\begin{align} \mathbb{P}\left(\not\exists u \in V \cap \mathcal{B}_0(R- 6\log\log n) : X_u = 0\right) \geq \mathbb{P}\left(\mathcal{E}\right)\cdot(1 - \sum_{u \in V \cap \mathcal{B}_0(R- 6\log\log n) }\mathbb{P}\left(X_u = 0|\mathcal{E}\right)) \in 1-n^{-\omega(1)}, \end{align}\] since \(|V \cap \mathcal{B}_0(R- 6\log\log n)| \in \mathcal{O}(n)\) w.e.h.p. That is, we showed that \(X_u \geq 1\) for all \(u \in V \cap \mathcal{B}_0(R- 6\log\log n)\) with probability \(1 - n^{-\omega(1)}\) as desired. This finishes the proofs since the desired runtime \(\mathcal{O}(1)\) is immediate. ◻
The following is an analogue statement of 6. Additionally, we show that any pair of vertices in the set \(S\) has a minimum angular distance ([item:min95angle]).
Lemma 9 (Step 2 of HRG-Shattering-MM). Let \(G\) be a threshold hyperbolic random graph and consider the set of vertices \(S\) in step 2 of HRG-Shattering-MM. Moreover, let \(c(C, \alpha)\) be a constant large enough and consider annulus \(\mathcal{A} \coloneq \mathcal{B}_0(R - {8\log\log n} + c(C, \alpha))\setminus \mathcal{B}_0(R - {8\log\log n} - c(C, \alpha))\). Then the following holds for \(S\) with probability \(1 - n^{-\omega(1)}\):
All vertices contained in the set \(S\) in the second step of HRG-Shattering-MM are in \(\mathcal{A}\); \(S \subseteq U_{(2)} \subseteq V \cap \mathcal{A}\).
In any sector with angle \(\frac{\log^{30\alpha}n}{n}\) there exists a vertex \(v\) such that \(v \in S\).
For any pair of vertices \(u,v \in S\) it holds that the angular distance is \(\delta_\varphi(u,v) \in \Omega(\log^8 n / n)\) and \(\{u,v\}\not\in E.\)
Proof. Items [item:ring2] and [item:active-sector2] follow directly from 6, since step 2 of HRG-Shattering-MM activates exactly the same set of vertices as step 1 of HRG-Shattering-MM and applies the same local selection rule. Therefore, the arguments used in the proof of 6 carry over verbatim: the set \(S\) is equivalent to \(I_{(1)}\).
To show that [item:min95angle] holds, suppose for contradiction that \(\{u,v\}\in E\). Since both \(u\) and \(v\) belong to \(S\), each has the largest ID in its closed neighbourhood. However, as \(u\in N(v)\) and \(v\in N(u)\), at most one of the two vertices can have the largest ID among the vertices of its closed neighbourhood, a contradiction. Hence \(u\) and \(v\) are not adjacent. Then, using that \(r(u), r(v) \in R - 8 \log\log n +\mathcal{O}(1)\) w.e.h.p. by [item:ring2], it holds w.e.h.p. using 1 in conjunction with [rmk:theta-monotonicity] that \[\delta_\varphi(u,v) > \theta_R(r(u),r(v)) \in \Omega(\log^8 n / n).\] Hence the claimed lower bound on the angular distance holds for every pair of vertices in \(S\) whenever the event of Item [item:ring2] occurs. Since this event holds with probability \(1-n^{-\omega(1)}\), the proof is complete. ◻
We now show that after the second step of HRG-Shattering-MM, a hyperbolic random graph is shattered into components of size at most \(\operatorname{\text{{\rm poly}}}\log n\) for maximal matching.
For a threshold hyperbolic random graph \(G\), it holds with probability \(1 - n^{-\omega(1)}\) that, after the second step of HRG-Shattering-MM, the largest connected component in \(G[V_{(2)}]\) is of size at most \(\log^{\mathcal{O}(1)}(n)\).16 The round complexity of HRG-Shattering-MM is \(\mathcal{O}(1)\) for \(\mathsf{LOCAL}\) and \(\mathcal{O}(\log^3\log n)\) for \(\mathsf{CONGEST}\).
Proof. \(\mathsf{LOCAL}\): We start by proving the \(\mathsf{LOCAL}\) part. Consider any pair of vertices \(s, s' \in S\) in the second step of HRG-Shattering-MM. We show for any of the two sets of vertices \(N_2(s) \cap V_{(1)} =: U(s)\) and \(N_2(s') \cap V_{(1)} =: U(s')\), that for any pair \(v \in U(s)\) and \(v' \in U(s')\) it holds \(\{v,v'\} \not \in E\) w.e.h.p. This ensures that for any pair \(s, s' \in S\), there are no conflicts when \(s\) and \(s'\) compute in parallel \(M(s)\) and \(M(s')\) in the \(\mathsf{LOCAL}\) computation of the second step in HRG-Shattering-MM. To see this, note first that \(\{s,s'\} \not\in E\) w.e.h.p. holds immediately by 9. Thus, it suffices to show that for any \(v \in U(s)\) and \(v' \in U(s')\setminus\{s'\}\) that \(\{v,v'\} \not \in E\) w.e.h.p. Now, recall that by 8 it holds w.e.h.p. for any \(v,u \in V_{(1)}\) that \(r(v),r(u) \geq R - 6\log\log n\). Moreover, by 9 it holds for some constant \(c\) that, \(r(s) \geq R - 8\log\log n -c\) w.e.h.p. Thus, for any vertices \(v, u \in U(s)\) it holds w.e.h.p. that \[\begin{align} \label{eq:two-hop-angle} \delta_\varphi(s,v) \leq \theta_R(r(s), r(u)) + \theta_R(r(v), r(u)) < \Theta(1) \log^{7} n / n, \end{align}\tag{4}\] by 1. Moreover, using 9 [item:min95angle], it holds w.e.h.p. for the pair \(s, s' \in S\) that \[\begin{align} \label{eq:separator-angle} \delta_\varphi(s,s') \in \Omega(\log^8 n / n). \end{align}\tag{5}\] Consequently, we obtain w.e.h.p. for any pair \(v \in U(s)\) and \(v' \in U(s')\setminus\{s'\}\) by combining 4 and 5 that \[\begin{align} \label{eq:two-hop-intersection} \delta_\varphi(v,v') = \delta_\varphi(s,s') - \delta_\varphi(s,v) - \delta_\varphi(s',v') \in \Omega(\log^8 n / n) \in \omega(\theta_R(r(v), r(v'))) \text{ w.e.h.p.,} \end{align}\tag{6}\] using 1 in conjunction with \(r(v') \geq R - 6\log\log n\) w.e.h.p. by 8 and \(r(v) \geq R - 8\log\log n -c\) w.e.h.p. by 9. It follows via 6 that no two-hop neighbourhoods \(U(s') \cap U(s)\) share any edge w.e.h.p. as desired, which, by union bound then also holds for any pair \(s, s' \in S\) w.e.h.p. We conclude that w.e.h.p. there are no conflicts for any matchings \(M(s)\) and \(M(s')\) and it is left to show that (a) we obtain the desired shattering and (b) the computation in HRG-Shattering-MM requires \(\mathcal{O}(1)\) rounds of the \(\mathsf{LOCAL}\) model.
(a) The shattering property follows by the same argument as in the proof of [pro:shattering-mis] after replacing the set \(I_{(1)}\) with \(S\). Indeed, by 8, every vertex in \(V \cap \mathcal{B}_0(R-6\log\log n)\) is already matched after Step 1 of HRG-Shattering-MM w.e.h.p. Moreover, for every \(s \in S\), the algorithm computes a maximal matching on the induced subgraph \(G[U(s)]\), where \(U(s)=N_2(s)\cap V_{(1)}\). Since \(N(s)\subseteq U(s)\), every neighbour of \(s\) is either matched itself or matched to another vertex in \(U(s)\). Consequently, \(N(s)\cap V_{(2)}=\emptyset\) for every \(s\in S\). Hence every unmatched vertex lies in the same outer annulus as in [pro:shattering-mis], and the remainder of the shattering argument applies verbatim.
(b) By 8 step 1 of HRG-Shattering-MM requires constant rounds. The communication in step 2 of HRG-Shattering-MM also requires constant rounds so that active vertices are informed if they participate in the set \(S\). Moreover, each vertex \(s \in S\) in parallel collects its two-hop neighbourhood of unmatched vertices in constant rounds and then, in further constant rounds, informs the vertices of the set \(U(s)\) about the matching \(M(s) \subseteq E(G_2)\). Hence, at most \(\mathcal{O}(1)\) communication rounds are required.
\(\mathsf{CONGEST}\): The correctness of the \(\mathsf{CONGEST}\) part follows from the same arguments as for the \(\mathsf{LOCAL}\) model. For the round complexity, we consider the following: first, for all \(s \in S\), in parallel, each induced subgraph \(G(U(s))\) can be informed by \(s\) to be active within \(2\) rounds in \(\mathsf{CONGEST}\). Next, fix a vertex \(s \in S\). Since there exists a constant \(c\) such that \(r(s) \geq R - 8\log\log n -c\) by 9, it follows by 3 that \(\mathbf{E}\left[\deg(s)\right] \in \Theta(\log^4 n)\). Thus, by a Chernoff and union bound it holds for any \(s \in S\) that \(\deg(s) \in \Theta(\log^4 n)\) w.e.h.p. Moreover, it holds for any \(v \in V_{(1)}\) that \(r(v) \geq R - 6\log\log n\) w.e.h.p. by 8. Subsequently, using 3 with a Chernoff and a union bound, we have for each \(v \in V_{(1)}\) that \(\deg(v) \in \mathcal{O}(\log^3 n)\). Hence, for all \(U(s)\) it holds \(|U(s)|\in \mathcal{O}(\log^7 n)\) w.e.h.p. Then, recall that for any pair of vertices \(s, s' \in S\) we showed that there are no edges between \(U(s)\) and \(U(s')\) w.e.h.p. Thus, using in parallel for each \(s \in S\) on the induced subgraph \(G[U(s)]\) [85], we obtain a maximal matching for all such induced subgraphs \(G[U(s)]\) in \(\mathcal{O}(\log^3\log n)\) rounds without conflicts w.e.h.p. ◻
Proof of [thm:mainpolylog] (MM part). \(\mathsf{LOCAL}\): For MM, we obtain a shattering where each connected component is of size at most \(\mathop{\mathrm{polylog}}n\) w.e.h.p. after \(\mathcal{O}(1)\) using [pro:shattering-mm]. Applying in parallel for each component [13], yields \(\tilde{\mathcal{O}}(\log^{5/3}\log n)\) rounds w.e.h.p..
\(\mathsf{CONGEST}\): The shattering into poly-logarithmic components in \(\mathcal{O}(\log^3\log n)\) rounds follows similarly from [pro:shattering-mm] w.e.h.p. Consequently, by applying in parallel for each connected component the maximal matching algorithm from [85], resulting in a round complexity of \(\mathcal{O}(\log^3\log n)\). ◻
In this section, we prove our lower bounds for MIS and MM on HRGs. We show the following statement.
To accomplish our lower bounds, we first establish that there exist many \(d\)-ary trees in HRGs with relatively large degree \(d\) and height \(h\) ([pro:d-regular-trees] in 6.1). We then use these substructures to obtain our lower bounds of [thm:lowerbound] in 6.2.
In this section, we prove the existence of (polynomially many) \(d\)-ary trees with \(n' \approx \sqrt{\log n}\) vertices, of any degree \(d \leq n'\) growing in \(n\) and height \(h \approx \log(n')/\log(d)\) in an HRG. Specifically, we show the following.
theoremtheoremregulartree Let \(G\) be a threshold hyperbolic random graph and let \(m:= m(n) \leq \sqrt{\log n}\) be a function growing in \(n\). Then, a.a.s., there exist \(n^{\Omega(1)}\) disjoint sectors \(\Phi_i\) such that for \(T_i := G[V \cap \Phi_i]\), \(T_i\) is a balanced \(d\)-ary tree17 with \(n' \in \Omega(\sqrt{\log n})\) vertices, degree \(d \in \Omega(m)\) and height \(h \in \Omega(\log_m\log n')\). In particular, for each such tree \(T_i\) with root \(u\)18, there exists a vertex \(v \in V(H)\) in the giant component \(H\) of \(G\) such that \(\{u,v\} \in E(H)\) is the unique edge between \(T_i\) and \(G \setminus V(T_i)\).
The parameter \(m\) of [pro:d-regular-trees] can be tuned to adjust the degree \(d\) and also implicitly the height \(h\) for a desired \(d\)-ary tree. For example, setting \(m = \sqrt{\log n}\) gives a "star graph" with degree \(\Omega(\sqrt{\log n})\) and height \(1\). For our purposes, we will set \(m = \log\log n\) later on in 6.2.
To prove [pro:d-regular-trees], we partition the hyperbolic disk into \(k\)-many sectors where \(k\) is polynomial in \(n\). For polynomially many such sectors, we reveal that they contain a \(d\)-ary tree with the parameters \(n', d,\) and \(h\). The rough idea of how we find this structure goes as follows. We fix a sector \(\Phi\) with angle \(\phi\) and aim to find a root node \(u\) located in layer \(\mathcal{L}_{\ell_0}\) intersecting sector \(\Phi\) (where \(\ell_0 \approx \log\log n\); see also 3a). Next we consider \(d\) equally sized boxes in a layer \(\mathcal{L}_{\ell_1}\) with \(\ell_1 < \ell_0\), where each box has width \(\phi_{\ell_1}\) and height \(1\). We choose these parameters such that they fulfil the following property: if a vertex \(v\) lies in a box in layer \(\mathcal{L}_{\ell_1}\), then it is adjacent to the root node \(u\) and \(v\) has distance at least \(R\) to any point in layer \(\mathcal{L}_{\ell_1}\) if the angular distance is larger than \(\phi_{\ell_1}\) (see 8 for a visualisation).
For each box in layer \(\mathcal{L}_{\ell_1}\), we find exactly one vertex such that the degree of \(u\) is \(d\) (\(u\) has \(d\) children). Moreover, for any pair of boxes in \(\mathcal{L}_{\ell_1}\), by our parameter choice \(\phi_{\ell_1}\), any pair of “children” of \(u\) has distance at least \(R\) and thus does not share an edge. Then we recursively build the tree by “assigning” any vertex \(v\) in layer \(\mathcal{L}_{\ell_1}\) a set of \(d\) boxes in a layer \(\mathcal{L}_{\ell_2}\), such that each box contains a vertex. We repeat this for each vertex in the boxes in \(\mathcal{L}_{\ell_2}\), obtaining boxes in layer \(\mathcal{L}_{\ell_3}\) and so on until we have boxes in layer \(\mathcal{L}_{\ell_h}\), yielding a tree of height \(h\). To prove that this indeed yields the sought-after \(d\)-ary tree structure, we first show that the graph resulting from the geometric embedding, given that we have exactly one vertex in each assigned box while the rest of the sector is empty, indeed results in a \(d\)-ary tree 10. Afterwards, we show that this event occurs with a small but not polynomial-decaying probability (11). This, in conjunction with the event that all vertices used to build the tree have no further neighbours, occurs for a fixed, slightly larger sector with probability \(n^{-o(1)}\). Since we have polynomial many sectors, there are, in expectation, \(n^{\Omega(1)}\) many sectors where the desired \(d\)-ary tree occurs. This, the bridge of the root to the giant component and the concentration result are addressed in the proof of [pro:d-regular-trees].
To make this formal, we first define for any vertex \(u\) the boxes where we wish to find exactly one vertex each, representing the children of \(u\) (see also 8b for a sketch of the definition).
Definition 1 (Boxes). Let \(\mathcal{D}_R\) be a hyperbolic disk with radius \(R = 2\log n + C\) and \(10000\leq m \leq \sqrt{\log n}\). Moreover, define
(root layer) \(\ell_0:= \lceil \log\log n \rceil\),
(layer of tree level \(i\)) \(\ell_i := \ell_0- \lceil 2i \cdot {\log m} \rceil\),
(box width in tree level \(i\)) \(\phi_{\ell_i} := \frac{100e^{\ell_{i} + C/2}}{n}\),
(height of tree) \(h := \lfloor{\log\log n}/(2\log m) \rfloor\) and
(degree in tree) \(d := \lfloor m/10000\rfloor\).
Consider any \(i \in [h]\) and a vertex \(u \in V_{\ell_i}\). Then for every \(j \in [d]\), the \(j\)-th box of vertex \(u\) is defined by the set of points \[\begin{align} \mathcal{B}_{u, j} := \left\{x \in \mathcal{L}_{\ell_{i+1}} : \varphi(u) + \frac{e^{\ell_i + C/2}}{10 n\cdot m} - (2j+1) \cdot \phi_{\ell_{i+1}} \leq\varphi(x) \leq \varphi(u) +\frac{e^{\ell_i + C/2}}{10 n\cdot m} - 2j \cdot \phi_{\ell_{i+1}} \right\}. \end{align}\]
Before embarking on the proof of our desired geometric properties, let us give some further intuition for our choice of parameters. Placing the "first" box of \(u\) at angular distance \(\approx \frac{e^{\ell_i}}{n\cdot m}\) ensures by 1 that \(u\) "barely" has an edge to a vertex in this box. Then, "moving" the subsequent boxes of \(u\) by angular difference \(\approx \frac{e^{\ell_i}}{n}\) ensures that the hyperbolic distance between \(u\) and any vertex contained in one of its boxes is at most \(R\) as well. Moreover, by “skipping” every second box and by our choice \(\phi_{\ell_i}\) we achieve that any pair of vertices in two different boxes on the same level does not have an edge. The multiplicative constant factors of \(\phi_{\ell_i}\) and \(d\) are chosen for convenience so that these properties are fulfilled. This gives us the desired property of a tree that \(u\) is a “parent” to all vertices in its “own” boxes \(\mathcal{B}_{u, j}\) while vertices that lie in boxes of \(u\) are the “children” of \(u\) and do not have an edge among each other.
To formalise our construction, we consider 9, which defines a nice sector. We shall prove that the construction indeed gives a \(d\)-ary tree in a sector \(\Phi\)19 (see also 8).
Lemma 10 (\(d\)-ary tree geometry). Let \(m \in \omega(1)\) and \(\Phi \subseteq \mathcal{D}_R\) be a sector with angular width \(\phi = \frac{5 \log n}{n\cdot \log m}\). Then, if \(\Phi\) is nice according to 9, the induced sub-graph \(G[V \cap \Phi]\) is a \(d\)-ary tree with degree \(d \in \Omega(m)\), height \(h \in \Omega(\log\log n/ \log m)\) and \(n' \in \Omega(\sqrt{\log n})\) vertices.
Proof. We show the desired statement as follows: first, we show that all boxes are within the sector \(\Phi\). Then, we establish for any vertex \(u \in V \cap \Phi\) that any point in a box of \(u\) has distance at most \(R\) to \(u\). Note that this already suffices to prove that all desired edges for a tree exist such that every vertex has \(d\) children (except for the leaves of the tree). Afterwards, we show that any vertex \(u\) has a distance larger than \(R\) to any point of a box that is not the box of \(u\) according to 1. Since the rest of the sector is empty, this ensures that we also have all non-edges. Using that we have \(h\) layers through which we iterate in line 4 of 9, this yields a tree of depth \(h\) where each vertex has \(d\) children because of the loop in line 8. For the desired properties of our tree, \(d\) and \(h\) follow directly from 1 and the number of vertices is given by \(n' \geq d^{h} \in \Omega(\sqrt{\log n})\).
By 1, a vertex \(u\) in layer \(\mathcal{L}_{\ell_i}\) has boxes with angular distance at most \(\frac{e^{\ell_i + C/2}}{10 n\cdot m}\) since \(d\cdot \phi_{\ell_{i+1}} \leq \exp(C/2 + \ell_i)/(100n\cdot m)\). Hence, the maximal angle that is spanned by any pair of points \(x,y\) in two different boxes is at most \[\delta_\varphi(x,y) \leq \sum_{i=0}^{h-1} \frac{e^{\ell_i + C/2}}{10 n\cdot m} \leq \Theta(1)\cdot \frac{e^{\ell_0}}{n\cdot m} \in o(\phi),\] where we used that \(\ell_i = \lceil \log\log n \rceil - \lceil 2i \cdot {\log m} \rceil\), \(\phi = 5\cdot \log n / (n\cdot \log m)\) and \(m \in \omega(1)\) in the last step. Thus, all boxes are contained in sector \(\Phi\) with angular width \(\phi\).
Recall that \(R = 2\log n + C\). To show that all desired edges exist, consider any vertex \(u \in V_{\ell_{i}}\). Note that any box \(\mathcal{B}_{u, j}\) is in layer \(\mathcal{L}_{\ell_{i+1}}\). Hence, using that \(d\cdot \phi_{\ell_{i+1}} \leq \exp(C/2 + \ell_i)/(100n\cdot m)\) it follows by the definition of \(\mathcal{B}_{u, j}\) that the angular distance between \(u\) and \(x \in \mathcal{B}_{u, j}\) is at most \[\begin{align} \delta_\varphi(u,x) \leq \frac{e^{\ell_i + C/2}}{10 n\cdot m} \leq e^{(R - r(u) - r(x))/2} \leq \theta_R(r(u), r(x)), \end{align}\] using that \(r(u) \leq R - \ell_i\), \(r(v) \leq R - \ell_i + 2\log m + 1\) and 1 in conjunction with \(\theta_R(r(u), r(x))\) being monotonically decreasing by [rmk:theta-monotonicity]. Thus, all desired edges exist.
Next, let \(u \in V_{\ell_i}\) and \(v \in V_{\ell_j}\). W.l.o.g. let \(\ell_i \geq \ell_j\). Then, if \(v\) is not in a box of \(u\), we have to show that \(d_h(u,v) > R\) to finish the proof. To this end, we distinguish two cases.
Case 1 [Vertex \(u\) is an ancestor of \(v\)]: Let \(w\) be the vertex that is in the \((d-1)\)-th box of \(u\), i.e., \(\mathcal{B}_{u, d-1}\) and note that, since \(d\cdot \phi_{\ell_{i+1}} \leq \exp(C/2 + \ell_i)/(100n\cdot m)\), \[\begin{align} \label{eq:lower-bound-angle} \delta_\varphi(u,v) \geq \delta_\varphi(u,w) > (1/10 -1/100)\frac{e^{C/2 + \ell_i}}{m\cdot n} , \end{align}\tag{7}\] using the definition for the box \(\mathcal{B}_{u, j}\) (1). On the other hand, using that \(u\) is an ancestor of \(v\) but \(v\) is not in a box of \(u\), we observe that \(\ell_i \geq \ell_j + 4\log m\). Thus, we obtain \[\begin{align} \label{eq:lower-bound-angle2} \theta_R(r(u), r(v)) \leq \Theta(1) \frac{e^{C/2 + \ell_i}}{m^2 \cdot n}, \end{align}\tag{8}\] by 1. Hence, by combining 7 and 8 we have by \(m \in \omega(1)\) \[\delta_\varphi(u,v) > \Theta(1)\frac{e^{C/2 + \ell_i}}{m\cdot n} > \Theta(1) \frac{e^{C/2 + \ell_i}}{m^2 \cdot n} \in \omega(\theta_R(r(u), r(v))),\] and there is no edge between \(u\) and \(v\) by 1. This concludes our first case.
Case 2 [Vertex \(u\) is not an ancestor of \(v\)]: Let \(w\) be the vertex that is the lowest common ancestor of \(u\) and \(v\). First, consider the case that \(w\) is the parent of both \(u\) and \(v\). In this case, \(w \in V \cap \mathcal{L}_{\ell_{i-1}}\) and as such, the minimal angular distance between \(u\) and \(v\), by 1, is \[\delta_\varphi(u,v) \geq \phi_{\ell_i} = \frac{100e^{\ell_{i} + C/2}}{n} > \theta_R(r(u), r(v)),\] where we applied 1 and \(r(u), r(v) \geq R - \ell_i -2\). Hence, there is no edge between \(u\) and \(v\).
To finalise the case, consider now that \(u\) and \(v\) do not share the same parent such that \(\ell_i > \ell_j\). Let the lowest common ancestor \(w\) be in layer \(\mathcal{L}_{\ell_k}\) and consider the two children \(u'\) and \(v'\) of \(w\), (where both \(u'\) and \(v'\) are in layer \(\mathcal{L}_{\ell_{k+1}}\)), that are ancestors of \(u\) and \(v\) respectively (possibly \(u' = u\)). Note that by \(d\cdot \phi_{\ell_{i+1}} \leq \exp(C/2 + \ell_i)/(100n\cdot m)\) and using 1 it follows that \[\begin{align} \label{eq:case2-part1} \delta_\varphi(u',v') \geq \phi_{\ell_{k+1}} = \frac{100 \cdot e^{\ell_k + C/2 }}{m^2 \cdot n}. \end{align}\tag{9}\] Then, using that \(\varphi(u)\geq \varphi(u')\) and that \(\varphi(v) \leq \varphi(v') + \sum_{t = j}^{k-1} \frac{e^{(C+ \ell_t + \ell_{t -1})/2}}{n}\) by 1, it follows \[\begin{align} \label{eq:case2-part2} \delta_\varphi(u,v) \geq \delta_\varphi(u',v') - \underbrace{\sum_{t = j}^{k-1} \frac{e^{(C+ \ell_t + \ell_{t -1})/2}}{n}}_{\mathcal{O}\left(\frac{e^{\ell_k}}{m^3 n}\right)}. \end{align}\tag{10}\] Plugging 9 into 10 and upper bounding the sum of exponentials by \(\mathcal{O}(\exp{(\ell_k - 3\log m)} /n)\) we obtain via \(m \in \omega(1)\) \[\begin{align} \delta_\varphi(u,v) \geq \frac{(100 - o(1)) e^{\ell_k + C/2}}{m^2 \cdot n} > \theta_R(r(u), r(v)), \end{align}\] using that \(\ell_k \geq \ell_i + 2\log m \geq \ell_j + 4\log m\) and \(r(u) \geq R - \ell_i - 2, r(v) \geq R - \ell_j - 2\) in conjunction with 1. Hence, there is no edge between \(u\) and \(v\) as desired. This concludes the case and the proof. ◻
We now prove that a fixed sector is nice with a small but not polynomial decaying probability.
Lemma 11 (\(d\)-ary tree probability). Let \(m \in \omega(1)\) and \(m \leq \sqrt{\log n}\). Moreover, let \(\Phi \subseteq \mathcal{D}_R\) be a sector with angular width \(\phi = \frac{5 \log n}{n\log m}\). Then, \(\Phi\) is nice according to 9 with probability \(n^{-o(1)}\).
Proof. We show that all boxes of a nice sector according to 9 have exactly one vertex with the desired probability \(n^{-o(1)}\). To this end, fix any \(i \in [h+1]\setminus\{0\}\) and consider any "non-root" layer \(\mathcal{L}_{\ell_i}\). Then, to fulfil the property of a nice sector according to 9, all \(d^i\) boxes include exactly one vertex. Let \(u\) be a vertex in layer \(\mathcal{L}_{\ell_{i-1}}\). Using 1 it holds for any single box in layer \(\mathcal{L}_{\ell_i}\) \[\begin{align} \mu(\mathcal{B}_{u, j}) = \frac{\phi_{\ell}}{2\pi} \cdot \mu(\mathcal{L}_{\ell_i}) = \Theta(1)\frac{e^{\ell_i(1-\alpha)}}{n} \in \Omega(1/n), \end{align}\] where we used ?? and for the last step that \(\ell_i \geq 0\). Thus, using the fact that the number of vertices in a box follows a Poisson point distribution, the probability that there is exactly one vertex in a box in layer \(\mathcal{L}_{\ell_i}\) is \[\begin{align} \label{eq:non-root-prob} p_i = n\cdot \mu(\mathcal{B}_{u, j}) \cdot e^{-n\cdot\mu(\mathcal{B}_{u, j})}\geq e^{\ell_i(1-\alpha)} \cdot \exp\left(-\Theta(1){e^{\ell_i(1-\alpha)}}\right). \end{align}\tag{11}\] Hence, using that any pair of boxes is disjoint, the probability that all \(d^i\) boxes in layer \(\mathcal{L}_{\ell_i}\) have exactly one vertex is given by \[\mathbf{P}\left(\mathcal{E}_{i}\right) = (p_i)^{d^i} \geq \exp\left(-\Theta(1){e^{\ell_i(1-\alpha)}}\right)^{d^i}.\] Next, we bound the probability that the desired "root" vertex in layer \(\mathcal{L}_{\ell_0}\) exists. Since the "box" of the "root" in 9 has angle \(\phi/5 = \log n /n\log(m)\) the probability that exactly one root vertex exists is by a Poisson distribution and ?? \[\begin{align} \label{eq:root-prob} p_0 = \frac{n\cdot\phi}{2\pi} \cdot \mu(\mathcal{L}_{\ell_0})\cdot \exp{\left(-\frac{n\cdot\phi}{2\pi} \cdot \mu(\mathcal{L}_{\ell_0})\right)} \geq \exp\left(-\Theta(1){e^{\ell_0(1-\alpha)}}\right). \end{align}\tag{12}\] since \(m \leq \sqrt{\log n}\) and \(\alpha < 1\). This implies that sector \(\Phi\) has vertices that are required to be nice, i.e., there is the root vertex and every box has exactly one vertex as required in 9, with probability \[\begin{align} \mathbf{P}\left(\mathcal{E}\right) \geq p_0 \cdot \prod_{i=1}^{h} \mathbb{P}\left(\mathcal{E}_{\ell_i}\right) \geq \prod_{i=0}^{h} \left((p_i)^{d^i}\right) \geq \prod_{i=0}^{h} \exp\left(-\Theta(1){e^{\ell_i(1-\alpha)}}\right)^{d^i}, \end{align}\] where we used 11 and 12 in the last step. Next, we use that \(\alpha > 1/2\) which yields \[\begin{align} \label{eq:probability-product} \mathbf{P}\left(\mathcal{E}\right) \geq \prod_{i=0}^{h} \exp\left(-\Theta(1){d^i}{e^{\ell_i/2}}\right). \end{align}\tag{13}\] Since \(\ell_i \leq \log\log n - 2i\cdot\log m + 1\), we obtain \[\begin{align} \; \mathbf{P}\left(\mathcal{E}\right) \geq \prod_{i=0}^{h} \exp\left(-\Theta(1){d^i}\cdot \sqrt{\log n}/m^i\right) \geq \prod_{i=0}^{h} \exp\left(-\Theta(1)\cdot \sqrt{\log n}/10000^i\right), \end{align}\] where we used \(d^i \leq (m/10000)^i\). Applying a geometric series, the above yields \[\begin{align} \; \mathbf{P}\left(\mathcal{E}\right) \geq \exp\left(-\Theta(1)\cdot\sum_{i=0}^{h} \sqrt{\log n}/10000^i\right) \geq e^{-\Theta(1)\sqrt{\log n}} \in n^{-o(1)}, \end{align}\] and we conclude that \(\Phi\) has all vertices that are required to be nice with a probability that is decaying more slowly than polynomial in \(n\).
Finally, we lift this to the statement that the sector is also nice with probability \(n^{-o(1)}\): we obtain a nice sector by showing that the sector, except for the boxes, is empty with essentially the same probability. To see this, note that by our choice of \(\phi \in \mathcal{O}(\log n /(m n))\) and using that the number of vertices in sector \(\Phi\) follow a Poisson distribution with expectation \(\mathbf{E}\left[|V \cap \Phi|\right] = n \phi /2\pi\), the probability that the sector is empty is given by \[\mathbf{P}\left(V \cap \Phi = \emptyset\right) = e^{-\mathbf{E}\left[|V \cap \Phi|\right]} \in e^{-\mathcal{O}(\log n /m)} \in n^{-o(1)},\] using \(m \in \omega(1)\). Hence, given any area \(\mathcal{S} \subseteq \mathcal{D}_R\), the probability that \(\Phi\setminus \mathcal{S}\) is empty is at least \(n^{-o(1)}\) (using the fact that the probability of an area being empty is monotonically increasing if the considered area is shrinking). Then, let \(\mathcal{S}\) be the area that we need to reveal for event \(\mathcal{E}\) and it follows that a sector \(\Phi\) is nice with probability \[\mathbf{P}\left(\Phi \text{ is nice}\right) \geq \mathbf{P}\left(\mathcal{E}\right)\cdot\mathbf{P}\left(V \cap (\Phi\setminus \mathcal{S}) = \emptyset | \mathcal{E}\right) \geq \mathbf{P}\left(\mathcal{E}\right)\cdot\mathbf{P}\left(V \cap \Phi = \emptyset\right) \in n^{-o(1)},\] as claimed. ◻
Next, we “embed” a nice sector in a larger “buffer” sector (see 3b for a sketch). This allows us, given that the “buffer sector” is large enough, to get independent probabilities among nice sectors to be disconnected from any other vertex in the disk. We use this to show that (polynomially) many \(d\)-ary trees exist and that they basically form their own component (except for a “cut edge” that goes from the root vertex of the tree to the giant, which is addressed in the “in particular” part of the statement ).
Proof. Throughout the proof, we partition the disk \(\mathcal{D}_R\) into \(k = \big\lceil \frac{n^{1-\frac{1}{2\alpha}}}{\log n}\rceil\) sectors such that for \(i \in [k]\), the angle \(\psi\) of any such sector \(\Psi_i\) is \(\psi \in \Theta\left(\frac{n^{\frac{1}{2\alpha}}\log n}{n}\right)\); see red sector in 3b for an illustration. Moreover, let \(r^* := R - \log n/\alpha - \log\log n /2\) and throughout the proof, we consider the area \(\mathcal{B}_0(r^*)\); see hatched area in 3.
We prove our theorem in two steps. In the first step of the proof, we show that any sector \(\Psi_i\setminus\mathcal{B}_0(r^*)\) contains a desired \(d\)-ary tree \(T_i\) with probability \(n^{-o(1)}\) if it contains a nice sector \(\Phi_i\); blue sector in 3. This is addressed in Claim [claim:indicator].
Then, in a second step, we show that we have \(n^{\Omega(1)}\) "buffer sectors" \(\Psi_i\) w.e.h.p. that contain a tree \(T_i\) and that if the area \(\mathcal{B}_0(r^*)\) is empty, then the neighbourhood of any \(d\)-ary tree \(T_i\) does not contain a vertex of any vertex of any other "buffer sector" \(\Psi_j\). Showing that \(\mathcal{B}_0(r^*)\) is empty a.a.s. then gives the desired probabilistic guarantee of any \(T_i\) not having any undesired edges. Finally, we finish the proof by showing that the edge of the "root" in \(T_i\) which exists due to Claim [claim:indicator] connects \(T_i\) to the giant in the desired fashion of the "in particular" statement of our theorem.
Now, for Claim [claim:indicator], we consider any \(\Psi_i\) and let \(\Phi_i \subset \Psi_i\) be the sector of width \(\phi = \frac{5\log n}{n\log m} \in o(\log n / n)\) as demanded in 11 that has the same bisector as \(\Psi_i\). In particular we consider a sector \(\Psi_i\) where the corresponding sector \(\Phi_i\) is nice. We introduce the following event \(\mathcal{E}_i\) which says that no vertex in a nice sector \(\Phi_i\) has an edge to a vertex in \(\Psi_i \setminus \Phi_i\) except for the “root” vertex \(u \in V_{\ell_0} \cap \Phi_i\) which has exactly one edge to a "special" vertex \(v\) in \(\Psi_i\setminus \Phi_i\) where \(v\) lies in layer \(\overline{\ell} := \lceil 2\log\log n/(1-\alpha)\rceil\). We later show that \(v\) is part of the giant component.
We define \(\mathcal{E}_i := \mathcal{E}_U \cap \mathcal{E}_{\text{root}}\) and we will show that this event holds with probability \(n^{-o(1)}\).
\(\mathbf{P}\left(\mathcal{E}_i \text{ }|\text{ } \Phi_i \text{ is nice}\right) \in n^{-o(1)}\).
Proof of claim. Event \(\mathcal{E}_U\): We start by showing the desired properties for all "non-root" vertices. Conditioning on the property that \(\Phi_i\) is nice according to 9, the probability for a vertex \(v \in V_{\ell_j}\) in a layer \(\ell_j < \ell_0\) to have no neighbour in \(\Psi_i \setminus \Phi_i\) is at least \[\begin{align} \label{eq:nice1} \mathbf{P}\left(N(v) \cap (\Psi_i \setminus \Phi_i) = \emptyset \text{ }|\text{ } v \in V_{\ell_j}, \Phi_i \text{ is nice}\right) \geq e^{-n\cdot\mu(\mathcal{B}_v(R) \cap \mathcal{B}_0(R))} \geq \exp{\left(-\Theta(1)e^{\ell_j /2}\right)}, \end{align}\tag{14}\] using that we have a Poisson distribution in conjunction with ?? .
Now, given that \(\Phi_i\) is nice, let \(\{v_1, v_2, \dots v_{n' - 1}\}\) be our set \(V \cap (\Phi_i \setminus \mathcal{L}_{\ell_0})=:U\), i.e., vertices in a nice sector except for the “root”, in no particular order. Moreover, let \(\mathcal{E}_t\) be the event that the set of vertices \(\{v_1, v_2,\dots v_{t-1}\}\) have no neighbour in \(\Psi_i \setminus \Phi_i\). Since we use a Poisson distribution, the probability of an area being empty is monotonically increasing if the considered area is shrinking and it holds for any \(t\) that \[\begin{align} \label{eq:nice2} \mathbf{P}\left(v_{t} \text{ has no neighbour in \Psi_i \setminus \Phi_i }| \mathcal{E}_t, \Phi_i \text{ is nice}\right) \geq \mathbf{P}\left(N(v_t) \cap (\Psi_i \setminus \Phi_i) = \emptyset \text{ }|\text{ } v \in V_{\ell_j}, \Phi_i \text{ is nice}\right). \end{align}\tag{15}\] Then, we obtain via 14 and 15 for the event \(\mathcal{E}_U= \{\)No vertex in \(U\) has any neighbour in \(\Psi_i \setminus (\Phi_i \cup \mathcal{B}_0(r^*))\}\) that \[\begin{align} \mathbf{P}\left(\mathcal{E}_U\right) &\geq \prod_{t = 1}^{n'-1}\mathbf{P}\left(v_{t} \text{ has no neighbour in \Psi_i \setminus \Phi_i }| \mathcal{E}_t, \Phi_i \text{ is nice}\right)\\ &\geq \prod_{j = 1}^{h}\prod_{v \in V_{\ell_j} \cap \Phi_i} \mathbf{P}\left(N(v) \cap (\Psi_i \setminus \Phi_i) = \emptyset \text{ }|\text{ } v \in V_{\ell_j}, \Phi_i \text{ is nice}\right) \geq \prod_{j = 1}^{h} \exp{(-\Theta(1)e^{\ell_j /2}})^{d^j}, \end{align}\] where we used in the last step that the number of vertices in layer \(\ell_j\) of a nice sector \(\Phi_i\) is \(d^j\). By the same arguments we applied for 13 , we then get that \[\begin{align} \label{eq:prob-of-U} \mathbf{P}\left(\mathcal{E}_U\right) \in n^{-o(1)}. \end{align}\tag{16}\] This concludes the part for "non-root" vertices.
Event \(\mathcal{E}_{\text{root}}\): Next, we consider the root vertex \(u\) and establish our “cut edge” \(\{u,v\}\) which we later show, connects \(u\) to the giant. Let \(u\) be the vertex in \(\Phi_i\) that lies in layer \(\mathcal{L}_{\ell_0}\), i.e., the root of our tree in the nice sector \(\Phi_i\). W.l.o.g. let \(\varphi(u) = 0\) and consider for \(\overline{\ell} = \lceil 2\log\log n/(1-\alpha)\rceil\) the set of points \[\mathcal{S}= \left\{x \in \mathcal{L}_{\overline{\ell}} : - \frac{e^{C/2}\log^{\frac{1}{2} + \frac{1}{1-\alpha}}(n)}{n} \leq \varphi(x) \leq - \frac{e^{C/2}\log^{\frac{1}{2} + \frac{1}{1-\alpha}}(n)}{2n}\right\} \subset \Psi_i.\] We proceed by setting out the following goals.
We show that any point in \(x\) has distance at most \(R\) to \(u\) while any other vertex in \(\Phi_i\), i.e., \(v \in (V \cap \Phi_i)\setminus \{u\}\) , has distance at least \(R\) to any point in \(\mathcal{S}\).
We show that there is exactly one vertex in \(\mathcal{S}\) with probability \(n^{-o(1)}\).
We show that, under the condition of event \(\mathcal{E}_U\), the probability of the event that \(u\) has no neighbour in \(\Psi_i \setminus (\Phi_i \cup \mathcal{S})\) is \(n^{-o(1)}\).
Note that if all goals are accomplished, then root \(u\) has all desired properties for event \(\mathcal{E}_{\text{root}}\). Moreover, if event \(\mathcal{E}_U\) occurs, there are no undesired edges of “non-root” vertices in \(\Psi_i \setminus \Phi_i\). Thus, since the events of (G2) and (G3) are holding independently with probability \(n^{-o(1)}\) while the event of (G1) is deterministic, we get that the intersection of the event of all goals combined occur with probability \(n^{-o(1)}\) such that \[\begin{align} \mathbf{P}\left(\mathcal{E}_i \text{ }|\text{ } \Phi_i \text{ is nice}\right) = \mathbf{P}\left(\mathcal{E}_{\text{root}} \cap \mathcal{E}_{U} \text{ }|\text{ } \Phi_i \text{ is nice}\right) \geq \mathbf{P}\left(\mathcal{E}_U\right)\cdot \mathbf{P}\left(\mathcal{E}_{\text{root}} | \Phi_i \text{ is nice and event } \mathcal{E}_U \text{ occurs.}\right) \in n^{-o(1)}, \end{align}\] as desired by using that event \(\mathcal{E}_U\) and \(\mathcal{E}_{\text{root}}\) are positively correlated and we previously established \(\mathbf{P}\left(\mathcal{E}_U\right) \in n^{-o(1)}\) in 16 . We continue by accomplishing each goal.
To tackle that \(u\) has distance at most \(R\) to \(x \in \mathcal{S}\), note that \[\delta_\varphi(u,x) \leq \frac{e^{C/2}\log^{\frac{1}{2} + \frac{1}{1-\alpha}}(n)}{n} \leq \frac{e^{(C + \ell_0 + \overline{\ell})/2}}{n}< \theta_R(r(u), r(x)),\] using that \(\ell_0 = \lceil \log\log n \rceil\) and \(\overline{\ell} = \lceil 2\log\log n/(1-\alpha) \rceil\) in conjunction with 1. Thus, \(u\) would have an edge to any potential vertex in \(\mathcal{S}\).
Next, to see that any other vertex in \(\Phi_i\) does not have an edge to a potential \(v \in V \cap \mathcal{S}\), note that any vertex \(w \in V \cap \Phi_i\) has angle \(\varphi(w) > \varphi(u) = 0\), since \(\Phi_i\) is nice. Hence, it holds for any \(x \in \mathcal{S}\), (by \(\ell_j < \ell_0\)), that for any \(w \in V \cap (\Phi_i \cap \mathcal{L}_{\ell_j})\) the angular distance is \[\delta_\varphi(w,x) \geq \frac{e^{C/2}\log^{\frac{1}{2} + \frac{1}{1-\alpha}}(n)}{2n} > \frac{1000 e^{(C + \ell_j + \overline{\ell})/2}}{n}> \theta_R(r(w), r(x)),\] since \(r(w) \geq R - \ell_j - 1 \geq R - \log\log n + 2\log m - 2\), \(m \in \omega(1)\) and \(r(x) \geq R - \overline{\ell} - 1 \geq R - 2\log\log/(1-\alpha) - 2\). Thus, any vertex in \(\Phi_i\), other than the “root” \(u \in V \cap \mathcal{L}_{\ell_0}\), would not have an edge to any vertex in \(\mathcal{S}\) by 1. This concludes the first goal.
For the second goal, we calculate the measure for area \(\mathcal{S}\) and obtain \[\mu(\mathcal{S}) = \mu(\mathcal{\mathcal{L}_{\overline{\ell}}})\cdot \frac{e^{C/2}\log^{\frac{1}{2} + \frac{1}{1-\alpha}}(n)}{4\pi n} \in 1/(n\log^{\Theta(1)} n),\] using the angle spanned by area \(\mathcal{S}\) and using ?? in conjunction with \(\overline{\ell} \in \Theta(\log\log n)\). Hence, the expected number of vertices in \(\mathcal{S}\) is \[\mathbf{E}\left[|V \cap \mathcal{S}|\right] \in 1/\log^{\Theta(1)} n.\] Then, we apply the Poisson distribution and obtain for the probability that the number of vertices in \(\mathcal{S}\) is exactly \(1\) that \[\mathbf{P}\left(|V \cap \mathcal{S}| = 1\right) = \mathbf{E}\left[|V \cap \mathcal{S}|\right] \cdot e^{-\mathbf{E}\left[|V \cap \mathcal{S}|\right]} \in n^{-o(1)},\] as set out by our second goal.
For our third and final goal, note that we have \[\begin{align} \mu(\mathcal{B}_u(R) \cap (\Psi_i \setminus (\Phi_i \cup \mathcal{S}))) \leq \mu(\mathcal{B}_u(R) \cap \mathcal{B}_0(R)) \leq \Theta(1)e^{\ell_0/2}/n, \end{align}\] using that \(u \in V \cap \mathcal{L}_{\ell_0}\) and ?? . Using the fact that the random variable of the number of vertices in \(\mathcal{B}_u(R) \cap \mathcal{B}_0(R)\) follows a Poisson distribution and that \(\ell_0 = \lceil \log\log n \rceil\), we then conclude that \[\mathbf{P}\left( V \cap \mathcal{B}_u(R) \cap (\Psi_i \setminus (\Phi_i \cup \mathcal{S})) = \emptyset|\mathcal{E}_U\right) \geq e^{-n\cdot \mu(\mathcal{B}_u(R) \cap \mathcal{B}_0(R))} \in n^{-o(1)},\] as we sought to show for the third goal, and in conclusion, it follows that \(\mathbf{P}\left(\mathcal{E}_i = 1 \text{ }|\text{ } \Phi_i \text{ is nice}\right) \in n^{-o(1)}\) as desired. ◻
To finish the proof of our statement, recall that \(\Psi_i\) has angle \(\psi \in \Theta\left(\frac{n^{\frac{1}{2\alpha}}\log n}{n}\right)\) and thus, we have \(k \in n^{\Omega(1)}\) such sectors in our entire disk. Now, let \(X_i\) be the indicator random variable that is \(1\) if the event \(\{\mathcal{E}_i \cap \Phi_i\) is nice\(\}\) occurs. Moreover, let \(X:= \sum_{i=1}^k X_i\), i.e., the number of "buffer sectors" \(\Psi_i\) that have the desired induced subgraph tree \(T_i = G[V \cap \Phi_i]\). In the following, we obtain a lower bound for \(X\). Recall that \(r^* = R - \log n / \alpha - \log\log n /2\). Then, using linearity of expectation, 11 and Claim [claim:indicator] we have \[\begin{align} \mathbf{E}\left[X\text{ }|\text{ }V \cap \mathcal{B}_0(r^*) = \emptyset\right] &= \sum_{i=1}^k \mathbf{P}\left(\mathcal{E}_i \cap \Phi_i \text{ is nice} \text{ }|\text{ }V \cap \mathcal{B}_0(r^*) = \emptyset\right) \\ &\geq k\cdot \mathbf{P}\left(\Phi_i \text{ is nice}\text{ }|\text{ }V \cap \mathcal{B}_0(r^*) = \emptyset\right)\cdot\mathbf{P}\left(\mathcal{E}_i \text{ }|\text{ } \Phi_i \text{ is nice}, V \cap \mathcal{B}_0(r^*) = \emptyset\right)\in n^{\Omega(1)} \end{align}\] using that a nice sector according to 9 has only vertices outside the area \(\mathcal{B}_0(r^*)\) and \(\mathcal{E}_i\) is independent of the event \(V \cap \mathcal{B}_0(r^*) = \emptyset\). To obtain concentration for \(X\), note that the sectors \(\Psi_i\) are disjoint and the events that they are nice are mutually independent: determining whether \(\Phi_i\) is nice only requires revealing the randomness within \(\Phi_i\). In similar fashion, event \(\mathcal{E}_i\) also only depends on the randomness in sector \(\Psi_i \setminus \mathcal{B}_0(r^*)\). Therefore, a Chernoff bound applies and we obtain \[\begin{align} \label{eq:X-concentration} \mathbf{P}\left(X \in n^{\Omega(1)} \text{ }|\text{ }V \cap \mathcal{B}_0(r^*) = \emptyset\right) \in 1 - n^{-\omega(1)}. \end{align}\tag{17}\] Next, we show for a sector \(\Psi_i\) where event \(\mathcal{E}_i\) occurs, that all neighbours of the nice sector \(\Phi_i\) are contained in \(\Psi_i\) if we condition on the event that \(V \cap \mathcal{B}_0(r^*) = \emptyset\). Note that if \(\Phi_i\) is nice, it holds for any vertex \(u \in V \cap \Phi_i\) that \(r(u) \geq R - \log\log n -2\). As such, the largest possible angle spanned among a vertex \(u \in V \cap \Phi_i\) and any vertex \(v \in V \cap \mathcal{D}_R\setminus\mathcal{B}_0(r^*)\) is \[\begin{align} \label{eq:angle-neighbour-bound} \theta_R(r(u), r(v)) \leq \theta_R(R - \log\log n - 2, \underbrace{R - \log n/\alpha - \log\log n /2}_{r^*}) \in \mathcal{O}\left(\frac{n^{\frac{1}{2\alpha}}\cdot \log^{3/4}n}{n}\right) \in o(\psi), \end{align}\tag{18}\] by 1. Hence, for any \(u \in V \cap \Phi_i\) it holds that \(N(u) \subseteq V \cap (\Psi_i\setminus \mathcal{B}_0(r^*))\).
In the next step, we bound the probability of the event that the area \(\mathcal{B}_0(r^*)\) is empty (see hatched area in 3b). Note that by 2 the expected number of vertices in this area is \(n\cdot\mu(B_0(r^*)) \in o(1)\). Hence, it holds via Poisson distribution that \[\begin{align} \label{eq:empty-inner-disk} \mathbf{P}\left(V \cap \mathcal{B}_0(r^*) = \emptyset\right) \in e^{-o(1)} \in 1 - o(1). \end{align}\tag{19}\] Thus, using 17 , 19 , and law of total probability we then get that \[\begin{align} \label{eq:X-concentration2} \mathbf{P}\left(X \in n^{o(1)}\right) \leq \mathbf{P}\left(V \cap \mathcal{B}_0(r^*) \neq \emptyset\right) + \mathbf{P}\left(V \cap \mathcal{B}_0(r^*) = \emptyset\right)\cdot \mathbf{P}\left(X \in n^{o(1)} \text{ }|\text{ }V \cap \mathcal{B}_0(r^*) = \emptyset\right) \in o(1). \end{align}\tag{20}\] It follows that, a.a.s., non of the induced subgraph tree \(T_i = G[V \cap \Phi_i]\) has any edge outside \(\Psi_i\) using 18 .
To wrap things up, consider the event that \(X \in n^{\Omega(1)}\) and let \(\mathcal{E}_{\text{giant}}\) be the event that any vertex in layer \(\overline{\ell} =\lceil 2\log\log n/(1-\alpha)\rceil\) is connected to the giant component of \(G\). Note that the intersection of the two events implies our theorem if also \(\mathcal{B}_0(r^*)\) contains no vertex, since event \(\mathcal{E}_i\) ensures that \(T_i\) is connected to the giant component in the desired way of the "in particular" statement of [pro:d-regular-trees].
To bound the probability of event \(\mathcal{E}_{\text{giant}}\) we observe the following: using [42] any vertex \(v \in V \cap \mathcal{L}_{\overline{\ell}}\) is connected to the set of vertices in \(V \cap \mathcal{B}_0(R/2)\) with probability \(1-o(1)\). Moreover, by [54], all vertices of the set \(V \cap \mathcal{B}_0(R/2)\) are part of the giant component with probability \(1-o(1)\). Thus, we conclude that also any \(v \in V \cap \mathcal{L}_{\overline{\ell}}\) is part of the giant component with probability \[\begin{align} \label{eq:giant} \mathbf{P}\left(\mathcal{E}_{\text{giant}}\right) \in 1-o(1). \end{align}\tag{21}\] Though, the three events \(X \in n^{\Omega(1)}\), \(\mathcal{E}_{\text{giant}}\) and \(V\cap \mathcal{B}_0(r^*) = \emptyset\) are not independent, a union bound of the complements yields via 19 , 20 and 21 \[\begin{align} \label{eq:wrap-up} \mathbf{P}\left(X \in n^{\Omega(1)} \cap \mathcal{E}_{\text{giant}} \cap V \cap \mathcal{B}_0(r^*) = \emptyset\right) \in 1 - o(1). \end{align}\tag{22}\] Using 22 then finishes the proof as the desired structural properties of a \(d\)-ary tree in a nice sector now follow from 10. ◻
In this section, we present a lower bound on the runtime of any distributed algorithm for solving MM and MIS on HRGs in the form of [thm:lowerbound]. Our proof strategy is as follows: first we show that an \(r\)-round randomised algorithm for MIS/MM on HRGs would imply an \(r\)-round randomised algorithm for MIS/MM on \(d\)-regular trees; this is addressed in 12. Then, we state a lemma for an established lower bound for MIS/MM on \(d\)-regular trees; see 13 for the statement. By stacking these two lemmas together, we obtain the desired lower bound.
To accomplish our goal, we make use of the following tree substructure in HRGs, which follows from [pro:d-regular-trees] by setting the parameter \(m(n) = \log\log n\).
Corollary 1 (MIS and MM obstruction). Let \(G\) be a threshold hyperbolic random graph. Then, a.a.s. there exists induced subgraph \(T_{\text{hrg}}\) that is a balanced \(d\)-ary tree with degree \(d_{\text{hrg}} \in \Theta(\log\log n)\), height \(h_{\text{hrg}} \in \Theta(\log_d\log n)\) and \(n'_{\text{hrg}} \in \Theta(\sqrt{\log n})\) vertices. Additionally, for each such tree \(T_{\text{hrg}}\) there exists exactly one vertex \(u \in V(T_{\text{hrg}})\) with a neighbour in \(v \in V(G)\setminus V(T_{\text{hrg}})\) where \(v\) is a vertex of the giant component of \(G\). In particular, vertex \(u\) has graph distance \(h_{\text{hrg}}\) to all leaves of \(T_{\text{hrg}}\).
The following lemma shows that if there existed a randomised algorithm that solves MIS/MM in \(o(\log\log n / \log\log\log n)\) rounds on hyperbolic random graphs, then this would imply a randomised algorithm that solves MIS/MM in \(o(\log\log n / \log\log\log n)\) rounds on balanced \(d\)-regular trees with degree \(d \approx \log\log n\) and height \(h \approx \log\log n /\log\log\log n\).
Lemma 12 (HRG-to-tree coupling). Let \(G\sim\mathcal{G}(n, \alpha, C)\) be a threshold hyperbolic random graph with the properties of induced trees \(T_{\text{hrg}}\) as stated in 1 and let \(A_{\text{hrg}}\) be a \(\mathsf{LOCAL}\) algorithm that solves MIS (MM) in \(r < h_{\text{hrg}}/200\) rounds with error probability \(p(n)\) on the giant component of \(G\).
Then, there exist an \(r\)-round \(\mathsf{LOCAL}\) algorithm \(A_{\text{tree}}\), which solves MIS (MM) on a balanced \(d\)-regular tree \(T\) with degree \(d=d_{\text{hrg}} + 1\) and height \(h = \lfloor h_{\text{hrg}}/100 \rfloor\) with error probability at most \(2p(n)\).
Proof. We aim to design an algorithm \(A_{\text{tree}}\) for our infinite family of balanced \(d\)-regular trees with degree \(d\) and height \(h\), such that \(A_{\text{tree}}\) has error probability at most \(2p(n)\). To do so, we run algorithm \(A_{\text{hrg}}\) on the vertices of \(T\), pretending that \(T\) is an HRG \(G\sim\mathcal{G}(n, \alpha, C)\) where \(n\) is suitably chosen such that \(n'_{\text{hrg}}\), \(d_{\text{hrg}}\), and \(h_{\text{hrg}}\) satisfy the conditions of 1. To prove that algorithm \(A_{\text{tree}}\) is well-defined, we consider its behaviour on a \(d\)-regular tree \(T\) with \(n'\) vertices, and also consider the giant component of an HRG \(G\) that contains an induced tree substructure \(T_{\text{hrg}}\) as stated in 1. Note that \(T_{\text{hrg}}\) contains \(T\) as an induced subgraph, and to show that the execution of \(A_{\text{tree}}\) is well-defined, we next consider a mapping of the vertices from \(T\) to \(T_{\text{hrg}}\).
Let \(U\subseteq V(T)\) be the set of vertices that are at a distance at most \(r\) from the leaves of \(T\), and let \(I=V(T)\setminus U\). Let \(N_{r}[U],N_{r}[I]\) be the \(r\)-hop neighbourhood of \(U,I\), respectively. For \(r < h\), we will map each of these subgraphs from \(T\) to \(T_{\text{hrg}}\) such that the \(r\)-hop local view of each vertex \(u\in V(T)\) is the same as the \(r\)-hop local view of a vertex \(v\in V(T_{\text{hrg}})\).
Since \(N_{r}[U]\) contains the vertices up to distance \(r\) from the leaves, if we assume that \(N_{r}[U]\) contains \(k\) leaves, we map the \(k\) leaves of \(N_{r}[U]\) to the \(k\) “leftmost” leaves of \(T_{\text{hrg}}\). After that, we map the rest of the vertices \(N_{r}[U]\) to vertices \(T_{\text{hrg}}\) so that the mapping preserves the \(r\)-hop neighbourhoods of the vertices in \(U\). Let \(M(U)\) be the set of vertices in \(T_{\text{hrg}}\) that \(N_{r}[U]\) was mapped to. With this partial mapping, we have guaranteed that if we run \(A_{\text{tree}}\) on \(T\), every vertex \(u\in U\) has a copy vertex \(v\in V(T_{\text{hrg}})\) with the same view. Additionally, by the choice of \(r<h/2\leq h_{\text{hrg}}/200\), no vertex in \(U\) can identify that they are not in an HRG \(G\).
We will now map the vertices from \(N_{r}[I]\). First, note that the diameter of the induced subgraph \(T':= N_{r}[I]\) is the same as \(T\), if we first remove the leaves of \(T\). Therefore, the induced subgraph of \(N_{r}[I]\) has height \(h-1\). Now let \(w\in V(T')\) be the unique vertex that is at distance \(h-1\) from every leaf of \(T'\). Also, among the vertices of \(T_{\text{hrg}}\) that are part of the partial mapping of \(N_{r}[U]\), let \(v\in M(U)\) be one of the vertices that are furthest away from the leaves in \(T_{\text{hrg}}\). Starting from \(v\), we traverse \(h-1\) hops towards the “root” of the \(T_{\text{hrg}}\), and let \(w'\in V(T_{\text{hrg})}\) be the vertex we arrive at. We map \(w\in V(T')\) to \(w'\in V(T_{\text{hrg}})\). After that, we map the rest of the vertices \(N_{r}[U]\) to vertices \(T_{\text{hrg}}\) so that the mapping preserves the \(r\)-hop neighbourhoods of the vertices in \(I\). Let \(M(I)\) be the set of vertices in \(T_{\text{hrg}}\) that \(N_{r}[I]\) was mapped to. Note that by our choice of \(w\), \(w'\), and \(r\), every vertex in \(M(I)\) is neither a leaf nor a root. This guarantees that the vertices in \(I \subset V(T)\) cannot distinguish that they are not on an HRG \(G\) after the execution of \(A_{\text{tree}}\). Additionally, every vertex in \(N_{r}[I]\) has an identical \(r\)-hop local view to its mapped vertex in \(M(I)\).
We now proceed by bounding the error probability of \(A_{\text{tree}}\). To this end, we utilise the mapping we laid out above. Since MIS and MM are locally checkable problems, if \(A_{\text{tree}}\) fails on \(T\), then there is a local witness of failure centred either at a vertex of \(U\) or at a vertex of \(I\). Denote these two events by \(F_U\) and \(F_I\). Thus, by union bound, we have \[\begin{align} \mathbb{P}\left(A_{\text{tree}}\text{ fails on } T\right) \le \mathbb{P}\left(F_U\right)+\mathbb{P}\left(F_I\right). \end{align}\] We first bound \(\mathbb{P}\left(F_U\right)\). To achieve this, we couple the randomness of the algorithm \(A_{\text{hrg}}\) on the vertices \(N_r[U]\) with the randomness of the algorithm \(A_{\text{tree}}\) on the mapped vertices in \(M(U)\). Since \(r\le h(n)/200\), no vertex in \(M(U)\) “sees” the unique attachment vertex of \(T_{\text{hrg}}\) to \(G\setminus T_{\text{hrg}}\). Hence, the radius-\(r\) views used by \(A_{\text{tree}}\) in \(T_{\text{hrg}}\) and the radius-\(r\) views seen by \(A_{\text{hrg}}\) in \(G\) are identical for all vertices relevant to a failure witness centred in \(U\). Therefore, under this coupling, every occurrence of \(F_U\) yields a local MIS, respectively, MM, violation in the output of \(A_{\text{hrg}}\) on \(G\). Consequently, by the hypothesis of our lemma statement, it holds \[\begin{align} \mathbb{P}\left(F_U\right)\le \mathbb{P}\left(A_{\text{hrg}}\text{ fails on } G\right)\le p(n). \end{align}\] It remains to bound \(\mathbb{P}\left(F_I\right)\). For vertices in \(I\), the relevant \(N_r[I]\) are far from the leaves. We couple the algorithm \(A_{\text{tree}}\) for \(N_r[I]\) with the algorithm \(A_{\text{hrg}}\) for the mapped vertices in \(M(I)\) that are also far from both the leaves and the unique attachment point to \(G\setminus T_{\text{hrg}}\). Thus, all radius-\(r\) views used by \(A_{\text{tree}}\) around a potential failure witness in \(I\) are identical to the corresponding views of \(A_{\text{hrg}}\) in an execution on \(G\). Hence every occurrence of \(F_I\) also yields a local violation of the output of \(A_{\text{hrg}}\) on \(G\), and therefore again \[\begin{align} \mathbb{P}\left(F_I\right)\le \mathbb{P}\left(A_{\text{hrg}}\text{ fails on } G\right)\le p(n). \end{align}\] Combining the two bounds gives \(\mathbb{P}\left(A_{\text{tree}}\text{ fails on } T\right) \le \mathbb{P}\left(F_U\right)+\mathbb{P}\left(F_I\right) \le 2p(n).\) Equivalently, if \(A_{\text{tree}}\) failed with probability greater than \(2p(n)\), then either \(F_U\) or \(F_I\) would occur with probability greater than \(p(n)\), contradicting the correctness guarantee of \(A_{\text{hrg}}\) under the corresponding coupling. ◻
Next, we show a statement which will imply that any randomised algorithm on a balanced \(d\)-regular tree with degree \(d \approx \log\log n\) and height \(h \approx \log\log n / \log\log\log n\) requires \(\Omega(\log\log n / \log\log\log n)\) rounds to solve MIS and MM respectively.
Lemma 13 ([52]). The randomised complexity with error probability \(p\) of MIS/MM on a balanced \(d\)-regular tree with \(n\) vertices is \(\Omega(\min\{d, \log_d n, \log_d \log {1/p}-\mathcal{O}(1)\})\).
Proof sketch:. The result directly follows with 1 (see Appendix 10 for a discussion), given that [52] show that the respective sequences to apply the theorem exist for MIS and MM. More detailed, [52] shows that such a sequence with \(t=\Theta(d)\) and \(f(d)=2^{d+1}\) exists for MIS. The main idea of the proof is that they construct a sequence and each problem in their sequence come with a fingerprint vector \(z\). For MIS the vector is of length \(2\) (formally \(len(z)=1\) in their terminology) and the fingerprint vector of the original MIS problem is \([1,0]\). Then, fingerprint vectors can be obtained iteratively. The \(j\)-entry of the fingerprint for the \(i\)-th problem in the sequence is the \(j\)-prefix sum of the \(i-1\)-st fingerprint. Thus, the fingerprint vectors evolve like \([1,0]\), \([1,1]\),\([1,2]\), until \([1,t]\). They show that the sequence can be extended by yet another problem whenever the \(l1\)-norm of the vector is \(\leq \Delta\). So, for MIS we can do \(\Delta-1\) such steps. The precise interpretation of the fingerprint vector and the actual description of the problems in the sequence is technical and we refer to their paper for more details. We obtain that [52] shows that the sequence for MIS is of length \(\Theta(d)\) with \(f(d)\leq 2^d(1+len(z))\) where \(len(z)=1\) for the MIS problem (in their paper \(len(z)=\beta\) which equals \(1\) for MIS, see Section 5.1 in the arxiv version of their work. Thus, via 1 the runtime bound for error probability \(p\) for MIS simplifies to \(\Omega\left(\min\left\{d,\log_{d}n, \log_{d}\log 1/p-\log_{d}\log 2^{d+1}\right\}\right)\). In particular, the term \(\log_{d}\log 2^{d+1} = \log_{d}({d+1}) \in \mathcal{O}(1)\).
For MM [90] (see e.g. Theorem 4.5) shows that the respective sequence with length \(t=2d-1\) and \(f(d)=5\) exists. For the latter, we cite from their work: Crucially, all the problems of the family are described using only \(5\) labels, and while the result of the round elimination technique may contain more than \(5\) labels, we will provide relaxations that allow us to map these problems back to this family exists, and hence, 1 yields a lower bound of \(\Omega\left(\min\left\{2d-1,\log_{d}n, \log_{d}\log 1/p-\log_{d}\log 5\right\}\right) \in \Omega\left(\min\left\{d,\log_{d}n, \log_{d}\log 1/p-\mathcal{O}(1)\right\}\right)\). ◻
We now put 1, 12 and 13 together to obtain our lower bounds of [thm:lowerbound].
Proof. Using 1 in conjunction with 12, it follows a.a.s. over the draw of \(G\sim\mathcal{G}(n, \alpha, C)\) that any \(o(\log\log n / \log\log\log n)\)-round randomised algorithm for MIS or MM on \(G\) that has an error probability of at most \(p\) would imply the existence of a \(o(\log\log n / \log\log\log n)\)-round randomised algorithm for MIS or MM on balanced \(d\)-regular trees, where the degree is \(d \in \Theta(\log\log n)\) and the tree has height \(h\in \Theta(\log\log n / \log\log\log n)\), with error probability at most \(2p\). Therefore, it suffices to show that no such algorithm exists on any such \(d\)-regular tree with running time \(o(\log\log n / \log\log\log n)\) and error probability at most \(2/n^{c}\), for every constant \(c>0\). Indeed, this immediately implies that a.a.s. over the draw of \(G\), there does not exist an \(o(\log\log n / \log\log\log n)\)-round randomised algorithm on \(G\) for MIS or MM with error probability \(p = 1/n^{c}\). Consequently, to establish our main result, it remains to prove that no such algorithm exists for balanced \(d\)-regular trees.
We show this by using 13. Fix any constants \(c, c_d, c_h > 0\) and consider a balanced \(d\)-regular tree \(T\) with degree \(d = c_d \cdot \log\log n\) and height \(h = c_h \cdot \log_d\log n\) giving a \(d\)-regular tree of size \(n' \geq (d-1)^{h} \geq \log^{c_n} n\) for some constant \(c_n > 0\). Since the runtime for any algorithm with error probability \(p'\) for MIS or MM on a balanced \(d\)-regular tree with \(n'\) vertices is \[\Omega(\min\{d, \log_d n', \log_d \log {1/p'}-\mathcal{O}(1)\})\] using 13, the runtime \(t\) for our tree \(T\) with parameters \(d = c_d \cdot \log\log n\), \(h = c_h \cdot \log_d\log n\) and \(n' \geq \log^{c_n} n\) is \[t\in \Omega\left(\min\left\{\log\log n, \frac{\log\log n}{\log\log\log n},\frac{\log\log n}{\log\log\log n}\right\}\right)\] setting \(p'=2/n^c\). Hence, \(t \in \Omega\left(\frac{\log\log n}{\log\log\log n}\right)\) and the theorem follows as this implies that a.a.s. for \(G\), an algorithm requires at least \(t\) rounds to obtain a probabilistic guarantee of \(1 - n^{-c}\) for MIS or MM for the giant component of \(G\). ◻
In this section, we assume that each vertex knows its geometric position. So a vertex \(v \in V(G)\) can also communicate its geometric coordinate \((r(v), \varphi(v))\) to any neighbour within one round. We prove the following. The Maximal Independent Set part we show in 7.1 and the Maximal Matching part in 7.2. Throughout the section we assume that radial and angular coordinates are unique to vertices; a property that holds almost surely. To describe our algorithms in both sections, we make use of a tiling which ensures that vertices in the same tile form a clique. Moreover, for any two vertices \(u\) and \(v\) in the same layer, if the tiles containing \(u\) and \(v\) are separated by sufficiently many intermediate tiles, then there is no edge between \(u\) and \(v\). This is useful since all vertices in a tile can communicate with each other within one round, finding a local solution that is globally valid if active tiles are far enough apart.
Throughout this section, let \(\ell \in [\lceil \log n / \log\log n \rceil]\). Recall the definition of a layer \[\mathcal{L}_{\ell}:= \mathcal{B}_0(R-\ell)\setminus \mathcal{B}_0(R-\ell-1),\] and for layer level \(\ell\), let \[N_\ell := 40\cdot \left\lceil \frac{n \cdot \pi e^{C/2}}{20e^{\ell}}\right\rceil,\] where the constant \(40\) is chosen for later convenience. Then, for any \(\ell\) and \(i \in [N_\ell]\) we define a tile \(\mathcal{T}_{\ell,i}\) by the set of points \[\begin{align} \label{eq:tile} \mathcal{T}_{\ell,i} := \{x \in \mathcal{L}_{\ell} : (i\cdot 2\pi)/N_{\ell} \leq \varphi(x) < ((i+1)\cdot 2\pi)/N_{\ell} \}. \end{align}\tag{23}\] For a sketch of the tiling, see also 4a. We remark that \(N_\ell\), the number of tiles per layer, is divisible by \(40\).
For any \(\ell\) it holds
\[N_{\ell} \equiv 0 \pmod{40}.\]
We sum up the key properties of our tiling in the following lemma. The first property tells us that vertices in the same tile form a clique, and the second property tells us that vertices in the same layer but with a constant number of tiles between them do not have an edge. Finally, the third property says that if a vertex \(u\) has a neighbour \(v\) which is located in a tile in a “smaller” layer, then \(v\) is not a neighbour to another vertex \(u'\) in the same layer as \(u\), if the number of tiles between \(u\) and \(u'\) is a large enough constant.
Lemma 14 (Tiling properties). The following holds for any \(\ell\) of our tiling:
for any \(i \in [N_\ell]\), it holds for any \(u \in V \cap \mathcal{T}_{\ell,i}\) and \(v \in V \cap \mathcal{T}_{\ell,i}\) that \(\{u,v\} \in E(G)\);
for any pair \(i, j\) such that \(|i-j| \geq 20\), it holds for any \(u \in V \cap \mathcal{T}_{\ell,i}\) and \(v \in V \cap \mathcal{T}_{\ell,j}\) that \(\{u,v\} \not\in E(G)\);
let \(\ell \geq \ell'\) and consider any pair \(i,j\) such that \(|i-j| \geq 40\). Moreover, let \(u \in V \cap \mathcal{T}_{\ell,i}\), \(u' \in V \cap \mathcal{T}_{\ell,j}\) and let \(v \in N(u) \cap \mathcal{L}_{\ell'}\). Then, it also holds \(\{u',v\} \not\in E(G)\).
Proof. [item:one]: W.l.o.g. let \(R - \ell - 1 < r(u) \leq r(v) \leq R - \ell\). Then, using that \(N_\ell = 40\cdot \left\lceil \frac{n \cdot \pi e^{C/2}}{20e^{\ell}}\right\rceil \geq \frac{2 n \cdot \pi e^{C/2}}{e^{\ell}}\), it holds for the angular distance between \(u, v\) that \[\begin{align} \delta_\varphi(u,v) &\leq 2\pi/ N_\ell \leq \frac{e^{\ell}}{ne^{C/2}} \leq \theta_R(R- \ell, R - \ell) \leq \theta_R(r(u), r(v)), \end{align}\] by applying the lower bound of 1 and \(\theta_R(\cdot, \cdot)\) is monotonically decreasing in both arguments (see e.g. [56]). Hence, \(u\) and \(v\) have an edge as desired.
[item:two]: W.l.o.g. let \(R - \ell - 1 < r(u) \leq r(v) \leq R - \ell\). Then, using that there are at least \(19\) tiles in between \(u\) and \(v\) and that \(N_\ell = 40\cdot \left\lceil \frac{n \cdot \pi e^{C/2}}{20e^{\ell}}\right\rceil \leq \frac{2 n \cdot \pi e^{C/2}}{e^{\ell}} + 1\), it holds for the angular distance between \(u, v\) that \[\begin{align} \delta_\varphi(u,v) &\geq 38\pi/ N_\ell \geq \frac{18 e^{\ell}}{ne^{C/2}} > \frac{\pi e^{\ell+1}}{ne^{C/2}} \geq \theta_R(R- \ell - 1, R - \ell - 1) \geq \theta_R(r(u), r(v)), \end{align}\] by applying the upper bound of 1 and the monotonicity of \(\theta_R(\cdot, \cdot)\). Hence, \(u\) and \(v\) have no edge as desired.
[item:three]: Since \(v \in N(u)\), it holds \[\theta_R(r(u), r(v)) \leq \pi e^{(R-r(v) -r(u))/2},\] using 1. W.l.o.g. let \(u\) have angular coordinate \(\varphi(u) = 0\) such that \[\varphi(v) \leq \frac{\pi \cdot e^{\ell+1-C/2}}{n},\] using that \(\theta_R(\cdot, \cdot)\) is monotonically decreasing in both arguments and \(r(u), r(v) \geq R - \ell - 1\).
On the other hand, we have \[\varphi(u') \geq 78\pi/ N_\ell \geq \frac{38 e^{\ell}}{ne^{C/2}} ,\] using that \(\varphi(u) = 0\) and between \(u\) and \(u'\) there are at least \(39\) tiles as \(|i-j|\geq 40\).
Combining these bounds on the angle, we get for the angular distance between \(u'\) and \(v\) that \[\delta_\varphi(u',v) = |\varphi(u') - \varphi(v)| \geq \frac{38 e^{\ell}}{ne^{C/2}} - \frac{\pi \cdot e^{\ell+1-C/2}}{n} > \frac{10 e^{\ell}}{ne^{C/2}} > \frac{\pi \cdot e^{\ell+1-C/2}}{n} > \theta_R(r(u'), r(v)),\] using the monotonocity of \(\theta_R(\cdot, \cdot)\), \(r(u'), r(v) \geq R - \ell - 1\) and 1. Subsequently, \(v \not\in N(u')\) which is what we sought to prove. ◻
We now exploit the knowledge of the coordinates, which implies the knowledge of which tile a vertex belongs to, in order to find an independent set in an HRG. We use the geometry and the coordinates in an algorithm we refer to as Embedding-aware-MIS, for which we first give an informal description.
We iterate through the layers in bottom-up fashion, starting with layer \(\mathcal{L}_{0}\) as follows: we “activate” tiles in the current layer such that “active tiles” are far enough apart from each other so that vertices in different active tiles do not have an edge. Thus, in parallel, we can then select one vertex in every active tile to participate in the independent set (as long as none of its neighbours already participates in the independent set). Since vertices in the same tile form a clique, this removes all vertices from a tile. We then iterate in this fashion through all tiles of a layer such that all vertices in a layer either participate in the independent set or have a neighbour that is part of the independent set. The key observation here is that we only need to iterate through \(\mathcal{O}(\log\log n)\) layers, since after iterating through layer \(\mathcal{L}_{0}\), every vertex with radius at most \(R - 4\log\log n\) has a neighbour that participates in the independent set (this is addressed in 15).
Indeed, we shall show in this section that this procedure, which is formalised in 10, gives a maximal independent set within \(\mathcal{O}(\log\log n)\) rounds.
We first observe that every 40-th tile in layer \(\mathcal{L}_{0}\) contains a vertex that is part of the independent set (if there is at least one vertex in this tile).
If \(|V \cap \mathcal{T}_{0,i}| \geq 1\) such that \(i \equiv 0 \pmod{40}\), then \(V \cap \mathcal{T}_{0,i}\) contains a vertex that is in the set \(I\).
Proof. Throughout the proof, we call a tile that contains at least one active vertex an active tile. Note that the tiles for which we wish to show the desired property are the tiles that are active in Embedding-aware-MIS if \(\ell, m = 0\), i.e., the first active tiles in the algorithm such that \(I = \emptyset\).
We fix any active tile, i.e., a tile \(\mathcal{T}_{0,i}\) such that \(i \equiv 0 \pmod{40}\). Since all vertices that share the same tile form a clique by 14 [item:one], in the Selection step, every active tile with at least one vertex has exactly one vertex \(u\) that tries to be part of the set \(I\). So it is sufficient to show that \(u\) has no active neighbour in any other active tile. By the Activation step and using [rem:divisible], it follows that between two active tiles there are at least \(39\) non-active tiles between two active tiles. Thus, since any other active tile is in the same layer as the tile \(\mathcal{T}_{0,i}\), any vertex \(v\) that is included in any other active tile is not a neighbour of \(u\) by 14 [item:two]. ◻
We make use of this observation to show that, after iterating through layer \(\mathcal{L}_{0}\), vertices with radial coordinate larger than \(R - 4\log\log n\) have a neighbour that is part of the independent set.
Lemma 15 (First iteration of Embedding-aware-MIS). When \(\ell \geq 1\) in Embedding-aware-MIS, any vertex \(u \in V\) with radial coordinate \(r(u) \geq R - 4\log\log n\) has a neighbour that is included in the set \(I\) with probability \(1 - n^{-\omega(1)}\).
Proof. Using Observation [obs:tile-layer0], it is sufficient to show that for any vertex \(u\) with radius \(r(u) \geq R - 4\log\log n\), there exists at least one neighbour \(v \in N(u) \cap \mathcal{T}_{0,i}\) such that \(i \equiv 0 \pmod{40}\).
To show this, fix any vertex \(u\) with radius \(r(u) \geq R - 4\log\log n\) and note that, since \(\theta_R(r(u), R) \in \Omega(\log^2(n)/n)\) (using 1) and the angular width of a tile \(\mathcal{T}_{0,i}\) is at most \(\Theta(1/n)\) (by 23 ), there are at least \(\Omega(\log^2 n)\) tiles \(\mathcal{T}_{0,i}\) such that \(i \equiv 0 \pmod{40}\) in the neighbourhood disk \(\mathcal{B}_u(R)\). Hence, it is sufficient to show that at least one of the \(\Omega(\log^2 n)\) tiles contains any vertex. Since tiles do not overlap and \(\mu(\mathcal{T}_{0,i}) = \mu(\mathcal{L}_{0})/{N_0} \in \Omega(1/n)\) (using ?? and angle \(1/N_0 \in \Theta(1/n)\) of a tile \(\mathcal{T}_{0,i}\)), we conclude that the expected number of vertices in distinct active tiles (tiles \(\mathcal{T}_{0,i}\) tiles such that \(i \equiv 0 \pmod{40}\)) in the neighbourhood disk \(\mathcal{B}_u(R)\) is \(\Omega(\log^2 n)\). Since the number of vertices in each tile follow a Poisson-distribution, a Chernoff-bound yields that there is at least one active tile in the neighbourhood radius of \(u\) with probability \(1 - n^{-\omega(1)}\). Hence,\(u\) has a neighbour \(v \in N(u) \cap \mathcal{T}_{0,i}\) that participates in the set \(I\) with probability \(1 - n^{-\omega(1)}\) and a union bound over all vertices gives the desired result. ◻
Next, we show that Embedding-aware-MIS produces an MIS on the vertices that were not removed by the previous lemma.
Lemma 16 (Embedding-aware-MIS gives MIS). After Embedding-aware-MIS has terminated, any vertex \(u \in V\) with radial coordinate \(r(u) \leq R - 4\log\log n\) is either in the set \(I\) or has a neighbour \(v \in N(u)\) that is in the set \(I\).
Proof. Note that if for all \(\ell \in [\lceil 4\log\log n\rceil + 1]\) it holds that after iteration \(\ell\) all vertices up to radius \(R - \ell\) are either part of the independent set or have a neighbour that is in the independent set, we are done.
We prove this by induction. For \(\ell = 0\), observe that for any \(m\) in 10 that for any pair of active vertices \(u,v\), either (a) \(u\) and \(v\) lie in the same tile, or (b) there are at least \(39\) tiles between \(u\) and \(v\). This follows from the Activation step of Embedding-aware-MIS together with [rem:divisible]. Hence, by 14, every active vertex belongs to a connected component that forms a clique among active vertices in the same tile. Thus, for any active tile, we can select any vertex \(u\) to be included in the \(I\), and all other vertices in the same tile will have \(u\) as a neighbour that is included in the set \(I\). This is achieved by the Selection step in Embedding-aware-MIS and as such, for \(\ell = 0\), all vertices \(V \cap \mathcal{L}_{\ell = 0}\) are either part of the set \(I\), or they have a neighbour that is in the set \(I\) since every tile gets processed exactly once. Thus, after iteration \(\ell = 0\), the set \(I\) is a maximal independent set for \(G[V \cap \mathcal{L}_{0}]\) since any tile in \(\mathcal{L}_{0}\) is exactly once active.
For the induction step, fix any \(\ell\) and suppose that, by the induction hypothesis, after processing layers \(0,\dots,\ell-1\), the set \(I\) is a maximal independent set of \(G\left[V \cap \left(\mathcal{D}_R\setminus \mathcal{B}_0(R-\ell-1)\right)\right].\) Consider iteration \(\ell\) of 10. Any vertex \(u \in V_\ell\) that already has a neighbour in \(I\) is not activated by the algorithm and is therefore already dominated. Hence, it remains to consider the active vertices, namely those in \(V_\ell \setminus N(I)\).
On the induced subgraph \(G[V_\ell \setminus N(I)]\), the same arguments as in the base case apply. For any fixed value of \(m\), active vertices are either contained in the same tile or are separated by at least \(39\) tiles. Thus, by 14, active vertices in different active tiles are non-adjacent, while active vertices within the same tile form a clique. Consequently, the Selection step chooses exactly one active vertex from every non-empty active tile, and every other active vertex in that tile has a neighbour that is added to \(I\).
Therefore, after iteration \(\ell\), every vertex in \(V_\ell\) either already had a neighbour in \(I\) before the iteration or is itself added to \(I\) or has a neighbour that is added to \(I\) during the iteration. Hence, after processing every tile in layer \(\ell\) exactly once, every vertex up to radius \(R-\ell\) is either contained in \(I\) or has a neighbour in \(I\), completing the induction. ◻
The following algorithm gives us a bound on the number of rounds for Embedding-aware-MIS
Lemma 17 (Embedding-aware-MIS Runtime). The algorithm Embedding-aware-MIS requires \(\mathcal{O}(\log \log n)\) rounds w.h.p. for \(\mathsf{CONGEST}\).
Proof. This follows from noting that Embedding-aware-MIS iterates through \(\mathcal{O}(\log\log n)\) layers and \(m \in \mathcal{O}(1)\) in line 4 of 10. Moreover, each iteration there is only one communication step where vertices communicate their angular coordinate in the Communication step. Though angular coordinates are a real number, it suffices to send \(\mathcal{O}(\log n)\) bits w.h.p. by 25 to obtain an ordering that is required for the Selection step. Hence \(\mathcal{O}(1)\) communication rounds per iteration are sufficient w.h.p. for \(\mathsf{CONGEST}\). Thus, we obtain that the number of rounds is \(\mathcal{O}(\log\log n)\) w.h.p. as claimed. ◻
We finish the section by using the established results to show that Embedding-aware-MIS produces a maximal independent set in \(\mathcal{O}(\log\log n)\) rounds.
There exists a deterministic \(\mathcal{O}(\log\log n)\)-round \(\mathsf{CONGEST}\) algorithm to find a Maximal Independent Set on HRGs w.h.p. if the vertices have access to their geometric coordinates.
Proof. Note that set \(I\) fulfils the properties of an MIS and thus, the result that Embedding-aware-MIS produces an MIS follows from 15 and 16. The runtime for \(\mathsf{CONGEST}\) follows from 17. ◻
In this section, we use geometry to find a maximal matching for a hyperbolic random graph in \(\mathcal{O}(\log\log\log n)\) rounds. We do this as follows:
First, in a preprocessing step, we match all vertices with degree at least \(\approx \log^3(n)\) with a vertex in layer \(\mathcal{L}_{0}\). This is accomplished by only marking vertices with degree \(\lceil\log^{3/2} n\rceil\) and letting them draw one random neighbour and mark this incident edge. Thereafter, we consider all vertices with radius at most \(R - 6\log\log n\) and show that all vertices such vertices are matched by this preprocessing step. This is shown in 8.
Thereafter, all that is left to do is to match all unmatched vertices which have radius at least \(R - 6\log\log n\). We partition the outer disk \(\mathcal{D}_R\setminus \mathcal{B}_0(R - 6\log\log n)\) into \(6\log\log n\) layers. We then apply our tiling with the properties shown in 14 to obtain for each induced subgraph of a layer a maximal matching in constant rounds, which can be done for each layer in parallel (see also 18).
Next, we use a divide-and-conquer approach to find a maximal matching for the outer disk. It can be shown that, if we have two annuli \(\mathcal{A}_1\) and \(\mathcal{A}_2\) for which we have each a maximal matching for the induced subgraphs, then we can extent these two matchings to a maximal matching for the induced subgraph of vertices in \(\mathcal{A}_1 \cup \mathcal{A}_2\) in constant rounds. We show this in 19.
Using 18 as a base case, and using that by 19 we can “merge” the matchings of two annuli in constant rounds, we then merge pairs of annuli in parallel in constant rounds. By this, we can double the size of the induced subgraph of an annulus that is matched in constant rounds, and we can do this for all disjoint annuli in parallel. Thus, since the outer disk consists of \(\mathcal{O}(\log\log n)\) layers, after \(\mathcal{O}(\log\log\log n)\) merging iterations, each of which takes constant rounds, we have found a maximal matching for the outer disk in \(\mathcal{O}(\log\log\log n)\) rounds.
The following shows that we can compute a maximal matching for an induced subgraph of a layer in constant rounds.
Lemma 18 (Maximal matching within a layer). Let \(G\) be a threshold hyperbolic random graph where vertices are given their geometric coordinates as an input. Let \(G_\ell := G[V_\ell]\). Then for all layers we can compute a maximal matching \(M_{\ell} \subseteq E(G_\ell)\) for \(G_\ell\) in \(\mathcal{O}(1)\) rounds w.h.p. for \(\mathsf{CONGEST}\).
Proof. We use the following claim, which says that in constant rounds we can compute a matching such that each tile20 contains at most one unmatched vertex21.
For any \(\ell\) let \(i \in [N_\ell]\) and let \(V'\subseteq V\) be a set of unmatched vertices. Moreover, let \(G' := G[V' \cap \mathcal{T}_{\ell,i}]\). Then we can compute a matching \(M' \subseteq E(G')\) with probability \(1 - n^{-10}\) for \(\mathsf{CONGEST}\) in \(\mathcal{O}(1)\) rounds, such that there exists at most one vertex in \(V'\) that is not incident to an edge of \(M'\).
Proof of claim. Consider any tile \(\mathcal{T}_{\ell,i}\) and let \(U:= V' \cap \mathcal{T}_{\ell,i}\). In constant communication rounds, any vertex \(u \in U\) can communicate the tile \(\mathcal{T}_{\ell,i}\) based on its radial and angular coordinate to all its neighbours \(N(u)\). Thus, after constant rounds, every vertex \(u\) has learned all vertices of \(U\) (the set of vertices located in the same tile as \(u\)) since the induced subgraph of a tile forms a clique (14). Next, every vertex \(u\) sends the first \(\lceil 13 \cdot \log n \rceil\) bits of \(\varphi(u)\) to all neighbour. Then, using 25, it holds with probability \(1-n^{-10}\) that every vertex is able to compute an ordering \(v_1 \prec v_2 \prec \dots\) which is the set \(U\), ordered by angular coordinates \(\varphi(v_i)\). Then, consider \(u = v_{j}\) in this ordering. If it has an odd index \(j\) based on our ordering, we add the edge \(\{\{u,v_{j+1}\}\}\) to the matching \(M'\). We now distinguish two cases which show that by this procedure, at most one vertex remains unmatched within a tile.
Case 1 [Number of vertices in a tile \(\mathcal{T}_{\ell,i}\) is even]: all vertices can be paired with any neighbour in the same tile \(\mathcal{T}_{\ell,i}\) using 14. By our procedure of pairing odd-even pairs, every vertex is matched.
Case 2 [Number of vertices in a tile \(\mathcal{T}_{\ell,i}\) is odd]: Observe that the vertex \(u\) with the largest angular coordinate is the only vertex in \(\mathcal{T}_{\ell,i}\) after our pairing procedure that has no matching partner. Thus, there is exactly one vertex which is unmatched.
Since our described procedure requires constant rounds and vertices communicate in any round at most \(\mathcal{O}(\log n)\) bits, the claim follows from the above two cases and that angular coordinates used are unique with probability \(1 - n^{-10}\). ◻
Now, applying Claim [claim:tile-single-vertex] to each tile of layer \(\ell\) in parallel, after constant rounds, there is at most \(1\) unmatched vertex in each tile \(\mathcal{T}_{\ell,i}\). Given this property, we iterate through \(\mathcal{L}_{\ell}\) as follows: for \(t \in [40]\) we attempt to match any possible unmatched vertex \(u \in V \cap \mathcal{T}_{\ell,i}\) where \(i \equiv t \;(\mathrm{mod}\;40)\), i.e., we “activate” every 40-th tile. We do so by first letting every unmatched vertex \(v\) send its layer \(\ell \in \mathcal{O}(\log n)\) to every neighbour \(N(v)\) in constant rounds. Thereafter, we activate any vertex \(u \in V \cap \mathcal{T}_{\ell,i}\) where \(i \equiv t \;(\mathrm{mod}\;40)\). Let \(N'(u)\) be the set of unmatched neighbours of \(u\) and \(u\) marks the set of unmatched vertices \(N'(u) \cap V_{\ell}\). Then \(u\) picks a vertex \(v\) uniform at random from \(N'(u) \cap V_{\ell}\) and {u,v} is added to \(M_{\ell}\). Since in each iteration \(t\), we have for any pair of active vertices \(u \in V \cap \mathcal{T}_{\ell,i}\) and \(u' \in V \cap \mathcal{T}_{\ell,j}\) that \(|i-j|\geq 40\), it follows by 14 that \(N(u) \cap N(u') = \emptyset\). Therefore, no conflict can occur: no unmatched vertex can receive matching requests from two distinct active vertices during the same iteration phase \(t\). Since after \(40\) iterations we tried to match any unmatched vertex in \(\mathcal{L}_{\ell}\) and each iteration requires \(\mathcal{O}(1)\) rounds while also \(t \in \mathcal{O}(1)\), the procedure terminates after constant rounds. Moreover, \(M_{\ell}\) is maximal for \(G_{\ell}\): when a tile \(\mathcal{T}_{\ell,i}\) is processed, the unique unmatched vertex of tile \(\mathcal{T}_{\ell,i}\) is matched whenever it has an unmatched neighbour. Therefore, after tile \(\mathcal{T}_{\ell,i}\) has been processed, either it has become matched or every neighbour remaining in later tiles is already matched (or non-existent). Hence, using that each tile gets processed exactly once, no edge can remain between two unmatched vertices and we obtain a maximal matching \(M_{\ell}\) for \(G_{\ell}\) as desired. ◻
In the following, let \(\mathcal{A}_{i,k} := \bigcup_{\ell=i}^{k+i-1}\mathcal{L}_{\ell}\). We refer to \(\mathcal{A}_{i,k}\) as an annulus of size \(k\) with starting layer \(i\). That is, an annulus of size \(k\) with starting layer \(i\) contains layer \(i\) up to layer \(k-i\) and in total \(k\) layers. Note that an annulus of size \(1\) is equivalent to a layer.
The following lemma says that if we are given two disjoint annuli for which we have a local maximal matching each, then we can enhance such matchings to a maximal matching for the induced subgraph of the two annuli combined in constant rounds (see also 11a for a sketch).
Lemma 19 (Merging two annuli). Let \(G\) be a threshold hyperbolic random graph where vertices are given their geometric coordinates as an input. Let \(i-j \geq k \geq k'\) such that \(U_1:= V \cap \mathcal{A}_{i,k}\), \(U_2:= V \cap \mathcal{A}_{j,k'}\), \(G_1:= G[U_1]\) and \(G_2:= G[U_2]\). If we are given matchings \(M_1 \subseteq E(G_1)\) and \(M_2 \subseteq E(G_2)\) such that \(M_1\) is a maximal matching for \(G_1\) and \(M_2\) a maximal matching for \(G_2\), then in \(\mathsf{CONGEST}\), we can compute in \(\mathcal{O}(1)\) rounds a matching \(M \subseteq E(G_1 \cup G_2)\setminus(M_1 \cup M_2)\) such that \(M_1 \cup M_2 \cup M\) is a maximal matching for \(G[U_1 \cup U_2]\).
Proof. We consider two cases. First, if the two annuli \(\mathcal{A}_{i,k}\) and \(\mathcal{A}_{j,k'}\) have distance larger than \(4\), i.e., for any pair of points \(x \in \mathcal{A}_{i,k}\) and \(y \in \mathcal{A}_{j,k'}\) it holds \(\mathop{\mathrm{d_h}}(x,y) \geq 4\).
Case 1 [For the two annuli it holds \(i-j \geq k + 4\)]: consider the given matchings \(M_1\) and \(M_2\). Under these matchings, all unmatched vertices \(v \in V \cap \mathcal{A}_{j,k'}\) send to any neighbour \(w \in N(v)\) their layer \(\ell\) in constant rounds. Then, if an unmatched vertex \(u \in V \cap \mathcal{A}_{i,k}\) has at least one unmatched neighbour in \(V \cap \mathcal{A}_{j,k'}\), let \(N'(u):= N(u) \cap \mathcal{A}_{j,k'}\) and \(u\) picks uniform at random a neighbour \(v\) from \(N'(u)\) and \(\{u,v\}\) is included in our matching \(M\) such that \(v \in V \cap \mathcal{A}_{j,k'}\). It is left to show that no unmatched vertex \(v \in V \cap \mathcal{A}_{j,k'}\) can receive matching requests from two distinct vertices \(u, u' \in V \cap \mathcal{A}_{i,k}\) in the same round. To see this, recall that under matching \(M_1\), for any unmatched pair \(u\) and \(u'\), there is no edge \(\{u,u'\} \in E\) since \(M_1\) is maximal such that two unmatched vertices cannot be adjacent. Hence, it holds \(\mathop{\mathrm{d_h}}(u,u') \geq R\). Now, for the sake of contradiction, assume that \(v \in N(u) \cap N(u')\). By our case assumption, it holds \(r(u), r(u') \leq r(v) - 4\). By an application of 5, \(v\) has distance at most \(R\) to both \(u\) and \(u'\) if and only if \(\mathop{\mathrm{d_h}}(u,u') \leq R\). A contradiction, and thus, we can obtain the desired matching \(M\) in constant rounds.
Case 2 [For the two annuli it holds \(i-j < k + 4\)]: let \(m := \min(i+3, k)\) and \(m' := \max(j, k'-3)\). We consider the two "border" annuli \[\mathcal{A}_1 := \bigcup_{\ell = i}^m (\mathcal{L}_{\ell}) \subseteq \mathcal{A}_{i,k} \text{ and } \mathcal{A}_2 := \bigcup_{\ell = m'}^{k'} (\mathcal{L}_{\ell}) \subseteq \mathcal{A}_{j,k'}.\]
For a visualisation, see the hatched area in 11b. Note that we can find a maximal matching \(M'\) for the induced subgraph \(G' = (V \cap (\mathcal{A}_{i,k} \cup (\mathcal{A}_{j,k'} \setminus \mathcal{A}_2)), E \setminus (M_1 \cup M_2))\) in constant rounds by the previous case when \(i-j \geq k + 4\). Analogously, we can find a matching \(M''\) for the induced subgraph \(G'' = (V \cap ( \mathcal{A}_{j,k'} \cup (\mathcal{A}_{i,k} \setminus \mathcal{A}_1)), E \setminus (M_1 \cup M_2 \cup M'))\) in constant rounds. Hence, after constant rounds, we can obtain a matching \(M_1 \cup M_2 \cup M' \cup M'' =: \tilde{M}\), such that for a maximal matching \(M \supseteq \tilde{M}\) where \(M \subseteq E(G_1 \cup G_2)\), all edges \(M\setminus \tilde{M}\) are formed among vertex pairs \(u,v\) where \(u \in V \cap \mathcal{A}_1\) and \(v \in V \cap \mathcal{A}_2\).
We then find such a set of edges \(M\setminus \tilde{M}\) in additional constant rounds as follows: we use \(m-i+1\) iterations (iterating layer by layer through \(\mathcal{A}_1\)), where we consider in iteration \(t \in [m-i+1]\) all unmatched vertices of the set \(V \cap \mathcal{L}_{m-t}\). That is, we iterate through the layers of \(\mathcal{A}_1\) in top-down fashion. By using Claim [claim:tile-single-vertex], we can assume that each tile in \(\mathcal{L}_{m-t}\) contains at most \(1\) unmatched vertex after constant rounds. Then, for \(t' \in [40]\), in step \((t, t')\), we activate the set of unmatched vertices where \(u \in V \cap \mathcal{T}_{\ell,b}\) fulfils \(b \equiv t' \;(\mathrm{mod}\;40)\). That is, we activate every \(40\)-th tile. We then match in parallel any such active vertex \(u\) with an unmatched vertex \(v \in V \cap \mathcal{A}_2\) uniform at random from the set of unmatched neighbours of \(u\) in \(\mathcal{A}_2\). This can be done in \(\mathcal{O}(1)\) rounds. Using 14, it holds for any pair of active vertices \(u,u'\) that \(N(u) \cap N(u') \cap \mathcal{A}_2 = \emptyset\) and thus, no conflict occurs by this procedure since any active tile contains at most one unmatched vertex. Since \(\mathcal{O}(40\cdot(m-i+1))\) is constant, this procedure requires \(\mathcal{O}(1)\) rounds. Using the obtained matching by this procedure in conjunction with our matching \(\tilde{M}\) (which we acquired in constant rounds), we obtain our desired matching \(M\) in \(\mathcal{O}(1)\) rounds, which is maximal by similar arguments as used in 18: in iteration \((t,t')\), after processing a tile, its unique unmatched vertex has either been matched or has no unmatched neighbour in \(\mathcal{A}_2\). Since every tile of \(\mathcal{A}_1\) is processed exactly once, no unmatched edge remains between the border annulli \(\mathcal{A}_1\) and \(\mathcal{A}_2\). Together with the maximality of \(\tilde{M}\) on all remaining vertex pairs, this implies that \(M\) is a maximal matching of \(G[U_1 \cup U_2]\) as desired. ◻
We now put everything together to obtain a maximal matching in \(\mathcal{O}(\log\log\log n)\) rounds if the embedding of the graph is given.
There exists a \(\mathcal{O}(\log\log\log n)\)-round \(\mathsf{CONGEST}\) algorithm to find a Maximal Matching on HRGs w.h.p., if the vertices have access to their geometric coordinates
Proof. Match inner-disk. By applying 8, in constant rounds we can find a matching \(M \subseteq E\), such that for any \(u \in V \cap \mathcal{B}_0(R - 6\log\log n)\), there exists an edge \(E(u) \in M\) w.e.h.p. For the rest of the proof we condition on this matching and focus on matching vertices outside \(\mathcal{B}_0(R - 6\log\log n)\) since all vertices in \(V \cap \mathcal{B}_0(R - 6\log\log n)\) are matched in constant rounds.
Match within layers. Next, we consider the \(\lceil 6\log\log n\rceil\) layers of the outer disk \(\mathcal{D}_R\setminus \mathcal{B}_0(R - \lceil 6\log\log n\rceil)\). Using 18, we compute in parallel, for every layer \(\mathcal{L}_{\ell}\) contained in \(\mathcal{D}_R\setminus \mathcal{B}_0(R - \lceil 6\log\log n\rceil)\), a maximal matching in \(G[V_\ell]\) in constant rounds w.h.p.
Match annuli. Now, we finish the proof by induction.
Base case: For each \(i \in [\lceil 3 \log\log n\rceil]\), consider the annuli \(\mathcal{A}_{2i,2} = \mathcal{L}_{2i} \cup \mathcal{L}_{2i+1}.\) That is, each annulus consists of two consecutive layers. Using 19, we compute in parallel, for every \(i \in [\lceil 3 \log\log n\rceil]\), a maximal matching for the annulus \(\mathcal{A}_{2i,2}\) in constant rounds.
Induction step: Consequently, we can repeat this process \(t\) times, doubling the annulus size in each step. More formally, in iteration \(t \in [\lceil \log_2\log\log n + 3 \rceil] \setminus\{0\}\), for \(i \in [\lceil 6\cdot 2^{-t} \log\log n\rceil]\) and \(k = 2^t\) we consider the annuli \(\mathcal{A}_{k \cdot i,k}\). By our induction hypothesis, for all \(i' \in [\lceil 6\cdot 2^{-(t-1)} \log\log n\rceil]\) and \(k' = 2^{t-1}\) we have a maximal matching for all the annuli \(\mathcal{A}_{k' \cdot i',k'}\) after \(\mathcal{O}(t-1)\) rounds. Thus, by an application of 19, we obtain, after an additional constant number of rounds, a maximal matching for all annuli \(\mathcal{A}_{k \cdot i,k}\) after iteration \(t\). That is, we “merge” two “neighbouring” annuli in constant rounds.
Setting \(t= \lceil \log_2\log\log n + 3 \rceil\) we obtain for an annulus of size \(k\geq \lceil 6\log\log n\rceil\) a maximal matching in \(\mathcal{O}(\log\log\log n)\) rounds. Thus, we can find a maximal matching for the outer disk \(\mathcal{D}_R\setminus \mathcal{B}_0(R - \lceil 6\log\log n\rceil)\) in \(\mathcal{O}(\log\log\log n)\) rounds. This, in conjunction with 8, yields a maximal matching for \(G\) in \(\mathcal{O}(\log\log\log n)\) rounds w.h.p. as claimed. ◻
In this section, we show that a standard Luby algorithm requires more than constant rounds so that the degree of every remaining vertex is at most \(n^{o(1)}\) (see [luby-aint-shatter-my-hrg] for a formal statement). Throughout the section, we call a vertex \(v \in V\) a leaf if \(\deg(v) = 1\). We use the following fact of HRGs ([23]).
Lemma 20 (Layer-leaves). Let \(G\sim\mathcal{G}(n, \alpha, C)\) be a threshold hyperbolic random graph and let \(u \in V(G)\) be a vertex in \(G\) with radius \(r \in R - \omega(1)\). Moreover, let \(N_0 := \mathcal{L}_{0} \cap N(u)\) and let \(L_0(u) := \{v \in N_0(u) : \deg(v) = 1 \}\). Then, with non-vanishing probability, \(|L_0(u)| \in \Theta(n\cdot e^{-r/2})\).
Let \(k, d \in \mathbb{Z}^+\) and let \(W = \{w_1, w_2, \dots, w_{k}\}\) be a set of vertices. We say that \(W\) is a similar degree path of length \(k\) if the induced subgraph \(G[W]\) is a path and if the degree of any vertex \(w \in W\) is \(\deg(w) \in \Theta(d)\). Moreover, we say that a vertex \(u \in V\setminus W\) with degree \(\deg(u) = d\) has a similar degree path \(W\) of length \(k\), if \(G[W\cup \{u\}]\) is a path. The following lemma says that most vertices have a similar degree path, where the upper bound on \(r\) guarantees expected degree at least polylogarithmic, while the lower bound ensures that the sectors contain polylogarithmically many vertices.
Lemma 21 (Similar degree path). Let \(G\sim\mathcal{G}(n, \alpha, C)\) be a threshold hyperbolic random graph and let \(u \in V(G)\) be a vertex in \(G\) with radius \(\frac{\log n - \log\log n}{\alpha} \leq r \leq 2\log n - \frac{2\log\log n}{1-\alpha}\). Then, for any constant \(c \in \mathcal{O}(1)\), \(u\) has a similar degree path \(W\) of length \(c\) and degree \(\Theta(n\cdot e^{-r/2})\) w.e.h.p. In particular, the similar degree path \(W\) is contained in the area \(\{x \in \mathcal{D}_R: r - 1/1000 \leq r(x) \leq r \text{ and } \varphi(u) \leq \varphi(x) \leq \varphi(u) + 6c\cdot (1 - 1/1000)e^{R/2-r}\}\).
Proof. Chosen with hindsight, let \(\psi := 2(1 - 1/1000)e^{R/2-r}\) and, w.l.o.g., let \(\varphi(u) = 0\). Then consider the sector \(\Phi := \{x \in \mathcal{D}_R: 0 \leq \varphi(x) \leq 3c \cdot\psi\}\), i.e., the sector that lies counter-clockwise to vertex \(u\), where one ray of \(\Phi\) intersects \(u\). We partition \(\Phi\) into \(3c\) equally sized sub-sectors, such that for all \(i \in [3c+1]\setminus\{0\}\), sector \(\Phi_i\) has angle \(\psi\). Now, let \(\mathcal{S}_i := \Phi_i \cap (\mathcal{B}_0(r)\setminus\mathcal{B}_0(r-1/1000))\). We show that for \(\mathcal{S}:= \bigcup_{i=1}^{3c}\), there exists a subset of vertices \(W \subseteq V \cap \mathcal{S}\), s.t. \(G[W]\) induces the desired similar degree path of \(u\) w.e.h.p. First, we note that by 3, every vertex \(w \in W\cup \{u\}\) has degree \(\mathbf{E}\left[\deg(w)\right] \in \Theta(n\cdot e^{-r/2})\). Since \(r \leq 2\log n - \frac{2\log\log n}{1-\alpha}\), and for every \(w \in W\) it holds \(r - 1/1000 \leq r(w) \leq r\) a Chernoff and union bound reveals for all \(w \in W\) that \(\deg(w) \in \Theta(\deg(u))\) with probability \(1 - n^{-\omega(1)}\) Next, we observe that the expected number of vertices in \(\mathcal{S}_i\) is \[\mathbf{E}\left[|V \cap \mathcal{S}_i|\right] = n \mu(\mathcal{S}_i) = \frac{n\psi}{2\pi}\mu(\mathcal{B}_0(r)\setminus\mathcal{B}_0(r-1/1000)) = \Theta(1)n\cdot e^{R/2-r}\cdot e^{-\alpha(R-r)},\] by our choice of \(\psi\) and 2. Using that hypothesis that \(r \leq 2\log n - \frac{2\log\log n}{1-\alpha}\) and \(R = 2\log n + C\) we then have \(\mathbf{E}\left[|V \cap \mathcal{S}_i|\right] \in \Omega(\log^2 n)\). Using Poisson distribution for each \(\mathcal{S}_i\) and a union bound for \(3c\) areas \(\mathcal{S}_i\), we have that each \(\mathcal{S}_i\) contains a vertex with probability \(1 - n^{-\omega(1)}\). We proceed by showing that, given that each \(\mathcal{S}_i\) contains at least one vertex, \(u\) has a similar degree path \(W\) of length \(c\), which finishes the proof. To this end, we use [89] and by our angular width \(\psi\), it follows for any \(x \in \mathcal{S}_{i}\) and \(y \in \mathcal{S}_{j}\) that \(\mathop{\mathrm{d_h}}(x,y) \leq R\) if \(|i - j| \leq 3\) while \(\mathop{\mathrm{d_h}}(x,y) > R\) if \(|i - j| \geq 5\) since \(r - 1/1000 \leq r(x), r(y) \leq r\). Consequently, if for every \(\mathcal{S}_i\) where \(i\) is divisible by \(3\) we include exactly one vertex \(w \in V \cap \mathcal{S}_i\) in the set \(W\), this ensures that \(G[W]\) induces a path of length \(c\). Moreover, \(u\) has an edge to the unique selected vertex \(w\in V\cap\mathcal{S}_3\), but not to any \(w \in V \cap \mathcal{S}_i\) if \(i\geq 6\). It follows that \(W\) is a similar degree path of length \(c\) of \(u\). Since every \(\mathcal{S}_i\) contains at least one vertex w.e.h.p. and all vertices \(w \in W\) have degree \(\deg(w) \in \Theta(\deg(u))\) w.e.h.p., a union bound over both events finishes the proof, where the ‘in particular’ statement directly follows from our construction of \(W\). ◻
We now show that an HRG has many vertices that have both polynomially many leaves and a similar degree path of constant length. This structure will be useful to show that Luby does not shatter an HRG in constant rounds.
Lemma 22 (Luby Obstruction). Let \(G \sim \mathcal{G}(n, \alpha, C)\) be a threshold hyperbolic random graph and let \(\varepsilon< 1 - 1/(2\alpha)\) be a constant larger than \(0\). Then, a.a.s. there exist a set of vertices \(U(\varepsilon) \subseteq V\) of size \(|U(\varepsilon)| \in \Omega\left(n^{1 - 1/2\alpha -\varepsilon}\right)\), where a vertex \(u \in U(\varepsilon)\) has the following properties.
Degree of \(u\): The degree of \(u\) is \(\deg(u) \in \Theta(n^{\varepsilon/ 2})\).
Leaves of \(u\): Vertex \(u\) has \(\Theta(n^{\varepsilon/ 2})\) neighbours that are leaves.
Similar degree path of \(u\): For any constant \(t \in \mathcal{O}(1)\), \(u\) has a similar degree path \(W(u)\) of length at least \(2t\).
Moreover, for any pair \(u, u' \in U(\varepsilon)\) it holds that for any pair \(v \in L(u) \cup W(u) \cup \{u\}\) and \(v' \in L(u') \cup W(u') \cup \{u'\}\) that \(\{v,v'\} \not\in E(G)\).
Proof. We partition the disk \(\mathcal{D}_R\) into \(\left\lfloor\frac{n\cdot 2\pi}{n^{1/(2\alpha) + \varepsilon}} \right\rfloor =: k\) sectors, (where \(k \geq n^{1 - 1/2\alpha -\varepsilon}\in n^{\Omega(1)}\) due to our choice of \(\varepsilon< 1 - 1/(2\alpha)\)), such that for \(i \in [k]\), a sector \(\Phi_i\) has angle \(\phi \in \Theta(n^{1/(2\alpha) + \varepsilon}/n)\). We consider one such sector \(\Phi_i\) and let \(X_i\) be the indicator random variable that is \(1\) if \(\Phi_i\) contains a vertex \(u\) with properties [item:property1], [item:property2] and [item:property3]. We first show that \(\mathbf{P}\left(X_i =1\right) >0\). Thereafter, we show for \(X:= \sum_{i=1}^k X_i\) that \(X \in \Omega(k)\) a.a.s. what proves our desired statement.
Without loss of generality, let \(\Phi_i\) have bisector \(0\). Then, for \(\hat{\ell} := \lceil \varepsilon\cdot \log n \rceil\), let \(\mathcal{S} := \{x \in \mathcal{L}_{\hat{\ell}}: - n^{\alpha \cdot \varepsilon}/n < \varphi(x) < n^{\alpha \cdot \varepsilon}/n\} \subset \Phi_i\). Hence, using ?? , we obtain \[\mathbf{E}\left[|V\cap \mathcal{S}|\right] = n\cdot\mu(\mathcal{S}) = \Theta(1)\cdot n^{\alpha\cdot \varepsilon}\cdot e^{-\alpha\hat{\ell}} \in \Theta(1).\] Thus, since \(|V\cap \mathcal{S}|\) follows a Poisson-distribution, we get \[\begin{align} \label{eq:Prob-event1} \mathbf{P}\left(\mathcal{E}\right) >0, \qquad \mathcal{E} := \{|V \cap \mathcal{S}|=1\}. \end{align}\tag{24}\] Next, condition on event \(\mathcal{E}\) and let let \(u\) denote the unique vertex in \(V\cap\mathcal{S}\). Then it follows \[\begin{align} \label{eq:Prob-event2} \mathbf{P}\left(\mathcal{E}_1 | \mathcal{E}\right) \in 1 - n^{-\omega(1)}, \qquad \mathcal{E}_1 := \{\deg(u) \in \Theta(n^{\varepsilon/ 2})\}, \end{align}\tag{25}\] by ?? and a Chernoff-bound since \(u \in V_{\hat{\ell}}\) where \(\hat{\ell} = \lceil \varepsilon\cdot \log n \rceil\).
Let \(|L_0(u)| := |\{v \in N_0(u) : \deg(v) = 1 \}|\) as in 20. Then, using 20 it follows \[\begin{align} \label{eq:Prob-event3} \mathbf{P}\left(\mathcal{E}_2 | \mathcal{E}\right) >0, \qquad \mathcal{E}_2 := \{|L_0(u)| \in \Theta(n^{\varepsilon/ 2})\}, \end{align}\tag{26}\] since \(u\) is in layer \(\hat{\ell} = \lceil\varepsilon\log n\rceil\) and thus, \(r(u) \leq R - \lceil\varepsilon\log n \rceil\).
Applying 21, it follows that \[\begin{align} \label{eq:Prob-event4} \mathbf{P}\left(\mathcal{E}_3 | \mathcal{E}\right) \in 1 - n^{-\omega(1)}, \qquad \mathcal{E}_3 := \{|W(u)| \geq 2t\}, \end{align}\tag{27}\] and we remark that for event \(\mathcal{E}_3\), the vertices of \(W(u)\) are within an area \(\mathcal{S}' := \{x \in \mathcal{D}_R: R- \hat{\ell} - (1 + 1/1000) \leq r(x) \leq R - \hat{\ell} \text{ and } \varphi(u) \leq \varphi(x) \leq \varphi(u) + 6t\cdot (1 - 1/1000)e^{-R/2 +\hat{\ell} + 1 + 1/1000}\}\) by the ‘in particular’ statement in 21. Thus, \(\mathcal{S} \subset \Phi_i\), since \(\Phi_i\) has an angle of \(\Theta(n^{1/(2\alpha) + \varepsilon}/n)\) while \(\mathcal{S}'\) spans an angle \(\Theta(n^{\varepsilon}/n) \in o(n^{1/(2\alpha) + \varepsilon}/n)\), using that \(R = 2\log n +C\), \(t \in \mathcal{O}(1)\) and \(\hat{\ell} = \lceil \varepsilon\log n \rceil\). Our desired ‘moreover’ statement then follows since for any pair of sectors \(\Phi_i\) and \(\Phi_j\) where \(\mathcal{E} \cap \mathcal{E}_1 \cap \mathcal{E}_2 \cap \mathcal{E}_3\) occur, we consider \(u \in V \cap \Phi_i\) and \(u' \in V \cap \Phi_j\) to be the vertices occurring due to event \(\mathcal{E}\) and it holds for any pair \(v \in L(u) \cup W(u) \cup \{u\}\) and \(v' \in L(u') \cup W(u') \cup \{u'\}\) that \(\{v,v'\} \not\in E(G)\). This is since every object associated with \(u\) that is contained in \(\Phi_i\), and distinct sector \(\Phi_j\) are separated by an angular distance exceeding the maximal adjacency angle of 1 for any pair of vertices \(v \in L(u) \cup W(u) \cup \{u\}\) and \(v' \in L(u') \cup W(u') \cup \{u'\}\); there are no edges between distinct structures belonging to different sectors. Subsequently, if \(X \in \Theta(k)\) a.a.s. this finishes the proof.
To obtain this result, note first that \[\mathbf{P}\left(X_i = 1\right) = \mathbf{P}\left(\mathcal{E} \cap \mathcal{E}_1 \cap \mathcal{E}_2 \cap \mathcal{E}_3\right).\]
Then, using 24 , 25 , 26 and 27 in conjunction with a union bound and conditional probabilities, we have that \[\mathbf{P}\left(\mathcal{E} \cap \mathcal{E}_1 \cap \mathcal{E}_2 \cap \mathcal{E}_3\right) \geq \left(1 - \mathbf{P}\left(\mathcal{E}_1^C | \mathcal{E}\right) - \mathbf{P}\left(\mathcal{E}_2^C | \mathcal{E}\right) -\mathbf{P}\left(\mathcal{E}_3^C | \mathcal{E}\right)\right)\mathbf{P}\left(\mathcal{E}\right) > 0,\] and we conclude that \(\mathbf{E}\left[X\right] \in \Theta(k) \in n^{\Omega(1)}\).
We wrap up the proof as follows: we show for any \(i \neq j\) that \(X_i\) and \(X_j\) are independent so that a Chernoff bound yields \(X \in \Omega(k)\). To this end, we show that events \(\mathcal{E}, \mathcal{E}_1, \mathcal{E}_2,\) and \(\mathcal{E}_3\) in a sector \(\Phi_i\) only depends on the randomness of the Poisson point process in \(\Phi_i\).
For event \(\mathcal{E}\) this is immediate since \(\mathcal{S} \subset \Phi_i\).
For event \(\mathcal{E}_1\), note that by 19 there is no vertex in \(\mathcal{B}_0(r^*)\) a.a.s. where \(r^* = R - \log n/\alpha - \log\log n /2\). Thus, since \(u \in V_{\hat{\ell}}\) and \(\hat{\ell} = \lceil \varepsilon\log n \rceil\), the sector that we need to reveal to obtain the degree of \(u\) has angle at most \[n^{\alpha \cdot \varepsilon}/n + 2\theta_R(r^*, R - \hat{\ell} -2) \in \mathcal{O}\left(e^{(\hat{\ell} - r^* )/2}\right) \in o(\phi) \text{ a.a.s.,}\] using \(\phi \in \Theta(n^{1/(2\alpha) + \varepsilon}/n)\) and 1 in conjunction with [rmk:theta-monotonicity]. That is, a smaller angle than the angle spanned by \(\Phi_i\).
By a similar argument we get for \(\mathcal{E}_2\), using that a leaf \(v \in L_0(u)\) has radius \(r(v) \geq R - 1\) and thus \(\delta_\varphi(u,v) \leq \theta_R(R- \hat{\ell} - 2, r(v)) \in \mathcal{O}(n^{-\varepsilon/2})\), a.a.s. that we do not require to reveal a sector with an angle that is larger than \(\mathcal{O}(n^{-\varepsilon/2}) + \theta_R(R-1, r^*) \in o(\phi)\) using 1.
For event \(\mathcal{E}_3\), recall that we remarked that the vertices of \(W(u)\) are within an area \[\mathcal{S}' := \{x \in \mathcal{D}_R: R- \hat{\ell} - (1 + 1/1000) \leq r(x) \leq R - \hat{\ell} \text{ and } \varphi(u) \leq \varphi(x) \leq \varphi(u) + 6t\cdot (1 - 1/1000)e^{-R/2 +\hat{\ell} + 1 + 1/1000}\}.\] Thus, a.a.s., we need to reveal a sector of angle at most \(\mathcal{O}\left(e^{-R/2 +\hat{\ell}}\right) + \theta_R(r^*, R - \hat{\ell} -2)\). Using 1, \(\hat{\ell} = \lceil \varepsilon\log n \rceil\) and \(r^* = R - \log n/\alpha - \log\log n /2\), this angle is again \(o(\phi)\) and thus, also completely contained within sector \(\Phi_i\).
Finally, using that \(\mathbf{E}\left[X\right] \in \Theta(k) \in n^{\Omega(1)}\), it follows by a Chernoff bound \[\mathbf{P}\left(\mathcal{B}_0(r^*) \cap V = \emptyset\right) \mathbf{P}\left(X \in \Omega(k) | \mathcal{B}_0(r^*) \cap V = \emptyset\right) \in (1-o(1)) \cdot (1-e^{-\Omega(k)}).\] That is, \(|U(\varepsilon)| \in \Omega(k) \in \Omega\left(n^{1 - 1/2\alpha -\varepsilon}\right)\) a.a.s. as desired. ◻
We now use 22 to show that for any constant \(t \in \mathcal{O}(1)\) number of iterations of Luby’s algorithm, there remains a vertex \(u\) with polynomial degree a.a.s. To formalise this, we write \(V_{(t)}\) for the set of vertices that are not removed by Luby after iteration \(t\), let \(G_t = [V_t]\) and \(\Delta(G_t) := \max_{v \in V_t}\deg_t(v)\) where for \(v \in V_t\) we used \(\deg_t(v):= |\{u \in N(v) : u \in V_t\}|\).
Let \(G \sim \mathcal{G}(n, \alpha, C)\) be a threshold hyperbolic random graph and let \(t \in \mathcal{O}(1)\) be any fixed constant. Then \(\Delta(G_t) \in n^{\Omega(1)}\) a.a.s.
Proof. With hindsight, we set \(\varepsilon:= \frac{1 - \frac{1}{2\alpha}}{2t(t+2)}\) and consider the set of vertices \(U(\varepsilon)\) of 22. We show that after \(t\) rounds there exists a.a.s. a vertex \(u \in U(\varepsilon)\) such that \(\deg_t(u) \in n^{\Omega(1)}\).
To this end, for \(i \in [t]\), let \(L_{(i)}(u) := \{v \in N(u) \cap V_{(i)} : \deg(v) = 1\}\) (the “remaining leaves” of \(u\) after round \(i\)), \(W_{(i)}(u):= \left(\bigcup_{j=1}^{2(t-i)}\{w_j\}\right)\cap V_{(i)}\)22 (the “remaining” similar degree path of \(u\) after round \(i\)) and we define the random variable \[\begin{align} K_{(i)} := |\{u \in U(\varepsilon) \cap V_{(i)} : |L_{(i)}(u)| \in n^{\Omega(1)} \text{ and } |W_{(i)}(u)| \geq 2(t-i)\}|. \end{align}\] Note that if \(K_{(t)} \geq 1\) a.a.s. this is sufficient to prove our desired statement and that \(K_{(0)} \in \Omega\left(n^{1 - 1/2\alpha - \frac{1 - \frac{1}{2\alpha}}{2t(t+2)}}\right)\) a.a.s. due to 22 and our choice of \(\varepsilon\). In order to prove \(K_{(t)} \geq 1\) a.a.s., we will use the following claim, which tells us that if \(K_{(i)}\) is polynomial in \(n\), then so is \(K_{(i+1)}\) a.a.s. given that \(i < t\).
Let \(\varepsilon:= \frac{1 - \frac{1}{2\alpha}}{2t(t+2)}\). Then, for \(i \in [t]\) it holds \[\mathbb{P}\left(K_{(i+1)} \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-(i+1)\varepsilon(2t+1)} \right) | K_{(i)} \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-i\varepsilon(2t+1)} \right)\right) \in 1 - o(1).\]
Proof of claim. Let \(U_{(i)} = \{u \in U(\varepsilon) \cap V_{(i)} : |L_{(i)}(u)| \in n^{\Omega(1)} \text{ and } |W_{(i)}(u)| \geq 2(t-1)\}\), i.e., \(|U_{(i)}| = K_{(i)}\). We analyse one round of Luby on \(U_{(i)}\). To this end, we fix a vertex \(u \in U_{(i)}\) and consider the following events:
Vertex \(u\) draws a random number for Luby that lies in the open interval \((1- \frac{2t + 1}{n^\varepsilon}, 1- \frac{2t}{n^\varepsilon})\).
For \(j \in [2(t-i)+1 ]\setminus\{0\}\), vertex \(w_j \in W(u)\)23 draws a random number for Luby that lies in the open interval \((1- \frac{2t + 1 -i}{n^\varepsilon}, 1- \frac{2t- i}{n^\varepsilon})\).
Any vertex \(v \in N((W(u) \cup \{u\}) \setminus (W(u)) \cup \{u\})\) draws a random number for Luby that lies in the open interval \((0, 1- \frac{2t + 1}{n^\varepsilon})\).
Observe that if event \(\mathcal{E}_u \cap \mathcal{E}_{\text{path}} \cap \mathcal{E}_{\text{leaves}} =: \mathcal{E}\) occurs, then \(u \in U_{(i+1)}\). To see this, note that all leaves of \(u\) draw a smaller number than \(u\) and \(u\) is not removed as it is the largest number among neighbours, except for \(w_1 \in N(u)\), which does not join the independent set since \(w_2 \in N(w_1)\) draws a larger number than \(w_1\). This “chain” of events continues for the entire similar degree path \(W_{(i)}(u)\) by the intersection of the two events \(\mathcal{E}_{\text{path}}\) and \(\mathcal{E}_{\text{leaves}}\). Hence, only the "terminal" vertex \(w_{2(t-i)}\) of \(W_{(i)}(u)\) that might join the independent set, removing all its neighbours including the neighbour \(w_{2(t-i-1)}\) but no other vertex of the similar degree path \(W_{(i)}(u)\). Since \(u\) has polynomially many leaves, all necessary conditions for \(u \in U_{(i+1)}\) are fulfilled. We continue by lower bounding the probability of event \(\mathcal{E}\).
Since the range of the open interval is \(n^{-\varepsilon}\) we have \[\begin{align} \label{eq:randomdraw1} \mathbb{P}\left(\mathcal{E}_u\right) = n^{-\varepsilon}. \end{align}\tag{28}\]
Similar the event \(\mathcal{E}_u\), the range of the interval of each random number is \(n^{-\varepsilon}\). Since each draw is independent and we draw at most \(2t\) random numbers, we obtain \[\begin{align} \label{eq:randomdraw2} \mathbb{P}\left(\mathcal{E}_{\text{path}}\right) \geq n^{-2t\varepsilon}. \end{align}\tag{29}\]
Since for any \(w \in W(u) \cup \{u\}\) we have \(\deg(w) \in \mathcal{O}(n^{\varepsilon/2})\), we can upper bound \(m := |N((W(u) \cup \{u\}) \setminus (W(u)) \cup \{u\})| \in \mathcal{O}(n^{\varepsilon/2})\) using that \(|W(u)|\) is constant since \(t \in \mathcal{O}(1)\). Applying this with the fact that the desired event for every vertex \(v\) of drawing a random number within the interval \((0, 1- \frac{2t + 1}{n^\varepsilon})\) has probability \(1- \frac{2t + 1}{n^\varepsilon}\), we obtain via independence among the events for each \(v\) \[\begin{align} \label{eq:randomdraw3} \mathbb{P}\left(\mathcal{E}_{\text{leaves}}\right) \geq \left(1- \frac{2t + 1}{n^\varepsilon}\right)^{m} \in \Omega(1), \end{align}\tag{30}\] where we used that \(t\) is constant and \(m \in \mathcal{O}(n^{\varepsilon/2})\).
Putting 28 , 29 and 30 together and using independence we then obtain \[\begin{align} \label{eq:final-prob} \mathbb{P}\left(\mathcal{E}\right) = \mathbb{P}\left(\mathcal{E}_u \cap \mathcal{E}_{\text{path}} \cap \mathcal{E}_{\text{leaves}}\right) \in \Omega(n^{-(2t +1)\varepsilon}). \end{align}\tag{31}\] To finish the proof of our claim, let \(X_u\) be the indicator random variable that event \(\mathcal{E}\) occurs for vertex \(u \in U_{(i)}\) and define the random variable \(X := \sum_{u\in U_{(i)}} X_u \geq K(i+1)\). We then get by linearity of expectation \(\mathbb{E}\left[X\right] \geq \mathbb{P}\left(\mathcal{E}\right)\cdot |U_{(i)}| \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-(i+1)\varepsilon(2t+1)} \right)\) by the assumption that \(K_{(i)} \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-i\varepsilon(2t+1)} \right)\) and 31 . Since the random variables \(X_u\) are independent by the ‘moreover’ statement of 22, a Chernoff bound gives the desired result that \(K_{(i+1)} \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-(i+1)\varepsilon(2t+1)} \right)\) a.a.s. ◻
Now, using that \(K_{(0)} \in \Omega\left(n^{1 - 1/2\alpha - \frac{1 - \frac{1}{2\alpha}}{2t(t+2)}}\right)\) a.a.s. by 22, and applying Claim [claim:appendix]iteratively for \(t \in \mathcal{O}(1)\) rounds, a union bound over the \(t\) induction steps yields \[\begin{align} \mathbb{P}\left(K_{(t)} \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-(t+1)\varepsilon(2t+1)} \right)\right) &\geq 1 - \mathbb{P}\left(K_{(0)} \not \in \Omega\left(n^{1 - 1/2\alpha - \frac{1 - \frac{1}{2\alpha}}{2t(t+2)}}\right)\right)\\ &- \mathbb{P}\left(\bigcup_{i=0}^{t-1} K_{(i+1)} \not \in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-(i+1)\varepsilon(2t+1)} \right) | K_{(i)}\in \Omega\left(n^{1 - 1/2\alpha - \varepsilon-i\varepsilon(2t+1)} \right)\right)\\ &\in 1- o(1). \end{align}\] This finishes the proof since by our choice \(\varepsilon:= \frac{1 - \frac{1}{2\alpha}}{2t(t+2)} > 0\), \(K_{(t)} \in n^{\Omega(1)}\) a.a.s. ◻
Lemma 23 (Chernoff bound). For \(i \in [k]\), let \(X_i \in \{0,1\}\) be independent random variables and \(X = \sum_i X_i\). Then \[\begin{align} \mathbb{P}\left(X\geq \frac{3}{2}\mathbb{E}\left[X\right]\right) \leq \exp{\left(-\frac{\mathbb{E}\left[X\right]}{12}\right)} \text{ and } \mathbb{P}\left(X\leq \frac{1}{2}\mathbb{E}\left[X\right]\right) \leq \exp{\left(- \frac{\mathbb{E}\left[X\right]}{8}\right)}. \end{align}\]
The following is a Chernoff24 type deviation bound (see e.g. [39]) which is necessary for the distribution of vertices as it follows a Poisson point distribution.
Lemma 24 (Poisson Chernoff bound). Let \(X\) have a Poisson distribution with mean \(\mathbf{E}\left[X\right]\). Then for \(a \geq e^{3/2}\), \[\begin{align} \mathbf{P}\left(X\leq \frac{1}{2}\mathbf{E}\left[X\right]\right) \leq \exp{(- \mathbf{E}\left[X\right]/8)} \text{ and } \mathbf{P}\left(X\geq a\cdot\mathbf{E}\left[X\right]\right) \leq \exp{(-a\cdot \mathbf{E}\left[X\right] /2)}. \end{align}\]
We extend Theorem 7.1 from [52] that can provide lower bounds for randomised algorithms in the \(\mathsf{LOCAL}\) model for problems that meet the prerequisites of the theorem. We emphasise that we do not claim any major contribution here; we only change their theorem to work with a general error probability instead of a hard-coded error probability of \(1/n\). Thus, we refer to [52] for more background and formal definitions, and here we focus on the differences in the proof. In their work, graph problems are modelled using the black-white formalism, where a locally checkable problem \(\Pi = (\Sigma, \mathcal{N}, \mathcal{E})\) is defined by a label space \(\Sigma\), node (white) constraints \(\mathcal{N}\), and hyperedge (black) constraints \(\mathcal{E}\) [52]. The goal in some computational models, like the \(\mathsf{LOCAL}\) model, is to compute a labelling that satisfies all constraints. MIS and MM can both easily be modelled in this formalism [52], [90]. Within this framework, Round Elimination (RE) acts as a mechanical operator, denoted by \(\overline{\mathcal{R}}(\mathcal{R}(\Pi))\), which transforms a given problem into a new variant that is exactly one round easier to solve [52]. In a nutshell, if one starts with a problem \(\Pi_0\) and applies this operator \(t\) times and obtains a problem \(\Pi_t\) that cannot be solved in \(0\) rounds, one obtains a \(t\) round lower bound for \(\Pi_0\). However, as analysed in [52], the allowed error probabilities also increase with each step of applying the operator and also the label space required to describe the problems increases doubly exponentially with each naïve application of the operator. Thus, to prevent the state space from expanding uncontrollably with each iteration, a relaxation is applied. A relaxation is a mapping to a structurally simpler problem \(\Pi'\) where any valid solution for \(\Pi\) remains a valid solution for \(\Pi'\) [52]. The label complexity function \(f(\Delta)\) serves as a strict upper bound on the size of the label alphabet (\(|\Sigma| \le f(\Delta)\)) across all intermediate problems and their relaxations throughout a sequence [52]. The main challenge is finding these relaxations such that the label complexity remains small and such that the sequence is long. For MIS and MM, one can limit the label complexity to \(2^{\Delta+1}\) and for MM to \(5\) while providing sequences of length \(\Theta(\Delta)\) [52], [90] (see 13). The following theorem implies that such round elimination lower bound sequences provide lower bounds for solving the respective problem in the randomised \(\mathsf{LOCAL}\) model with a certain error probability \(p\).
Theorem 1 (Extension of Theorem 7.1 in [52]). Let \(\Pi_0 \rightarrow \Pi_1 \rightarrow \cdots \rightarrow \Pi_t\) be a sequence of problems. Assume that, for all \(0 \le i < t\), and for some function \(f\), the following holds:
There exists a problem \(\Pi_i'\) that is a relaxation of \(\mathcal{R}(\Pi_i)\);
\(\Pi_{i+1}\) is a relaxation of \(\overline{\mathcal{R}}(\Pi_i')\);
The number of labels of \(\Pi_i\), and the ones of \(\Pi_i'\), are upper bounded by \(f(\Delta)\).
Also, assume that \(\Pi_t\) has at most \(f(\Delta)\) labels and is not \(0\)-round solvable in the deterministic port numbering model, even if the port numbering assignment satisfies some local constraints \(C^{\mathrm{port}}\).
Then, \(\Pi_0\) requires \[\Omega\left(\min\left\{t,\log_{\Delta}n, \log_{\Delta}\log 1/p-\log_{\Delta}\log f(\Delta)\right\}\right)\] rounds in the randomised \(\mathsf{LOCAL}\) model with local error probability \(p\), even if the port numbering satisfies some local constraints \(C^{\mathrm{port}}\).
Proof. The proof is almost verbatim along the proof of the almost identical claim in [52]. Thus, we focus on the changes to their proof that are limited to the following two aspects:
Change 1: In their proof, the error probability is hardcoded to \(p=1/n\). This is done in [52]. Instead, doing similar calculations with a general parameter \(p\) for the error probability leads to the following claim:
Modified A.5: Let \(\Pi_0 \rightarrow \Pi_1 \rightarrow \cdots \rightarrow \Pi_t\) be a sequence of problems satisfying the conditions of 1. Any randomised algorithm in the port numbering model25 running in strictly less than \(\min\{t, \frac{1}{10} (\log_{\Delta}(\log(1/p)) - \log_{\Delta} \log f(\Delta))\}\) rounds must fail with probability at least \(p\).
Proof. Applying [52], after \(t\) rounds the error probability of any randomised algorithm is \[p_{\text{Error}} \geq \frac{1}{{f(\Delta)^{\Delta}}^{10 t}} \geq \frac{1}{{f(\Delta)^{\log (1/ p) / \log f(\Delta)}}} \geq p,\] where in the second step we used our hypothesis \(t < \frac{1}{10} (\log_{\Delta}( \log( 1/p)) - \log_{\Delta} \log f(\Delta))\). So we obtain \(p_{\text{Error}} \geq p\) as desired. ◻
Change 2: In the proof of [52], the previous claim is lifted to the \(\mathsf{LOCAL}\) model using an indistinguishability argument. This argument only works if nodes cannot see the whole graph, requiring to cap the lower bound at \(\Omega(\log_{\Delta}n)\), which is the diameter of an \(n\)-node \(\Delta\)-regular tree. In their work, this does not appear in their theorem as the cap is always larger than the bound \(\log_{\Delta}( \log( 1/p)) - \log_{\Delta} \log f(\Delta))\) imposed when used with a hardcoded \(p=1/n\). Thus, for our theorem, the lower bound cannot exceed \(\log_{\Delta}n\), to which we adjust by adding it as a third argument in the min function.
Modified A.6: Let \(\Pi_0 \rightarrow \Pi_1 \rightarrow \cdots \rightarrow \Pi_t\) be a sequence of problems satisfying the conditions of1. Any randomised algorithm running in the \(\mathsf{LOCAL}\) model, in \(\Delta\)-regular balanced trees of \(n\) nodes, that fails with probability at most \(p\), requires \(\Omega\left(\min\left\{t,\log_{\Delta}(\log( 1/p)) -\log_{\Delta}\log f(\Delta), \log_{\Delta} n \right\}\right)\) rounds.
Note that this is exactly our desired statement 1. ◻
By the following lemma, \(\mathcal{O}(\log n)\) bits per angular coordinate suffice so that each angular coordinate is unique w.h.p.
Lemma 25. Let \(G\) be a threshold hyperbolic random graph and for any pair of vertices \(u,v \in V(G)\), let \(\varphi_b(u)\) and \(\varphi_b(v)\) be the first \(b\) random bits of the angular coordinates \(\varphi(u)\) and \(\varphi(v)\) respectively. Then, for any pair \(u,v\), any constant \(c > 0\) and \(b = \lceil c\cdot \log n \rceil\), it holds \(\varphi_b(u) \neq \varphi_b(v)\) with probability \(1 - \mathcal{O}(n^{-c+5/2})\).
Proof. Throughout the proof, we work with the binomial model of hyperbolic random graphs, in which the number of vertices is fixed to be \(n\). At the end of the proof, the result is transferred to the Poisson model with \(\mathbb{E}[N]=n\), incurring only a minor additional error probability (see also [43] for a detailed discussion).
Let \(b = \lceil c\cdot \log n \rceil\) as given by our hypothesis and let \(\mathcal{E}\) be the event that there exists a pair of vertices \(u,v \in [n]\) where \(\varphi_b(u) = \varphi_b(v)\). Since the angular coordinates are sampled independently and uniformly from \([0,2\pi)\), the first \(b\) bits of every angular coordinate form an independent uniformly distributed bit string in \(\{0,1\}^b\). Hence, for every fixed pair of distinct vertices \(u,v \in [n]\) we have \(\mathbf{P}\left(\varphi_b(u)=\varphi_b(v)\right) = 2^{-b}\). Subsequently, it holds for our "collision" event \(\mathcal{E}\) by a union bound \[\mathbf{P}\left(\mathcal{E}\right) \leq {n\choose 2}\cdot 2^{-b}\in\mathcal{O}(n^{-c+2}),\] where in the last step we plugged in \(b = \lceil c\cdot \log n \rceil\). We then obtain the desired probability bound by the transfer from the binomial model to Poisson [43] which yields \(\mathbf{P}\left(\mathcal{E}\right)\cdot\mathcal{O}(n^{1/2}) \in \mathcal{O}(n^{-c+5/2})\). ◻
This research was funded in whole or in part by the Austrian Science Fund (FWF) https://doi.org/10.55776/P36280, https://doi.org/10.55776/I6915. For open access purposes, the author has applied a CC BY public copyright license to any author-accepted manuscript version arising from this submission.↩︎
This research was partially funded by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) – project number 390859508.↩︎
That is, the number of vertices that have degree \(d\) is roughly \(\approx n\cdot d^{-\tau}\) (where \(\tau \in (2, 3)\) is the power-law exponent).↩︎
We say an event \(\mathcal{E}\) occurs with extremely high probability (w.e.h.p.) if \(\mathbb{P}\left(\mathcal{E}\right) \in 1 - n^{-\omega(1)}\).↩︎
For MM, we obtain the slightly stronger bound of \(\mathcal{O}(\log^3\log n)\).↩︎
With probability \(1-o(1)\) over the draw of the hyperbolic random graph \(G\).↩︎
There exists an algorithm where the error probability is at most \(1/n^{c}\) for some constant \(c>0\).↩︎
A rooted tree where each vertex except for the leaves has branching factor \(d\).↩︎
We say an event \(\mathcal{E}\) occurs with high probability (w.h.p.) if \(\mathbb{P}\left(\mathcal{E}\right) \in 1 - \mathcal{O}(1/n)\).↩︎
A unique connected component of size \(\Theta(n)\).↩︎
The largest distance between vertices is \(\Theta(\log n)\).↩︎
This can be seen by taking the line graph \(H\) of graph \(G\) such that an MM on the line graph \(H\) is equivalent to an MIS on \(G\).↩︎
This work was honoured with the FOCS’19 best paper award.↩︎
By \(V_{(2)}\) we refer to the set of vertices that are neither in the independent set nor have a neighbour in the independent set after step 2 of HRG-Shattering-MIS.↩︎
We write \(N_2(u)\) for the two-hop neighbourhood \(N_2(u) := \{v \in V : d_G(u,v) \leq 2\}\).↩︎
By \(V_{(2)}\) we refer to the set of unmatched vertices after step 2 of HRG-Shattering-MM.↩︎
A rooted tree where each vertex, except for the leaves, has \(d\) children.↩︎
The root node \(u\) is the vertex with graph distance \(d_G(u,w) = h\) to each leaf \(w \in V(T_i)\).↩︎
The value for \(m\) in 9 is here the same as the parameter \(m\) in [pro:d-regular-trees].↩︎
We call a vertex \(u\) unmatched if \(\exists v \in N(u)\) such that \(E(v) \cap M = \emptyset\) and \(E(u) \cap M = \emptyset\), i.e., under current matching \(M\), it is still possible to add an edge \(e\) incident to \(u\) such that \(M \cup e\) is a valid matching.↩︎
Recall that \(W(u) := \{w_1, w_2, \dots w_k\}\) is the similar degree path of \(u\) such that for \(j \in [k]\) it holds \(\{w_{j}, w_{j+1}\} \in E\).↩︎
Where \(w_j \in W(u)\) is the \(j\)-th vertex of the similar degree path \(W(u)\) such that \(\{u, w_1\} \in E\) and \(\{w_j, w_{j+1}\} \in E\).↩︎
For convenience, we refer to both as a Chernoff bound whenever using either of the two.↩︎
See [52] for the precise definition of the port numbering model.↩︎