Color-coding

In computer science and graph theory, the term color-coding refers to an algorithmic technique which is useful in the discovery of network motifs. For example, it can be used to detect a simple path of length k in a given graph. The traditional color-coding algorithm is probabilistic, but it can be derandomized without much overhead in the running time.

Color-coding also applies to the detection of cycles of a given length, and more generally it applies to the subgraph isomorphism problem (an NP-complete problem), where it yields polynomial time algorithms when the subgraph pattern that it is trying to detect has bounded treewidth.

The color-coding method was proposed and analyzed in 1994 by Noga Alon, Raphael Yuster, and Uri Zwick.

Results

The following results can be obtained through the method of color-coding:

• For every fixed constant k, if a graph G = (V, E) contains a simple cycle of size k, then such a cycle can be found in:
  • [math]\displaystyle{ O(|V|^\omega) }[/math] expected time, or
  • [math]\displaystyle{ O(|V|^\omega \log |V|) }[/math] worst-case time, where ω is the exponent of matrix multiplication.
• For every fixed constant k, and every graph G = (V, E) that is in any nontrivial minor-closed graph family (e.g., a planar graph), if G contains a simple cycle of size k, then such cycle can be found in:
  • O(V) expected time, or
  • O(V log V) worst-case time.
• If a graph G = (V, E) contains a subgraph isomorphic to a bounded treewidth graph which has O(log V) vertices, then such a subgraph can be found in polynomial time.

The method

To solve the problem of finding a subgraph [math]\displaystyle{ H = (V_H, E_H) }[/math] in a given graph G = (V, E), where H can be a path, a cycle, or any bounded treewidth graph where [math]\displaystyle{ |V_H| = O(\log |V|) }[/math], the method of color-coding begins by randomly coloring each vertex of G with [math]\displaystyle{ k = |V_H| }[/math] colors, and then tries to find a colorful copy of H in colored G. Here, a graph is colorful if every vertex in it is colored with a distinct color. This method works by repeating (1) random coloring a graph and (2) finding colorful copy of the target subgraph, and eventually the target subgraph can be found if the process is repeated a sufficient number of times.

Suppose a copy of H in G becomes colorful with some non-zero probability p. It immediately follows that if the random coloring is repeated 1/p times, then this copy is expected to become colorful once. Note that though p is small, it is shown that if [math]\displaystyle{ |V_H| = O(\log |V|) }[/math], p is only polynomially small. Suppose again there exists an algorithm such that, given a graph G and a coloring which maps each vertex of G to one of the k colors, it finds a copy of colorful H, if one exists, within some runtime O(r). Then the expected time to find a copy of H in G, if one exists, is [math]\displaystyle{ O(\tfrac{r}{p}) }[/math].

Sometimes it is also desirable to use a more restricted version of colorfulness. For example, in the context of finding cycles in planar graphs, it is possible to develop an algorithm that finds well-colored cycles. Here, a cycle is well-colored if its vertices are colored by consecutive colors.

Example

An example would be finding a simple cycle of length k in graph G = (V, E).

By applying random coloring method, each simple cycle has a probability of [math]\displaystyle{ k!/k^k \gt e^{-k} }[/math] to become colorful, since there are [math]\displaystyle{ k^k }[/math] ways of coloring the k vertices on the cycle, among which there are [math]\displaystyle{ k! }[/math] colorful occurrences. Then an algorithm (described next) can be used to find colorful cycles in the randomly colored graph G in time [math]\displaystyle{ O(V^\omega) }[/math], where [math]\displaystyle{ \omega }[/math] is the matrix multiplication constant. Therefore, it takes [math]\displaystyle{ e^k\cdot O(V^\omega) }[/math] overall time to find a simple cycle of length k in G.

The colorful cycle-finding algorithm works by first finding all pairs of vertices in V that are connected by a simple path of length k − 1, and then checking whether the two vertices in each pair are connected. Given a coloring function c : V → {1, ..., k} to color graph G, enumerate all partitions of the color set {1, ..., k} into two subsets C₁, C₂ of size [math]\displaystyle{ k/2 }[/math] each. Note that V can be divided into V₁ and V₂ accordingly, and let G₁ and G₂ denote the subgraphs induced by V₁ and V₂ respectively. Then, recursively find colorful paths of length [math]\displaystyle{ k/2 - 1 }[/math] in each of G₁ and G₂. Suppose the boolean matrix A₁ and A₂ represent the connectivity of each pair of vertices in G₁ and G₂ by a colorful path, respectively, and let B be the matrix describing the adjacency relations between vertices of V₁ and those of V₂, the boolean product [math]\displaystyle{ A_1BA_2 }[/math] gives all pairs of vertices in V that are connected by a colorful path of length k − 1. Thus, the recursive relation of matrix multiplications is [math]\displaystyle{ t(k) \le 2^k\cdot t(k/2) }[/math], which yields a runtime of [math]\displaystyle{ 2^{O(k)}\cdot V^\omega }[/math]. Although this algorithm finds only the end points of the colorful path, another algorithm by Alon and Naor that finds colorful paths themselves can be incorporated into it.

