Abstract 1 Introduction 2 Delayed decomposition 3 Delayed rank 4 Approximating the complete interval minor number References

A Polynomial-Time Approximation Algorithm for Complete Interval Minors

Romain Bourneuf ORCID Univ. Bordeaux, CNRS, Bordeaux INP, LaBRI, UMR 5800, F-33400 Talence, France Julien Cocquet ORCID ENS de Lyon, Université Claude Bernard Lyon 1, CNRS, Inria, LIP UMR 5668, Lyon, France Chaoliang Tang ORCID Shanghai Center for Mathematical Science, Fudan University, Shanghai, China
ENS de Lyon, Université Claude Bernard Lyon 1, CNRS, Inria, LIP UMR 5668, Lyon, France
Stéphan Thomassé ORCID ENS de Lyon, Université Claude Bernard Lyon 1, CNRS, Inria, LIP UMR 5668, Lyon, France
Abstract

As shown by Robertson and Seymour, deciding whether the complete graph Kt is a minor of an input graph G is a fixed parameter tractable problem when parameterized by t. From the approximation viewpoint, a substantial gap remains: there is no PTAS for finding the largest complete minor unless 𝖯=𝖭𝖯, whereas the best known result is a polytime O(n)-approximation algorithm by Alon, Lingas and Wahlén.

We investigate the complexity of finding Kt as interval minor in ordered graphs (i.e. graphs with a linear order on the vertices, in which intervals are contracted to form minors). Our main result is a polytime f(t)-approximation algorithm, where f is triply exponential in t but independent of n. The algorithm is based on delayed decompositions and shows that ordered graphs without a Kt interval minor can be constructed via a bounded number of three operations: closure under substitutions, edge union, and concatenation of a stable set. As a byproduct, graphs avoiding Kt as an interval minor have bounded chromatic number.

Keywords and phrases:
Approximation algorithm, Ordered graphs, Interval minors, Delayed decompositions
Category:
APPROX
Copyright and License:
[Uncaptioned image] © Romain Bourneuf, Julien Cocquet, Chaoliang Tang, and Stéphan Thomassé; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Approximation algorithms analysis
; Theory of computation Graph algorithms analysis ; Mathematics of computing Graph theory
Related Version:
Full Version: https://arxiv.org/abs/2505.05997
Acknowledgements:
We thank Julien Duron for his help with the tedious calculations.
Funding:
The authors are partially supported by the French National Research Agency under research grants ANR GODASse ANR-24-CE48-4377 and ANR Twin-width ANR-21-CE48-0014-01. The third author is also supported by Chinese Scholarship Council (CSC).
Editors:
Alina Ene and Eshan Chattopadhyay

1 Introduction

Complete minors in graphs form an extensively studied subject, notably featuring the fundamental result of Robertson and Seymour [27], which asserts that testing whether the complete graph Kt is a minor of an input graph G on n vertices can be done in time f(t)n3. This was proved in the series of Graph Minors papers, which also provided a decomposition theorem of Kt-minor-free graphs [28] (whose bounds were recently improved significantly by Gorsky, Seweryn and Wiederrecht [16]). One of the key basic facts allowing such a result is that Kt-minor-free graphs are sparse, i.e. they have a linear number of edges. Recently, Korhonen, Pilipczuk and Stamoulis [22] provided an algorithm running in time f(t)n1+o(1) for the same problem, improving on a quadratic algorithm from Kawarabayashi, Kobayashi and Reed [19]. From the approximation point of view, the landscape is much less understood. The first hardness result is due to Eppstein [11]: finding the size of a largest complete minor is 𝖭𝖯-hard. This was extended by Wahlén [29], who showed that there is no polynomial time approximation scheme (PTAS) for the size of a largest complete minor, unless 𝖯=𝖭𝖯. The current best known approximation factor achievable in polynomial time is O(n), as shown by Alon, Lingas and Wahlén [1]. Up to this day, it is still open whether this problem admits a polytime f(OPT) approximation algorithm. The goal of this paper is to investigate the corresponding problem for complete interval minors in ordered graphs.

An ordered graph (G,<) is a graph G=(V,E) equipped with a linear order < on its vertices. We will consistently denote the number of vertices by n, and the number of edges by m. We usually denote the vertices as v1,,vn, enumerated according to <. For simplicity of notation, we will sometimes omit the order < and talk about an ordered graph G. An interval minor of (G,<) is obtained by deleting some edges of G as well as iteratively contracting pairs of vertices which are consecutive in <. An example is displayed in Figure 1.

Figure 1: A K4 interval minor: the red zones represent the intervals we contract and the red edges are the remaining edges. Observe that this interval minor model is not a minor model since we contracted non-connected subsets of vertices.

The central computational problem on interval minors is the following:

Interval Minor Detection

Input: Two ordered graphs G and H.

Question: Is H an interval minor of G?

The complete study of the computational complexity of this question is probably very challenging, and a good first step is to limit ourselves (as for usual minors) to the case where H is a complete graph. The crucial fact to observe is that Kt-interval-minor-free ordered graphs can be very complex. Say that a bipartite ordered graph (G,<) with edges between two parts X and Y is monotone if either X<Y or Y<X. Then, observe that the monotone Kt,t does not contain K4 as an interval minor. In particular, K4-interval-minor-free ordered graphs can have a quadratic number of edges. They also form a large class of graphs (with growth at least 2n2/4) as they contain all subgraphs of the monotone Kn/2,n/2. Therefore, any attempt to effectively construct Kt-interval-minor-free ordered graphs must involve an operation allowing the creation of arbitrary monotone bipartite graphs. The landscape of Kt-interval-minor-free ordered graphs thus seems very different from the one of Kt-minor-free graphs. However, a strong connection exists between these two worlds: the right analogy with graph minors involves monotone-Kt,t-interval-minor-free ordered graphs. To see this, let us first recall the celebrated Marcus-Tardos theorem [24]:

Theorem 1 ([24]).

There exists a function f such that every ordered graph on n vertices with at least f(t)n edges contains a monotone Kt,t as an interval minor.

The notion of interval minors was formally introduced by Fox [12] to prove bounds on the function f of Theorem 1. The upper bound was later improved by Cibulka and Kynčl [10], again using interval minors. Theorem 1 gave a positive answer to a conjecture of Füredi and Hajnal [13], stating that for every permutation matrix P, there exists a constant cP such that every n×n binary matrix with at least cPn entries 1 contains P as a pattern. This result was later generalized to matrices in higher dimension by Klazar and Marcus [21]. The corresponding bound was later improved by Geneson and Tian [15], once more using interval minors. In [12], Fox suggested to conduct a thorough study of ordered graphs avoiding a given interval minor, with the hope of developing a theory analogue to the graph minor theory of Robertson and Seymour. There have been several papers going in that direction, see for instance [25, 23, 20]. In this paper, we continue this line of research by providing a decomposition theorem for Kt-interval-minor-free graphs.

The Marcus-Tardos theorem therefore implies that monotone-Kt,t-interval-minor-free ordered graphs are sparse, as are Kt-minor-free graphs (in the non-ordered case). A second analogy between these two classes appears from the computational angle:

Theorem 2 ([17, 5]).

There exists a function g such that testing whether an ordered graph (G,<) contains a monotone Kt,t as an interval minor can be done in time g(t)|V(G)|.

The cornerstone of this algorithm is the 𝖥𝖯𝖳 algorithm of Guillemot and Marx [17] for detecting a pattern in a permutation. Their algorithm amounts to detecting a monotone Kt,t as an interval minor in a monotone matching. The generalization to arbitrary ordered graphs follows from the notion of twin-width introduced in [6]. More specifically, approximating the twin-width of an ordered graph can be done in 𝖥𝖯𝖳 time, see [5]. In a nutshell, if an ordered graph (G,<) does not contain a monotone complete bipartite graph as an interval minor, then its twin-width is bounded, and therefore dynamic programming can test any first-order formula in 𝖥𝖯𝖳 time. This gives a win/win algorithm: if the twin-width is large compared to t, simply return Yes since (G,<) must contain Kt,t as interval minor, and if not, one can test if Kt,t is an interval minor of (G,<), as first-order logic can encode interval minors. Indeed, (G,<) contains an ordered graph (H,<) on vertex set u1<<uh if and only if there exist vertices x1,,xh1V(G) such that x1<<xh1, and for every edge uiujE(H), there exists yi,yjV(G) such that xi1<yixi, xj1<yjxj and yiyjE(G) (with the obvious adaptation for u1 and uh). Let us end this detour about Kt,t interval minors with the central open question in this topic: Given a graph G which admits an order < avoiding Kt,t as an interval minor, can we efficiently compute an order < such that (G,<) avoids Kf(t),f(t) as an interval minor. This question is equivalent to ask for a f(OPT)-approximation algorithm of the twin-width for sparse graphs (see Theorem 2.12 in [4]).

