Extensive Simulations for Longest Common Subsequences: Finite Size Scaling, a Cavity Solution, and Configuration Space properties

Abstract

The Longest Common Subsequence (LCS) Problem asks for the longest sequence of (non-contiguous) matches between two given strings of characters. Using extensive Monte Carlo simulations, we find a finite size scaling law of the form E(L)/N =C + A/(N^1/2 ln N)+... for the mean LCS length of two random strings of size N over S letters. We provide precise estimates of C for S between 2 and 15. We consider also a related Bernoulli Matching model where the different entries of an N times M array are independently occupied with probability 1/S. In that case we find the expression of the limit of L(N,M)/N as N grows to infinity, as a function of r=M/N. This expression provides a very good approximation for the Random String model, which gets more and more accurate as S increases. The question of the ``universality class'' of the LCS problem is also considered. For the Bernoulli Matching model we find very good agreement with recent scaling predictions of Hwa and Lassig for Needleman-Wunsch sequence alignment. We find however that the variance of the LCS length has a different scaling different in the Random String model, suggesting that long-ranged correlations among the matches are relevant in this model. We finally study the ``ground state'' properties of this problem. We find in particular that the number of solutions typically grows exponentially with N, i.e. this system has a residual entropy at T=0. Also the overlap between two LCSs chosen at random is found to be self averaging and to aproach a definite value q(S)<1 as N grows.

17 pages, 12 figures. Accepted for publication in EPJ B

Cited by in corpus (11)

• Exact Asymptotic Results for a Model of Sequence Alignment
• Standard deviation of the longest common subsequence
• A Central Limit Theorem for the Length of the Longest Common Subsequences in Random Words
• Mean-Field Approximations to the Longest Common Subsequence Problem
• Bethe Ansatz in the Bernoulli Matching Model of Random Sequence Alignment
• Exact solution of the Bernoulli matching model of sequence alignment
• New alphabet-dependent morphological transition in a random RNA alignment
• Sparse long blocks and the variance of the longest common subsequences in random words
• Unzipping of two random heteropolymers: Ground state energy and finite size effects
• Sparse Long Blocks and the Micro-Structure of the Longest Common Subsequences
• A novel Recurrence-Transience transition and Tracy-Widom growth in a cellular automaton with quenched noise