Abstract 1 Introduction 2 Preliminaries 3 Degeneracy of Hyperbolic Random Graphs 4 Clique Number of Hyperbolic Random Graphs 5 Geometric Inhomogeneous Random Graphs 6 Conclusion References

Hyperbolic Random Graphs: Clique Number and Degeneracy with Implications for Colouring

Samuel Baguley ORCID Hasso Plattner Institute, University of Potsdam, Germany Yannic Maus ORCID TU Graz, Austria Janosch Ruff ORCID Hasso Plattner Institute, University of Potsdam, Germany George Skretas ORCID Hasso Plattner Institute, University of Potsdam, Germany
Abstract

Hyperbolic random graphs inherit many properties that are present in real-world networks. The hyperbolic geometry imposes a scale-free network with a strong clustering coefficient. Other properties like a giant component, the small world phenomena and others follow. This motivates the design of simple algorithms for hyperbolic random graphs.

In this paper we consider threshold hyperbolic random graphs (HRGs). Greedy heuristics are commonly used in practice as they deliver a good approximations to the optimal solution even though their theoretical analysis would suggest otherwise. A typical example for HRGs are degeneracy-based greedy algorithms [Bläsius, Fischbeck; Transactions of Algorithms ’24]. In an attempt to bridge this theory-practice gap we characterise the parameter of degeneracy yielding a simple approximation algorithm for colouring HRGs. The approximation ratio of our algorithm ranges from (2/3) to 4/3 depending on the power-law exponent of the model. We complement our findings for the degeneracy with new insights on the clique number of hyperbolic random graphs. We show that degeneracy and clique number are substantially different and derive an improved upper bound on the clique number. Additionally, we show that the core of HRGs does not constitute the largest clique.

Lastly we demonstrate that the degeneracy of the closely related standard model of geometric inhomogeneous random graphs behaves inherently different compared to the one of hyperbolic random graphs.

Keywords and phrases:
hyperbolic random graphs, scale-free networks, power-law graphs, cliques, degeneracy, vertex colouring, chromatic number
Copyright and License:
[Uncaptioned image] © Samuel Baguley, Yannic Maus, Janosch Ruff, and George Skretas; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Random network models
Related Version:
Full Version: https://arxiv.org/abs/2410.11549 [1]
Acknowledgements:
The authors would like to thank Thomas Bläsius for enriching discussions and Maximilian Katzmann for running experiments indicating a gap between core and maximum inner-neighbourhood during a research visit of JR in Karlsruhe.
Funding:
This research was partially funded by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) – project number 390859508, and by Austrian Science Fund (FWF) https://doi.org/10.55776/P36280, https://doi.org/10.55776/I6915, and https://doi.org/10.55776/DOC183.
Editors:
Olaf Beyersdorff, Michał Pilipczuk, Elaine Pimentel, and Nguyễn Kim Thắng

1 Introduction

Many real-world networks have a heterogeneous degree distribution, close to a power-law, as well as a constant clustering coefficient. The hyperbolic random graph model (HRG) introduced by Krioukov et. al. [29] combines both properties [19], which has led to considerable interest in recent years. Various aspects of HRGs have been studied, including clique size [8, 9, 16], treewidth [10], minimum vertex cover size [4, 5], and diameter [17, 24, 33]. A hyperbolic random graph is a graph embedded in the hyperbolic plane where pairs of vertices have an edge if they are close according to the hyperbolic distance. It is generated by randomly throwing n vertices on a disk of radius R (dependent on n). Since hyperbolic space grows exponentially, most vertices are of small degree and located close to the boundary of the disk, while few vertices have large degree and lie close to the centre. This distribution of vertices leads to the power-law degree distribution of HRGs.

In this paper, we study the vertex colouring problem on HRGs, along with the related concepts of clique number and degeneracy. The k-colouring problem asks to colour the vertices of a graph with k colours, while assigning different colours to adjacent vertices. The chromatic number χ(G) is the minimum number of colours needed to colour a graph in such a way. For k3 the problem is one of the original NP-hard problems [13, 21]. On general graphs, even loosely approximating the chromatic number is particularly hard [36].

The answer to the colouring problem is closely related to the clique number ω(G) and the degeneracy of a graph. The clique number is the number of vertices of the largest clique of the graph and it serves as a natural lower bound to the chromatic number. The degeneracy κ(G) is the minimum integer k for which there exists an ordering of the vertex set of G, V=(v1,v2,,vn), such that for every index i[n1], vi has at most k neighbours with greater index. Any graph G can be easily coloured with κ(G)+1 colours by iterating through the vertices in reverse order and simply colouring a vertex with a colour not used by any of its higher ranked neighbours.

We aim to study these structural parameters that are not only fundamental to the model but also in general for algorithm design in various models of computation [2, 12, 18]. The most prominent large clique in an HRG is formed by the vertices in the graph’s core [8]. Simply put, the core emerges among polynomially many vertices of distance at most R/2 from the centre of the disk, which due to the triangle inequality form a clique. We denote the size of this clique by σ(G). At this point one may wonder whether the core forms the largest clique of the graph, a statement that we disprove in Proposition 16. Nevertheless we show that the largest clique can at most be small constant factor larger than the core.

Theorem (Simplified version of Theorem 20).

There exists a constant δ>0 such that for any threshold HRG G, σ(G)+1ω(G)4/3δσ(G) holds w.e.h.p. 111An event holds with extremely high probability (w.e.h.p.), if for every c>1, there exists an n0 such that for every nn0 the event holds with probability at least 1nc.

This upper bound improves on prior work [8], which showed that there exists some constant c>1 such that ω(G)cσ(G) holds w.e.h.p., but without providing any upper bound on c.

We can now see that the core and clique are not the same. Nonetheless, (i) this theorem shows that the largest clique size and the core size are closely related, and (ii) the core of HRGs is a very well understood object, whereas the largest clique is not. Thus the natural approach to bound the degeneracy is to use the core. We show the following theorem.

Theorem (Simplified version of Theorems 9 and 9).

There exist constants δ1,δ2>0 such that for any threshold HRG G, (1+δ1)σ(G)κ(G)(4/3δ2)σ(G) holds w.e.h.p.

The main surprise of this theorem is that the degeneracy is bounded away from the coresize by a constant factor. As the chromatic number is lower bounded by the core size, the upper bound on the degeneracy in this theorem implies a simple algorithm colouring with at most (4/3δ2)χ(G) colours. The approximation guarantee of this algorithm ranges from 2/3 to 4/3 depending on further model parameters, see Section 2 and Theorem 11 for details. In any case, this improves on the previously best approximation ratio of 2 [7, Lemma 7]. The algorithm iteratively removes vertices of degree at most (4/3δ2)σ(G)1 and then colours them in the reverse order. The next thing one would hope is to be able colour the graph with ω(G) colours using the same process. In [3], the authors conducted experiments where they were iteratively removing the vertex with the smallest degree of the graph, up to vertices with residual degree equal to ω(G). In their findings, this process did not remove every vertex, implying ω(G) ¡ κ(G) for their generated graphs. We substantiate their findings by providing a rigorous proof demonstrating that the clique number is, in fact, a constant factor smaller than the degeneracy.

Theorem (Simplified version of Theorem 13).

There exists a constant ε>0 such that for any threshold HRG G, ω(G)(1ε)κ(G) holds w.e.h.p.

Our final contribution is to study the degeneracy of geometric inhomogeneous random graphs (GIRGs) [23], a sibling to HRGs. The GIRGs also combine heterogeneity and high clustering. For most properties GIRGs and HRGs exhibit the same behaviour. Perhaps the first paper to find a difference between them is [9], where the authors show that the minimum number of maximal cliques in the two models differ. We show a significant discrepancy for the degeneracy of GIRGs compared to that of HRGs, see Figure 1 and Corollary 24.

Outline.