Let us go back to our original goal of approximately detecting Kt as an interval minor. We follow the classical strategy of providing an effective decomposition of Kt-interval-minor-free ordered graphs using a bounded number of tractable operations. The key-tool for this is the notion of delayed decomposition, used in [26] and formally defined in [7] to show that graphs of bounded twin-width are polynomially χ-bounded. It was later used as the central tool to show that pattern-avoiding permutations are the product of a bounded number of separable permutations [3]. We first introduce three operations on classes of ordered graphs. Given two classes 𝒞,𝒞 of ordered graphs, we define:

  • The substitution closure of 𝒞 is the class 𝒞¯ of ordered graphs which can be obtained from 𝒞 by iterating an arbitrary number of substitutions of elements of 𝒞. The substitution of a vertex v of an ordered graph (G,<) by an ordered graph (H,<) consists of replacing v by a copy of (H,<) and joining all neighbors of v to all vertices of H. The vertices of H remain ordered according to <, at the position of v in <.

  • The edge union of 𝒞 and 𝒞 is the class 𝒞𝒞 of all ordered graphs of the form (GG,<) where (G,<)𝒞 and (G,<)𝒞 are two graphs on the same set of vertices ordered by the same linear order <, and GG is their edge union.

  • The independent concatenation of 𝒞 is the class 𝒞+ of all ordered graphs (G+,<+) which can be obtained from a graph (G,<)𝒞 by adding an independent set I at the beginning or at the end of <, which can be connected arbitrarily to the vertices of G. Formally, either V(G)<+I or I<+V(G).

We define the rank of an ordered graph (G,<) inductively, as follows:

  • The empty graph has rank 0.

  • The rank of (G,<) is the smallest integer r such that (G,<) can be built from graphs of rank at most r1, using one of the following operations: substitution closure, edge union or independent concatenation.

Let us illustrate this parameter. The class of ordered edgeless graphs has rank 1, as they can be obtained by an independent concatenation from the empty graph. The class of monotone bipartite graphs has rank 2, as we can perform another round of independent concatenation. In particular, this means that the class of ordered graphs with rank at most 2 has unbounded twin-width, and even unbounded VC-dimension. The class of ordered complete graphs has rank at most 3, since the edge graph K2 is obtained at rank 2, and its closure under substitution gives all complete graphs. A slightly more involved example is that of the ordered path (Pn,<), in which each vertex is joined to its successor in <.

Note that the subgraph Mo of Pn consisting of the odd indexed edges is a matching consisting of edges e1<e3<, hence can obtained by substituting the vertices of an independent set by edges. In particular Mo has rank 3, and the even indexed edges subgraph Me also has rank 3. Finally, Pn=MoMe has rank 4. This is illustrated in Figure 2.

Figure 2: Construction showing that the ordered P6 has rank at most 4. The edge graph K2 has rank 2, which is why we go from rank 1 to rank 3.

The central result of this paper is the following:

Theorem 3.

Every Kt-interval-minor-free ordered graph (G,<) has bounded rank.

The cornerstone of Theorem 3 is the interplay between the two operations of substitution closure and edge union (as is the case for the construction of a path from two matchings). Before introducing the notion of delayed decomposition, let us first formalize what it means exactly that a graph G belongs to the substitution closure of a class 𝒞. The adequate representation of G must take into account that some vertices have been repeatedly substituted by a graph of 𝒞, hence G can be expressed as a tree of substitutions. Formally, a 𝒞-substitution tree of G is a rooted tree T whose leaves are the vertices of G. Moreover, every internal node x of T is labeled by a graph Gx of 𝒞, whose vertices are the children of x in T. Finally, two vertices u,v of G form an edge if and only if given their closest common ancestor x, the two children u,v of x which are the respective ancestors of u,v form an edge in Gx. By construction, the graphs in the substitution closure of 𝒞 are precisely the ones representable by a 𝒞-substitution tree. These structured trees were introduced by Gallai to study partial orders [14]. The most popular examples of substitution trees are those used to decompose cographs (P4-free graphs) using binary trees in which the graphs Gx are the edge and the non edge. These trees are well-suited for ordered graphs, as the order < can be represented as the left-to-right order on the leaves.

However, most graphs G do not admit a non-trivial substitution tree, i.e. one in which the root is not labeled by G. For instance paths of length at least 3 are indecomposable (or prime) with respect to substitutions. A simple way to strengthen substitution trees is to delay the effect of the graphs Gx: instead of creating edges between children, they act on their grandchildren. Let us formalize this:

A delayed structured tree is a rooted tree T whose leaves are the vertices of a graph G (the realization of T). Moreover, every internal node x of T is labeled by a quotient graph Gx whose vertices are the grandchildren of x in T. Finally, two vertices u,v of G are connected if and only if given their closest common ancestor x, the two grandchildren u,v of x which are the respective ancestors of u,v are connected in Gx. For an example, see Figure 4. Authorizing this delay results in a tool which is much more expressive than substitution trees. In particular, given a class of graphs 𝒞, the class 𝒞 of realizations of delayed structured trees whose quotient graphs Gx belong to 𝒞 is much harder to grasp than the mere substitution closure of 𝒞.

However, the class 𝒞 is not too complex compared to 𝒞: given the delayed structured tree T, consider the quotient graphs Gx of the nodes x with even depth, and of the ones with odd depth. Now form two trees Te and To from T, by setting in Te all quotient graphs Gx at odd depth as edgeless graphs, and setting in To all quotient graphs Gx at even depth as edgeless graphs. The key observation is that the realization Ge of Te belongs to the substitution closure of 𝒞. The same holds for Go, and thus G is the edge union of two graphs in the substitution closure of 𝒞.

In a nutshell, the proof of Theorem 3 is now straightforward: we just have to show that every ordered graph (G,<) without Kt-interval minor is the realization of a delayed structured tree whose quotient graphs are simpler than G, and that these simpler graphs are again decomposable into simpler objects, and that this process has a bounded height of recursion. This is the main technical part of the paper. Note that this process is too tame to create high-entropy objects such as monotone bipartite graphs. Here a choice has to be made: either declare that the basic class is indeed all monotone bipartite graphs, and just consider substitution closure and edge union, or set the empty graph as our basic class, and authorize independent set concatenation. We chose the latter convention as it fits more in our algorithmic purpose, but we feel that the former is more in the spirit of classical graph decompositions, with only two tame operations, and the basic class consisting of the indecomposable (or prime) monotone bipartite graphs.

The proof of Theorem 3 is algorithmic, and a sequence of operations for constructing (G,<) can be effectively computed in polynomial time by iterating delayed decompositions. This is the first phase of our algorithm to detect Kt as an interval minor: either we fail to achieve bounded rank and then find Kt as an interval minor, or we compute a sequence of operations achieving bounded rank for (G,<). Unfortunately, a graph can have bounded rank and still contain a Kt interval minor. Indeed, as we saw before, complete graphs have rank at most 3. The nice point is that large complete subgraphs are the only reason why an ordered graph of bounded rank can contain arbitrarily large complete interval minors. The second phase of the algorithm uses the decomposition of (G,<) and either outputs a Kt subgraph, or certifies that G has no Kf(t) interval minor. As a result, we show:

Theorem 4.

There is a triply exponential function f and a decision algorithm which, given as input an ordered graph (G,<) with n vertices and m edges and an integer t, satisfies the following:

  • If the algorithm returns Yes then (G,<) contains Kt as an interval minor.

  • If the algorithm returns No then (G,<) does not contain Kf(t) as an interval minor.

  • The algorithm runs in time O(tmn2).

Furthermore, when the algorithm returns Yes, it can also output a model of a Kt interval minor within the same running time.

The bound on the approximation factor is primarily due to the second phase of the algorithm. Specifically, we need to argue that if G has a huge complete interval minor and small rank, then we can detect a Kt subgraph in G.

