题意描述
求两个字符串的最长公共子序列的长度
动态规划
用二维数组 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
以上两份代码皆可通过此题。