Derandomization

The derandomization of color-coding involves enumerating possible colorings of a graph G, such that the randomness of coloring G is no longer required. For the target subgraph H in G to be discoverable, the enumeration has to include at least one instance where the H is colorful. To achieve this, enumerating a k-perfect family F of hash functions from {1, ..., |V|} to {1, ..., k} is sufficient. By definition, F is k-perfect if for every subset S of {1, ..., |V|} where [math]\displaystyle{ |S| = k }[/math], there exists a hash function h in F such that h : S → {1, ..., k} is perfect. In other words, there must exist a hash function in F that colors any given k vertices with k distinct colors.

There are several approaches to construct such a k-perfect hash family:

1. The best explicit construction is by Moni Naor, Leonard J. Schulman, and Aravind Srinivasan, where a family of size [math]\displaystyle{ e^k k^{O(\log k)} \log |V| }[/math] can be obtained. This construction does not require the target subgraph to exist in the original subgraph finding problem.
2. Another explicit construction by Jeanette P. Schmidt and Alan Siegel yields a family of size [math]\displaystyle{ 2^{O(k)}\log^2 |V| }[/math].
3. Another construction that appears in the original paper of Noga Alon et al. can be obtained by first building a k-perfect family that maps {1, ..., |V|} to {1, ..., k²}, followed by building another k-perfect family that maps {1, ..., k²} to {1, ..., k}. In the first step, it is possible to construct such a family with 2nlog k random bits that are almost 2log k-wise independent, and the sample space needed for generating those random bits can be as small as [math]\displaystyle{ k^{O(1)}\log |V| }[/math]. In the second step, it has been shown by Jeanette P. Schmidt and Alan Siegel that the size of such k-perfect family can be [math]\displaystyle{ 2^{O(k)} }[/math]. Consequently, by composing the k-perfect families from both steps, a k-perfect family of size [math]\displaystyle{ 2^{O(k)}\log |V| }[/math] that maps from {1, ..., |V|} to {1, ..., k} can be obtained.

In the case of derandomizing well-coloring, where each vertex on the subgraph is colored consecutively, a k-perfect family of hash functions from {1, ..., |V|} to {1, ..., k!} is needed. A sufficient k-perfect family which maps from {1, ..., |V|} to {1, ..., k^k} can be constructed in a way similar to the approach 3 above (the first step). In particular, it is done by using nklog k random bits that are almost klog k independent, and the size of the resulting k-perfect family will be [math]\displaystyle{ k^{O(k)}\log |V| }[/math].

The derandomization of color-coding method can be easily parallelized, yielding efficient NC algorithms.

Applications

Recently, color-coding has attracted much attention in the field of bioinformatics. One example is the detection of signaling pathways in protein-protein interaction (PPI) networks. Another example is to discover and to count the number of motifs in PPI networks. Studying both signaling pathways and motifs allows a deeper understanding of the similarities and differences of many biological functions, processes, and structures among organisms.

Due to the huge amount of gene data that can be collected, searching for pathways or motifs can be highly time consuming. However, by exploiting the color-coding method, the motifs or signaling pathways with [math]\displaystyle{ k=O(\log n) }[/math] vertices in a network G with n vertices can be found very efficiently in polynomial time. Thus, this enables us to explore more complex or larger structures in PPI networks.

Further reading

• Alon, N.; Dao, P.; Hajirasouliha, I.; Hormozdiari, F.; Sahinalp, S. C. (2008). "Biomolecular network motif counting and discovery by color coding". Bioinformatics 24 (13): i241–i249. doi:10.1093/bioinformatics/btn163. PMID 18586721. 
• Hüffner, F.; Wernicke, S.; Zichner, T. (2008). "Algorithm Engineering for Color-Coding with Applications to Signaling Pathway Detection". Algorithmica 52 (2): 114–132. doi:10.1007/s00453-007-9008-7.