See Figure 1 for a table with our results, as well as a plot comparing the bounds of our theorems for various model parameters. In Section 1.1 we provide a detailed discussion of our results and techniques. Section 3 contains bounds on the degeneracy of HRGs (Theorem 9). In Section 4, we show the gap between clique number and degeneracy (Theorem 13), as well as bounds on the clique number (Theorem 20). Finally, Section 5 contains results about the degeneracy of GIRGs. Statements where proofs or details are omitted due to space constraints can be found in the full version [1].

Figure 1: Results on the degeneracy κ(G) and the size of the largest clique ω(G) in hyperbolic random graphs (HRG) and geometric inhomogeneous random graphs (GIRG). The bounds hold w.e.h.p. and are stated in comparison to the core size σ(G). Each curve represents the multiplicative factor in front of σ(G) for κ(G) and ω(G) (y-axis) depending on the parameter α(1/2,1) (x-axis). Prior work is listed with white background, whereas our results are listed with blue background.

1.1 Discussion of our Results and Techniques

HRGs have a power-law degree distribution [19, 34], that is, the probability that a vertex has degree k is given by k(2α+1) . The model parameter α(1/2,1) controls the power-law exponent and all our results, particularly the size of the aforementioned constant factor gap depends on the choice of α. For the ease of presentation this overview largely omits this dependence, but the summary of our results in Figure 1 plots it in detail.

Upper bound on degeneracy (Theorem 9).

One consequence of generating a graph in hyperbolic space is that vertices tend to have fewer neighbours with increasing radius, i.e., the expected number of neighbours of a vertex decreases with the distance from the centre of the hyperbolic disc. This produces the power-law degree distribution that is valuable in modelling real-world networks. It also leads to a simple approach for upper bounding degeneracy: instead of removing vertices ordered by (increasing) degree, we remove them by (decreasing) radius. If k is such that each vertex has at most k neighbours of smaller radius, then k is an upper bound on the degeneracy.

The notion governing this approach is the inner-neighbourhood of a vertex, see Figure 2. The inner-neighbourhood of a vertex u with radius r, denoted Γ(u), is the set of vertices of distance at most R from u, i.e., they are neighbours of u, and with radius at most r, i.e., they are closer to the centre of the disc than u. The size |Γ(u)| of the inner-neighbourhood is called the inner-degree of u. The expected value of |Γ(u)| scales with the area of u’s inner-ball (r)=u(R)0(r). More precisely, 𝔼[|Γ(u)|]=(n1)μ((r)). See Figure 2 a for a visual representation. The vertex of largest inner-degree is denoted u; the choice of u and the value of |Γ(u)| depend on the random distribution of the vertices.

Figure 2: Illustration of the inner-neighbourhood. (a) The pink area is the ball 0(r). The hatched area u=u(R)0(r) is the inner-ball of u. The vertices Vu form the inner-neighbourhood Γ(u). (b) Sketch of proof of Lemma 6. Vertex u is the vertex with the largest inner-degree. Each vertex vU=V0(r) has at least nearly as many neighbours within 0(r) as |Γ(u)| (w.e.h.p.).

If vertices are removed from outside to inside, then each vertex at the time of its removal will have degree less or equal to the inner-degree of u. To bound the degeneracy, we derive a probabilistic upper bound for |Γ(u)|. We do this by finding the radius that maximises μ((r)), which upper bounds 𝔼[|Γ(u)|] (for every vertex u). Since |Γ(u)| is concentrated we can apply a Chernoff and a union bound to obtain a high probability upper bound for |Γ(u)|. Despite the simplicity of the inner-neighbourhood, we elaborate on this key concept as it is not only crucial to our upper bound on degeneracy but also for most of our other results discussed later on.

Lower bound on degeneracy (Theorem 9).

In Lemma 6 we show that the maximum inner-degree also produces a lower bound on the degeneracy κ(G), in the sense that κ(G)(1o(1))|Γ(u)| w.e.h.p. This yields asymptotically tight bounds on κ(G). We prove the lower bound of κ(G) by considering the subgraph G induced by the vertices that have smaller radius than u (see Figure 2 b). Note that G contains vertices that are not neighbours of u. We show that every vertex u in G has, w.e.h.p., at least (1o(1))|Γ(u)| neighbours in G; this statement uses the choice of u and is not true if u was an arbitrary vertex. Now, in any ordering of the vertices of G, the first vertex has at least (1o(1))|Γ(u)| neighbours of greater index, implying (1o(1))|Γ(u)|κ(G)κ(G). While this provides a lower bound that asymptotically matches our upper bound, a lower order gap remains.

Gap between clique number and degeneracy (Theorem 13).

The most immediate lower bound for degeneracy is the clique number [2], because in every ordering of the vertices of G, the vertex of the clique with the lowest index has at least ω(G)1 neighbours of higher index. We prove that the lower bounds on the degeneracy that are derived from the clique number are strictly worse than the bounds discussed above obtained via analysing the inner-degree. Our approach to show this is the following: we take an arbitrary clique K and the vertex u with the largest radius in K. Let G′′ be the subgraph induced by Γ(u). Next, we partition G′′ intro three sets of vertices, each containing at least a constant fraction of the vertices of G′′ w.e.h.p. See Figure 3 b for an illustration of this partition. Lastly, we use purely geometric arguments to show that K cannot contain vertices from all three sets. The gap then follows as the left out set contains a constant fraction of u’s inner neighbourhood. This gap between ω(G) and κ(G) makes approaches for computing the clique number via degeneracy, like that of Walteros and Buchanan [35], unsuitable for HRGs.

Clique number and the core (Theorem 20).

The upper bound on the clique number is proven via a geometric approach. Let u,v,w be any three vertices. If the vertices are far apart they cannot be contained in a clique. Otherwise their pairwise distance is at most R. Using the hyperbolic of version Jung’s theorem [20, 14, 15] implies that they are contained in a ball of small radius and we show that also has a small area. Hence the expected number of vertices in is small as well. As this expectation is well concentrated whenever the area is significantly large, this bound holds with extremely high probability when introducing lower order deviations from the expectation. Via a union bound over all possible (n3) triples of vertices we rule out that any clique contained in such a covering ball is large. The main claim now follows because any clique has to be contained in one of these coverings balls. The question of whether ω(G)/σ(G)1 as n remains open.

2 Preliminaries

Hyperbolic Random Graphs.

We follow the formalisation of hyperbolic random graphs introduced in [34], which is known as the native representation. We denote by 2=[0,)×[0,2π) the hyperbolic plane in the polar coordinate system, where a point x2 is parametrised by a radius r(x) and an angle φ(x). We equip 2 with a metric dh(x,y) characterised by

cosh(dh(x,y))=cosh(r(x))cosh(r(y))sinh(r(y))sinh(r(x))cos(φ(x)φ(y)). (1)

This metric is what gives 2 a hyperbolic geometry, of curvature -1, as opposed to the Euclidean metric. We equip 2 with the topology induced by dh.

The geometric space of most importance in this work is the bounded disk in 2 defined by 𝒟R=[0,R]×[0,2π), where R=2log(n)+C with CΘ(1). We refer to point (0,0) as the centre of this disk. The space 𝒟R inherits the topology of 2, and from now on we shall only consider subsets of this space – thus, for example, a ball around a point x𝒟R is defined by the set x(ε)={y𝒟R:dh(x,y)ε}𝒟R.

We now introduce a probability measure μ on 𝒟R, which is parametrised by the model parameter α(1/2,1), and was first defined by Papadopoulos et. al. [34]. For measurable 𝒮𝒟R, define

μ(𝒮)=Sρ(x)dx,ρ(x)=αsinh(αx)2π(cosh(αR)1),

where ρ is the density of μ with respect to the Lebesgue measure on 𝒟R. This measure differs from the uniform probability measure on 𝒟R in that it puts more mass at the centre of the disk; both measures coincide at α=1. The benefit of μ lies in the properties it induces in our central object of study, the hyperbolic random graph.

Threshold hyperbolic random graph (HRG).

