Chromatic polynomial

Short description: Function in algebraic graph theory

All non-isomorphic graphs on 3 vertices and their chromatic polynomials, clockwise from the top. The independent 3-set: k³. An edge and a single vertex: k²(k – 1). The 3-path: k(k – 1)². The 3-clique: k(k – 1)(k – 2).

The chromatic polynomial is a graph polynomial studied in algebraic graph theory, a branch of mathematics. It counts the number of graph colorings as a function of the number of colors and was originally defined by George David Birkhoff to study the four color problem. It was generalised to the Tutte polynomial by Hassler Whitney and W. T. Tutte, linking it to the Potts model of statistical physics.

History

George David Birkhoff introduced the chromatic polynomial in 1912, defining it only for planar graphs, in an attempt to prove the four color theorem. If [math]\displaystyle{ P(G, k) }[/math] denotes the number of proper colorings of G with k colors then one could establish the four color theorem by showing [math]\displaystyle{ P(G, 4)\gt 0 }[/math] for all planar graphs G. In this way he hoped to apply the powerful tools of analysis and algebra for studying the roots of polynomials to the combinatorial coloring problem.

Hassler Whitney generalised Birkhoff’s polynomial from the planar case to general graphs in 1932. In 1968, Ronald C. Read asked which polynomials are the chromatic polynomials of some graph, a question that remains open, and introduced the concept of chromatically equivalent graphs. Today, chromatic polynomials are one of the central objects of algebraic graph theory.

Definition

All proper vertex colorings of vertex graphs with 3 vertices using k colors for [math]\displaystyle{ k=0,1,2,3 }[/math]. The chromatic polynomial of each graph interpolates through the number of proper colorings.

For a graph G, [math]\displaystyle{ P(G,k) }[/math] counts the number of its (proper) vertex k-colorings. Other commonly used notations include [math]\displaystyle{ P_G(k) }[/math], [math]\displaystyle{ \chi_G(k) }[/math], or [math]\displaystyle{ \pi_G(k) }[/math]. There is a unique polynomial [math]\displaystyle{ P(G,x) }[/math] which evaluated at any integer k ≥ 0 coincides with [math]\displaystyle{ P(G,k) }[/math]; it is called the chromatic polynomial of G.

For example, to color the path graph [math]\displaystyle{ P_3 }[/math] on 3 vertices with k colors, one may choose any of the k colors for the first vertex, any of the [math]\displaystyle{ k - 1 }[/math] remaining colors for the second vertex, and lastly for the third vertex, any of the [math]\displaystyle{ k - 1 }[/math] colors that are different from the second vertex's choice. Therefore, [math]\displaystyle{ P(P_3,k) = k \cdot (k-1) \cdot (k-1) }[/math] is the number of k-colorings of [math]\displaystyle{ P_3 }[/math]. For a variable x (not necessarily integer), we thus have [math]\displaystyle{ P(P_3,x)=x(x-1)^2=x^3-2x^2+x }[/math]. (Colorings which differ only by permuting colors or by automorphisms of G are still counted as different.)

Deletion–contraction

The fact that the number of k-colorings is a polynomial in k follows from a recurrence relation called the deletion–contraction recurrence or Fundamental Reduction Theorem. It is based on edge contraction: for a pair of vertices [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math] the graph [math]\displaystyle{ G/uv }[/math] is obtained by merging the two vertices and removing any edges between them. If [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math] are adjacent in G, let [math]\displaystyle{ G-uv }[/math] denote the graph obtained by removing the edge [math]\displaystyle{ uv }[/math]. Then the numbers of k-colorings of these graphs satisfy:

[math]\displaystyle{ P(G,k)=P(G-uv, k)- P(G/uv,k) }[/math]

Equivalently, if [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math] are not adjacent in G and [math]\displaystyle{ G+uv }[/math] is the graph with the edge [math]\displaystyle{ uv }[/math] added, then

