思路:
动态规划
dp[i][j] = s1的前i个和s2的前j个的公共子序列长度
递推公式:
if s1[i-1] = s2[j-1]:
dp[i][j] = dp[i-1][j-1]+1;
else:
dp[i][j] = max(dp[i][j-1],dp[i-1][j]);
(证明:1、dp[i][j]显然不小于dp[i][j-1]和dp[i-1][j]
2、用反证法可得dp[i][j]不同时大于dp[i][j-1]和dp[i-1][j]
故等式成立)
初始情况:
dp[i][0] = 0;
dp[0][j] = 0;
确定dp的顺序:
从上到下、从左到右,
或从左到右、从上到下
#include <cstdio>
#include <iostream>
using namespace std;
int dp[10000][10000];
//公共子序列
int main(){
string s1,s2;
while(cin>>s1>>s2){
int n1 = s1.size();
int n2 = s2.size();
int i,j;
for(i=0;i<=n1;i++){
dp[i][0] = 0;
}
for(i=0;i<=n2;i++){
dp[0][i] = 0;
}
for(i=1;i<=n1;i++){
for(j=1;j<=n2;j++){
if(s1[i-1]==s2[j-1]){
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[n1][n2]);
}
return 0;
}