给出两个字符串A B,求A与B的最长公共子序列(子序列不要求是连续的)。
比如两个串为:
abcicba
abdkscab
ab是两个串的子序列,abc也是,abca也是,其中abca是这两个字符串最长的子序列。
Input
第1行:字符串A
第2行:字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列,如果有多个,随意输出1个。
Input示例
abcicba
abdkscab
Output示例
abca
解析:
最长公共子序列是典型的动态规划问题,把一个大问题拆分成多个子问题,每个子问题都存储起来方便下次调用,每次都只处理当前状态,处理时可根据已知情况得出结果,也就是按照一个顺序处理,每次处理时都根据之前的处理结果得出当前方案。最长公共子序列需要把两个串中任意两个数都要对比,所以我们用二维数组来存储对比结果,这样对比下去时间复杂度是O(m*n)递推公式为:
if(s1[i]==s2[j])
m[i][j]=m[i-1][j-1]+1;
else
m[i][j]=ma(m[i-1][j],m[i][j-1]);
这样递推下来m[i-1][j-1]就是最长公共子序列的长度。
关于如何输出他们可以根据之前建立的表和递推关系来回溯。
(抄的一张图)
我的回溯方式是:
while(i>0&&j>0)
{
if(s1[i]==s2[j])
{
s[n]=s1[i];
i--;
j--;
n--;
}
else
if(m[i-1][j]>m[i][j-1])
i--;
else
j--;
}
详细代码:
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define N 1010
int ma(int a,int b)
{
if(a>b)
{
return a;
}
return b;
}
char s1[N],s2[N],s[N];
int n1,n2,n,m[N][N];
int main()
{
int i,j,k;
scanf("%s %s",s1+1,s2+1);
memset(m,0,sizeof(s));
for(i=1;s1[i];i++){
for(j=1;s2[j];j++){
if(s1[i]==s2[j])
{
m[i][j]=m[i-1][j-1]+1;
}
else
{
m[i][j]=ma(m[i-1][j],m[i][j-1]);
}
}
}
i--;
j--;
n=N-1;
while(i>0&&j>0)
{
if(s1[i]==s2[j])
{
s[n]=s1[i];
i--;
j--;
n--;
}
else
if(m[i-1][j]>m[i][j-1])
i--;
else
j--;
}
printf("%s",s+n+1);
return 0;
}