First, observe that if G can be obtained by independent concatenation from G then G still contains a huge complete interval minor. Similarly, if G is the edge union of G1 and G2 then it follows from Ramsey’s theorem that one of G1 and G2 still contains a large complete interval minor. However, since we will apply this argument repeatedly, the standard version of Ramsey’s theorem would yield a tower-type bound. To avoid this, we establish a Ramsey-type result for interval minors (see Section 4.2), which eventually allows us to reduce the bound to exponential-type. Finally, if G is obtained by substitution closure then either G contains a clique which can be easily detected from the delayed decomposition tree, or one of the quotient graphs still contains a large complete interval minor (see Section 2.3).

Before diving into the details, we now discuss two hardness results in the context of Kt-minors.

Hardness of ordering a graph

A natural question regarding interval minors in ordered graphs is to find the minimum t, such that an input graph G admits an ordering < with no Kt interval minor. Denote this value by kim(G). This question is particularly interesting since the same problem for Kt,t instead of Kt amounts to approximating the twin-width of a sparse graph. Unfortunately the answer for Kt is as hard as approximating the chromatic number χ of a graph:

Lemma 5.

The parameters kim(G) and χ(G) are functionally equivalent.

Proof.

If a graph has chromatic number t, ordering its vertices according to the color classes directly gives an order without K2t interval minor. To see it, consider any partition of the vertices into 2t intervals. At most t1 intervals can contain vertices from two different color classes, so at least t+1 intervals must be entirely contained in a color class. By the pigeonhole principle, two distinct intervals must be contained inside the same color class, and therefore cannot be connected by an edge. Conversely, if a graph G has an ordering < such that (G,<) does not contain a Kt interval minor, Theorem 3 implies that (G,<) has bounded rank. Say that a class 𝒞 of graphs is χ-bounded if there exists a function h such that every graph G in 𝒞 satisfies χ(G)h(ω(G)), where ω(G) is the maximum size of a clique of G. We speak of polynomial χ-boundedness if h is a polynomial.

Claim.

For every r0, the class 𝒞r of ordered graphs with rank at most r is polynomially χ-bounded.

Proof of the Claim.

We prove it by induction on r. The statement is trivial for r=0 since the empty graph has clique number 0 and chromatic number 0. For the induction step, we just have to show that the three operations preserve polynomial χ-boundedness. Independent concatenation increases the chromatic number by at most one and cannot decrease the clique number. Similarly, the chromatic number of the edge union of two graphs is at most the product of their chromatic numbers, and the clique number of the edge union of two graphs is at least the maximum of the two clique numbers. The fact that the substitution closure preserves polynomial χ-boundedness was shown by Chudnovsky, Penev, Scott and Trotignon [9].

Then, if (G,<) does not contain a Kt interval minor then G does not contain a clique of size t, and therefore its chromatic number is bounded.

The fact that Kt-interval-minor-free ordered graphs have bounded chromatic number directly implies that the class of ordered graphs which do not contain some (arbitrary but fixed) ordered matching as a subgraph has bounded chromatic number (the two statements are in fact easily equivalent). This problem was introduced in [2] by Axenovich, Rollin and Ueckerdt, where they ask the question for some specific matchings. A far-reaching generalization was announced by Briański, Davies and Walczak [8]: for any ordered matching M, the class of ordered graphs which do not contain M as an induced ordered subgraph is χ-bounded. A study of the list chromatic number of ordered graphs avoiding an induced subgraph was also conducted by Hajebi, Li and Spirkl [18].

Hardness of detecting complete interval minors in ordered graphs

In our main result, we give a polynomial approximate algorithm to detect whether an ordered graph contains a Kt interval minor. For the hardness part, we show that deciding whether the complete interval minor number is at least t is 𝖭𝖯-complete in general. Formally, we consider the following problem.

Complete Interval Minor

Input: An ordered graph G and an integer t.

Question: Is Kt an interval minor of G?

Theorem 6.

Complete Interval Minor is 𝖭𝖯-complete.

Proof.

We first show that the problem is in 𝖭𝖯. To do so, observe that we can describe a Kt interval minor model of G by giving the intervals. Then, it can easily be checked in polynomial time that these intervals indeed form a model of G.

We show that the problem is 𝖭𝖯-hard by reduction from Clique. Consider an instance (G,k) of Clique: we are asked whether G contains a clique of size at least k. We build an instance (G^,t) of Complete Interval Minor as follows. Write V(G)={v1,,vn}. Consider the ordered graph (G^,<) defined as

  • V(G^)={v1,,vn}{u1,,un1}, where the ui are fresh vertices,

  • vivjE(G^) if and only if vivjE(G), viujE(G^) for every i[n] and j[n1], and uiujE(G^) for every ij[n1],

  • v1<u1<v2<u2<<vn1<un1<vn.

Informally, the ordered graph (G^,<) is obtained from G by arbitrarily ordering the vertices of G, and then adding a universal vertex between any two vertices of G. See Figure 3 for an illutration. Finally, we set t=n1+k. Note that ((G^,<),t) can be built from (G,k) in polynomial time.

We now prove that G has a clique of size k if and only if Kt is an interval minor of (G^,<). Suppose first that G has a clique of size k induced by the vertices vi1,,vik. Then, the vertices vi1,,vik,u1,,un1 induce a clique of size k+n1=t in G^, hence (G^,<) contains Kt as an interval minor. Conversely, suppose that Kt is an interval minor of (G^,<). Since there are only n1 vertices uj, there are at least t(n1)=k intervals in the interval model of Kt which do not contain a vertex uj. By definition of the order <, any interval that does not contain a vertex uj is of the form {vi}. Let {vi1},,{vik} be these k intervals. Then, the vertices vi1,,vik induce a clique of size k in G.

Figure 3: A graph G drawn with an arbitrary order and the corresponding ordered graph G^. The fresh vertices u1 and u2 are depicted in red, as well as the edges incident to them.

Perspectives

The obvious next step is of course to obtain a more digestible approximation factor than the triply exponential function f. We did not try very hard to show lower bounds, but failed to rule out constant factor approximation. A very natural problem to ask is the existence of an 𝖥𝖯𝖳 algorithm to find a Kt interval minor. This looks really challenging, as the only strategy seems to obtain a much more precise description of Kt-interval-minor-free ordered graphs which would allow dynamic programming. One of the main conclusions of this study is the versatility of delayed decompositions. It really suffices to apply them in a canonical way, and wonder whether the quotient graphs are simpler than the original one. For which (not necessarily ordered) graph classes does this machinery provide a bounded rank decomposition?

Organization of the paper

In Section 2, we formally define delayed decompositions, show that they can be computed efficiently, and review some of their properties. Then, in Section 3, we introduce a variant of rank tailored for algorithmic use, the delayed rank, which is based on delayed decompositions. We study some of its basic properties, before proving that ordered graphs with large delayed rank contain large complete interval minors. Finally, in Section 4, we present the algorithm of Theorem 4. To bound its approximation factor, we prove a Ramsey-type result for interval minors.

2 Delayed decomposition

This section focuses on the notion of delayed decomposition, a key tool for our approximation algorithm. We start by introducing delayed decompositions in Section 2.1. Then, in Section 2.2, we prove that all (ordered) graphs admit so-called distinguishing delayed decompositions, and that they can be computed in linear time. Finally, in Section 2.3, we establish a variant of Kőnig’s lemma on trees, and explore its consequences related to delayed decompositions.

2.1 Definition

We consider rooted trees T, and call the vertices of T nodes. The ancestors of a node x are the nodes in the unique path from x to the root r of T, and the parent of x is the first node on this path (other than x). The root has no parent. We also speak of grandparents, descendants, children, grandchildren, siblings (nodes with same parent), and cousins (non-sibling nodes with the same grandparent). The set of leaves of a tree T is denoted by L(T). For any node xV(T), we denote by L(x)L(T) the set of leaves which are descendants of x.

An ordered tree (T,<) is a rooted tree T equipped with a linear order < on L(T), such that for each node xV(T), the leaves L(x) form an interval of <. It is natural to think of < as a left-to-right order on the leaves of T. If (T,<) is an ordered tree and x,yV(T) are not in an ancestor–descendant relationship, then L(x) and L(y) are disjoint, and each of them is an interval for <, hence either L(x)<L(y), or L(y)<L(x). We then naturally extend < to x,y by x<y in the former case, and y<x in the latter. In particular, < induces a linear order <x on the children of any internal node x. Therefore, given a child y of a node x, we can speak of the predecessor and the successor of y, and of consecutive children of x. We can also consider the first child of x, which is the <x-minimum child of x, and the last child of x, defined analogously.

