Shortest common supersequence problem

From HandWiki

In computer science, the shortest common supersequence of two sequences X and Y is the shortest sequence which has X and Y as subsequences. This is a problem closely related to the longest common subsequence problem. Given two sequences X = < x1,...,xm > and Y = < y1,...,yn >, a sequence U = < u1,...,uk > is a common supersequence of X and Y if items can be removed from U to produce X and Y. A shortest common supersequence (SCS) is a common supersequence of minimal length. In the shortest common supersequence problem, two sequences X and Y are given, and the task is to find a shortest possible common supersequence of these sequences. In general, an SCS is not unique.

For two input sequences, an SCS can be formed from a longest common subsequence (LCS) easily. For example, the longest common subsequence of X[math]\displaystyle{ [1..m] = abcbdab }[/math] and Y[math]\displaystyle{ [1..n] = bdcaba }[/math] is Z[math]\displaystyle{ [1..L] = bcba }[/math]. By inserting the non-LCS symbols into Z while preserving their original order, we obtain a shortest common supersequence U[math]\displaystyle{ [1..S] = abdcabdab }[/math]. In particular, the equation [math]\displaystyle{ L + S = m + n }[/math] holds for any two input sequences.

There is no similar relationship between shortest common supersequences and longest common subsequences of three or more input sequences. (In particular, LCS and SCS are not dual problems.) However, both problems can be solved in [math]\displaystyle{ O(n^k) }[/math] time using dynamic programming, where [math]\displaystyle{ k }[/math] is the number of sequences, and [math]\displaystyle{ n }[/math] is their maximum length. For the general case of an arbitrary number of input sequences, the problem is NP-hard.[1]

Shortest common superstring

The closely related problem of finding a minimum-length string which is a superstring of a finite set of strings S = { s1,s2,...,sn } is also NP-hard.[2] Several constant factor approximations have been proposed throughout the years, and the current best known algorithm has an approximation factor of 2.475.[3] However, perhaps the simplest solution is to reformulate the problem as an instance of weighted set cover in such a way that the weight of the optimal solution to the set cover instance is less than twice the length of the shortest superstring S. One can then use the O(log(n))-approximation for weighted set-cover to obtain an O(log(n))-approximation for the shortest superstring (note that this is not a constant factor approximation).

For any string x in this alphabet, define P(x) to be the set of all strings which are substrings of x. The instance I of set cover is formulated as follows:

  • Let M be empty.
  • For each pair of strings si and sj, if the last k symbols of si are the same as the first k symbols of sj, then add a string to M that consists of the concatenation with maximal overlap of si with sj.
  • Define the universe [math]\displaystyle{ \mathcal U }[/math] of the set cover instance to be S
  • Define the set of subsets of the universe to be { P(x) | xSM }
  • Define the cost of each subset P(x) to be |x|, the length of x.

The instance I can then be solved using an algorithm for weighted set cover, and the algorithm can output an arbitrary concatenation of the strings x for which the weighted set cover algorithm outputs P(x).[4]

Example

Consider the set S = { abc, cde, fab }, which becomes the universe of the weighted set cover instance. In this case, M = { abcde, fabc }. Then the set of subsets of the universe is

[math]\displaystyle{ \begin{align} \{ P(x) | x\in S\cup M \} &= \{ P(x) | x\in \{ abc, cde, fab, abcde, fabc \} \} \\ &= \{ P(abc), P(cde), P(fab), P(abcde), P(fabc) \} \} \\ &= \{ \{a,b,c,ab,bc,abc\}, \{c,d,e,cd,de,cde\},\ldots, \{f,a,b,c,fa,ab,bc,fab,abc,fabc\} \} \} \\ \end{align} }[/math]

which have costs 3, 3, 3, 5, and 4, respectively.

References

  1. David Maier (1978). "The Complexity of Some Problems on Subsequences and Supersequences". J. ACM (ACM Press) 25 (2): 322–336. doi:10.1145/322063.322075. 
  2. Kari-Jouko Räihä, Esko Ukkonen (1981). "The shortest common supersequence problem over binary alphabet is NP-complete". Theoretical Computer Science 16 (2): 187–198. doi:10.1016/0304-3975(81)90075-x. 
  3. Matthias Englert and Nicolaos Matsakis and Pavel Vesel (2022). "Improved Approximation Guarantees for Shortest Superstrings using Cycle Classification by Overlap to Length Ratios". Proceedings of the 54th Annual ACM SIGACT Symposium on Theory of Computing: 317–330. doi:10.1145/3519935.3520001. ISBN 9781450392648. https://dl.acm.org/doi/10.1145/3519935.3520001. 
  4. Vazirani, p. 20.

External links