输出最长公共子序列

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示例

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>