A delayed structured tree (T,<,{Gx}xV(T)) is an ordered tree (T,<), equipped with, for each node xV(T), a graph Gx on the grandchildren of x. We will often refer to these graphs Gx as the quotient graphs. This is analogous to the trees describing substitutions (module decomposition trees), except that the graphs Gx are defined on the grandchildren instead of the children, hence “delayed”. We add the technical requirement that each leaf is a single child (with no siblings), so that whenever xy are leaves, their closest ancestor is at distance at least 2.

The realization of a delayed structured tree (T,<,{Gx}xV(T)) is the ordered graph GT=((L(T),ET),<), where for two leaves x,y of T, we have the edge xyET if and only if xyE(Gz), where z is the closest ancestor of x,y, and x,y are the grandchildren of z which are ancestors of x,y respectively. If (G,<) is the realization (T,<,{Gx}xV(T)), we say that (T,<,{Gx}xV(T)) is a delayed decomposition of (G,<). See Figure 4 for an example.

Figure 4: A delayed structured tree and its realization. The edges in the realization are colored according to their corresponding quotient graph.

A delayed decomposition (T,<,{Gx}xV(T)) of an ordered graph (G,<) is distinguishing if it satisfies the following property: For every node x which has at least 3 children, for every two consecutive children x1,x2 of x, there exists some vertex vV(G)L(x) which is adjacent to all vertices in one of L(x1),L(x2) and to no vertex in the other.

As we shall see in Theorem 8, every graph has a distinguishing delayed decomposition, which can be computed in linear time. When we talk about the distinguishing delayed decomposition of a graph, we mean the one computed by Theorem 8.

The following result was already observed in [7].

Lemma 7 ([7, Lemma 2.1]).

Let (G,<) be the realization of a delayed structured tree (T,<,{Gx}xV(T)). Then, (G,<) is the edge union of two ordered graphs, each of which can be obtained by substitutions from the quotient graphs {Gx}xV(T).

2.2 Computing a distinguishing delayed decomposition

In this section, we prove that every ordered graph has a distinguishing delayed decomposition, which can be computed in linear time.

To discuss the running time of the algorithm, we first detail how we store ordered graphs. Throughout this article, we assume that the edge set is stored as a list of edges, and the vertex set as a list of vertices. There are two natural ways to store the order on the vertex set. The first one is to store the ordered set of vertices explicitly, for instance by assuming that the list of vertices is sorted according to the order. In particular, after a O(n)-time precomputation, all vertices can store their rank in the order. We call this representation explicit. The second one is to assume that the order can be determined from the labels of the vertices, which is the case for instance if the vertices are ordered by increasing labels. We call this representation implicit. Observe that an explicit representation can be computed in time O(nlog(n)) from an implicit representation by sorting the vertex set according to the order, and then storing the sorted list explicitly. Given an explicit representation of an ordered graph, we can relabel all the vertices in time O(n) so that the vertex set is [n], equipped with its natural linear order.

An algorithm to compute a distinguishing delayed decomposition of an ordered graph was described in [7], where delayed decompositions were formally introduced, and it was shown in [3] how to implement this algorithm in linear time in the context of permutations. A straightforward adaptation to the context of ordered graphs yields the following result.

Theorem 8 ([3]).

There is an algorithm which, given as input an explicit ordered graph (G,<) with n vertices and m edges, computes in time O(m+n) a distinguishing delayed decomposition (T,<,{Gx}xV(T)) of (G,<).

 Remark 9.

If (G,<) is an ordered graph with n vertices and m edges, the distinguishing delayed decomposition computed by Theorem 8 has the property that the total number of edges over all quotient graphs is at most m. Furthermore, the number of quotient graphs is equal to the number of nodes of T which have grandchildren, which is at most n1. By construction, for every node x of T and every child y of x, the children of y form an independent set in Gx.

2.3 A result on trees and its consequences

A classical result on trees states that every tree with a very large number of leaves contains either a node with large degree, or a path containing many nodes of degree at least 3. The goal of this section is to prove a generalization of this statement in the context of ordered trees. More precisely, consider an ordered tree whose leaves are grouped into disjoint intervals. We prove that if there is a very large number of intervals, either there is a node which “splits” a large number of intervals between its children, or there is a long path which progressively “peels” a large number of intervals.

Given an ordered tree (T,<), an interval family of (T,<) is a set of disjoint intervals of the linear order (L(T),<). A node xV(T) is b-branching if there is a subset of b intervals such that every interval of is included in L(x) and there is no child y of x such that L(y) intersects two intervals of . See Figure 5 for an illustration. An -interval path is a sequence of intervals I1,,I such that there exists nodes x1,,x of T such that:

  • L(xj) contains IjI for all j=1,,.

  • L(xj)Ij1= for all j=2,,.

Note that the nodes xj are pairwise distinct, and xj is a descendant of xi whenever i<j. Also, when j3, the parent p(xj) of xj satisfies L(p(xj))Ij2= since p(xj) is a (not necessarily proper) descendant of xj1.

Figure 5: A 4-branching vertex x, with being the set of red intervals, and a 4-interval path I1,I2,I3,I4.
Lemma 10.

Given an interval family of size 2(b+2)t of an ordered tree (T,<), there is either a b-branching node or a t-interval path in T.

Proof.

For every node x of T, let w(x) be the number of intervals of which are entirely contained inside L(x). Let B(x) be the set of children y of x such that w(y)0. Note that x is (at least) |B(x)|-branching, and so we can conclude if |B(x)|b. Therefore, we can assume that |B(x)|<b for every xV(T).

Consider a node x with w(x)2b2. Let 1 be the set of intervals which are contained inside some L(y) for yB(x) and 2 be the set of intervals which are contained inside L(x) but are not in 1. Note that if |2|2b1, the node x is b-branching. Indeed, in that case it suffices to order 2 with respect to <, and to select every other interval to form a b-branching subfamily .

So there is a child y of x such that (using that w(x)2b2 in the last inequality):

w(y)|1|/|B(x)|(w(x)2b+1)/(b1)w(x)/b.

Starting from the root r, we define a sequence of vertices (r=x1,,xt) forming a descending chain in T, such that w(xi+1)w(xi)3 for every i[t1] and w(xi)2(b+2)t+1i for every i[t].

Start by setting x1 to be the root r, and note that w(x1)=2(b+2)t. Suppose that we already defined x1,,xi1 satisfying the desired property, with it. Consider a descendant xi of xi1 with maximum w(xi), such that w(xi)w(xi1)3, and among them one which is as close to the root as possible in T. By definition of xi, we have:

w(p(xi))w(xi1)22(b+2)t+2i22(b+2)222b2.

Thus, p(xi) has a child y with

w(y)w(p(xi))/b(w(xi1)2)/b.

Therefore, by maximality of w(xi), we have:

w(xi)(w(xi1)2)/b(2(b+2)t+2i2)/b2(b+2)t+1i.

For every i[t1], since w(xi+1)w(xi)3, there are at least three intervals which are entirely contained in L(xi) and which are not entirely contained in L(xi+1). Therefore, there is an interval Ii which is entirely contained in L(xi) and such that IiL(xi+1)=. Finally, w(xt)2(b+2)>0 so there exists an interval ItL(xt).

If (T,<,{Gx}xV(T)) is a delayed structured tree, a leaf y of T is h-heavy if there are at least h ancestors x of y which are not isolated in the graph Gp2(x), where p2(x) is the grandparent of x in T.

Lemma 11.

Let (T,<,{Gx}xV(T)) be a delayed structured tree and (GT,<) be its realization. Let be an interval family of (T,<) forming a complete interval minor in (GT,<). If has a (2t1)-interval path in T, then there is a (2t3)-heavy leaf in (T,<,{Gx}xV(T)).

Proof.

Let I1,,I2t1 be the intervals of the (2t1)-interval path of in T, and x1,,x2t1 be nodes of T such that for all j[2t1], L(xj) contains IjI2t1, and for all j[[2,2t1]], L(xj)Ij1=.

