算法第四周 H : 最长公共子序列

Description

求两个序列的最每组测试样例都为一行，两组字符串，每组不超过1000，用空格隔开。求最长公共子序列，都为小写字母。

Input

每组测试样例都为一行，两组字符串，每组不超过1000，用空格隔开。

Output

对于每个测试实例，输出最长公共子序列的长度，每个实例的输出占一行。

Sample Input

abcfbc abfcab
programming contest 
abcd mnp

Sample Output

4
2
0

首先要明白什么叫最长公共子序列，比如两个序列s1={1,3,4,5,6,7,7,8},s2={3,5,7,4,8,6,7,8,2}，

则最长公共子序列为{3,4,6,7,8} ，只需要求出两个序列的公共部分，不要求连续，但顺序必须一致

解题思路：

使用动态规划，定义一个二维数组dp[][]。首先我们需要理解dp代表的含义是什么

dp[i][j] 的含义是：对于 s1[1..i] 和 s2[1..j]，它们的 LCS 长度是 dp[i][j]

比如dp[2][4] ：对于s1[1...2]（即''ac''）和s2[1...4]（即"babc"），它们的LCS长度为2，所以我们想得到dp[2][4]的答案为2。

然后我们就要根据字符串的特征，来寻找状态转移方程，来求取dp 

例一个例子（S1={1,3,4,5,6,7,7,8}和S2={3,5,7,4,8,6,7,8,2}），并结合上图来说：

       假如S1的最后一个元素 与 S2的最后一个元素相等，那么S1和S2的LCS就等于 {S1减去最后一个元素} 与 {S2减去最后一个元素} 的 LCS  再加上 S1和S2相等的最后一个元素。

if (str1[i - 1] == str2[j - 1]) dp[i][j] = dp[i - 1][j - 1] + 1;

       假如S1的最后一个元素 与 S2的最后一个元素不等（本例子就是属于这种情况），那么S1和S2的LCS就等于 ： {S1减去最后一个元素} 与 S2 的LCS， {S2减去最后一个元素} 与 S1 的LCS 中的最大的那个序列。
 

else
dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);

 完整代码：

#include<iostream>
#include<cstring>
#include<vector>
#include<math.h>
using namespace std;
int dp[1005][1005];
int main()
{
    string str1, str2;
	while(cin>>str1>>str2)
	{
		memset(dp, 0, sizeof(dp));
		int len1 = str1.length();
		int len2 = str2.length();
		for (int i = 1; i <=len1; i++)
		{
			for (int j = 1; j <=len2; j++)
			{
                //S1的最后一个元素 与 S2的最后一个元素相等
				if (str1[i - 1] == str2[j - 1]) dp[i][j] = dp[i - 1][j - 1] + 1;
				else
                //不相等，则选择较大的那一个
					dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
			}
		}
		cout << dp[len1][len2] << endl;
	}
}