最长公共子序列LCS简介

算法简介

一个数列 ，如果分别是两个或多个已知数列的子序列，且是所有符合此条件序列中最长的，则称为已知序列的最长公共子序列。

算法思想

貌似不难吧。。。动动脑筋。
运用DP的思想，对于当前枚举到的两个数a[i],b[j]，如果a[i]==b[j]，那么f[i][j]=f[i-1][j-1]+1，否则就取f[i-1][j]和f[i][j-1]的最大值。

可以写出状态转移方程式：

if (a[i]==b[j])    f[i][j]=f[i-1][j-1]+1;else    f[i][j]=max(f[i-1][j],f[i][j-1]);

划水

模板

#include<cstdio>#include<cstring>#include<algorithm>#define MAXN 1000using namespace std;char s1[MAXN+5],s2[MAXN+5];int f[MAXN+5][MAXN+5];int main(){    scanf("%s%s",s1,s2);    int n1=strlen(s1),n2=strlen(s2);    memset(f,0,sizeof(f));    for (int i=1;i<=n1;i++)        for (int j=1;j<=n2;j++)            if (s1[i-1]==s2[j-1])                f[i][j]=f[i-1][j-1]+1;            else                f[i][j]=max(f[i-1][j],f[i][j-1]);    printf("%d\n",f[n1][n2]);    return 0;}

算法拓展

某些题目可能要求你求出LCS。 此时我们可以把它往回倒推，比如像这样：

int x=n,y=m,node=f[n][m];    while (x>0&&y>0){        if (s1[x-1]==s2[y-1]){//如果这两个相等            x--; y--; s[node-1]=s1[x]; node--;//往回退        }        else//不相等就不取并调整            if (f[x-1][y]>f[x][y-1]) x--;            else y--;    }

还有些可能有三个数列，这时加一维即可。

 for (int i=1;i<=n;i++)        for (int j=1;j<=m;j++)            for (int k=1;k<=l;k++)                if (s1[i-1]==s2[j-1]&&s2[j-1]==s3[k-1])                    f[i][j][k]=f[i-1][j-1][k-1]+1;                else                    f[i][j][k]=max(max(f[i-1][j][k],f[i][j-1][k]),f[i][j][k-1]);

求LCS也一样：

int x=n,y=m,z=l,node=f[n][m][l];    while (x>0&&y>0&&z>0)        if (s1[x-1]==s2[y-1]&&s2[y-1]==s3[z-1]){            x--; y--; z--; node--; s[node]=s1[x];        }        else{            if (f[x][y][z]==f[x-1][y][z]) x--;            if (f[x][y][z]==f[x][y-1][z]) y--;            if (f[x][y][z]==f[x][y][z-1]) z--;        }