Let yL(T) be any leaf in I2t1. We show that y is (2t3)-heavy. Let i[2t3]. Since forms a complete interval minor in (GT,<), there is an edge between some vertex yiIi and I2t1. Let zi be the deepest node of T such that I2t1{yi}L(zi). Then, zi is a descendant of xi and a proper ancestor of xi+1 (hence all zi are distinct). Furthermore, since zi is a proper ancestor of xi+1 then the grandchild wi of zi such that yL(wi) is an ancestor of xi+2. In particular, I2t1L(x2t1)L(xi+2)L(wi). Let wi be the grandchild of zi such that yiL(wi). By maximality of the depth of zi, wi and wi are cousins in T. Since there is an edge between yi and I2t1 in GT, there is an edge between wi and wi in Gzi. Thus, we found 2t3 distinct ancestors w of y which are not isolated in Gp2(w), i.e. y is (2t3)-heavy.

Lemma 12.

Let (T,<,{Gx}xV(T)) be a delayed structured tree containing a (2t3)-heavy leaf. Then, its realization (GT,<) contains a clique of size t as a subgraph.

Proof.

Let y be a (2t3)-heavy leaf of T. Let x1,,x2t3 be ancestors of y such that each xi is not isolated in the graph Gp2(xi). Up to renaming them, we can assume that xi is a proper ancestor of xj whenever i<j. We build a sequence (y0,,yt1) of vertices of GT with the following properties:

  • For every i[[0,t2]], yi is adjacent to every vertex in L(x2i+1), and

  • For every i[t1], yiL(x2i1).

Note that these properties immediately imply that y0,,yt1 induce a clique of size t in GT.

Let i[[0,t2]] and let xi be a neighbor of x2i+1 in the graph Gp2(x2i+1). Let yi be an arbitrary vertex in L(xi). Then, yi is adjacent to all vertices in L(x2i+1). If i>0 then p2(x2i+1) is a descendant (not necessarily proper) of x2i1 so yiL(x2i1). Finally, we choose yt1 to be any vertex of L(x2t3).

3 Delayed rank

The notion of rank is natural and convenient for structural analysis, but not so much for algorithmic purposes, since there does not seem to be a simple way of computing it, or even approximating it. To overcome this issue, we introduce an alternative to the rank, based on delayed decompositions, which we call the delayed rank. This delayed rank can easily be computed, and shares some key structural properties with the rank. It is the key notion for our approximation algorithm.

3.1 Definition and first properties

We add some information to the delayed structured trees. Consider a delayed structured tree (T,<,{Gx}xV(T)). We label the nodes x of T as follows (see Figure 6):

  • If x does not have a grand-parent, we label it with (the first two layers of the tree).

  • Else, if there is a cousin x of x such that x<x and there is an edge xx in Gp2(x) then we label x with R.

  • Otherwise, if there is a cousin x of x such that x<x and there is an edge xx in Gp2(x) then we label x with L.

  • Else, we label x with O.

Figure 6: Labeling of a delayed structured tree.
Lemma 13.

Let (T,<,{Gx}xV(T)) be a distinguishing delayed decomposition of an ordered graph (G,<). If x is a node of T with at least 3 children, there are no consecutive children y1,y2 of x both labelled with O.

Proof.

Let x be such a node, and y1,y2 be consecutive children of x. Since the delayed decomposition (T,<,{Gx}xV(T)) is distinguishing, there is some vertex vV(G)L(x) which is adjacent to all vertices in one of L(y1),L(y2) and none in the other. Let u1 be any vertex in L(y1) and u2 any vertex in L(y2). Since vL(x) then u1 and v have the same last common ancestor z as u2 and v, and z is a proper ancestor of x. Suppose by contradiction that z is a proper ancestor of p(x). Then, the grandchild of z which is an ancestor of u1 is also the grandchild of z which is an ancestor of u2, call it x. Let v be the ancestor of v which is a grandchild of z. If xvE(Gz) then u1vE(G) and u2vE(G), a contradiction. Similarly, if xvE(Gz) then u1vE(G) and u2vE(G), again a contradiction. Therefore, z is the parent of x. Let v be be the ancestor of v which is a grandchild of z. Then, v is a cousin of y1 and y2. If i{1,2} is such that v is adjacent to all vertices in L(yi) then vyiE(Gp2(yi)) so one of y1 and y2 is not labelled with O.

In view of Lemma 13, if x is a node of T with at least 3 children, and y is a child of x labelled O which is not the first child of x, then the predecessor y of y has either label L or label R. If y has label L, we refine the label of y to OL and if y has label R, we refine the label of y to OR.

Given a node x of a delayed structured tree, if y is a vertex of Gx whose parent is labelled with R then there exists some vertex y>V(Gx) which is adjacent to y. This observation will be our main tool to prove that graphs with large delayed rank have large complete interval minors (see Theorem 19). For this reason, we define the type of a node xV(T) as the label of its parent in T.

The delayed rank of an ordered graph (G,<) is defined as follows:

  • If (G,<) is monotone bipartite, (G,<) has delayed rank 0,

  • Otherwise, we compute the distinguishing delayed decomposition of (G,<). For each quotient graph (Gx,<), if (Gx,<) is monotone bipartite we say that (Gx,<) is a refined quotient graph of (G,<). Otherwise, note that x has at least 3 children (since there are only edges between cousins in Gx). In that case, if the first child of x is labelled with O, we remove all its children from Gx. Then, we partition the vertices into types R,L,OR and OL. Set R={R,OR} and L={L,OL}. We partition the edges into four types, depending on whether the type of their left endpoint is in R or L, and similarly for their right endpoint, denote these four types by RR,RL,LR and LL. The refined quotient graphs of (Gx,<) are then defined as follows:

    • The graph induced by the edges RR, to which we remove all vertices before the first vertex of type R and after the last vertex of type R.

    • The graph induced by the edges RL, to which we remove all vertices before the first vertex of type L and after the last vertex of type L.

    • The graph induced by the edges LR, to which we remove all vertices before the first vertex of type R and after the last vertex of type R.

    • The graph induced by the edges LL, to which we remove all vertices before the first vertex of type L and after the last vertex of type L.

    Then, the delayed rank of (G,<) is 1 more than the maximum delayed rank of all its refined quotient graphs.

Observe first that siblings form an independent set in Gx, so removing all the children of the first child of x is just removing an independent set, at the beginning of the order of Gx. Similarly, the vertices before the first vertex of type R are at the beginning of the order of Gx and all have type either L or OL, so they cannot induce any edge of type RR or LR. As for the vertices after the last vertex of type R, they are at the end of the order of Gx and all have type either L,OL or OR, and in the latter case they are siblings and come before the vertices of type L and OL. Thus, all these vertices cannot induce any edge of type RR or LR. Similar observations hold for the vertices before the first vertex of type L and after the last vertex of type L.

We now give some simple properties of the delayed rank. The first one gives a simple way of computing the delayed rank of an ordered graph. We first need some definitions.

If G is an ordered graph, we define a sequence (𝒢i(G))i0 of sets of graphs as follows. Set 𝒢0(G)={G}. Suppose that 𝒢i(G) is defined. Then, for every H𝒢i(G), add all the refined quotient graphs of H to 𝒢i+1(G).

The following results are immediate from the definitions, and can be proved easily by induction on r.

Lemma 14.

An ordered graph G has delayed rank at least r if and only if 𝒢r(G).

Lemma 15.

If G,H are ordered graphs such that H𝒢r(G) for some r0 then H is a subgraph of G.

The following result bounds the size of each 𝒢r(G). It will be important in the analysis of the running time of the algorithm of Theorem 4.

Lemma 16.

If G is an ordered graph with n vertices and m edges, for every r1, the total number of edges of all graphs in 𝒢r(G) is at most m. Thus, for every r1, there are at most 4mn graphs in 𝒢r(G).

Proof.

We prove the first property by induction on r. For r=1, Remark 9 implies that the total number of edges of all quotient graphs of G is at most m. Since the refined quotient graphs of G are obtained from the quotient graphs by partitioning the edges and removing some edges, the total number of edges of all graphs in 𝒢1(G) is at most m. Suppose now that the property holds for some r1. By definition, 𝒢r+1(G)=H𝒢r(G)𝒢1(H). By the induction hypothesis at rank 1, if H has m edges then the total number of edges of all graphs in 𝒢1(H) is at most m. By the induction hypothesis at rank r, the total number of edges over all H𝒢r(G) is at most m. Thus, the total number of edges of all graphs in 𝒢r+1(G) is at most m. This proves the first part of the statement.