[math]\displaystyle{ P(G,k)= P(G+uv, k) + P(G/uv,k) }[/math]

This follows from the observation that every k-coloring of G either gives different colors to [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math], or the same colors. In the first case this gives a (proper) k-coloring of [math]\displaystyle{ G+uv }[/math], while in the second case it gives a coloring of [math]\displaystyle{ G/uv }[/math]. Conversely, every k-coloring of G can be uniquely obtained from a k-coloring of [math]\displaystyle{ G+uv }[/math] or [math]\displaystyle{ G/uv }[/math] (if [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math] are not adjacent in G).

The chromatic polynomial can hence be recursively defined as

[math]\displaystyle{ P(G,x)=x^n }[/math] for the edgeless graph on n vertices, and

[math]\displaystyle{ P(G,x)=P(G-uv, x)- P(G/uv,x) }[/math] for a graph G with an edge [math]\displaystyle{ uv }[/math] (arbitrarily chosen).

Since the number of k-colorings of the edgeless graph is indeed [math]\displaystyle{ k^n }[/math], it follows by induction on the number of edges that for all G, the polynomial [math]\displaystyle{ P(G,x) }[/math] coincides with the number of k-colorings at every integer point x = k. In particular, the chromatic polynomial is the unique interpolating polynomial of degree at most n through the points

[math]\displaystyle{ \left \{ (0, P(G, 0)), (1, P(G, 1)), \ldots, (n, P(G, n)) \right \}. }[/math]

Tutte’s curiosity about which other graph invariants satisfied such recurrences led him to discover a bivariate generalization of the chromatic polynomial, the Tutte polynomial [math]\displaystyle{ T_G(x,y) }[/math].

Examples

Chromatic polynomials for certain graphs

Triangle [math]\displaystyle{ K_3 }[/math]	[math]\displaystyle{ x(x-1)(x-2) }[/math]
Complete graph [math]\displaystyle{ K_n }[/math]	[math]\displaystyle{ x(x-1)(x-2)\cdots(x-(n-1)) }[/math]
Edgeless graph [math]\displaystyle{ \overline K_n }[/math]	[math]\displaystyle{ x^n }[/math]
Path graph [math]\displaystyle{ P_n }[/math]	[math]\displaystyle{ x(x-1)^{n-1} }[/math]
Any tree on n vertices	[math]\displaystyle{ x(x-1)^{n-1} }[/math]
Cycle [math]\displaystyle{ C_n }[/math]	[math]\displaystyle{ (x-1)^n+(-1)^n(x-1) }[/math]
Petersen graph	[math]\displaystyle{ x(x-1)(x-2) \left (x^7-12x^6+67x^5-230x^4+529x^3-814x^2+775x-352 \right) }[/math]

Properties

For fixed G on n vertices, the chromatic polynomial [math]\displaystyle{ P(G, x) }[/math] is a monic polynomial of degree exactly n, with integer coefficients.

The chromatic polynomial includes at least as much information about the colorability of G as does the chromatic number. Indeed, the chromatic number is the smallest positive integer that is not a zero of the chromatic polynomial,

[math]\displaystyle{ \chi (G)=\min\{ k\in\mathbb{N} : P(G, k) \gt 0 \}. }[/math]

The polynomial evaluated at [math]\displaystyle{ -1 }[/math], that is [math]\displaystyle{ P(G,-1) }[/math], yields [math]\displaystyle{ (-1)^{|V(G)|} }[/math] times the number of acyclic orientations of G.

