动态规划之最长公共子序列

题目描述：

给出两个字符串，求出这样的一个最长的公共子序列的长度：子序列中的每个字符都能在两个原串中找到，而且每个字符的先后顺序和原串中的先后顺序一致。

范例输入：

abcfbc  abfcab
programming  contest
abcd  mnp

范例输出：

4
2
0

---

示例代码：

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

---

讲解：
1.

设置maxlen数组来存放 s1以i-1为终点，s2以j-1为终点 的最长公共子序列的长度(maxlen[2][2]即表示s1[1]为终点，s2[1]为终点所对应的最长公共子序列的长度).

由以上maxlen数组的定义可以看出来，maxlen[i][0]均为0，maxlen[0][j]均为0.
我们以此为初始条件进行下一步的递推.

当s1[i-1]==s2[j-1]时，有maxlen[i][j]=maxlen[i-1][j-1]+1,这是易于理解的.

当s1[i-1]!=s2[j-1]时，有maxlen[i][j]=max(maxlen[i-1][j],maxlen[i][j-1]);此处不易理解，作以下解释：

1.举一个栗子： 现在有 abccfc 与 abcfba 两个字符串，当我们求maxlen[4][4]时，也就是s1以f为终点，s2以b为终点的最长公共子序列的长度，很显然，f！=b，所以调用maxlen[i][j]=max(maxlen[i-1][j],maxlen[i][j-1])
所以也就是求maxlen[3][4]与maxlen[4][3]的较大的那一个，而maxlen[4][3]为4，maxlen[3][4]为3，所以maxlen[4][4]为4.
2.更加浅显的想 ：maxlen[4][4]肯定不可能比maxlen[4][3]和maxlen[3][4]小，因为其中一个序列多了一个元素，自然不可能比短的串还要小.
那么为什么不能都大于这两者呢？
我们分开来讨论，如果maxlen[4][4]>maxlen[4][3]，那么相比之下s2中新增加的那个元素一定是在最长子序列中的，且处于末尾。
另外如果maxlen[4][4]也同时大于maxlen[3][4]，那么同理s1中新增加的那个元素一定在最长子序列中，且处于末尾。
继续想，如果maxlen[4][4]同时满足以上两个条件，那么，一定会有s1新增的元素与s2新增的元素相等，这就与前提条件 s1[i-1]!=s2[j-1] 不相符.

最长公共子序列就完成啦！