If m=0 then G is monotone bipartite so 𝒢1(G)= and |𝒢1(G)|4mn. Otherwise, by Remark 9, G has at most n1 quotient graphs, each of which gives rise to at most 4 refined quotient graphs, so |𝒢1(G)|4(n1)4mn. Now, suppose r>1. By the first part of the statement, there are at most m graphs H𝒢r1(G) with at least one edge. All other graphs in 𝒢r1(G) are monotone bipartite, hence they have no refined quotient graphs. If H𝒢r1(G) then H has at most n vertices by Lemma 15 so, by Remark 9, H has at most n quotient graphs. Each of these quotient graphs gives rise to at most 4 refined quotient graphs, so |𝒢1(H)|4n. Taking the union over all H𝒢r1(G) which have at least one edge, we get |𝒢r(G)|4mn.

By Lemma 7, delayed substitution closure can alternatively be seen as substitution closure followed by edge union. Therefore, a simple inductive proof shows that the notion of delayed rank is simply a refinement of the notion of rank.

Lemma 17.

If G has delayed rank r then G has rank at most 7r+2.

3.2 Delayed rank and complete interval minors

We now prove the key result about graphs of large delayed rank: they contain large complete interval minors. The next lemma is the engine of our proof, it allows us to measure the progress we do in building the large complete interval minor when going from delayed rank r to delayed rank r+1.

We consider looped interval minors: an ordered graph (H,<) (possibly with a loop on each vertex) with vertex set v1<<vh is an interval minor of an ordered graph (G,<) if there exists a partition of V(G) into intervals I1,,Ih such that whenever vivjE(H), there is an edge in G between Ii and Ij. A looped Kt is an ordered clique of size t, with a loop on every vertex. A left-lazy looped Kt is an ordered clique of size t, with a loop on every vertex, except the first one. A right-lazy looped Kt is an ordered clique of size t, with a loop on every vertex, except the last one. A lazy looped Kt is an ordered clique of size t, with a loop on every vertex, except the first one and the last one. See Figure 7 for an illustration of all these graphs.

Figure 7: The various looped interval minors we consider. The arrows indicate the possible outcomes of Lemma 18.
Lemma 18.

Let r1 and let 𝒞r be the class of ordered graphs with delayed rank r. Then, the following assertions hold.

  1. 1.

    If every ordered graph in 𝒞r contains a looped Kt interval minor, then every ordered graph in 𝒞r+1 contains a left-lazy or a right-lazy looped Kt+1 interval minor.

  2. 2.

    If every ordered graph in 𝒞r contains a left-lazy or a right-lazy looped Kt interval minor, then every ordered graph in 𝒞r+1 contains either a lazy looped Kt+1 interval minor or a looped Kt interval minor.

  3. 3.

    If every ordered graph in 𝒞r contains a lazy looped Kt interval minor, then every ordered graph in 𝒞r+1 contains a left-lazy or a right-lazy looped Kt interval minor.

Proof.

Consider an ordered graph (G,<)𝒞r+1. Compute the delayed decomposition of (G,<). By definition, there is a refined quotient graph H of G some type (RR,RL,LR or LL) which is in 𝒞r. We consider several cases depending on the type of H.

  • If H has type RR then H is the graph induced by the edges RR, to which we remove all vertices before the first vertex of type R and after the last vertex of type R. Suppose that H contains a lazy looped Kt interval minor, and let I1<<It be the intervals corresponding to this interval minor. Recall that (I1,,It) is a partition of V(H). We argue that each interval Ij contains a vertex of type R.

    Since the first and the last vertex of H are vertices of type R, this is true for I1 and It. Consider now any Ij with 1<j<t. Since the Kt interval minor is looped, there is an edge between two vertices of Ij, which means that there are two vertices in Ij which are of type R but don’t have the same parent (since siblings form an independent set in the quotient graphs). Therefore, there is a vertex of type R in Ij. Let It+1={xV(G):V(H)<x}. If yjIj is a vertex of type R, it follows from the definition of the type R that yj has a neighbor in It+1.

    Using It+1, we can then extend any looped Kt to a right lazy looped Kt+1, any left-lazy looped Kt to a lazy looped Kt+1, any right-lazy looped Kt to a looped Kt and any lazy looped Kt to a right-lazy looped Kt.

  • If H has type RL then H is the graph induced by the edges RL, to which we remove all vertices before the first vertex of type L and after the last vertex of type L. Suppose that H contains a lazy looped Kt interval minor, and let I1<<It be the intervals corresponding to this interval minor. Recall that (I1,,It) is a partition of V(H). We argue that each interval Ij contains a vertex of type L.

    Since the first and the last vertex of H are vertices of type L, this is true for I1 and It. Consider now any Ij with 1<j<t. Since the Kt interval minor is looped, there is an edge between two vertices of Ij, which means that there exist u<vIj such that the type of u is in R={R,OR}, and the type of v is in L={L,OL}. Thus, there exist consecutive vertices u<vIj such that the type of u is in R and the type of v is in L. This implies that the type of v is L. Therefore, there is a vertex with type L in Ij. Let I0={xV(G):x<V(H)}. If yjIj is a vertex of type L, it follows from the definition of the type L that yj has a neighbor in I0.

    Using I0, we can extend any looped Kt to a left lazy looped Kt+1, any left-lazy looped Kt to a looped Kt, any right-lazy looped Kt to a lazy looped Kt+1 and any lazy looped Kt to a left-lazy looped Kt.

  • The case LR is similar to the case RL (except that we prove that each Ij contains a vertex of type R), and the case LL is similar to the case RR (except that we prove that each Ij contains a vertex of type L).

The result then follows immediately since we can extend any lazy looped Kt interval minor of H as desired, and since H𝒞r.

We easily deduce the main result of this section from Lemma 18.

Theorem 19.

Every ordered graph with delayed rank at least 3r2 contains a Kr interval minor.

Proof.

We prove by induction on r that every ordered graph with delayed rank at least 3r2 has a looped Kr interval minor. For r=1, if G has delayed rank at least 1 then G contains at least one edge so G has a looped K1 interval minor. Suppose that the property holds for 3r2. By Lemma 18, every graph of delayed rank at least 3r1 contains a left-lazy or a right-lazy looped Kt+1 interval minor. Applying Lemma 18 again, every graph of delayed rank at least 3r contains either a lazy looped Kt+2 interval minor or a looped Kt+1 interval minor. Using Lemma 18 once more, every graph of delayed rank at least 3r+1 contains either a left-lazy or a right-lazy looped Kt+2 interval minor, which itself contains a looped Kt+1 interval minor.

Now, Theorem 3 follows immediately from combining Theorem 19 with Lemma 17.

4 Approximating the complete interval minor number

In this section, we present an algorithm to approximate the size of a largest complete interval minor in an ordered graph G. In Section 4.1, we give some implementation details on parts of the algorithm. In Section 4.2, we give a Ramsey-type theorem in the context of interval minors, which is crucial for bounding the approximation factor of the algorithm. In Section 4.3, we describe and analyse the algorithm.

4.1 More algorithmic tools

By Theorem 8, the delayed decomposition of an explicit ordered graph with n vertices and m edges can be computed in time O(n+m). From there, it is simple to compute the refined quotients in the same running time.

Lemma 20.

There is an algorithm which, given as input an explicit ordered graph G with n vertices and m edges, computes all the refined quotients of G in time O(n+m).

Proof.

First, we start by checking whether G is monotone bipartite. To do so, we iterate over all edges to find the largest left endpoint of an edge, call it , and the smallest right endpoint of an edge, call it r. If <r then G is monotone bipartite and has no refined quotients, so we return the empty set. Otherwise, G is not monotone bipartite.

In that case, we compute the delayed decomposition (T,<,{Gx}xV(T)) of G in time O(n+m) using Theorem 8, with all the Gx stored explicitly. Then, we compute the label of every node x (L,R or O) by looking at its neighborhood in the graph Gp2(x). Thus, in time O(n+m), we can compute the label of all nodes xV(T). From this, we can compute the type (L,R,OL,OR or O) of every node of T, in time O(n+m) overall. Then, for each quotient graph Gx, we check whether it is monotone bipartite. If so, we add it to the set of refined quotient graphs. Otherwise, we continue. Using the same method as above, this can be done in time linear in the size of Gx, so in time O(n+m) overall.

