最长公共子序列针对小字符集的算法

Posted ACMsong

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了最长公共子序列针对小字符集的算法相关的知识,希望对你有一定的参考价值。

一般对于两个字符串,长度分别为n和m,其时间复杂度为O(nm)。

但是针对小字符集的情况,可以把复杂度降低到O(n^2+km),其中n为两个字符串较短的长度。这种方法对于两个字符串长度相差很大的情况比O(nm)要优化很多。

 

就假设所有的字符都是小写字母,这样就符合小字符集的前提了。设较短的字符串为S1,较长的字符串为S2。字符串下标从1开始。

S2字符串每个位置右边第一个字符是可以通过O(km)预处理得到的。其中k为小字符集的字符个数,m为较长的那个字符串的长度。

用next[i][j]表示S2[i]右边第一个(char)(‘a‘+j)的位置。

设dp[i][j]表示S1匹配了前i位,长度为j的最长公共子序列与S2匹配到的最靠左的位置。如果不存在则为length(S2)+1。

dp[i][0] = 0

如果S2的dp[i-1][j-1]这个位置右边第一个与S1[i]相等的位置存在,那么dp[i][j] = min{ dp[i-1][j], next[dp[i-1][j-1]][S1[i]] }。

否则dp[i][j] = dp[i-1][j]。

对于每一个不是-1的dp,记录下j,最后取一个最大就是最长公共子序列了。

因此,总的复杂度为O(n^2+km)。

 

具体实现:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int INF=0x3f3f3f3f;
const int maxn=1005;
const int maxm=1000005;
char s1[maxn],s2[maxm];
int dp[maxn][maxn];
int next[maxm][26]; 

int main()
{
    scanf("%s%s",s1+1,s2+1);  
    int l1=strlen(s1+1);
    int l2=strlen(s2+1);
    for (int i=0;i<maxm;i++)
        for (int j=0;j<26;j++)
            next[i][j]=l2+1;
    for (int i=0;i<maxn;i++)
        for (int j=0;j<maxn;j++)
            dp[i][j]=l2+1;
    for (int i=l2-1;i>=0;i--)
    {
        for (int j=0;j<26;j++)
        {
            char cc=‘a‘+j;
            if (s2[i+1]==cc) next[i][j]=i+1;
            else next[i][j]=next[i+1][j];
        }
    }
    for (int i=1;i<=l1;i++) dp[i][0]=0;
    int ans=0;
    for (int i=1;i<=l1;i++)
    {
        for (int j=1;j<=i;j++)
        {
            if (next[dp[i-1][j-1]][s1[i]-‘a‘]!=l2+1) dp[i][j]=min(dp[i-1][j],next[dp[i-1][j-1]][s1[i]-‘a‘]);
            else dp[i][j]=dp[i-1][j];
            if (dp[i][j]!=l2+1) ans=max(ans,j);
        }
    }
    printf("%d\n",ans);
    return 0;
}

 

以上是关于最长公共子序列针对小字符集的算法的主要内容,如果未能解决你的问题,请参考以下文章

最长公共子序列算法求破!杭电1159

算法设计之最长公共子序列

算法设计之最长公共子序列

算法设计之最长公共子序列

5.23算法练习最长公共子序列

ACM/ICPC 之 最长公共子序列计数及其回溯算法(51Nod-1006(最长公共子序列))