A (threshold) hyperbolic random graph or HRG is a pair G=(V,E) defined by the following procedure. First, n vertices are sampled independently at random in 𝒟R according to μ. Then any two vertices u,vV are connected by an edge if and only if their distance dh(u,v) is at most R. We write G𝒢(n,α,C) to denote a graph generated in this way. A vertex uV is identified by its point coordinates in 2 and we write V𝒜 to denote the set of vertices that are located in an area 𝒜𝒟R.

The use of μ to distribute vertices in 𝒟R has the effect of giving G a power-law degree distribution, as was shown in [19, 34]. It is sometimes convenient to characterise connection of vertices in terms of their angular distance, and to that end we define

θR(r1,r2)=arccos(cosh(r1)cosh(r2)cosh(R)sinh(r1sinh(r2)),

which per (1) yields the following observation.

Observation 1.

Two vertices u and v are connected if and only if their angular distance is less than θR(r(u),r(v)).

We also make use of the following expression of the distribution of the radius of a vertex u.

(r(u)r)=μ(0(r))=0rππρ(x)dθdx=cosh(αr)1cosh(αR)1=(1o(1))eα(Rr) (2)

We briefly note that a variant of HRGs exists in which vertices are not connected purely according to whether their distance is below a threshold, but rather with probability p(u,v)=(1+exp(12T(dh(u,v)R)))1 determined by both distance and a “temperature” parameter T (see e.g. [30, §3.1]).

Degeneracy, clique number, chromatic number and core.

For a graph G=(V,E), the degeneracy κ(G) is the minimum integer k for which there exists an ordering of the vertex set of G, V=(v1,v2,,vn), such that for every index i[n1], vi has at most k neighbours with greater index. The clique number ω(G) is the size of the largest clique of G. The chromatic number χ(G) is the smallest number of colours required so that a conflict-free vertex colouring is possible for G. The core of a hyperbolic random graph is the set of vertices with radius at most R/2 and we denote its size by σ(G). Since for any points u,v0(R/2) the distance is at most dh(u,v)R, the core forms a clique. Finally, since the core is a clique, any vertex of a clique needs a different colour in a conflict-free colouring, and χ(G)κ(G)+1 (see e.g. [31, Lemma 4]), we have the following chain of inequalities.

Observation 2.

Let G𝒢(n,α,C) be a threshold HRG. Then,

σ(G)ω(G)χ(G)κ(G)+1.

Concentration bounds.

We use the following Chernoff bounds [32, Theorem 4.4].

Theorem 3 (Chernoff bound).

For i[k], let Xi{0,1} be independent random variables and X=iXi. Then for ε(0,1),

(X(1+ε)𝔼[X])eε2𝔼[X]/3 and (X(1ε)𝔼[X])eε2𝔼[X]/2.

3 Degeneracy of Hyperbolic Random Graphs

A tool we make use of several times is the inner-ball of a point x𝒟R, defined by x=x(R)0(r(x)). The inner-ball of a vertex is the inner-ball of the point using the vertex’ coordinates. The inner-neighbourhood of a vertex u is the set of vertices (excluding u) contained in its inner-ball, that is, the neighbours of u that have a smaller radius than u (see Figure 2 a), and is denoted Γ(u).

We upper bound the degeneracy κ(G) via the inner-neighbourhood by using the following informal process. Consider a graph G and order its vertices (v1,v2,,vn) by decreasing radius, so r(vi)r(vi+1), and iteratively remove vertices from G one-by-one, from lower to higher index. Note that the set of neighbours of vi that have greater index than i coincide with Γ(vi). This implies that the degree of each vertex vi at the time of its removal is |Γ(vi)|. Let u be the vertex vk that maximises |Γ(vk)|. As the largest degree of a vertex at the time of its removal is given by |Γ(u)|, we get the following upper bound for the degeneracy.

Observation 4.

Let G𝒢(n,α,C) be a threshold HRG and let u be the vertex of G with the largest inner-degree in G. Then κ(G)|Γ(u)|.

We will now show that the largest inner-degree does not only yield an immediate upper bound on the degeneracy, but also a lower bound. Informally, this lower bound follows from the following argument. Order the vertices of the graph in descending order of their radius, that is, (v1,v2,,vn) such that r(vi)r(vi+1). For i[n], let Gi=G{v1,v2,,vi}. Let vk be the node with the maximum inner neighbourhood in G. We show (with a probabilistic guarantee) that the graph Gk has minimum degree (1o(1))|Γ(vk)|. Before we make this formal, we introduce a slightly modified version of [19, Lemma 3.3], that implies the following: the closer a vertex is to the origin, the more neighbours it has in expectation. This can be derived via the fact that the angle θR(r,x) is monotonically decreasing in x.

Corollary 5.

Let r,s,t[0,R) with s<t. Then μ(s(R)0(r))μ(t(R)0(r)).

Corollary 5 tells us that any vertex with radius at most r has in expectation at least as many neighbours up to radius r, as the expected inner-degree of a vertex with radius exactly r. In order to derive a high-probability bound on |Γ(u)| it now suffices to show concentration around the expectation of all considered neighbourhoods.

Since the size of every clique K is upper bounded by the inner-degree of the vertex in uK that has the largest radius, ω(G) is a lower bound for the maximum inner-degree. We can now lower bound the degeneracy κ(G) based on the largest inner-degree. Note that, in contrast to the upper bound of Observation 4, the lower bound is not deterministic.

Lemma 6.

Let G𝒢(n,α,C) be a threshold HRG. Then κ(G)(1o(1))|Γ(u)| w.e.h.p.

Proof.

For any subgraph HG, it is clear that κ(G)minvV(H)deg|H(v), where deg|H(v)=wV{v}𝟙{{v,w}E(H)}. We let H be the (random) subgraph of G created by restricting G to vertices that land in 0(r), that is, that have radius at most r, and keeping the same edges. Then for any vV,

𝔼[deg|H(v)|r(v)] =wV{v}(w0(r)r(v)(R)|r(v))𝟙{r(v)r}
=(n1)μ(0(r)r(v)(R))𝟙{r(v)r}
(n1)μ(0(r)r(R))𝟙{r(v)r}. (Corollary 5)

We write γ(n1)μ(0(r)r(R)). Since deg|H(v) is a sum of Bernoulli random variables that are all independent under (|r(v)), a Chernoff bound (Theorem 3) gives

(minvV(H)deg|H(v) <(1ε)γ)=(vV{deg|H(v)<(1ε)γ,r(v)r})
n(deg|H(v)<(1ε)γ,r(v)r) (Union bound)
=n𝔼[(deg|H(v)<(1ε)γ|r(v))𝟙{r(v)r}] (Tower rule)
n𝔼[eε2𝔼[deg|H(v)|r(v)]/2𝟙{r(v)r}] (Chernoff bound)
neγε2/2(r(v)r)neγε2/2.

Taking r to be the argmax of γ yields γ(r)𝔼[Γ(u)]. By Observation 2 and applying (2) at R/2 we have that γ(r)𝔼[Γ(u)]𝔼[σ(G)]nΩ(1). Thus, choosing ε=1/log(n), we obtain for any constant c as long as n is large enough

(κ(G)<(1ε)|Γ(u)|)(κ(G)<(1ε)γ(r))neγ(r)ε2/2nenΩ(1)/log(n)n1c,

that is, κ(G)(1o(1))|Γ(u)| w.e.h.p.

In the rest of this section we derive bounds for the largest inner-degree on HRGs that hold w.e.h.p. which, by Observations 4 and 6, translate into results for the degeneracy. Since a vertex v belongs to the inner-neighbourhood of a vertex u, if and only if it resides in the inner-ball of u, we can use the measure of an inner-ball u to bound the maximum inner-degree of a graph. Moreover, the measure of the inner-ball is invariant under rotation around the origin, which is why we write μ((r)) instead of μ(u) for r=r(u). We sum up our results for the area of the inner-ball in the following technical lemma.

Lemma 7 (Volume of the inner-ball).

Let ΔΘ(1) and let u𝒟R with radius r=R/2+Δ. Then, depending on Δ, there exist constants γ,η, such that

μ((r)) (1+Θ(eαR))αeαr(α1/2)(2πe12(2α1)(2rR)(2π(α1/2)α)) and
μ((r)) (1+Θ(eαR))αeαrα1/2(γe12(2α1)(2rR)η),

where

γ,η={1,12α for Δ0,433,12α(1433)(43)(α1/2) for Δlog(4/3),12,12α(1433)(43)(α1/2)(43312)2(α1/2) for Δlog(2),23,12α(1433)(43)(α1/2)(43312)2(α1/2)(1223)2(2α1) for Δlog(2).
Lemma 8.

Let r be the radial coordinate of the point in 𝒟R that maximises the measure of the inner-ball. Then r=R/2+log(αηγ(1α))/(2α1).

Proof.

Using the expression of μ((r)) in Lemma 7,

ddrμ((r)) =(1+Θ(eαR))(ηα2eαrα1/2γαeαrα1/2((1α)e12(2α1)(2rR))).

Setting this equal to 0 and solving for r yields r=R/2+log(αηγ(1α))/(2α1).

We use Lemma 7 to upper and lower bound the degeneracy. The lower bound tells us that there exists a constant δα>0 such that κ(G)(1+δα)σ(G) w.e.h.p. The constant δα is increasing with increasing 1/2<α<1, see Figure 1.

Theorem 9 (Bounds on the degeneracy).

Let G𝒢(n,α,C) be a threshold HRG. Then w.e.h.p. its degeneracy κ(G) satisfies

(4o(1))π(2(1α)π2α(π2))1α2α1σ(G)κ(G)((4/3)α+o(1))σ(G).

Proof sketch.

The upper bound is obtained by using r as stated in Lemma 8 for the upper bound of the inner-ball with (n1)μ((r)), and using a Chernoff bound along with a union bound which gives an upper bound for |Γ(u)| w.e.h.p. Observation 4 then yields the upper bound for the degeneracy. For the lower bound, we first show that there exists a vertex with a radius r~ close in value to r w.e.h.p. Then (1o(1))(n1)μ((r~)) lower bounds |Γ(u)| w.e.h.p. This gives a lower bound for the degeneracy, due to Lemma 6.

Applying Observations 2 and 9, the following is immediate.

Corollary 10 (Bounds on the chromatic number).

Let G𝒢(n,α,C) be a threshold HRG. Then w.e.h.p. its chromatic number is σ(G)χ(G)((4/3)α+o(1))σ(G).

Our structural results directly produce algorithmic applications. The small gap between degeneracy and core translates into an efficient approximation algorithm to colour a HRG.

Theorem 11.

Let G𝒢(n,α,C) be a threshold HRG. Then an approximate vertex colouring of G can be computed in time 𝒪(n) with approximation ratio ((4/3)α+o(1)) w.e.h.p.

Proof sketch..

Using a smallest-last vertex ordering [31] the number of colours required is upper-bounded by κ(G). Computing the smallest-last vertex ordering, and then using the ordering to colour the graph, both takes linear as the giant component is sparse w.e.h.p. by [25, Corollary 17]. The approximation ratio is achieved by comparing the lower bound of the chromatic number in Corollary 10 to the upper bound of the degeneracy in Theorem 9.

4 Clique Number of Hyperbolic Random Graphs

Recall that for any graph G, the clique number ω(G), chromatic number χ(G), and degeneracy κ(G) are related via the inequalities ω(G)χ(G)κ(G)+1. For this reason we are interested in the relationship between clique number and degeneracy for hyperbolic random graphs. In Section 4.1 we show that the two differ and that for HRGs, the degeneracy is strictly larger than the clique number by a constant multiplicative constant. In Section 4.2 we give new insights about where in the hyperbolic disk the largest clique is formed. We then conclude the section by providing a new upper bound for the clique number in Section 4.3 that states a leading constant in front of the size of the core, and that is increasing in α.

4.1 The gap between Clique Number and Degeneracy

Because of the centralising effect of hyperbolic geometry, one might hope to show that the clique contained in the core of the disk is the largest, and that ω(G)=(1o(1))κ(G). This would achieve two things on HRGs: first, it would imply a tight bound for the chromatic number χ(G), sandwiching it between clique number and degeneracy. Moreover, it would also imply a linear time (1+o(1))-approximation algorithm for the two NP-complete problems clique number and chromatic number using a smallest-last vertex ordering (see [31]).

In this section we disprove these claims. We show that there exists a constant gap between clique number and degeneracy; this is the content of Theorem 13. Before embarking on the details of the proof, we first sketch its idea. For any clique K, its size is bounded by the inner-degree |Γ(u)|, where u is the vertex with largest radius among the vertices of K. For r(u)R/2+o(1) and r(u)R/2+ω(1), |Γ(u)| is already smaller by a multiplicative constant than the lower bound for the degeneracy given in Theorem 9. What remains is to extend the result to r(u)R/2+Θ(1), which requires more intricate arguments and is addressed in the following lemma.

Lemma 12.

Let G𝒢(n,α,C) be a threshold HRG, let ΔΘ(1) and let K be any clique where uK is the vertex with largest radius rK=R/2+Δ. Then w.e.h.p. there exists a constant ε(0,1) such that |K|(1ε)|Γ(u)|.

Proof.

Figure 3: Illustration of a separator. (a) The line aφ partitions disk 𝒟R into two halfdisks 𝒟1 (blue) and 𝒟2 (pink). The hypercycle φ (hatched area) is defined by the line aφ. Points located in different halfdisks and outside the hatched area have distance at least R. (b) The separator (hatched area) separates the inner-neighbourhood (grey area) of a vertex u of radius rK into three sub-areas 𝒮0, 𝒮1 and 𝒮2. Any vertex located in 𝒮1 has no edge to a vertex in 𝒮2.

Let G be the inner-neighbourhood of u, i.e., the induced subgraph G=G[Γ(u)]. Then KV(G), and thus ω(G)|K|. We show that ω(G)<(1ε)|Γ(u)| w.e.h.p.

We accomplish this by proving that there exists a colouring for G with (1ε)|V(G)|=(1ε)|Γ(u)| many colours. This implies an upper bound for the chromatic number χ(G), which also serves as an upper bound for ω(G). We do this by partitioning the inner-neighbourhood u into three disjoint sub-regions 𝒮0𝒮1𝒮2 such that no vertex in 𝒮1 is adjacent to any vertex in 𝒮2. Thus, 𝒮0 separates 𝒮1 from 𝒮2 and we can colour the set of vertices V𝒮1 with the same colours as V𝒮2. Thus if min{|V𝒮1|,|V𝒮2|}ε|Γ(u)| w.e.h.p., our desired statement will be proven.

To find such a separator 𝒮0, we follow the lines drawn by Bläsius et. al. [10] where hypercycles for hyperbolic random graphs were introduced (see Figure 3 a). A hypercycle φ (of radius R/2) is defined as follows: let aφ denote the line whose points have angle φ and φ+π. Then φ:={u𝒟R:dh(u,aφ)R/2}, i.e. the set of points with distance at most R/2 to line aφ. Consider the point u=(rK,φ) and let 𝒮0:=uφ. To define 𝒮1 and 𝒮2 separate 𝒟R into the two disjoint halfdisks 𝒟1={x𝒟R:φ(x)φπ} and 𝒟2={x𝒟R:φ(x)φπ} (see Figure 3 a). Then we define 𝒮1:=(u𝒟1)φ and symmetrically 𝒮2:=(u𝒟2)φ.

We observe that any point w𝒮1 has distance at least R to any point v𝒮2. This can be shown for example by observing that the geodesic between w and v must pass through some point xaφ, and so dh(w,v)=dh(w,x)+dh(x,v)dh(w,aφ)+dh(aφ,v)>R. This ensures our objective of separating the two regions via 𝒮0.

Next, we show that there exists a constant ε>0 small enough, such that μ(𝒮1)=μ(𝒮2)εnα (a sketch of the idea is given in Figure 3 b). Setting ϕ(x)=max(0,θR(rK,x)θR(x,x)/2) we derive by symmetry

μ(𝒮1)=μ(𝒮2)=R/2r0ϕ(x)ρ(x)dϕdx=R/2R/2+Δϕ(x)ρ(x)dx.

Now we choose another constant Δ<Δ that fulfills 2θR(R/2+Δ,R/2+Δ)θR(R/2+Δ,R/2+Δ)cΘ(1). This is possible because ΔΩ(1). By our choice of Δ we then obtain

μ(𝒮1) R/2+ΔR/2+Δϕ(x)ρ(x)dx=R/2+ΔR/2+Δ(θR(R/2+Δ,x)θR(x,x)2)ρ(x)dx
cR/2+ΔR/2+Δρ(x)dx=c(cosh(α(R/2+Δ))cosh(α(R/2+Δ)))cosh(αR)1Ω(nα),

where the last line follows since θR(,) is monotonically decreasing and since ρ(x)=αsinh(αx)cosh(αR)1, |ΔΔ|>0 and R=2log(n)+C. Therefore μ(𝒮1)Ω(𝔼[|Γ(u)|]/n) w.e.h.p., since

lim infnnμ(𝒮1)𝔼[|Γ(u)|]=lim infnμ(𝒮1)nαn1α𝔼[|Γ(u)|]>0,

where 𝔼[|Γ(u)|]nμ(0(R/2+Δ))𝒪(n1α) by Equation 2, because Δ𝒪(1). Thus there exists some ε>0 for which, for n large enough,

𝔼[|V𝒮2|]=𝔼[|V𝒮1|]=nμ(𝒮1)(11/log(n))2ε𝔼[|Γ(u)|];

since |V𝒮1||Γ(u)| a.s. and is strictly smaller with positive probability, then ε<1.

Using a Chernoff bound for both |V𝒮1| and |V𝒮2|, we obtain that neither random variable is smaller than (11/log(n))𝔼[|V𝒮1|](11/log(n))1ε𝔼[|Γ(u)|] w.e.h.p. On the other hand, another application of a Chernoff bound reveals |Γ(u)|(1+1/log(n))𝔼[|Γ(u)|] w.e.h.p., since for r(u)R/2+Θ(1) we have 𝔼[|Γ(u)|]Θ(n(1α)) using Lemma 7. A union bound then shows that w.e.h.p., min{|V𝒮1|,|V𝒮2|}ε𝔼[|Γ(u)|]11/log(n)ε|Γ(u)|. As argued above, a naïve colouring that colours the vertices of 𝒮1 and 𝒮2 with the same set of colours yields the upper bound.

Theorem 13 (Clique-degeneracy-gap).

Let G𝒢(n,α,C) be a threshold HRG. Then w.e.h.p. there exists a constant ε(0,1) such that κ(G)/ω(G)>1+ε.

Proof.

Let K be the largest clique of G, and let u be the vertex of K with maximal radius rKr(u). We assume the following cases for rK which cover all possibilities:
Case 1 [rKR/2+ω(1)]: Observe that KV(rK). Hence, |K||V(rK)| and

μ((rK))(1+Θ(eαR))αeαrKα1/2(γe12(2α1)(2rKR)η)Θ(1)nαe(1α)ω(1),

by Lemma 7 since γ and η are both constants. Taking the expectation and using a Chernoff bound we have |K|o(n1α) w.e.h.p. Recall that κ(G)>(1+δ)eαC/2n1α w.e.h.p. by Theorem 9. It follows κ(G)/|K|ω(1) w.e.h.p.
Case 2 [rKR/2+o(1)]: Observe that the size of any clique KV0(r) is upper bounded by X=|V0(r)|. By Equation 2 we have μ(0(rK))(1+o(1))nαeαC/2. Hence, we get 𝔼[X](1+o(1))n1αeαC/2. Since X is a binomial random variable we can apply a Chernoff bound, which yields |K|X(1+o(1))n1αeαC/2 w.e.h.p. Taking ε=δ/(1+δ) with δ as in Theorem 9, we have

|K|1ε(1+o(1))(1+δ)n1αeαC/2<(1+o(1))κ(G) w.e.h.p.

Case 3 [rKR/2+Θ(1)]: Let ξo(1). Using Lemmas 12 and 6, we get w.e.h.p. that

ω(G)=|K|1ε1ξ|Γ(u)|1ε1ξ|Γ(u)|(1ε)κ(G),

for an adequate choice of ξ.

Corollary 14.

Let G𝒢(n,α,C) be a threshold HRG. Then there exists a constant ε(0,1) such that w.e.h.p., ω(G)(4(1ε)/3)ασ(G).

Proof.

This follows from Theorems 9 and 13.

4.2 Cliques larger than the Core

In this section, we show that there exist a super-constant number of unique cliques that contain the core, but are strictly larger than it. The overall argument goes as follows. We consider a set of vertices with radial coordinates slightly outside the core. We show that any vertex of this set has a constant probability to be adjacent to all vertices belonging to the core, and thus induces a clique that is larger than the core itself. To this end, the following lemma concerned about points close to the core proves useful.

Lemma 15.

Let k{0} and ξk=log(1+logk(n)n1α)o(1). Consider two points with radial coordinates r=R/2+ξk and x=R/2. Then θR(r,x)π2logk(n)nα1.

Lemma 15 lower bounds the angle distance such that two points have distance at most R.

Proposition 16.

Let G𝒢(n,α,C) be a threshold HRG. Then w.e.h.p. there exist ω(log(n)) cliques that are larger than σ(G).

Proof.

Figure 4: Sketch of the proof idea for Proposition 16. (a) Illustration of the set of points 𝒜 (blue) of width ξ slightly outside of the core (pink). The intersection with a sector of angle ψ and the band 𝒜 forms a box 𝒜i that contains a vertex w.e.h.p. The number of non-intersecting boxes is k=2π/ψω(log(n)). (b) A vertex u located on the boundary of the area 𝒜. The hatched area u with angle at most ϕ is the corresponding forbidden area of u. Any point in u has distance at least R to u. An adequate choice for the width ξ yields that this area is empty with constant probability.

For the proof we consider the Poissonized version of the HRG model (see e.g. [26, 27]). The upshot of this model is that it allows us to analyse disjoint areas in the hyperbolic disk independently. Since the final result holds w.e.h.p. this directly carries over to the uniform model w.e.h.p. ([22, Lemma 3.9]). We start by defining an area 𝒜 close to the core. Let ξ=log(1+log3(n)n1α)o(1) and consider a band of points 𝒜:=0(R/2+ξ)0(R/2). We see via R=2log(n)+C that

μ(𝒜)=R/2R/2+ξsinh(αx)cosh(αR)1=cosh(α(R/2+ξ))cosh(αR/2)cosh(αR)1=(1+Θ(eαR))(eα(R/2ξ)eαR/2)=(1+Θ(eαR))(nαeαC/2(1+log3(n)nα1)nαeαC/2)=(1+Θ(eαR))(eαC/2log3(n)n1)Θ(log3(n)/n). (3)

Let A=V𝒜. Then 𝔼[|A|]Θ(log3(n)). We further partition the band 𝒜 into k=log3/2(n) sectors 𝒜1,,𝒜k, each of equal size (see Figure 4 a), and Ai=V𝒜i. Since 𝔼[|A|]Θ(log3(n)) and k=log3/2(n), we have for any i[k] that 𝔼[|Ai|]=𝔼[|A|]kω(log(n)). Since we are in the Poissonized model we get (|Ai|=0)=e𝔼[|Ai|]nω(1). Subsequently, a union bound yields that there is no sector 𝒜i that is empty w.e.h.p.

In the second step, we show that for any vertex uA, the probability that there exists a vertex in its forbidden area u:={x0(R/2):dh(u,x)>R} (see Figure 4 b) is strictly less than 1. Since r(u)R/2+ξ, we have that uA has distance at most R (and thus, an edge) to any vertex in the area 0(R/2ξ). Hence, we have u0(R/2)0(R/2ξ). Now, for R/2xR/2ξ we seek to find the angle size ϕ=2(πθR(r(u),x)) in order to upper bound μ(u). Since θR(r,x) is monotonically decreasing in x we have ϕ2(πθR(R/2+ξ,R/2)) and we obtain

μ(u)=R/2ξR/20ϕρ(x)dϕdx 2(πθR(R/2+ξ,R/2))R/2R/2+ξρ(x)dx
2(πθR(R/2+ξ,R/2))μ(𝒜),

where we used R/2ξR/2ρ(x)dxR/2R/2+ξρ(x)dx=μ(𝒜). By our choice of ξ=log(1+log3(n)nα1)o(1), we can apply Lemma 15 and obtain θR(R/2+ξ,R/2)π2log3(n)n1α. In Equation 3 we established that μ(𝒜)Θ(log3(n)/n). Combining the two leads to μ(u)𝒪(log9(n)nα3).

Notice that for α(1/2,1), the measure of the forbidden area of a vertex uA is μ(u)o(1/n). Hence, writing F=Vu, we get 𝔼[|F|]o(1), i.e., the expected number of vertices in u is vanishing. Applying Markov’s inequality then gives us (|F|1)o(1). Thus p(|F|<1)Ω(1), so the forbidden area is empty with constant probability.

We now establish our third and final desired property. Namely, we construct a subset UA, consisting of vertices whose forbidden areas are empty and disjoint and for which 𝔼[|U|]ω(log(n)).

Recall that we partitioned 𝒜 into k=log3/2(n)ω(log(n)) sectors. Let Xi be the indicator that there exists a vertex u𝒜i whose forbidden area u is empty. Then X=i=1k/2X2i is a loose lower bound on the number of sectors with this property. By linearity of expectation, p constant and k=log3/2(n) we then obtain

𝔼[X]𝔼[i=1k/2X2i]pΘ(1)log3/2(n)ω(log(n)).

We proceed by showing independence among the random variables Xi and Xj for |ji|>1. To this end, observe that the angle ψ (see Figure 4 b) spanned by any sector 𝒜i is ψ=2π/k2πlog3/2(n). In contrast, we recall ϕ2(πθR(R/2+ξ,R/2))4log3(n)nα1. Since log3(n)nα1o(log3/2(n)) we conclude ϕo(ψ). This implies that the forbidden areas u and v of any u𝒜i, v𝒜j are disjoint. Thus Xi and Xj are independent.

To wrap things up, recall that in the first step we established that each sector 𝒜i contains a vertex w.e.h.p. Moreover, since the Xi are independent, we have by a Chernoff bound that Xω(log(n)) w.e.h.p. Though these two events are not independent, we can apply the union bound to their complements to obtain that w.e.h.p. ω(log(n)) vertices outside of 0(R/2) are adjacent to all vertices in 0(R/2), which finishes the proof.

4.3 Upper Bound on the Clique Number

Recall that two vertices u,v are adjacent if and only if dh(u,v)R and that we call the clique formed in 0(R/2) the core whose size is σ(G)=(1o(1))eαC/2n1α w.e.h.p. The core size is a lower bound for the clique number ω(G). We have established that the largest clique is smaller than the degeneracy w.e.h.p. (Theorem 13), and in this section we further investigate an upper bound for ω(G). We note that the upper bound we derive in this section implies Theorem 13 for α large enough (see Figure 1). However, for smaller α, the upper bound for ω(G) is larger than the lower bound (Theorem 9) for κ(G) in the HRG model, and thus does not directly imply Theorem 13 for these values for α.

Before going into details, we lay out our proof strategy. We aim to bound the region where a clique can be located. Since vertices are adjacent if and only if their (hyperbolic) distance is at most R, this can be done by characterising a shape that covers any hyperbolic region of diameter R. A classic result by Jung [20] answers the question of how large the radius of a ball in Euclidean space needs to be at most, so that its interior can contain an entire set of points of fixed diameter. The hyperbolic version of this result was discovered by Dekster [14, 15] nearly a century later. He extended Jung’s result to (among other geometries) hyperbolic space and we apply it as follows: we identify 𝒪(n3) many balls where one of these balls contains the clique of largest size ω(G). This clique (and all the other identified cliques) needs to be located in a ball Bx(r) with radius r large enough. We use the hyperbolic variant of Jung’s theorem to upper bound r which, in turn, allows us to upper bound the area of this ball. This yields an upper bound for the amount of vertices one such ball Bx(r) could contain w.e.h.p., leading to an upper bound for ω(G). Since we only need to consider at most 𝒪(n3) balls, a union bound is sufficient to derive the same bound for the worst case. We work with the following version of Jung’s theorem for hyperbolic geometry.

Theorem 17.

[14, Theorem 2] Let 𝒦d be compact and suppose that for any y,z𝒦, dh(y,z)D. Then there exists xd such that 𝒦x(r) for r satisfying

D2sinh1(d+12dsinh(r)).

In the hyperbolic plane 2, this simplifies to the following.

Corollary 18.

Let 𝒦2 be compact and suppose that for any y,z𝒦, dh(y,z)R. Then there exists x2 such that 𝒦x(r) for r satisfying rR/2+log(2/3).

Proof.

Using that d=2, we directly get from Theorem 17 for diameter R that R2sinh1(3/4sinh(r)). Rearranging and using for x that sinh(x)=12ex(1e2x) yields

R/2rlog(3/4)+log((1e2r)(1eR)).

Solving for r and using that rR/2log(n)+C/2 in conjunction with recalling that CΘ(1), it follows that rR/2+log(2/3).

Our next observation follows from the definition of μ, and formalises the intuition that μ puts more mass at the centre of the disk.

Observation 19.

Let 0<rR and u,v𝒟R with r(u)r(v). Then μ(u(r))μ(v(r)).

Recall that σ(G) denotes the core size |V0(R/2)| which is a lower bound for the clique number ω(G) (see Observation 2), and that σ(G)=(1o(1))eαC/2n1α w.e.h.p. We state our upper bound relative to this lower bound.

Theorem 20 (Clique upper bound).

Let G𝒢(n,α,C) be a threshold HRG with α(1/2,1). Then w.e.h.p.

ω(G)((4/3)α/2+o(1))σ(G).

Proof.

Consider any triplet of vertices u,v,wV with pairwise distance at most R, so that they are pairwise adjacent. Since u,v,w are a.s. in general position, there is a unique ball x(r) such that u,v,w lie on the boundary x(r). By Corollary 18, rR/2+log(4/3) since max(dh(u,v),dh(v,w),dh(u,w))R. Over all possible triplets u,v,wV this gives us a set B of at most (n3) closed balls. Any clique must be contained in one of these balls, and therefore so is the largest clique. Thus upper bounding the number of vertices for each individual ball B yields an upper bound on the size of the largest clique.

We now upper bound the expected number of vertices in one ball B. To this end we fix a ball =x(r)B and let Y be the random variable counting the number of vertices in . The balls in B are identically (though clearly not independently) distributed. Since vertices are thrown independently according to μ, we have that Y3Bin(n3,μ()), and so

𝔼[Y]=3+(n3)μ() 3+(n3)μ(0(r)) (Observation 19)
3+(n3)μ(0(R/2+log(2/3))) (Corollary 18)
((2/3)α+o(1))eαC/2n1α. (Equation 2)

Thus ((4/3)α/2+o(1))σ(G)(1+1/log(n))𝔼[Y] w.e.h.p. This is relevant to the bound in the theorem statement because it implies that

(ω(G)>((4/3)α/2+o(1))σ(G)) (maxBY>((4/3)α/2+o(1))σ(G))
(maxBY>(1+1/log(n))𝔼[Y])+nc

for arbitrary c. Thus to finish the proof we need to show concentration, which via a union bound over all triplets u,v,w will yield the result. To show concentration we apply a Chernoff bound Theorem 3. Using ε=(1/log(n)) we obtain

(Y>(1+1/log(n))𝔼[Y]) e𝔼[Y]/(3(1/log(n))2)eΘ(1)(n1α)/log2(n)nc

for any choice of c, since for α<1, Θ(1)(n1α)nΘ(1) and lim infnnΘ(1)log2(n)ω(log(n)). Finally, to show that this holds w.e.h.p. for all balls in B, we use that |B|(n3)<n3, so that

(maxBY(1+1/log(n))𝔼[Y]) B(Y(1+1/log(n))𝔼[Y])
n3(Y(1+1/log(n))𝔼[Y])nc+3.

A further refinement of the “clique covering” argument of Theorem 20 should be possible. Any clique has by definition a diameter of at most R, and so the shape in 2 of diameter R with maximal area would provide an improved upper bound via a similar covering argument. It is not clear what a tight bound would be, and ω(G)(1+o(1))σ(G) may be possible.

5 Geometric Inhomogeneous Random Graphs

Geometric Inhomogeneous Random Graphs or GIRGs were introduced in [11] as an alternative model to HRGs that capture many of the same properties, in particular the power-law degree distribution. In their most general form, GIRGs strictly generalise HRGs, but they are more often studied in a slightly restricted form; comparisons are made in [6, 28]. In this restricted form, called the standard GIRG model by [16], any HRG G can be coupled with two GIRGs H1 and H2 such that H1GH2, where denotes graph inclusion.

Because of this relationship, GIRGs are used as proxies for HRGs in some theoretical and experimental works. This is partly done because GIRGs are (by design) far more tractable than HRGs. It is therefore valuable to understand differences between the two models. In [6] experimental evidence was given to suggest that the “sandwiching” of an HRG by two standard GIRGs is not tight. In Corollary 24 we provide a theoretical result demonstrating a difference between the two models.

Definition 21 (Standard GIRG model).

Let β(2,3), λΘ(1), and n. A geometric inhomogeneous random graph G𝒢(n,β,λ) is a random graph with vertex set V={v1,,vn} satisfying the following properties.

  1. 1.

    Every uV is equipped with a random tuple (wu,xu), where weight wu[1,) has density f(y)=(β1)yβ and coordinate xu is drawn uniformly at random from [0,1];

  2. 2.

    Any pair of vertices u,vV are connected if and only if min{|xuxv|,1|xuxv|}t(u,v), where t(u,v)=12(λwuwvn).

One way of thinking of a GIRG is that vertices are being thrown uniformly at random onto the 1-dimensional torus 𝕋1, and connected according to whether their distance is below their threshold t(u,v). The weights are drawn according to a Pareto distribution. Analogously to HRGs, for a vertex u of a GIRG we define the inner-degree of u to be |Γ(u)|=|{vV|u and v are connected and wvwu}|. The proofs of Observations 4 and 6 can be adapted to the GIRG model to characterise the degeneracy via the largest inner-degree.

Corollary 22.

Let G𝒢(n,2α+1,λ) be a standard GIRG. Consider the vertex u with the largest inner-degree in G. Then w.e.h.p. κ(G)=(1o(1))|Γ(u)|.

Corollary 22 allows us to state a tight bound for the degeneracy in comparison to the core of the GIRG, which is defined to contain all vertices of weight w^n/λ, and has size σ(G)=(1±o(1))λαn1α w.e.h.p. This is analagous to the core of an HRG, which is the clique formed by vertices of radius at most R/2, regardless of their angular coordinates.

Theorem 23.

Let G𝒢(n,2α+1,λ) be a threshold GIRG. The degeneracy is w.e.h.p.

κ(G)=(2±o(1))(2(1α))(1α)/(2α1)σ(G).

Proof.

We bound the maximal inner-degree |Γ(u)|; the statement then follows from Corollary 22. Notice that, independent of the geometric distance, a vertex with weight w is adjacent to any vertex with weight w if wnwλ since t(w,w)=λww2nλwnwλ2n=12 which is the maximal distance between two points in the unit torus. Thus, using β=2α+1, the probability that a vertex v is in the inner-neighbourhood of a vertex u with weight w is

(vΓ(u)) =({{u,v}E}{Wvw})
=(Wvnwλ)+2wnwλt(y,w)2αy(2α+1)𝑑y (t(w,n/wλ)=1/2)
=(nwλ)2α+2αλwnwnwλy2α𝑑y (Pareto and threshold)
=(nwλ)2α+2αλwn(2α1)(w12α(nwλ)12α) (y2α𝑑y=[y12α12α])
=(nwλ)2α+αα1/2(λw2(1α)n(nwλ)2α). (4)

Next, we calculate the value w, which maximises the expected inner-degree. To this end, we consider the probability measure of the inner-neighbourhood, take its derivative with respect to w and set it equal to 0. Differentiating yields

ddw(vΓ(u)) =2(nw)(2α+1)(1α)α(n2αw2λnw4αλ2α)2α1,

and solving for the maximum reveals w=(2(1α))14α2nλΘ(n).

We plug in w for the weight of u denoted by u into (vΓ(u)) and get by (5) that

(vΓ(u))=2(2(1α))(1α)/(2α1)(nλ)α.

Recalling that σ(G)=(1±o(1))λαn1α w.e.h.p., the upper bound now follows from the expectation of |Γ(u)| and applying a Chernoff bound in conjunction with a union bound. The lower bound is established by showing that there exists a vertex within the range of weights w~=[w(1+nα1log2(n))1/(2α),w] w.e.h.p. and lower bound the inner-degree of such vertex. Using the Pareto distribution and wΘ(n), we calculate the probability for a vertex u to belong to the range of weights w~. We obtain

(Wuw~) =(w(1+nα1log2(n))1/(2α)Wuw) (Range of w~)
=(Wuw)(Wuw(1+nα1log2(n))1/(2α))
=1(w)2α(1(w(1+nα1log2(n))1/(2α))2α) (Pareto)
=(w(1+nα1log2(n))1/(2α))2α(w)2α
=(w)2αnα1log2(n)
=Θ(1)log2(n)n. (wΘ(n))

By this we have 𝔼[|Vw~|]ω(log(n)). Using a Chernoff bound there exists a vertex within the desired weight range w~ w.e.h.p. To conclude the proof we lower bound the inner-degree of a vertex u~ included in the weight range w~. Note that w~=(1o(1))w=(2o(1)(1α))14α2nλ. We then have via Section 5

𝔼[|Γ(u~)|]=(n1)(vΓ(u~))(2o(1))(2(1α))(1α)/(2α1)λαn1α.

A final application of a Chernoff bound then ensures the concentration to finish the proof.

Comparing the lower bound of the degeneracy for GIRGs given in Theorem 23 to the upper bound of a HRG we obtained in Theorem 9 we draw the conclusion that the degeneracy-to-core ratio between the two models is fundamentally different.

Corollary 24 (GIRG-HRG degeneracy difference).

Fix an α(1/2,1). Let G𝒢(n,2α+1,λ) be a standard GIRG and H𝒢(n,α,C) be a threshold HRG. Then w.e.h.p.

|κ(G)σ(G)κ(H)σ(H)|Θ(1).

6 Conclusion

We have shown that the clique number, degeneracy, and chromatic number of HRGs are asymptotically (with small differences in the leading O-notation constants) as large as the core, though the clique number and degeneracy differ significantly. Our upper bound on the degeneracy provides a constant factor approximation algorithm for the graph colouring problem. The approximation ratio ranges from 2/3 to 4/3 and depends on the model parameter α. This raises several open questions and future research directions.

  • Is the chromatic number bounded away from the degeneracy, the clique number, or both?

  • Can HRGs be coloured optimally in polynomial time or is it NP-complete?

  • What are the asymptotics of ω(G)/σ(G)? Is the clique number a constant factor larger than the core and has similar behaviour as the degeneracy?

There are further directions of research such as determining other differences between HRGs and GIRGs or designing colouring algorithms for HRGs in various models of computation.

References

  • [1] Samuel Baguley, Yannic Maus, Janosch Ruff, and George Skretas. Hyperbolic random graphs: Clique number and degeneracy with implications for colouring. CoRR, 2024. doi:10.48550/arXiv.2410.11549.
  • [2] Leonid Barenboim and Michael Elkin. Distributed Graph Coloring: Fundamentals and Recent Developments. Morgan & Claypool Publishers, 2013. doi:10.1007/978-3-031-02009-4.
  • [3] Thomas Bläsius and Philipp Fischbeck. On the external validity of average-case analyses of graph algorithms. ACM Trans. Algorithms, 2024. doi:10.1145/3633778.
  • [4] Thomas Bläsius, Philipp Fischbeck, Tobias Friedrich, and Maximilian Katzmann. Solving vertex cover in polynomial time on hyperbolic random graphs. Theory Comput. Syst., 2023. doi:10.1007/S00224-021-10062-9.
  • [5] Thomas Bläsius, Tobias Friedrich, and Maximilian Katzmann. Efficiently approximating vertex cover on scale-free networks with underlying hyperbolic geometry. Algorithmica, 2023. doi:10.1007/S00453-023-01143-X.
  • [6] Thomas Bläsius, Tobias Friedrich, Maximilian Katzmann, Ulrich Meyer, Manuel Penschuck, and Christopher Weyand. Efficiently Generating Geometric Inhomogeneous and Hyperbolic Random Graphs. In ESA, 2019. doi:10.4230/LIPIcs.ESA.2019.21.
  • [7] Thomas Bläsius, Tobias Friedrich, Maximilian Katzmann, and Daniel Stephan. Strongly Hyperbolic Unit Disk Graphs. In STACS, 2023. doi:10.4230/LIPIcs.STACS.2023.13.
  • [8] Thomas Bläsius, Tobias Friedrich, and Anton Krohmer. Cliques in hyperbolic random graphs. Algorithmica, 2017. doi:10.1007/s00453-017-0323-3.
  • [9] Thomas Bläsius, Maximillian Katzmann, and Clara Stegehuis. Maximal cliques in scale-free random graphs. Network Science, 2024. doi:10.1017/nws.2024.13.
  • [10] Thomas Bläsius, Tobias Friedrich, and Anton Krohmer. Hyperbolic random graphs: Separators and treewidth. In ESA, 2016. doi:10.4230/LIPIcs.ESA.2016.15.
  • [11] Karl Bringmann, Ralph Keusch, and Johannes Lengler. Geometric inhomogeneous random graphs. Theoretical Computer Science, 2019. doi:10.1016/j.tcs.2018.08.014.
  • [12] Aleksander Bjørn Grodt Christiansen, Krzysztof Nowicki, and Eva Rotenberg. Improved dynamic colouring of sparse graphs. In STOC, 2023. doi:10.1145/3564246.3585111.
  • [13] Stephen A. Cook. The complexity of theorem-proving procedures. In STOC, 1971. doi:10.1145/800157.805047.
  • [14] B. V. Dekster. The Jung theorem for spherical and hyperbolic spaces. Acta Mathematica Hungarica, 1995. doi:10.1007/bf01874495.
  • [15] B. V. Dekster. The Jung theorem in metric spaces of curvature bounded above. Proceedings of the American Mathematical Society, 1997. doi:10.1090/s0002-9939-97-03842-2.
  • [16] Tobias Friedrich, Andreas Göbel, Maximilian Katzmann, and Leon Schiller. Cliques in high-dimensional geometric inhomogeneous random graphs. SIAM Journal on Discrete Mathematics, 2024. doi:10.1137/23m157394x.
  • [17] Tobias Friedrich and Anton Krohmer. On the Diameter of Hyperbolic Random Graphs. SIAM Journal on Discrete Mathematics, 2018. doi:10.1137/17M1123961.
  • [18] Mohsen Ghaffari and Christoph Grunau. Dynamic o(arboricity) coloring in polylogarithmic worst-case time. In STOC, 2024. doi:10.1145/3618260.3649782.
  • [19] Luca Gugelmann, Konstantinos Panagiotou, and Ueli Peter. Random Hyperbolic Graphs: Degree Sequence and Clustering. In ICALP, 2012. doi:10.1007/978-3-642-31585-5_51.
  • [20] Heinrich Jung. Ueber die kleinste Kugel, die eine räumliche Figur einschliesst. Journal für die reine und angewandte Mathematik, 1901. URL: http://eudml.org/doc/149122.
  • [21] Richard M. Karp. Reducibility among Combinatorial Problems, pages 85–103. Springer US, 1972. doi:10.1007/978-1-4684-2001-2_9.
  • [22] Maximilian Katzmann. About the analysis of algorithms on networks with underlying hyperbolic geometry. doctoralthesis, Universität Potsdam, 2023. doi:10.25932/publishup-58296.
  • [23] Ralph Keusch. Geometric Inhomogeneous Random Graphs and Graph Coloring Games. PhD thesis, ETH Zurich, 2018. doi:10.3929/ethz-b-000269658.
  • [24] Marcos Kiwi and Dieter Mitsche. A bound for the diameter of random hyperbolic graphs. In ANALCO, 2015.
  • [25] Marcos Kiwi and Dieter Mitsche. Spectral gap of random hyperbolic graphs and related parameters. The Annals of Applied Probability, 2018. doi:10.1214/17-aap1323.
  • [26] Marcos Kiwi and Dieter Mitsche. On the second largest component of random hyperbolic graphs. SIAM Journal on Discrete Mathematics, 2019. doi:10.1137/18m121201x.
  • [27] Marcos Kiwi, Markus Schepers, and John Sylvester. Cover and hitting times of hyperbolic random graphs. Random Structures & Algorithms, 2024. doi:10.1002/rsa.21249.
  • [28] Júlia Komjáthy and Bas Lodewijks. Explosion in weighted hyperbolic random graphs and geometric inhomogeneous random graphs. Stochastic Processes and their Applications, 2020. doi:10.1016/j.spa.2019.04.014.
  • [29] Dmitri Krioukov, Fragkiskos Papadopoulos, Maksim Kitsak, Amin Vahdat, and Marián Boguñá. Hyperbolic geometry of complex networks. Phys. Rev. E, 2010. doi:10.1103/PhysRevE.82.036106.
  • [30] Anton Krohmer. Structures & algorithms in hyperbolic random graphs. doctoralthesis, Universität Potsdam, 2016. URL: https://publishup.uni-potsdam.de/frontdoor/index/index/docId/39597.
  • [31] David W. Matula and Leland L. Beck. Smallest-last ordering and clustering and graph coloring algorithms. Journal of the ACM, 1983. doi:10.1145/2402.322385.
  • [32] Michael Mitzenmacher and Eli Upfal. Probability and Computing: Randomized Algorithms and Probabilistic Analysis. Cambridge University Press, 2005. doi:10.1017/CBO9780511813603.
  • [33] Tobias Müller and Merlijn Staps. The Diameter of KPKVB Random Graphs. Advances in Applied Probability, 2019. doi:10.1017/apr.2019.23.
  • [34] Fragkiskos Papadopoulos, Dmitri Krioukov, Marian Boguna, and Amin Vahdat. Greedy forwarding in dynamic scale-free networks embedded in hyperbolic metric spaces. In INFOCOM, 2010. doi:10.1109/infcom.2010.5462131.
  • [35] Jose L. Walteros and Austin Buchanan. Why is maximum clique often easy in practice? Operations Research, 2020. doi:10.1287/opre.2019.1970.
  • [36] David Zuckerman. Linear degree extractors and the inapproximability of max clique and chromatic number. Theory of Computing, 2007. doi:10.4086/toc.2007.v003a006.