LCS 最大公共子序列算法详解

LCS问题简述：最长公共子序列

一个数列 S，如果分别是两个或多个已知数列的子序列，且是所有符合此条件序列中最长的，则S 称为已知序列的最长公共子序列。

LCS问题的分支：最长公共子串与最长公共子序列

子串（Substring）是串的一个连续的部分，子序列（Subsequence）则是从不改变序列的顺序，而从序列中去掉任意的元素而获得的新序列；更简略地说，前者（子串）的字符的位置必须连续，后者（子序列LCS）则不必。比如字符串acdfg同akdfc的最长公共子串为df，而他们的最长公共子序列是adf。

给出str1串=“abcd”, str2串= “becd”；

当str1用前一个的时候，str2分别用前1234个扫一遍出来第一行的值，当第i个与第j个不相等的时候， dp[i][j] 等于他前面扫过的所有面积中最大的那个数（因为要找到最大的那个么）而最大的数，一行当中目前最大的就是dp[j][i-1](这里j是行数，i是列数)因为每一个行的值，都是前面所有扫过的最大值）而前j-1行最大的数就是dp[j-1][i]（因为只能选前i个，所以面积有所缩小，面积就是以当前的i和j为长宽包括的面积）然后当第i个与第j个相等的时候，他就等于dp[i-1][j-1]+1(因为dp[i-1][j-1]是当前最大的那个公共串，前面能够证明)。下图就是这个str1和str2的具体示意图。

而hdu 1159 和 poj 1458 都是这种题

模板代码如下:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

const int maxn=1010;
int dp[maxn][maxn];
char a[maxn],b[maxn];

using namespace std;

int main()
{
	while(~scanf("%s%s",a+1,b+1))
	{
		memset(dp,0,sizeof(dp));
		int i,j;
		for( i=1;a[i];i++)
		{
			for(j=1;b[j];j++)
			{
				if(a[i]==b[j])dp[i][j]=dp[i-1][j-1]+1;
				else dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
			}
		}
		printf("%d\n",dp[i-1][j-1]);
	}
	return 0;
}