Faster Edge Coloring by Partition Sieving
Abstract
In the Edge Coloring problem, we are given an undirected graph with vertices and edges, and are tasked with finding the smallest positive integer so that the edges of can be assigned colors in such a way that no two edges incident to the same vertex are assigned the same color. Edge Coloring is a classic NP-hard problem, and so significant research has gone into designing fast exponential-time algorithms for solving Edge Coloring and its variants exactly. Prior work showed that Edge Coloring can be solved in time and polynomial space, and in graphs with average degree in time and exponential space, where .
We present an algorithm that solves Edge Coloring in time and polynomial space. Our result is the first algorithm for this problem which simultaneously runs in faster than time and uses only polynomial space. In graphs of average degree , our algorithm runs in time, which has far better dependence in than previous results. We also consider a generalization of Edge Coloring called List Edge Coloring, where each edge in the input graph comes with a list of colors, and we must determine whether we can assign each edge a color from its list so that no two edges incident to the same vertex receive the same color. We show that this problem can be solved in time and polynomial space. The previous best algorithm for List Edge Coloring took time and space.
Our algorithms are algebraic, and work by constructing a special polynomial based off the input graph that contains a multilinear monomial (i.e., a monomial where every variable has degree at most one) if and only if the answer to the List Edge Coloring problem on the input graph is YES. We then solve the problem by detecting multilinear monomials in . Previous work also employed such monomial detection techniques to solve Edge Coloring. We obtain faster algorithms both by carefully constructing our polynomial , and by improving the runtimes for certain structured monomial detection problems using a technique we call partition sieving.
Keywords and phrases:
Coloring, Edge coloring, Chromatic index, Matroid, Pfaffian, Algebraic algorithmFunding:
Shyan Akmal: Partially funded by the Ministry of Education and Science of Bulgaria (support for INSAIT, part of the Bulgarian National Roadmap for Research Infrastructure).Copyright and License:
![[Uncaptioned image]](x1.png)
2012 ACM Subject Classification:
Theory of computation Graph algorithms analysis ; Theory of computation Parameterized complexity and exact algorithmsAcknowledgements:
We thank the anonymous reviewers for helpful feedback on this work.Editors:
Olaf Beyersdorff, Michał Pilipczuk, Elaine Pimentel, and Nguyễn Kim ThắngSeries and Publisher:

1 Introduction
Coloring graphs is a rich area of research in graph theory and computer science. Given a graph , a proper vertex coloring of is an assignment of colors to its nodes such that no two adjacent vertices receive the same color. In the Vertex Coloring problem, we are given an undirected graph on vertices, and are tasked with computing the smallest positive integer such that admits a proper vertex coloring using colors. Vertex Coloring is a classic combinatorial problem, NP-hard to solve even approximately [27]. An influential line of research has worked on designing faster and faster exponential-time algorithms for this problem [24, 17, 13, 8], culminating in an algorithm solving Vertex Coloring in time and space [11], where throughout we write as shorthand for .
In this paper, we study exact algorithms for the closely related Edge Coloring problem. Given a graph , a proper edge coloring of is an assignment of colors to its edges such that no two edges incident to a common vertex receive the same color. The smallest positive integer such that admits a proper edge coloring using colors is the chromatic index of , denoted by . In the Edge Coloring problem, we are given an undirected graph on vertices and edges, and are tasked with computing .
Let denote the maximum degree of the input graph . For each vertex in , a proper edge coloring of must assign distinct colors to the edges incident to . Thus . A classic result in graph theory known as Vizing’s theorem proves that in fact , so that the simple lower bound is close to the truth. Moreover, a proper edge coloring of using colors can be found in near-linear time [2]. Consequently, solving Edge Coloring amounts to distinguishing between the cases and . Despite this powerful structural result, the Edge Coloring problem is NP-hard just like Vertex Coloring, even for graphs where all vertices have degree exactly [21].
In recent years, researchers have made major strides in our knowledge of algorithms for approximate variants of Edge Coloring in the distributed [19, 5, 3] and dynamic settings [15, 14, 6]. In comparison, much less is known about the exponential-time complexity for solving the Edge Coloring problem exactly.
A simple way to solve Edge Coloring is by reduction to Vertex Coloring. Given the input graph , we can construct in polynomial time the line graph whose nodes are edges of , and whose nodes are adjacent if the corresponding edges in are incident to a common vertex. By construction, the solution to Vertex Coloring on gives the solution to Edge Coloring on . If has edges and maximum degree , then has vertices and maximum degree . Using the aforementioned time and space algorithm for Vertex Coloring, this immediately implies a time and space algorithm for Edge Coloring. In graphs with maximum degree , Vertex Coloring can be solved in time and space, where [9]. By applying the above reduction, Edge Coloring can similarly be solved in time and space, for .
Line graphs are much more structured objects than generic undirected graphs (for example, there are small patterns that line graphs cannot contain as induced subgraphs [4, Theorem 3]), so one might hope to solve Edge Coloring faster without relying on a black-box reduction to Vertex Coloring.
Nonetheless, there is only one result from previous work which solves Edge Coloring on general undirected graphs without using this reduction.
Specifically, [10, Section 5.6] shows that Edge Coloring can be solved in time and polynomial space.
In comparison, it remains open whether Vertex Coloring can be solved in time using only polynomial space.
In summary, algorithms from previous work
can solve Edge Coloring in
-
time using polynomial space, or in
-
faster than time in bounded degree graphs, using exponential space.
Given this state of affairs, it is natural to ask:
Can Edge Coloring be solved in faster than time, using only polynomial space?
We answer this question affirmatively by proving the following result.
Theorem 1.1.
There is a randomized algorithm which solves Edge Coloring with high probability and one-sided error in time and polynomial space.
If the input graph has average degree , then , so Theorem 1.1 equivalently states that Edge Coloring can be solved in time for . In comparison, the current fastest algorithm for Vertex Coloring on graphs with average degree takes time, where [18, Lemma 4.4 and Theorem 5.4]. Prior to our work, the fastest algorithm for Edge Coloring in the special case of regular graphs (i.e., graphs where all vertices have the same degree) took time [10, Theorem 6]. The regularity assumption simplifies the Edge Coloring problem significantly, as in this case the problem reduces to finding a collection of mutually disjoint perfect matchings, which together consist of only edges. Theorem 1.1 improves upon this runtime with an algorithm that succeeds on all undirected graphs, not just regular graphs.
In the case of regular graphs, we obtain improvements beyond Theorem 1.1 when the graphs have high degree. We say a graph is -regular if all its nodes have degree . Given a positive integer , let
(1) |
denote the Harmonic number. We prove the following result.
Theorem 1.2.
For each integer , there is a randomized algorithm that solves Edge Coloring with high probability and one-sided error on -regular graphs in time and polynomial space, for .
For example for , Theorem 1.2 shows that we can solve Edge Coloring in 6-regular graphs in time, slightly faster than the runtime presented in Theorem 1.1. The savings in the exponent are better for larger , approaching a runtime of as the degree grows since we have .
We also consider a generalization of Edge Coloring called List Edge Coloring. In this problem, we are given an undirected graph with edges as before, together with a positive integer and a list of possible colors for each edge . We are tasked with determining whether admits a proper edge coloring, which assigns each edge a color from the list . If we set to be the maximum degree of and set for each edge , we recover the standard Edge Coloring problem. Using more sophisticated arguments, we generalize Theorem 1.1 to the List Edge Coloring problem.
Theorem 1.3.
There is a randomized algorithm that solves List Edge Coloring with high probability and one-sided error in time and polynomial space.
Prior to our work, no algorithm was known for List Edge Coloring which simultaneously ran in faster than time and used polynomial space.
1.1 Our Techniques
Our algorithms for Edge Coloring and List Edge Coloring are algebraic, and involve reducing these problems to certain monomial detection tasks, by now a common paradigm in graph algorithms. To achieve the runtime and polynomial space bounds for these problems, we combine and improve techniques in polynomial sieving, matroid constructions, and enumeration using matrices.
The main technical ingredient in our work builds off the determinantal sieving technique introduced in [16]. Given a homogeneous multivariate polynomial of degree and a linear matroid on of rank , this technique allows one to test in time whether contains a multilinear monomial (i.e., a monomial where each variable has degree at most one) that forms a basis in . We recall the formal definition of matroids in Section 2. For the purpose of this overview, a matroid is a family of subsets of variables of satisfying certain properties, and a basis in is a set in this family whose size equals the rank . In our applications, we construct to enumerate certain structures in the input graph. Finding a multilinear monomial in corresponds to identifying a particularly nice structure in the graph (e.g., a solution to the Edge Coloring problem). Determinantal sieving is useful because one can pick in such a way that identifying a monomial in that forms a basis in recovers extra information about the relevant structures in the graph.
We build off determinantal sieving and introduce the partition sieving method. This technique involves a partition matroid , which is a simple type of matroid defined using an underlying partition of the variable set of . Given a partition matroid whose partition has parts, and a polynomial compatible with (“compatible” is a technical condition defined in Section 3 – intuitively, it requires that monomials in have degrees respecting the partition defining ), we can test whether contains a multilinear monomial forming a basis in in time and polynomial space. This offers an exponential speed-up over determinantal sieving. To apply partition sieving, we need to carefully design both the polynomial and the partition matroid to be compatible in our technical framework.
We design the polynomial using the notion of the Pfaffian of a matrix. If is a symbolic skew-symmetric adjacency matrix of , then its Pfaffian is a polynomial that enumerates the perfect matchings in . The Edge Coloring problem may be rephrased as asking whether the edge set of can be partitioned into a collection of matchings. One can design a polynomial that enumerates matchings by computing a product of Pfaffians of adjacency matrices. However, this simple construction does not meet the technical conditions we need to employ partition sieving. Instead we utilize the Ishikawa-Wakayama formula [22], a convolution identity for Pfaffians and determinants, to design a polynomial that enumerates matchings which satisfy certain degree constraints. Specifically, the formula lets us construct a polynomial whose monomials correspond to -tuples of matchings in , where each vertex in appears as an endpoint in exactly of the matchings in the tuple.
We design the partition matroid in terms of a dominating set of the input graph. A subset of vertices is a dominating set in if every vertex in is adjacent to a node in . Given a dominating set in , we show that it is always possible to construct a partition matroid with parts, compatible with our enumerating polynomial (in our proofs, is actually only compatible with a polynomial obtained from by the inclusion-exclusion principle, but we ignore that distinction in this overview). We achieve compatability using the degree condition imposed on , mentioned at the end of the previous paragraph. Partition sieving then lets us detect a multilinear monomial in , and thus solve Edge Coloring, in time.
Lemma 1.4.
There is a randomized algorithm that, given a graph on vertices and edges and a dominating set of , solves Edge Coloring on with high probability and one-sided error in time and polynomial space.
In other words, by finding smaller dominating sets in , we can solve Edge Coloring faster. In polynomial time, it is easy to construct a dominating set of size at most in any connected graph with vertices. Consequently, the partition sieving framework lets us solve Edge Coloring in time. We obtain faster algorithms using the fact that if a graph with nodes has minimum degree two, then it contains a dominating set of size at most [29]. We make this result effective, showing how to find such a dominating set in time and polynomial space. When solving Edge Coloring, it turns out one can assume without loss of generality that the graph has minimum degree two, and so this framework yields a time algorithm for the problem. The situation is far more complex for List Edge Coloring, where unit-degree vertices cannot be removed freely. For this problem, we construct a more complicated polynomial , by enforcing additional matroid conditions in the Ishikawa-Wakayama formula. Intuitively, we identify a subgraph of minimum degree two in , and use to determine whether admits a list edge coloring which can be extended to all of .
Organization
In Section 2 we review notation, basic assumptions, and useful facts about polynomials, matrices, and matroids. In Section 3 we introduce our partition sieving method for monomial detection. In Section 4 we present a generic algorithmic template for solving coloring problems using polynomials. In Section 5 we apply this framework to design our algorithms for Edge Coloring and prove Theorems 1.1 and 1.2. We conclude in Section 6 with a summary of our work and some open problems. Our List Edge Coloring algorithm and the proof of Theorem 1.3 is deferred to the full version of this paper.
2 Preliminaries
General Notation
Given a positive integer , we let denote the set of the first consecutive positive integers. Given a set and positive integer , we let denote the family of subsets of of size . Given a positive integer , we let denote the Harmonic number, defined in Equation 1.
Graph Notation and Assumptions
Throughout, we consider graphs which are undirected. We let and denote the number of vertices and edges in the input graphs for the problems we study. We assume that the input graphs are connected, so that . This is without loss of generality, since when solving Edge Coloring and List Edge Coloring on disconnected graphs, it suffices to solve the problems separately on each connected component. Given a graph and vertex , we let denote the number of neighbors of (i.e., the degree of in ). For the List Edge Coloring problem where each edge is assigned a list of colors from a subset of , we assume that for all vertices in , since if some vertex has degree greater than , its incident edges must all be assigned distinct colors in a proper edge coloring, and the answer to the List Edge Coloring problem is trivially no.
We let denote the set of perfect matchings of – partitions of the vertex set of into parts of size two, such that each part consists of two adjacent vertices.
Dominating Sets
A dominating set in a graph is a subset of vertices of such that every vertex in is either in or adjacent to a node in . We collect some useful results relating to finding small dominating sets in graphs.
Lemma 2.1 (Simple Dominating Set).
There is a polynomial time algorithm which takes in a connected graph on nodes and outputs a dominating set for the graph of size at most .
Lemma 2.2 (Dominating Sets in Dense Graphs [12, 29]).
If has nodes, and every vertex in has degree at least two, then has a dominating set of size at most .
Lemma 2.3 (Dominating Sets in Regular Graphs [1, 33]).
Any -regular graph on vertices has a dominating set of size at most .
Lemma 2.4 (Dominating Set Algorithm).
Let be a graph obtained by starting with , and repeatedly deleting vertices of degree one until no vertices of degree one remain. Suppose at most unit-degree vertices were deleted from to produce . Then a minimum-size dominating set of can be found in time and polynomial space.
Lemma 2.4 is proved in the full version of this paper.
Finite Field Arithmetic
Throughout, we work with polynomials and matrices over a field of characteristic two, where is sufficiently large. Arithmetic operations over take time. All elements satisfy . Consequently, given any element , we can compute its unique square root as by repeated squaring in time.
Polynomials
Let be a polynomial over a set of variables . A monomial is a product for some nonnegative integers . A monomial is multilinear if for each . We say appears in the polynomial if the coefficient of in is nonzero. We define the support of to be the set of indices
of variables which appear with positive degree in .
Similarly, we define the odd support of to be the set
of indices of variables which appear with odd degree in . By definition, for all monomials .
The degree of a monomial is defined to be Given , we define the degree of restricted to the variables indexed by to be
The total degree (or just degree) of is defined to be
where the maximum is taken over all monomials with nonzero coefficient in . We say that a polynomial is homogeneous if for all monomials in .
We recall the Lagrange Interpolation formula, which allows one to recover coefficients of a low degree univariate polynomial from a small number of evaluations of that polynomial.
Proposition 2.5 (Lagrange Interpolation).
Let be a univariate polynomial of degree less than . Suppose that for distinct . Then
So given evaluations of at distinct points, we can compute all the coefficients of in polynomial time.
We also record the following observation, proven for example in [30, Theorem 7.2], which shows that to test whether a low degree polynomial is nonzero, it suffices to check whether a random evaluation of over a sufficiently large field is nonzero.
Proposition 2.6 (Schwartz-Zippel Lemma).
Let be a nonzero polynomial over a finite field of degree at most . If each variable of is assigned an independent, uniform random value from , then the corresponding evaluation of is nonzero with probability .
Matrices
Given a matrix and sets of rows and columns , we let denote the submatrix of restricted to rows in and columns in . We let denote the submatrix of restricted to rows in but using all columns, and denote the submatrix of restricted to columns in but using all rows. If the rows and column sets of are identical, we write as shorthand for . We let denote the all-zeros matrix, whose dimensions will always be clear from context. We let denote the transpose of . A matrix is symmetric if , and skew-symmetric if it is symmetric and has zeros along its main diagonal (we employ this definition because we work over characteristic two).
If is a square matrix, we let denote the determinant of .
Given a set , a perfect matching on is a partition of into sets of two elements each. We let denote the set of perfect matchings on . Given a skew-symmetric matrix with rows and columns indexed by a set , we define the Pfaffian of to be
(2) |
The standard definition for Pfaffians involves a sign term (see e.g., [31, Section 7.3.2]), which we ignore here because we work over a field of characteristic two. It is well known that for all skew-symmetric , we have (see e.g., [31, Proposition 7.3.3] for a proof). Consequently, the Pfaffian of matrix over the field can be computed in time by using the equation .
We record the following useful facts for computation with Pfaffians and determinants.
Proposition 2.7 (Direct Sums).
For any skew-symmetric matrices and , we have
over a field of characteristic two.
Proposition 2.7 is proved in the full version of this paper.
Proposition 2.8 (Ishikawa-Wakayama Formula).
For a skew-symmetric -matrix and a -matrix with , we have
Proposition 2.8 is due to [22] (see [23, Section 4.3] for a recent alternate proof).
Matroids
A matroid is a pair , consisting of a ground set and a collection of subsets of called independent sets, satisfying the following axioms:
-
1.
.
-
2.
If and , then .
-
3.
For any with , there exists an element such that .
A maximal independent set in a matroid is called a basis. It is well known that all bases of have the same size , which is referred to as the rank of . We say spans if is a superset of a basis of . If there exists a matrix with column set such that for every , the set is independent in if and only if has full rank, then we say that is a linear matroid represented by . Provided we work over a large enough field of size , it is known that if is linear, then it can be represented by a matrix whose number of rows equals the rank of (see e.g., [26] or [28, Section 3.7]).
Given matroids over disjoint ground sets respectively, we define the direct sum
of these matroids to be the ground set together with the family of subsets of with the property that is independent in for all . It is well known that is a matroid as well, whose rank is equal to the sum of the ranks of the .
The following result lets us represent direct sums of linear matroids [28, Section 3.4].
Proposition 2.9 (Matroid Sum).
Given linear matroids represented by matrices respectively, their direct sum is represented by the block-diagonal matrix
In our Edge Coloring algorithm, we make use of two simple types of matroids: uniform and partition matroids. Given a ground set and a nonnegative integer , the uniform matroid over of rank has as its independent sets all subsets of of size at most .
Proposition 2.10 (Uniform Matroid Representation [28, Section 3.5]).
Given a set , a nonnegative integer , and a field with , we can construct an matrix representing the uniform matroid over of rank in time.
Given a set partitioned into parts , and a list of nonnegative integer capacities , the partition matroid over the ground set for this partition and list of capacities has as its independent sets all subsets satisfying for all . Equivalently, is the direct sum of uniform matroids of rank over for all . By definition, has rank . By Propositions 2.9 and 2.10 we see that we can efficiently construct representations of partition matroids [28, Proposition 3.5].
Proposition 2.11 (Partition Matroid Representation).
Given a partition matroid over of rank , and a field with , we can construct an matrix representing in time.
3 Polynomial Sieving
In this section, we develop faster algorithms for a certain monomial detection problem. This is our main technical result, which allows us to obtain faster algorithms for Edge Coloring and List Edge Coloring. We recall the following well-known sieving technique:
Proposition 3.1 (Inclusion-Exclusion).
Let be a polynomial over a field of characteristic two, with variable set . Fix . Let be the polynomial consisting of all monomials (with corresponding coefficients) in that are divisible by . Then
where is the polynomial obtained from by setting for all .
See [35, Lemma 2] for a proof of Proposition 3.1.
Combining Propositions 2.5 and 3.1, we obtain the following result.
Lemma 3.2 (Coefficient Extraction).
Let be a polynomial over a field of characteristic two, with variable set . For , let be the coefficient of in , viewed as a polynomial over the variable set . Then we can evaluate at a point as a linear combination of evaluations of .
Proof.
Fix a subset . Let be the polynomial with variable set obtained by substituting each variable with in . Let be the coefficient of in . By definition of , is the polynomial obtained by taking the monomials in whose degree restricted to equals . By Proposition 2.5, we can compute at any point as a linear combination of evaluations of .
By definition, is the polynomial obtained by taking monomials of that are divisible by , and then setting all variables in equal to . So by Proposition 3.1 we can evaluate as the sum of evaluations of . Then by the conclusion of the previous paragraph, we get that we can evaluate at any point as a linear combination of evaluations of as claimed.
The inclusion-exclusion principle allows us to sieve for the multilinear term in a polynomial that is the product of all its variables. This was extended to the parameterized setting, where the goal is to sieve for multilinear terms of degree in time [7, 10, 36], and then further extended to multilinear monomial detection in the matroid setting [16].
Lemma 3.3 (Basis Sieving [16, Theorem 1.1]).
Let be a homogeneous polynomial of degree over a field of characteristic 2, and let be a matroid on of rank , with a representation over given. Then there is a randomized algorithm running in time and polynomial space, which determines with high probability if contains a multilinear monomial such that is a basis of .
Moreover, [16] introduced a variant of sieving that takes the odd support into account. This tool is key to our exponential speed-up for the Edge Coloring problem.
Lemma 3.4 (Odd Sieving [16, Theorem 1.2]).
Let be a polynomial of degree over a field of characteristic 2, and let be a matroid on of rank , represented by a matrix over . Then using evaluations of and polynomial space, we can determine if contains a monomial such that spans .
Let be a homogeneous polynomial of degree in variables . Consider a partition into parts . Let be integers with
Let be the partition matroid defined by the partition of into the parts and the capacities for . We say is compatible with if for every and every monomial in , the degree of restricted to equals .
The following result gives an improvement over Lemma 3.3 for polynomials that are compatible with partition matroids.
Theorem 3.5 (Partition Sieving).
Let be a polynomial of degree over a field of characteristic 2 with , and let be a partition matroid on of rank with partition classes. If is compatible with , then there is a randomized algorithm which runs in time and polynomial space that determines with high probability if contains a monomial such that is a basis of .
Proof.
Let be the partition and be the capacities defining the partition matroid . By assumption, for all . Let be the partition matroid over defined by the same partition as , but with capacities for . From this construction, has rank .
Our algorithm to determine whether contains a monomial whose support spans works as follows. By Proposition 2.11, a linear representation of over can be computed in polynomial time. We run the odd sieving algorithm of Lemma 3.4 with and to determine in time whether contains a monomial whose odd support spans . We return YES if such a monomial exists, and return NO otherwise. We now prove that this algorithm is correct.
First, suppose that has a monomial whose support is a basis of . By the compatibility condition, is multilinear, which implies that . Since is a basis of , and thus spans , it follows that spans .
Conversely, suppose has a monomial such that spans . We claim that spans . Indeed, since spans , the set contains at least variables of for each . Hence its superset contains at least variables of each part as well.
Claim 3.6.
The set contains exactly variables of each part .
Proof.
Suppose to the contrary that does not contain variables of for some index . Then contains exactly variables of . Since is compatible with , the monomial has degree exactly when restricted to . The only way this is possible is if , and has exactly variables in of degree one and a single variable in with degree two. The variable of degree two does not show up in the odd support of , which implies that has at most elements, which contradicts the assumption that spans .
By 3.6, the set has exactly variables from for each . Since is compatible with , this implies that has degree one in every variable in . Hence spans , as claimed.
This shows that the algorithm is correct, and proves the desired result.
4 Coloring Algorithm Template
Recall that in the List Edge Coloring problem we are given a graph , an integer , and lists of colors for every edge , and are tasked with determining if contains an assignment of colors to each edge such that no two edges incident to the same vertex are assigned the same color. We call such an assignment a proper edge coloring of with respect to the lists , or just a proper list edge coloring if the are clear from context.
For each , let be the subgraph of , where we define to be the set of all edges in the graph which have available as a color. Let denote the collection of -tuples of matchings such that
-
1.
for all , is a matching in , and
-
2.
for all vertices , appears as the endpoint of exactly of the matchings.
Observe that admits a proper list edge coloring if and only if contains a tuple of edge-disjoint matchings. Indeed, suppose admits a proper list edge coloring. For each , the set of edges assigned color must be a matching in . Moreover, every edge is assigned a color, so each vertex appears as an endpoint in of these matchings, hence the matchings form an edge-disjoint tuple in . Conversely, given an edge-disjoint tuple , assigning the edges in color gives a color to every edge in by condition 2 above, and thus yields a proper list edge coloring by condition 1.
The basic idea of our algorithmic template is to construct a polynomial which enumerates a certain subfamily . We then argue that contains a multilinear monomial satisfying certain matroid constraints if and only if contains a tuple of edge-disjoint matchings. Using similar reasoning to the previous paragraph, this then lets us solve List Edge Coloring by running monomial detection algorithms on .
4.1 Enumerating Matchings
In this section, we design a polynomial whose monomials enumerate tuples of matchings satisfying special conditions. To construct , we first define certain polynomial matrices.
For each edge we introduce a variable , and for each pair we introduce a variable . Let be the set of variables and be the set of variables.
For each , define to be the matrix with rows and columns indexed by with
For each , define to be a copy of . Let
be the union of these copies. Let be the symmetric matrix with rows and columns indexed by , defined by taking
(3) |
and identifying as the row and column set indexing .
Note that is , and over a field of characteristic two is skew-symmetric.
For each vertex , let be the set of copies of .
For each vertex , let be a linear matroid of rank over . These matroids will be specified later on, in our algorithms for Edge Coloring and List Edge Coloring. Let be the direct sum
(4) |
of these matroids.
By Proposition 2.7 the matroid is linear and has rank equal to the sum of the degrees of vertices in , which is . Let be a matrix representing .
Let be the collection of -tuples of matchings in such that
-
1.
all edges in have lists containing the color ,
-
2.
every vertex is incident to exactly edges across the matchings , and
-
3.
for each vertex , the set of copies taken over all such that appears as an endpoint of an edge in forms an independent set in .
We define the polynomial
(5) |
The next result shows that enumerates tuples of matchings from .
Lemma 4.1 (Enumerating Matchings).
We have
where each is a constant depending only on .
Proof.
Let be a subset of of size . By Proposition 2.8, we have
(6) |
By Equation 3 and Proposition 2.7, for every subset of size we have
Substituting the above equation into Equation 6 yields
Recall that given a graph , we let denote the set of perfect matchings of . By expanding out the definition of in the above equation, we see that
where is defined as a copy of , i.e., . By interchanging the order of the product and the summation, we see that the above is equivalent to summing over all collections , where is a matching in for each :
Let be a tuple such that the summand corresponding to in the above equation appears with nonzeo coefficient in . Since is nonzero if and only if , the tuple meets condition 1 in the definition of . By the definition of , we have if and only if is independent in for every vertex . Thus meets condition 3 the definition of . Since has rank for each vertex , has size , and the sum of degrees of vertices in is , we see that the only way for to be independent in for all is if in fact always. So satisfies condition 2 from the definition of .
Thus, the nonzero terms in the above sum correspond precisely to tuples , which proves the desired result.
4.2 Partition Sieving With Dominating Sets
In this section, we describe a generic way of going from a dominating set in to a certain partition matroid . After defining this , we will apply Theorem 3.5 to and to achieve a fast algorithm for detecting a multilinear monomial in .
Recall that a dominating set of a graph is a subset of vertices such that every vertex in is either in or adjacent to a vertex in .
Fix a dominating set of the input graph . Let . Let be the set of edges in with one endpoint in and one endpoint in , and be the corresponding set of variables. For each , let
be the set of edges incident to which connect to a vertex in . By definition, we have a partition
Let be the partition matroid over with respect to the above partition, and with capacities for each , where is viewed as a subgraph of . Since is a dominating set, every vertex is adjacent to some node in , and thus the capacities are all positive.
We now use this partition matroid to sieve over the polynomial , defined in Equation 5.
Lemma 4.2 (Accelerated Multilinear Monomial Detection).
There is a randomized algorithm that determines with high probability whether has a monomial divisible by , running in time and polynomial space.
Proof.
Let be the coefficient of
in . By definition, contains a monomial divisible by if and only if contains a monomial divisible by .
Let be the coefficient of in , viewed as a polynomial with variable set . Note that the total degree of is . Take a uniform random assignment over to the variables in . By Proposition 2.6, if is nonzero, then it remains nonzero with high probability after this evaluation, provided we take to be sufficiently large. For the rest of this proof, we work with polynomials after this random evaluation, so that is a polynomial in , is a polynomial in , and is a field element.
By the discussion in the first paragraph of this proof, it suffices to check whether contains the monomial in time and polynomial space.
We observe the following helpful property of .
Claim 4.3.
The polynomial is compatible with the matroid .
Proof.
Take an arbitrary monomial in . By definition, there is a corresponding monomial
in . By Lemma 4.1, for every , has degree exactly when restricted to the variables corresponding to edges containing . Consequently, for all , has degree exactly
when restricted to variables with . Since this holds for arbitrary monomials in , the polynomial is compatible with the partition matroid as claimed.
By Lemma 4.1, , viewed as a polynomial in , is homogeneous of degree . Hence , viewed as a polynomial in , is homogeneous of degree . By 4.3, is compatible with . Since is a partition matroid with positive capacities and parts in its underlying partition, by Theorem 3.5 we can determine whether contains a monomial in the variables such that is a basis of using
(7) |
evaluations of . Since is the unique basis of , this means we can test whether contains the monomial using evaluations of .
By Lemma 3.2, can be evaluated at any point by computing a linear combination of evaluations of . Combining this observation with Equation 7, we get that we can determine whether contains the monomial by computing a linear combination of
evaluations of . Each evaluation of takes polynomial time, so the desired result follows.
5 Edge Coloring Algorithm
In this section, we present our algorithm for Edge Coloring and prove Theorems 1.1 and 1.2. Recall that in this problem, we are given an input graph with vertices, edges, and maximum degree . To solve the problem, it suffices to determine whether admits a proper edge coloring using at most colors (i.e., whether ).
5.1 Removing Low Degree Vertices
Lemma 5.1 (Unit-Degree Deletion).
Let be a graph with a vertex of degree one. Let be the graph obtained by deleting from . Then if and only if .
Proof.
If admits a proper edge coloring with colors, then this coloring is also a proper edge coloring for with colors.
Conversely, suppose admits a proper edge coloring with colors. Since , there is a unique vertex in that is adjacent to. Then
because has maximum degree . So vertex has edges using at most distinct colors incident to it. Taking the coloring of and assigning edge a color not used by the edges incident to in recovers a proper edge coloring of with at most colors.
5.2 Instantiating the Template
Recall the definitions from Section 4. We view the Edge Coloring problem as a special case of List Edge Coloring, where and every edge is assigned the list . For each vertex , we set to be the uniform matroid over of rank . Then by Equation 4, is the partition matroid over defined by the partition
and the capacities for parts . By Proposition 2.11, we can construct the matrix representing in polynomial time. In this case, condition 3 is identical to condition 2 in the definition of the collection .
Lemma 5.2 (Characterization by Polynomials).
We have if and only if the polynomial has a monomial divisible by .
Proof.
Suppose contains a monomial divisible by . Since is homogeneous of degree in the variables, by Lemma 4.1 this means there exists a -tuple of edge-disjoint matchings , such that every edge of appears in some matching . Consequently, assigning the edges in color yields a proper edge coloring of using at most colors.
Conversely, suppose has a proper edge coloring with colors. Without loss of generality, let the set of colors used be . For each , let be the set of edges assigned color . Since this is a proper edge coloring, each is a matching in . Since every edge is assigned a unique color, each vertex has exactly edges across the matchings . Hence . Then by Lemma 4.1, contains a monomial divisible by .
See 1.4
Proof.
Let have vertex set and edge set , and let . By Lemma 5.2, we have if and only if has a monomial divisible by . By Lemma 4.2, we can determine whether such a monomial exists in
time and polynomial space, as desired.
Having established Lemma 1.4, we are ready to present out Edge Coloring algorithms.
See 1.1
Proof.
To solve Edge Coloring, it suffices to determine whether .
If has a vertex of degree one, delete it. Keep performing such deletions, until we are left with a graph containing no vertices of degree one. By repeated application of Lemma 5.1, the resulting graph admits a proper edge coloring using colors if and only if . If the resulting graph has constant size, we can determine this trivially. Each such deletion decreases the number of edges and vertices in the graph by one, and thus also decreases the value of . Consequently, to show the desired result, we may assume without loss of generality that has minimum degree two, and at least eight vertices.
Since has minimum degree two and nodes, by Lemma 2.2 the graph has a dominating set of size at most . By Lemma 2.4, we can then find a dominating set in with in time and polynomial space. Hence by Lemma 1.4 we can solve Edge Coloring on in time and polynomial space, as desired.
See 1.2
Proof.
Let be the input graph. By exhaustive search over all subsets of vertices, we can find a minimum-size dominating set of in time and polynomial space. Being -regular, has edges, so we find in time.
6 Conclusion
In this paper, we presented the first algorithms for Edge Coloring and List Edge Coloring which run in faster than time while using only polynomial space. Our algorithm is based off a partition sieving procedure that works over polynomials of degree and partition matroids with parts. We showed how to implement this sieving routine in time when the polynomial is in some sense compatible with the matroid, an improvement over the runtime that arose in previous techniques.
Overall, we solve Edge Coloring in faster than time by
- Step 1
-
designing a polynomial enumerating matchings meeting certain degree constraints,
- Step 2
-
finding a dominating set of size at most in the input graph, and
- Step 3
-
applying partition sieving with a partition matroid on parts.
Implementing step 2 above was possible for Edge Coloring because we could assume that the input graph had minimum degree two without loss of generality, and such graphs on nodes always have dominating sets of size at most . For List Edge Coloring we cannot make this assumption, but nonetheless implement a variant of step 2 by identifying a subgraph with minimum degree two and using extension matroids to handle unit-degree vertices for free. For all sufficiently large , it is also known that graphs on nodes with minimum degree three have dominating sets of size at most [34]. If degree-two vertices could somehow be “freely removed” from the input graph just like unit-degree vertices can, one could use this bound on the dominating set to solve coloring problems faster. More general connections between the minimum degree of a graph and the minimum size of its dominating set have also been studied [20].
Understanding the precise time complexity of Edge Coloring and its variants remains an interesting open research direction. Is Edge Coloring solvable in time, or at least in time? Establishing any nontrivial lower bounds for Edge Coloring would also be very interesting. The reductions from [21] imply that solving Edge Coloring in graphs with maximum degree requires time, assuming the Exponential Time Hypothesis (ETH), a standard hypothesis from fine-grained complexity. Does ETH similarly imply a or even a time lower bound for Edge Coloring in dense graphs (this question was previously raised in [25, Open problem 4.25])?
References
- [1] Vladimir I Arnautov. Estimation of the exterior stability number of a graph by means of the minimal degree of the vertices. Prikl. Mat. i Programmirovanie, 11(3-8):126, 1974.
- [2] Sepehr Assadi, Soheil Behnezhad, Sayan Bhattacharya, Martín Costa, Shay Solomon, and Tianyi Zhang. Vizing’s theorem in near-linear time, 2024. arXiv:2410.05240.
- [3] Alkida Balliu, Sebastian Brandt, Fabian Kuhn, and Dennis Olivetti. Distributed edge coloring in time polylogarithmic in . In Proceedings of the 2022 ACM Symposium on Principles of Distributed Computing, PODC ’22, pages 15–25. ACM, July 2022. doi:10.1145/3519270.3538440.
- [4] Lowell W Beineke. Derived graphs and digraphs. Beiträge zur graphentheorie, pages 17–33, 1968.
- [5] Anton Bernshteyn. A fast distributed algorithm for -edge-coloring. Journal of Combinatorial Theory, Series B, 152:319–352, January 2022. doi:10.1016/j.jctb.2021.10.004.
- [6] Sayan Bhattacharya, Martín Costa, Nadav Panski, and Shay Solomon. Arboricity-Dependent Algorithms for Edge Coloring. In Proceedings of the 19th Scandinavian Symposium and Workshops on Algorithm Theory (SWAT 2024), volume 294 of Leibniz International Proceedings in Informatics (LIPIcs), pages 12:1–12:15, Dagstuhl, Germany, 2024. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.SWAT.2024.12.
- [7] Andreas Björklund. Determinant sums for undirected Hamiltonicity. SIAM Journal on Computing, 43(1):280–299, 2014. doi:10.1137/110839229.
- [8] Andreas Björklund and Thore Husfeldt. Exact algorithms for exact satisfiability and number of perfect matchings. Algorithmica, 52(2):226–249, December 2007. doi:10.1007/s00453-007-9149-8.
- [9] Andreas Björklund, Thore Husfeldt, Petteri Kaski, and Mikko Koivisto. Trimmed moebius inversion and graphs of bounded degree. Theory of Computing Systems, 47(3):637–654, January 2009. doi:10.1007/s00224-009-9185-7.
- [10] Andreas Björklund, Thore Husfeldt, Petteri Kaski, and Mikko Koivisto. Narrow sieves for parameterized paths and packings. Journal of Computer and System Sciences, 87:119–139, 2017. doi:10.1016/j.jcss.2017.03.003.
- [11] Andreas Björklund, Thore Husfeldt, and Mikko Koivisto. Set partitioning via inclusion-exclusion. SIAM Journal on Computing, 39(2):546–563, 2009. doi:10.1137/070683933.
- [12] M Blank. An estimate of the external stability number of a graph without suspended vertices. Prikl. Mat. i Programmirovanie, 10:3–11, 1973.
- [13] Jesper Makholm Byskov. Enumerating maximal independent sets with applications to graph colouring. Operations Research Letters, 32(6):547–556, November 2004. doi:10.1016/j.orl.2004.03.002.
- [14] Aleksander B. G. Christiansen, Eva Rotenberg, and Juliette Vlieghe. Sparsity-Parameterised Dynamic Edge Colouring. In Proceedings of the 19th Scandinavian Symposium and Workshops on Algorithm Theory (SWAT 2024), volume 294 of Leibniz International Proceedings in Informatics (LIPIcs), pages 20:1–20:18, Dagstuhl, Germany, 2024. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.SWAT.2024.20.
- [15] Ran Duan, Haoqing He, and Tianyi Zhang. Dynamic Edge Coloring with Improved Approximation, pages 1937–1945. Society for Industrial and Applied Mathematics, January 2019. doi:10.1137/1.9781611975482.117.
- [16] Eduard Eiben, Tomohiro Koana, and Magnus Wahlström. Determinantal sieving. In David P. Woodruff, editor, Proceedings of the 2024 ACM-SIAM Symposium on Discrete Algorithms (SODA 2024), pages 377–423. SIAM, 2024. doi:10.1137/1.9781611977912.16.
- [17] David Eppstein. Small Maximal Independent Sets and Faster Exact Graph Coloring, pages 462–470. Springer Berlin Heidelberg, 2001. doi:10.1007/3-540-44634-6_42.
- [18] Alexander Golovnev, Alexander S. Kulikov, and Ivan Mihajlin. Families with infants: Speeding up algorithms for NP-hard problems using FFT. ACM Trans. Algorithms, 12(3):35:1–35:17, 2016. doi:10.1145/2847419.
- [19] David G. Harris. Distributed local approximation algorithms for maximum matching in graphs and hypergraphs. In Proceedings of the 60th Annual Symposium on Foundations of Computer Science (FOCS 2019), pages 700–724, 2019. doi:10.1109/FOCS.2019.00048.
- [20] Michael A. Henning. Bounds on domination parameters in graphs: a brief survey. Discuss. Math. Graph Theory, 42(3):665–708, 2022. doi:10.7151/DMGT.2454.
- [21] Ian Holyer. The NP-completeness of edge-coloring. SIAM Journal on Computing, 10(4):718–720, 1981. doi:10.1137/0210055.
- [22] Masao Ishikawa and Masato Wakayama. Minor summation formula of pfaffians. Linear and Multilinear algebra, 39(3):285–305, 1995.
- [23] Tomohiro Koana and Magnus Wahlström. Faster algorithms on linear delta-matroids. arXiv preprint arXiv:2402.11596, 2024. doi:10.48550/arXiv.2402.11596.
- [24] Eugene L. Lawler. A note on the complexity of the chromatic number problem. Inf. Process. Lett., 5(3):66–67, 1976. doi:10.1016/0020-0190(76)90065-X.
- [25] Moshe Lewenstein, Seth Pettie, and Virginia Vassilevska Williams. Structure and hardness in P (dagstuhl seminar 16451). Dagstuhl Reports, 6(11):1–34, 2016. doi:10.4230/DAGREP.6.11.1.
- [26] Daniel Lokshtanov, Pranabendu Misra, Fahad Panolan, and Saket Saurabh. Deterministic truncation of linear matroids. ACM Transactions on Algorithms, 14(2):14:1–14:20, 2018. doi:10.1145/3170444.
- [27] Carsten Lund and Mihalis Yannakakis. On the hardness of approximating minimization problems. Journal of the ACM, 41(5):960–981, September 1994. doi:10.1145/185675.306789.
- [28] Dániel Marx. A parameterized view on matroid optimization problems. Theoretical Computer Science, 410(44):4471–4479, 2009. doi:10.1016/j.tcs.2009.07.027.
- [29] William McCuaig and Bruce Shepherd. Domination in graphs with minimum degree two. Journal of Graph Theory, 13(6):749–762, 1989. doi:10.1002/JGT.3190130610.
- [30] Rajeev Motwani and Prabhakar Raghavan. Randomized Algorithms. Cambridge University Press, August 1995. doi:10.1017/cbo9780511814075.
- [31] Kazuo Murota. Matrices and matroids for systems analysis, volume 20. Springer Science & Business Media, 1999.
- [32] Oystein Ore. Theory of graphs. In Colloquium Publications. American Mathematical Society, 1962.
- [33] Charles Payan. Sur le nombre d’absorption d’un graphe simple, 1975.
- [34] Bruce A. Reed. Paths, stars and the number three. Comb. Probab. Comput., 5:277–295, 1996. doi:10.1017/S0963548300002042.
- [35] Magnus Wahlström. Abusing the Tutte matrix: An algebraic instance compression for the K-set-cycle problem. In Proceedings of the 30th International Symposium on Theoretical Aspects of Computer Science (STACS 2013), volume 20 of LIPIcs, pages 341–352. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2013. doi:10.4230/LIPICS.STACS.2013.341.
- [36] Ryan Williams. Finding paths of length in time. Information Processing Letter, 109(6):315–318, 2009. doi:10.1016/J.IPL.2008.11.004.