We then remove from each Gx the first vertex as long as it is of type O. Then, for every edge of every Gx, we can compute its type by looking at the types of its endpoints. Once this is done, it is simple to compute the graphs induced by each of the four edge types. Finally, removing all vertices up to the first vertex of some type and after the last vertex of some type can also be done in time O(n+m) overall. Note that the refined quotients are all stored explicitly. The total size of a delayed structured tree (T,<,{Gx}xV(T)) is |V(T)|+xV(T)|E(Gx)|. The next result follows from a simple top-down dynamic programming algorithm.

Lemma 21.

There is an algorithm which, given as input a delayed structured tree (T,<,{Gx}xV(T)) of total size s, decides in time O(s) whether there is a h-heavy leaf in T.

4.2 A Ramsey-type result for complete interval minors

Before we move on to the algorithm, we present a Ramsey-type result for complete interval minors, which will be crucial for the bound on the “approximation factor” of our algorithm. Ramsey’s theorem states that in every red/blue coloring of the edges of Kn, there is a monochromatic clique of size log(n)/2, and this bound is essentially tight up to constant factors. We consider the analog question in the context of interval minors for ordered graphs. A red/blue coloring of the edges of the ordered graph Kn contains a red (resp. blue) complete interval minor of size t if there exists a partition of V(G) into t intervals such that there is a red (resp. blue) edge between any two of these intervals. A monochromatic complete interval minor of size t is a red or a blue complete interval minor of size t.

We now prove that the bounds are much better for monochromatic complete interval minors than for monochromatic cliques.

Lemma 22.

For every red/blue coloring of the edges of the ordered graph Kn, there is a monochromatic complete interval minor of size 2log(n)1.

Proof.

Suppose that there does not exist a red complete interval minor of size 2log(n)1. We prove that for every ilog(n)1, we can find 2i intervals, each of size at least n/2ilog(n), with only blue edges between any two of these intervals. The property is trivial for i=0. Suppose that the property holds for some i<log(n)1. Then, there exist 2i intervals, each of size at least n/2ilog(n), with only blue edges between any two of them. Consider one such interval I, and cut it into 2log(n)1 subintervals, each of size at least |I|/2log(n)1. Then, each subinterval has size at least |I|/2log(n)1n/2(i+1)log(n)11n/2(i+1)log(n). Indeed:

n2(i+1)log(n)11n2(i+1)log(n) 2n2(i+1)log(n)2(i+1)log(n)n2(i+1)log(n)
2n2(i+1)log(n)n
n2(i+1)log(n)
log(n)(i+1)log(n)
ilog(n)1.

Since there does not exist a red complete interval minor of size 2log(n)1, there are two subintervals with no red edge between them, hence only blue edges between them. Doing this in each of the 2i intervals, we find 2i+1 subintervals, each of size at least n/2(i+1)log(n), with only blue edges between any two of them. This result for i=log(n)1 proves the existence of a blue complete interval minor of size 2log(n)1.

The previous bound is almost sharp. To see it, set q=2log(n)loglog(n) and consider a random red/blue edge coloring of the edges of the ordered graph Kq. Then, taking lexicographic products of this graph, it can be shown that for n large enough, with high probability, this yields a red/blue coloring of the edges of the ordered graph Kn where the largest monochromatic complete interval minor has size 22log(n)loglog(n).

4.3 The algorithm

We now move to the description of the algorithm. We will need the following technical lemma to bound the “approximation factor” of our algorithm.

Lemma 23.

For every t1, there exists a function ht:{0,1,,3t2}+, such that the following holds:

  1. 1.

    f:tht(0) satisfies f(t)=222O(t),

  2. 2.

    For every r[3t2], ht(r)2log((ht(r1)/2)1/(2t1)4)114,

  3. 3.

    For every r{0,,3t2}, ht(r)4.

Proof.

For every integer t1 and every integer r{0,,3t2}, set

gt(r)=43t2r+43t2r13log(512t)

and ht(r)=22g(r). It can be verified that the functions ht have the desired properties.

We are now ready to prove Theorem 4. In the following statement, we assume that the ordered graph (G,<) is given explicitly.

Theorem 4. [Restated, see original statement.]

There is a triply exponential function f and a decision algorithm which, given as input an ordered graph (G,<) with n vertices and m edges and an integer t, satisfies the following:

  • If the algorithm returns Yes then (G,<) contains Kt as an interval minor.

  • If the algorithm returns No then (G,<) does not contain Kf(t) as an interval minor.

  • The algorithm runs in time O(tmn2).

Furthermore, when the algorithm returns Yes, it can also output a model of a Kt interval minor within the same running time.

Proof.

The high-level description of the algorithm is extremely simple: we compute the delayed rank of (G,<). If this rank is at least 3t2, we return Yes. Otherwise, we look whether there is a (2t3)-heavy leaf in one of the delayed structured trees that we computed. If so, we return Yes. Otherwise, we return No.

More formally, we compute the sets 𝒢0(G),,𝒢3t2(G). If 𝒢3t2(G), we return Yes. Otherwise, if there exists some H𝒢0(G)𝒢3t2(G) whose delayed decomposition tree contains a (2t3)-heavy leaf then we return Yes. Otherwise, we return No.

For every t1, let ht be the function provided by Lemma 23. Then, define f:tht(0), and note that f(t)=222O(t).

Claim.

If the algorithm returns Yes then (G,<) has a Kt interval minor.

Proof of the Claim.

By Lemma 14, if 𝒢3t2(G) then G has delayed rank at least 3t2 so Theorem 19 implies that G has a Kt interval minor. If there exists some H𝒢0(G)𝒢3t2(G) whose delayed decomposition tree contains a (2t3)-heavy leaf then by Lemma 15, H is a subgraph of G and H is the realization of a delayed structured tree that contains a (2t3)-heavy leaf so by Lemma 12, H contains a clique of size t. Therefore, G itself contains a clique of size t, hence a Kt interval minor.

Claim.

If G has a Kf(t) interval minor then the algorithm returns Yes.

Proof of the Claim.

Suppose that there is no graph H𝒢0(G)𝒢3t2(G) whose delayed decomposition tree contains a (2t3)-heavy leaf (otherwise we return Yes). We prove by induction that for every 0r3t2, there exists a graph Hr𝒢r(G) which has a complete interval minor of size ht(r). In particular, this will imply that 𝒢3t2(G), so the algorithm returns Yes. For r=0, set H0=G𝒢0(G), which has a complete interval minor of size f(t)=ht(0) by assumption. Suppose that the property holds for some r{0,1,,3t3}. Consider an interval family which realizes the complete interval minor of size ht(r) in Hr. In particular, since ht(r)4 then Hr is not bipartite so the refined quotient graphs of Hr are in 𝒢r+1(G). Let (T,<,{Gx}xV(T)) be the delayed decomposition of Hr. Set br=(ht(r)/2)1/(2t1)2, so that ht(r)=2(br+2)2t1. By Lemma 11, since Hr doesn’t contain a (2t3)heavy leaf, there is no (2t1)-interval path in Hr. Then, Lemma 10 implies that there is a br-branching node x in T, which means that the corresponding quotient graph Gx has a complete interval minor of size br. After possibly removing the first vertices of type O, the resulting subgraph of Gx still has a complete interval minor of size br2 since we only removed an independent set. Applying Lemma 22 twice, we get that one of the types RR,RL,LR and LR satisfies that the subgraph induced by the edges of this type has a complete interval minor of size at least 2log(br2)11. Removing the first vertices and the last vertices, which both form a stable set, decreases the size of the complete interval minor by at most 4, so one of the refined quotient graphs Hr+1 of Hr has a complete interval minor of size at least 2log(br2)114=2log((ht(r)/2)1/2t14)114ht(r+1) by definition of ht. Thus, Hr+1𝒢r+1(G) has a complete interval minor of size at least ht(r+1), as desired.

Claim.

This algorithm can be implemented to run in time O(tmn2).

Proof of the Claim.

By Lemma 20, the set of refined quotient graphs (stored explicitly) of an explicit ordered graph with v vertices and e edges can be computed in time O(v+e). Thus, 𝒢1(G) can be computed in time O(n+m), with all graphs being stored explicitly.

