LCS(最长公共子序列)

题意描述

求两个字符串的最长公共子序列的长度

动态规划

用二维数组 C[i][j] 记录串x1x2⋯xi与y1y2⋯yj 的 LCS长度，则可得到状态转移方程:

代码实现：

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#define MAXN 1000
using namespace std;
typedef long long ll;
int dp[MAXN+5][MAXN+5];
int main()
{
	ios::sync_with_stdio(false);
	string a,b;
	while(cin>>a>>b){
		memset(dp,0,sizeof(dp));
		int len1=a.size();
		int len2=b.size();
		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;
				else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
			}
		}
		cout<<dp[len1][len2]<<endl;

滚动数组优化

上述解法空间复杂度为O(n^2)，当n特别大时，空间耗费极高。我们可以发现，C[i][j]的转移只与
C[i-1][j-1]、C[i][j-1]、C[i-1][j]有关。所以二维数组第一维最大容量为2，我们只需要将第一维模2即可。
代码实现：

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#define MAXN 1000
using namespace std;
typedef long long ll;
int dp[5][MAXN+5];
int main()
{
	ios::sync_with_stdio(false);
	string a,b;
	while(cin>>a>>b){
		memset(dp,0,sizeof(dp));
		int len1=a.size();
		int len2=b.size();
		for(int i=1;i<=len1;++i){
			for(int j=1;j<=len2;++j){
				if(a[i-1]==b[j-1]) dp[i%2][j]=dp[(i-1)%2][j-1]+1;
				else dp[i%2][j]=max(dp[(i-1)%2][j],dp[i%2][j-1]);
			}
		}
		cout<<dp[len1%2][len2]<<endl;
	}
 } 

推荐LCS习题 HDU Common Subsequence
题目链接：http://poj.org/problem?id=1458
以上两份代码皆可通过此题。