The derivative evaluated at 1, [math]\displaystyle{ P'(G, 1) }[/math] equals the chromatic invariant [math]\displaystyle{ \theta(G) }[/math] up to sign.

If G has n vertices and c components [math]\displaystyle{ G_1, \ldots, G_c }[/math], then

• The coefficients of [math]\displaystyle{ x^0, \ldots, x^{c-1} }[/math] are zeros.
• The coefficients of [math]\displaystyle{ x^c, \ldots, x^n }[/math] are all non-zero and alternate in signs.
• The coefficient of [math]\displaystyle{ x^n }[/math] is 1 (the polynomial is monic).
• The coefficient of [math]\displaystyle{ x^{n-1} }[/math] is [math]\displaystyle{ -|E(G)| }[/math].

We prove this via induction on the number of edges on a simple graph G with [math]\displaystyle{ n }[/math] vertices and [math]\displaystyle{ k }[/math] edges. When [math]\displaystyle{ k = 0 }[/math], G is an empty graph. Hence per definition [math]\displaystyle{ P(G, x)= x^n }[/math]. So the coefficient of [math]\displaystyle{ x^{n-1} }[/math] is [math]\displaystyle{ 0 }[/math], which implies the statement is true for an empty graph. When [math]\displaystyle{ k = 1 }[/math], as in G has just a single edge, [math]\displaystyle{ P(G, x) = x^n - x^{n-1} }[/math]. Thus coefficient of [math]\displaystyle{ x^{n-1} }[/math] is [math]\displaystyle{ -1 = -|E(G)| }[/math]. So the statement holds for k = 1. Using strong induction assume the statement is true for [math]\displaystyle{ k = 0,1,2,\ldots,(k-1) }[/math]. Let G have [math]\displaystyle{ k }[/math] edges. By the contraction-deletion principle,

[math]\displaystyle{ P(G, x) = P(G-e, x) - P(G/e, x) }[/math],

Let [math]\displaystyle{ P(G-e, x) = x^n - a_{n-1}x^{n-1} + a_{n-2}x^{n-2}-\cdots }[/math], and [math]\displaystyle{ P(G/e, x) = x^{n-1} - b_{n-2} x^{n-2} + b_{n-3}x^{n-3}-\cdots }[/math].
Hence [math]\displaystyle{ P(G, x) = x^n - (a_{n-1} +1)x^{n-1} +\cdots }[/math].
Since [math]\displaystyle{ G-e }[/math] is obtained from G by removal of just one edge e, [math]\displaystyle{ a_{n-1} = k - 1 }[/math], so [math]\displaystyle{ a_{n-1} + 1 = k }[/math] and thus the statement is true for k.

• The coefficient of [math]\displaystyle{ x^1 }[/math] is [math]\displaystyle{ (-1)^{n-1} }[/math] times the number of acyclic orientations that have a unique sink, at a specified, arbitrarily chosen vertex.
• The absolute values of coefficients of every chromatic polynomial form a log-concave sequence.
• [math]\displaystyle{ \scriptstyle P(G, x) = P(G_1, x)P(G_2,x) \cdots P(G_c,x) }[/math]

The last property is generalized by the fact that if G is a k-clique-sum of [math]\displaystyle{ G_1 }[/math] and [math]\displaystyle{ G_2 }[/math] (i.e., a graph obtained by gluing the two at a clique on k vertices), then

[math]\displaystyle{ P(G, x) = \frac{P(G_1,x)P(G_2,x)}{x(x-1)\cdots(x-k+1)}. }[/math]

A graph G with n vertices is a tree if and only if

[math]\displaystyle{ P(G, x) = x(x-1)^{n-1}. }[/math]

Chromatic equivalence

The three graphs with a chromatic polynomial equal to [math]\displaystyle{ (x-2)(x-1)^3x }[/math].

Two graphs are said to be chromatically equivalent if they have the same chromatic polynomial. Isomorphic graphs have the same chromatic polynomial, but non-isomorphic graphs can be chromatically equivalent. For example, all trees on n vertices have the same chromatic polynomial. In particular, [math]\displaystyle{ (x-1)^3x }[/math] is the chromatic polynomial of both the claw graph and the path graph on 4 vertices.

A graph is chromatically unique if it is determined by its chromatic polynomial, up to isomorphism. In other words, G is chromatically unique, then [math]\displaystyle{ P(G, x) = P(H, x) }[/math] would imply that G and H are isomorphic. All cycle graphs are chromatically unique.

Chromatic roots

A root (or zero) of a chromatic polynomial, called a “chromatic root”, is a value x where [math]\displaystyle{ P(G, x)=0 }[/math]. Chromatic roots have been very well studied, in fact, Birkhoff’s original motivation for defining the chromatic polynomial was to show that for planar graphs, [math]\displaystyle{ P(G, x)\gt 0 }[/math] for x ≥ 4. This would have established the four color theorem.

No graph can be 0-colored, so 0 is always a chromatic root. Only edgeless graphs can be 1-colored, so 1 is a chromatic root of every graph with at least one edge. On the other hand, except for these two points, no graph can have a chromatic root at a real number smaller than or equal to 32/27. A result of Tutte connects the golden ratio [math]\displaystyle{ \phi }[/math] with the study of chromatic roots, showing that chromatic roots exist very close to [math]\displaystyle{ \phi^2 }[/math]: If [math]\displaystyle{ G_n }[/math] is a planar triangulation of a sphere then

[math]\displaystyle{ P(G_n,\phi^2) \leq \phi^{5-n}. }[/math]

While the real line thus has large parts that contain no chromatic roots for any graph, every point in the complex plane is arbitrarily close to a chromatic root in the sense that there exists an infinite family of graphs whose chromatic roots are dense in the complex plane.

Colorings using all colors

For a graph G on n vertices, let [math]\displaystyle{ e_k }[/math] denote the number of colorings using exactly k colors up to renaming colors (so colorings that can be obtained from one another by permuting colors are counted as one; colorings obtained by automorphisms of G are still counted separately). In other words, [math]\displaystyle{ e_k }[/math] counts the number of partitions of the vertex set into k (non-empty) independent sets. Then [math]\displaystyle{ k! \cdot e_k }[/math] counts the number of colorings using exactly k colors (with distinguishable colors). For an integer x, all x-colorings of G can be uniquely obtained by choosing an integer k ≤ x, choosing k colors to be used out of x available, and a coloring using exactly those k (distinguishable) colors. Therefore:

[math]\displaystyle{ P(G,x) = \sum_{k=0}^x \binom{x}{k} k! \cdot e_k = \sum_{k=0}^x (x)_k \cdot e_k }[/math],

where [math]\displaystyle{ (x)_k = x(x-1)(x-2)\cdots(x-k+1) }[/math] denotes the falling factorial. Thus the numbers [math]\displaystyle{ e_k }[/math] are the coefficients of the polynomial [math]\displaystyle{ P(G,x) }[/math] in the basis [math]\displaystyle{ 1,(x)_1,(x)_2,(x)_3,\ldots }[/math] of falling factorials.

Let [math]\displaystyle{ a_k }[/math] be the k-th coefficient of [math]\displaystyle{ P(G,x) }[/math] in the standard basis [math]\displaystyle{ 1,x,x^2,x^3,\ldots }[/math], that is:

[math]\displaystyle{ P(G,x) = \sum_{k=0}^n a_k x^k }[/math]

Stirling numbers give a change of basis between the standard basis and the basis of falling factorials. This implies:

[math]\displaystyle{ a_k = \sum_{j=0}^{n} (-1)^{j-k} \begin{bmatrix}j\\k\end{bmatrix} e_j }[/math]   and   [math]\displaystyle{ e_k = \sum_{j=0}^n \begin{Bmatrix}j\\k\end{Bmatrix} a_j }[/math].

Categorification

The chromatic polynomial is categorified by a homology theory closely related to Khovanov homology.

Algorithms

Chromatic polynomial
Input	Graph G with n vertices.
Output	Coefficients of [math]\displaystyle{ P(G, x) }[/math]
Running time	[math]\displaystyle{ O(2^nn^r) }[/math] for some constant [math]\displaystyle{ r }[/math]
Complexity	P-hard
Reduction from	#3SAT

#k-colorings
Input	Graph G with n vertices.
Output	[math]\displaystyle{ P(G, k) }[/math]
Running time	In P for [math]\displaystyle{ k=0,1,2 }[/math]. [math]\displaystyle{ O(1.6262^n) }[/math] for [math]\displaystyle{ k=3 }[/math]. Otherwise [math]\displaystyle{ O(2^nn^r) }[/math] for some constant [math]\displaystyle{ r }[/math]
Complexity	P-hard unless [math]\displaystyle{ k=0,1,2 }[/math]
Approximability	No FPRAS for [math]\displaystyle{ k\gt 2 }[/math]

Computational problems associated with the chromatic polynomial include

• finding the chromatic polynomial [math]\displaystyle{ P(G, x) }[/math] of a given graph G;
• evaluating [math]\displaystyle{ P(G, x) }[/math] at a fixed x for given G.

The first problem is more general because if we knew the coefficients of [math]\displaystyle{ P(G, x) }[/math] we could evaluate it at any point in polynomial time because the degree is n. The difficulty of the second type of problem depends strongly on the value of x and has been intensively studied in computational complexity. When x is a natural number, this problem is normally viewed as computing the number of x-colorings of a given graph. For example, this includes the problem #3-coloring of counting the number of 3-colorings, a canonical problem in the study of complexity of counting, complete for the counting class #P.

Efficient algorithms

For some basic graph classes, closed formulas for the chromatic polynomial are known. For instance this is true for trees and cliques, as listed in the table above.

Polynomial time algorithms are known for computing the chromatic polynomial for wider classes of graphs, including chordal graphs and graphs of bounded clique-width. The latter class includes cographs and graphs of bounded tree-width, such as outerplanar graphs.

Deletion–contraction

The deletion-contraction recurrence gives a way of computing the chromatic polynomial, called the deletion–contraction algorithm. In the first form (with a minus), the recurrence terminates in a collection of empty graphs. In the second form (with a plus), it terminates in a collection of complete graphs. This forms the basis of many algorithms for graph coloring. The ChromaticPolynomial function in the Combinatorica package of the computer algebra system Mathematica uses the second recurrence if the graph is dense, and the first recurrence if the graph is sparse. The worst case running time of either formula satisfies the same recurrence relation as the Fibonacci numbers, so in the worst case, the algorithm runs in time within a polynomial factor of

[math]\displaystyle{ \phi^{n+m}=\left (\frac{1+\sqrt{5}}{2} \right)^{n+m}\in O\left(1.62^{n+m}\right), }[/math]

on a graph with n vertices and m edges. The analysis can be improved to within a polynomial factor of the number [math]\displaystyle{ t(G) }[/math] of spanning trees of the input graph. In practice, branch and bound strategies and graph isomorphism rejection are employed to avoid some recursive calls, the running time depends on the heuristic used to pick the vertex pair.

Cube method

There is a natural geometric perspective on graph colorings by observing that, as an assignment of natural numbers to each vertex, a graph coloring is a vector in the integer lattice. Since two vertices [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math] being given the same color is equivalent to the [math]\displaystyle{ i }[/math]’th and [math]\displaystyle{ j }[/math]’th coordinate in the coloring vector being equal, each edge can be associated with a hyperplane of the form [math]\displaystyle{ \{x\in R^d:x_i=x_j\} }[/math]. The collection of such hyperplanes for a given graph is called its graphic arrangement. The proper colorings of a graph are those lattice points which avoid forbidden hyperplanes. Restricting to a set of [math]\displaystyle{ k }[/math] colors, the lattice points are contained in the cube [math]\displaystyle{ [0,k]^n }[/math]. In this context the chromatic polynomial counts the number of lattice points in the [math]\displaystyle{ [0,k] }[/math]-cube that avoid the graphic arrangement.

Computational complexity

The problem of computing the number of 3-colorings of a given graph is a canonical example of a #P-complete problem, so the problem of computing the coefficients of the chromatic polynomial is #P-hard. Similarly, evaluating [math]\displaystyle{ P(G, 3) }[/math] for given G is #P-complete. On the other hand, for [math]\displaystyle{ k=0,1,2 }[/math] it is easy to compute [math]\displaystyle{ P(G, k) }[/math], so the corresponding problems are polynomial-time computable. For integers [math]\displaystyle{ k\gt 3 }[/math] the problem is #P-hard, which is established similar to the case [math]\displaystyle{ k=3 }[/math]. In fact, it is known that [math]\displaystyle{ P(G, x) }[/math] is #P-hard for all x (including negative integers and even all complex numbers) except for the three “easy points”. Thus, from the perspective of #P-hardness, the complexity of computing the chromatic polynomial is completely understood.

In the expansion

[math]\displaystyle{ P(G, x)= a_1 x + a_2x^2+\cdots +a_nx^n, }[/math]

the coefficient [math]\displaystyle{ a_n }[/math] is always equal to 1, and several other properties of the coefficients are known. This raises the question if some of the coefficients are easy to compute. However the computational problem of computing a_r for a fixed r ≥ 1 and a given graph G is #P-hard, even for bipartite planar graphs.

No approximation algorithms for computing [math]\displaystyle{ P(G, x) }[/math] are known for any x except for the three easy points. At the integer points [math]\displaystyle{ k=3,4,\ldots }[/math], the corresponding decision problem of deciding if a given graph can be k-colored is NP-hard. Such problems cannot be approximated to any multiplicative factor by a bounded-error probabilistic algorithm unless NP = RP, because any multiplicative approximation would distinguish the values 0 and 1, effectively solving the decision version in bounded-error probabilistic polynomial time. In particular, under the same assumption, this rules out the possibility of a fully polynomial time randomised approximation scheme (FPRAS). For other points, more complicated arguments are needed, and the question is the focus of active research. (As of 2008), it is known that there is no FPRAS for computing [math]\displaystyle{ P(G, x) }[/math] for any x > 2, unless NP = RP holds.

References

• Biggs, N. (1993), Algebraic Graph Theory, Cambridge University Press, ISBN 978-0-521-45897-9 
• Chao, C.-Y.; Whitehead, E. G. (1978), "On chromatic equivalence of graphs", Theory and Applications of Graphs, Lecture Notes in Mathematics, 642, Springer, pp. 121–131, ISBN 978-3-540-08666-6 
• Dong, F. M.; Koh, K. M.; Teo, K. L. (2005), Chromatic polynomials and chromaticity of graphs, World Scientific Publishing Company, ISBN 978-981-256-317-0 
• Dehmer, Matthias, ed. (2011), "Graph polynomials and their applications I: The Tutte polynomial" (in en-gb), Structural Analysis of Complex Networks, doi:10.1007/978-0-8176-4789-6_9, ISBN 978-0-8176-4789-6 
• Giménez, O.; Hliněný, P.; Noy, M. (2005), "Computing the Tutte polynomial on graphs of bounded clique-width", Proc. 31st Int. Worksh. Graph-Theoretic Concepts in Computer Science (WG 2005), Lecture Notes in Computer Science, 3787, Springer-Verlag, pp. 59–68, doi:10.1007/11604686_6, ISBN 978-3-540-31000-6 
• "Inapproximability of the Tutte polynomial", Information and Computation 206 (7): 908, 2008, doi:10.1016/j.ic.2008.04.003 
• Helme-Guizon, Laure; Rong, Yongwu (2005), "A categorification of the chromatic polynomial", Algebraic & Geometric Topology 5 (4): 1365–1388, doi:10.2140/agt.2005.5.1365 
• Huh, June (2012), "Milnor numbers of projective hypersurfaces and the chromatic polynomial of graphs", Journal of the American Mathematical Society 25 (3): 907–927, doi:10.1090/S0894-0347-2012-00731-0 
• Jackson, B. (1993), "A Zero-Free Interval for Chromatic Polynomials of Graphs", Combinatorics, Probability and Computing 2 (3): 325–336, doi:10.1017/S0963548300000705 
• Jaeger, F.; Vertigan, D. L.; Welsh, D. J. A. (1990), "On the computational complexity of the Jones and Tutte polynomials", Mathematical Proceedings of the Cambridge Philosophical Society 108 (1): 35–53, doi:10.1017/S0305004100068936, Bibcode: 1990MPCPS.108...35J 
• "Hard enumeration problems in geometry and combinatorics", SIAM J. Algebr. Discrete Methods 7 (2): 331–335, 1986, doi:10.1137/0607036 
• Makowsky, J. A.; Rotics, U.; Averbouch, I.; Godlin, B. (2006), "Computing graph polynomials on graphs of bounded clique-width", Proc. 32nd Int. Worksh. Graph-Theoretic Concepts in Computer Science (WG 2006), Lecture Notes in Computer Science, 4271, Springer-Verlag, pp. 191–204, doi:10.1007/11917496_18, ISBN 978-3-540-48381-6 
• Naor, J.; Naor, M.; Schaffer, A. (1987), "Fast parallel algorithms for chordal graphs", Proc. 19th ACM Symp. Theory of Computing (STOC '87), pp. 355–364, doi:10.1145/28395.28433, ISBN 978-0897912211 .
• Oxley, J. G.; Welsh, D. J. A. (2002), "Chromatic, flow and reliability polynomials: The complexity of their coefficients.", Combinatorics, Probability and Computing 11 (4): 403–426, doi:10.1017/S0963548302005175 
• Pemmaraju, S.; Skiena, S. (2003), Computational Discrete Mathematics: Combinatorics and Graph Theory with Mathematica, Cambridge University Press, section 7.4.2, ISBN 978-0-521-80686-2 
• "An introduction to chromatic polynomials", Journal of Combinatorial Theory 4: 52–71, 1968, doi:10.1016/S0021-9800(68)80087-0, http://dml.cz/bitstream/handle/10338.dmlcz/128421/CzechMathJ_43-1993-3_10.pdf 
• Sekine, Kyoko; Imai, Hiroshi; Tani, Seiichiro (1995), "Computing the Tutte polynomial of a graph of moderate size", in Staples, John; Eades, Peter; Katoh, Naoki et al., Algorithms and Computation, 6th International Symposium, ISAAC '95, Cairns, Australia, December 4–6, 1995, Proceedings, Lecture Notes in Computer Science, 1004, Springer, pp. 224–233, doi:10.1007/BFb0015427 
• "Chromatic Roots are Dense in the Whole Complex Plane", Combinatorics, Probability and Computing 13 (2): 221–261, 2004, doi:10.1017/S0963548303006023 
• Stanley, R. P. (1973), "Acyclic orientations of graphs", Discrete Math. 5 (2): 171–178, doi:10.1016/0012-365X(73)90108-8, http://math.mit.edu/~rstan/pubs/pubfiles/18.pdf 
• Voloshin, Vitaly I. (2002), Coloring Mixed Hypergraphs: Theory, Algorithms and Applications., American Mathematical Society, ISBN 978-0-8218-2812-0 
• Algorithms and Complexity, Prentice–Hall, 1986, ISBN 978-0-13-021973-2 

External links

• Weisstein, Eric W.. "Chromatic polynomial". http://mathworld.wolfram.com/ChromaticPolynomial.html. 
• PlanetMath Chromatic polynomial
• Code for computing Tutte, Chromatic and Flow Polynomials by Gary Haggard, David J. Pearce and Gordon Royle: [1]