题目链接这里呀
1006 最长公共子序列Lcs
基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注
给出两个字符串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
分析:注意哦,子序列不是子段,子序列可以是不连续的。
求最长公共子序列(LCS)用动态规划。
两字符串str1,str2,用dp[i][j]表示str1[1]~str1[i]和str2[1]~str2[j]两个子段最长公共子序列的长度。
若str1[i]==str2[j],那么此时的LCS长度就相当于在(i-1,j-1)的LCS长度上加1,即dp[i][j]=dp[i-1][j-1]+1;
若str1[i]!=str2[j],此时的LCS长度应取(i-1,j)和(i,j-1)的较大者,因为要求最长嘛!即dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
所以,dp[i][j]={ dp[i-1][j-1]+1; (str1[i]==str2[j])
max(dp[i-1][j],dp[i][j-1];(str1[i]!=str2[j])
}
dp[len1][len2]即为LCS长度;
若要求出LCS,还需要用到栈的数据结构哦。
用两个‘指针’ant1,ant2倒序指向两字符串中的字符,
若str1[ant1]==str1[ant2],则入栈(任一个就行哦),同时ant1–,ant2–;
若不相等,肯定首先字符不能入栈了哦,那么两个‘指针’怎么办呢?
两种情况:
如果dp[ant1-1][ant2]>dp[ant1][ant2-1],说明(ant1-1,ant2)这一部分的LCS比(ant1,ant2-1)的要长哦,那么当然就要研究看str1[ant1-1]和str1[ant2]是否相等咯,也就是ant1–,前移;
否则,ant2–。
LCS的全部字符入栈后,再依次输出。
代码如下:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX_N 1000
#include<stack>
int dp[MAX_N+11][MAX_N+11];
char str1[MAX_N+11],str2[MAX_N+11];
int main()
{
scanf("%s %s",str1+1,str2+1);
memset(dp,0,sizeof(dp));
str1[0]='0';
str2[0]='0';
int len1=strlen(str1)-1;
int len2=strlen(str2)-1;
int ans=0;
int i,j;
for(i=1;i<=len1;i++)
{
for(j=1;j<=len2;j++)
{
if(str1[i]==str2[j])
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
}
ans=max(ans,dp[i][len2]);//更新为当前最长,最后得出全部最长
}
int ant1=len1;
int ant2=len2;
stack<char> stk;
while(ant1>0&&ant2>0)
{
if(str1[ant1]==str2[ant2])//若对应两字符相等,说明是LCS组成部分,入栈
{
stk.push(str1[ant1]);
ant1--;
ant2--;
}
else if(dp[ant1-1][ant2]>dp[ant1][ant2-1])
ant1--;
else
ant2--;
}
while(!stk.empty() )
{
printf("%c",stk.top());
stk.pop() ;
}
printf("\n");
return 0;
}