基准时间限制: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示例
abcicbaabdkscab
Output示例
abca
相信大家在刚刚接触动态规划时候,基本都接触过求最长公共子序列的题目,而它的状态转移方程也并不难写:
这里的主要目的并不在此,因此我假设大家都能理解上述的状态转移方程。现在我们来思考一下上面那道题目,也即输出求得的最长公共子序列。
不难想到,我们可以通过在状态转移时记录每一步的前一步状态,然后最后通过层层回溯,找到每一步的单个序列字母,看起来问题似乎解决了?
然而事实上这样并不太行(或者说我不会…)。因为你并不知道上一步的状态是否是由 s[i]==t[j] 这个条件推导出来的,因为只有这个状态你才可以输出序列,否则会有重复(您可以试试)。
因此,我们可以换种思路,通过层层递归来输出每一步的解。表面上这种想法跟前面的一样,事实上在实现的时候并不相同。我们给每一步状态是由哪个条件得到分别做个标记,这样你可以通过标记来判断下一步递归是(i-1, j-1)还是(i-1, j)和(i, j-1),也就可以正确地输出序列。
感觉代码更容易理解,So…
<span style="font-size:14px;">//TC_1006.cpp -- 最长公共子序列Lcs
/*
if( s[i]==t[j] )
dp[i+1][j+1] = dp[i][j] + 1;
else
dp[i+1][j+1] = max(dp[i+1][j], dp[i][j+1]);
*/
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <queue>
#include <map>
using namespace std;
typedef long long ll;
const int maxn = 1000 + 10;
char a[maxn], b[maxn], k[maxn];
int dp[maxn][maxn];
int stal[maxn][maxn];
void Print(int i, int j)
{
if( i==0 || j==0 ) return ;
if( !stal[i][j] )
{
Print(i-1, j-1);
printf("%c", a[i-1]);
}
else if( stal[i][j]==1 )
Print(i-1, j);
else
Print(i, j-1);
}
int main()
{
scanf("%s %s", a, b);
int len1 = strlen(a), len2 = strlen(b);
for( int i=1; i<=len1; i++ )
{
for( int j=1; j<=len2; j++ )
{
if( a[i-1]==b[j-1] )
{
dp[i][j] = dp[i-1][j-1] + 1;
stal[i][j] = 0;
}
else
{
if( dp[i-1][j]>dp[i][j-1] )
{
dp[i][j] = dp[i-1][j];
stal[i][j] = 1;
}
else
{
dp[i][j] = dp[i][j-1];
stal[i][j] = -1;
}
}
}
}
Print(len1, len2);
return 0;
}</span>