Suppose that we already computed 𝒢r(G) for some 1r<3t2, with all graphs being stored explicitly. Then, Lemmas 15 and 16 imply that it contains at most 4mn graphs, each with at most n vertices, and with at most m edges in total over all graphs in 𝒢r(G). Thus, 𝒢r+1(G) can be computed in time H𝒢r(G)O(|V(H)|+|E(H)|)=O(4mnn+m)=O(mn2), with all graphs being stored explicitly. Therefore, computing 𝒢0(G)𝒢3t2(G) can be done in time O(tmn2).

Finally, given an explicit ordered graph with v vertices and e edges, its delayed decomposition can be computed in time O(v+e), hence the corresponding delayed structured tree has total size O(v+e). Then, by Lemma 21, the algorithm can compute in time O(v+e) whether there is a (2t3)-heavy leaf in it. Thus, by iterating over all graphs in 𝒢0(G)𝒢3t2(G), the algorithm can check whether there exists some H𝒢0(G)𝒢3t2(G) whose delayed structured tree contains a (2t3)-heavy leaf. Since by Remark 9 there are O(tmn) such graphs and together they contain at most O(tm) edges, this can be done in time

H𝒢0(G)𝒢3t2(G)O(|V(H)|+|E(H)|)=O(tmnn+tm)=O(tmn2).

For the last statement, observe that the proofs of Lemmas 12, 18, and 19 are all algorithmic and can be implemented efficiently. Note also that in the course of the algorithm, we compute all the graphs in 𝒢0(G)𝒢3r2(G) and their delayed decompositions.

References

  • [1] Noga Alon, Andrzej Lingas, and Martin Wahlen. Approximating the maximum clique minor and some subgraph homeomorphism problems. Theoretical Computer Science, 374(1):149–158, 2007. doi:10.1016/j.tcs.2006.12.021.
  • [2] Maria Axenovich, Jonathan Rollin, and Torsten Ueckerdt. Chromatic number of ordered graphs with forbidden ordered subgraphs. Combinatorica, 38(5):1021–1043, October 2018. doi:10.1007/s00493-017-3593-0.
  • [3] Édouard Bonnet, Romain Bourneuf, Colin Geniet, and Stéphan Thomassé. Factoring pattern-free permutations into separable ones. In Proceedings of the 2024 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 752–779, 2024. doi:10.1137/1.9781611977912.30.
  • [4] Édouard Bonnet, Colin Geniet, Eun Jung Kim, Stéphan Thomassé, and Rémi Watrigant. Twin-width II: Small classes. Combinatorial Theory, June 2022. doi:10.5070/C62257876.
  • [5] Édouard Bonnet, Ugo Giocanti, Patrice Ossona de Mendez, Pierre Simon, Stéphan Thomassé, and Szymon Toruńczyk. Twin-width IV: Ordered graphs and matrices. Journal of the ACM, 71(3), June 2024. doi:10.1145/3651151.
  • [6] Édouard Bonnet, Eun Jung Kim, Stéphan Thomassé, and Rémi Watrigant. Twin-width I: Tractable FO model checking. Journal of the ACM, 69(1), November 2021. doi:10.1145/3486655.
  • [7] Romain Bourneuf and Stéphan Thomassé. Bounded twin-width graphs are polynomially χ-bounded. Advances in Combinatorics, 2025(2):19p, 2025. doi:10.19086/aic.2025.2.
  • [8] Marcin Briański, James Davies, and Bartosz Walczak. Colouring graphs without an induced ordered matching, In preparation.
  • [9] Maria Chudnovsky, Irena Penev, Alex Scott, and Nicolas Trotignon. Substitution and χ-boundedness. Journal of Combinatorial Theory, Series B, 103(5):567–586, 2013. doi:10.1016/j.jctb.2013.02.004.
  • [10] Josef Cibulka and Jan Kynčl. Better upper bounds on the Füredi-Hajnal limits of permutations. In Proceedings of the 2017 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 2280–2293, 2017. doi:10.1137/1.9781611974782.150.
  • [11] David Eppstein. Finding large clique minors is hard. Journal of Graph Algorithms and Applications, 13(2):197–204, February 2009. doi:10.7155/jgaa.00183.
  • [12] Jacob Fox. Stanley-Wilf limits are typically exponential, 2013. arXiv:1310.8378.
  • [13] Zoltán Füredi and Péter Hajnal. Davenport-Schinzel theory of matrices. Discrete Mathematics, 103(3):233–251, 1992. doi:10.1016/0012-365X(92)90316-8.
  • [14] Tibor Gallai. Transitiv orientierbare graphen. Acta Mathematica Academiae Scientiarum Hungarica, 18(1):25–66, March 1967. doi:10.1007/BF02020961.
  • [15] Jesse T. Geneson and Peter M. Tian. Extremal functions of forbidden multidimensional matrices. Discrete Mathematics, 340(12):2769–2781, 2017. doi:10.1016/j.disc.2017.08.017.
  • [16] Maximilian Gorsky, Michał T. Seweryn, and Sebastian Wiederrecht. Polynomial bounds for the graph minor structure theorem, 2025. doi:10.48550/arXiv.2504.02532.
  • [17] Sylvain Guillemot and Daniel Marx. Finding small patterns in permutations in linear time. In Proceedings of the 2014 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 82–101, 2014. doi:10.1137/1.9781611973402.7.
  • [18] Sepehr Hajebi, Yanjia Li, and Sophie Spirkl. List-3-coloring ordered graphs with a forbidden induced subgraph. SIAM Journal on Discrete Mathematics, 38(1):1158–1190, 2024. doi:10.1137/22M1515768.
  • [19] Ken ichi Kawarabayashi, Yusuke Kobayashi, and Bruce Reed. The disjoint paths problem in quadratic time. Journal of Combinatorial Theory, Series B, 102(2):424–435, 2012. doi:10.1016/j.jctb.2011.07.004.
  • [20] Vít Jelínek and Stanislav Kučera. On the structure of matrices avoiding interval-minor patterns. Advances in Applied Mathematics, 101:70–99, 2018. doi:10.1016/j.aam.2018.07.005.
  • [21] Martin Klazar and Adam Marcus. Extensions of the linear bound in the Füredi–Hajnal conjecture. Advances in Applied Mathematics, 38(2):258–266, 2007. doi:10.1016/j.aam.2006.05.002.
  • [22] Tuukka Korhonen, Michał Pilipczuk, and Giannos Stamoulis. Minor containment and disjoint paths in almost-linear time. In 2024 IEEE 65th Annual Symposium on Foundations of Computer Science (FOCS), pages 53–61, Los Alamitos, CA, USA, October 2024. IEEE Computer Society. doi:10.1109/FOCS61266.2024.00014.
  • [23] Yaping Mao, Hongjian Lai, Zhao Wang, and Zhiwei Guo. Interval minors of complete multipartite graphs, 2015. arXiv:1508.01263.
  • [24] Adam Marcus and Gábor Tardos. Excluded permutation matrices and the Stanley–Wilf conjecture. Journal of Combinatorial Theory, Series A, 107(1):153–160, 2004. doi:10.1016/j.jcta.2004.04.002.
  • [25] Bojan Mohar, Arash Rafiey, Behruz Tayfeh-Rezaie, and Hehui Wu. Interval minors of complete bipartite graphs. Journal of Graph Theory, 82(3):312–321, 2016. doi:10.1002/jgt.21903.
  • [26] Michał Pilipczuk and Marek Sokołowski. Graphs of bounded twin-width are quasi-polynomially χ-bounded. Journal of Combinatorial Theory, Series B, 161:382–406, 2023. doi:10.1016/j.jctb.2023.02.006.
  • [27] N. Robertson and Paul D. Seymour. Graph Minors. XIII. The disjoint paths problem. Journal of Combinatorial Theory, Series B, 63(1):65–110, 1995. doi:10.1006/jctb.1995.1006.
  • [28] Neil Robertson and Paul D. Seymour. Graph Minors. XVI. Excluding a non-planar graph. Journal of Combinatorial Theory, Series B, 89(1):43–76, 2003. doi:10.1016/S0095-8956(03)00042-X.
  • [29] Martin Wahlén. On the complexity of approximating the Hadwiger number. Theoretical Computer Science, 410(8):994–996, 2009. doi:10.1016/j.tcs.2008.12.020.