我们定义状态dp(i,j)为第一个字符串的前i位于第二个字符串的前j位的最长公共序列的长度,那么很容易求出状态转移方程为:
if(a[i]==b[j]) dp(i,j)=dp(i-1,j-1)+1;
else dp(i,j)=Max{dp(i-1,j),dp(i,j-1)};
相等情况很好理解,不相等的情况其实也很简单,a串的前i个字符和b串的前j个字符的最长公共序列要么是a的前i个和b的前j-1个组成的,要么是a的前i-1个和b的前j个组成的,那么此时要不要比较一下dp(i-1,j-1)呢,其实没有必要,因为a[i]!=b[j],所以我们可以“删掉”一个来形成一个子问题,也就是说,a[i]可能和b的前j-1个中的一个字符相同导致组成了一个更长的公共序列,那么dp(i-1,j-1)<=Max{dp(i-1,j),dp(i,j-1)},所以我们没必要再去比较。另外,注意dp(i,0)=0,目标是求出dp(lena,lenb),并且数据所给的字符串可能含有空格。
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
char a[1010],b[1010];
int dp[1010][1010];
int main()
{
while(cin.getline(a,1010)){
cin.getline(b,1010);
int lena=strlen(a),lenb=strlen(b);
memset(dp,0,sizeof(dp));
for(int i=1;i<=lena;i++)
for(int j=1;j<=lenb;j++){
if(a[i-1]==b[j-1]) dp[i][j]=dp[i-1][j-1]+1;
else dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
}
cout<<dp[lena][lenb]<<endl;
}
return 0;
}