25-最长公共子序列

问题描述：

给出两个字符串，求出这样的一个最长的公共子序列的长度：子序列中的每个字符都能在两个原串中找到，而且每个字符的先后顺序和原串中的先后顺序一致。
Sample Input：
abcfbc abfcab
programming contest
abcd mnp
Sample Output：
4
2
0

问题分析：

第一个案例中最长的公共子序列为abcb，第二个为on，第三个没有。

前面说过动态规划的四个步骤为：1、将原问题分解为子问题；2、确定状态；3、确定一些初始状态（边界状态）的值；4. 确定状态转移方程。这四步中最难的是如何将问题划分，确定状态。而且没有什么固定的套路可寻，只能具体问题具体分析，多做多练。

1、划分子问题，寻找状态
输入两个串 s1,s2，设 MaxLen (i,j) 表示 : s1 的左边 i 个字符形成的子串与 s2 左边的 j 个字符形成的子串的最长公共子序列的长度 (i,j 从0开始算）。

MaxLen(i,j) 就是本题的“状态”。假定 len1 = strlen(s1)，len2 = strlen(s2），那么题目就是要求 MaxLen(len1,len2)

2、确定初始状态
其中一个字符长度为0时，最长公共子序列为0。
MaxLen(n,0) = 0 ( n=0…len1）
MaxLen(0,n) = 0 ( n=0…len2）

3、确定状态转移方程

if(s1[i-1] == s2[j-1]) //s1 的最左边字符是 s1[0]
	MaxLen(i,j) = MaxLen(i-1,j-1) + 1;
else
	MaxLen(i,j) = max(MaxLen(i,j-1), MaxLen(i-1,j));

下面来讲解一下状态转移方程。

1. 当s1[i-1] == s2[j-1]时，自然s1[i-1]与s2[j-1]匹配，使得最长公共子序列长度加1。
2. 当s1[i-1] != s2[j-1]时，为什么MaxLen(i,j) = max(MaxLen(i,j-1), MaxLen(i-1,j));成立。下面用反证法来证明一下。需要证明 MaxLen(S1,S2) 不会比 MaxLen (S1,S2_j-1) 和 MaxLen(S1_i-1,S2) 两着之中任何一个小，也不会比两者都大，则等式必然成立。

首先 MaxLen(s1,s2) 不会比MaxLen (s1,s2_j-1) 和 MaxLen(s1_i-1,s2) 任何一个小，因为字符串长度增加了，公共子序列只可能增加或者不变。其次，若MaxLen(s1,s2) 比MaxLen (s1,s2_j-1) 大，说明s2_j 是匹配序列的最后一个字符，若MaxLen(s1,s2) 比 MaxLen(s1_i-1,s2) 大，说明s1_i是匹配序列的最后一个字符，则一定有s1[i-1] == s2[j-1]，这一结果与s1[i-1] != s2[j-1]相矛盾。

从图中看的比较清楚，当s1[i-1] == s2[j-1]时，maxlen[i][j]等于左上角元素加1，当s1[i-1] != s2[j-1]时，maxlen[i][j]等于正上方与左侧元素的较大值。

复杂度分析：
m,n 是两个字符串长度，所以一共有 m*n 种状态，每种状态求解的时间复杂度为常数，所以总的时间复杂度为 O(mn)。

完整C++代码

#include<iostream>
#include<algorithm>
using namespace std;
char s1[1000];
char s2[1000];
int maxlen[1000][1000];
int main()
{
	while (cin >> s1 >> s2)
	{
		int len1 = strlen(s1);
		int len2 = strlen(s2);
		for (int i = 0; i <= len1; i++)
			maxlen[i][0] = 0;
		for (int j = 0; j <= len2; j++)
			maxlen[0][j] = 0;
		for (int i = 1; i <= len1; i++)
		{
			for (int j = 1; j <= len2; j++)
			{
				if (s1[i-1] == s2[j-1])
					maxlen[i][j] = maxlen[i - 1][j - 1] + 1;
				else
					maxlen[i][j] = max(maxlen[i-1][j],maxlen[i][j-1]);
			}
		}
		cout << maxlen[len1][len2]<<endl;
	}
	return 0;
}