给出两个字符串A B,求A与B的最长公共子序列(子序列不要求是连续的)。


比如两个串为:




abcicba


abdkscab




ab是两个串的子序列,abc也是,abca也是,其中abca是这两个字符串最长的子序列。




Input

第1行:字符串A
第2行:字符串B
(A,B的长度 <= 1000)


Output

输出最长的子序列,如果有多个,随意输出1个。


Input示例

abcicba
abdkscab


Output示例

abca



分析:动态规划

 

dp[i][j] 表示字符串A以第i个位置,字符串B以第j个位置的最长公共子序列的长度

dp[i][j] = dp[i - 1][j - 1] + 1 if a[i] == a[j]

else dp[i][j] == max(dp[i - 1][j] , dp[i][j - 1]);

最大长度就是 dp[n][m] ,n 为A的长度 ,m为B的长度

还原字符串 ,只需要回到 dp[i][j] 刚开始发生改变的地方即可

 

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn = 1E3 + 10;
char a[maxn],b[maxn],ans[maxn];
int dp[maxn][maxn];
int main(){
scanf("%s%s",a + 1,b + 1);
int n = strlen(a+1),m = strlen(b+1);
memset(dp,0,sizeof(dp));
for(int i = 1 ; i <= n ; ++i){
for(int j = 1 ; j <= m ; ++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]);
}
}
int cur = 0;
for(int i = n,j = m;dp[i][j];--i,--j){//返回到第一次更新值的地方
while(dp[i][j] == dp[i - 1][j]) --i;
while(dp[i][j] == dp[i][j - 1]) --j;
ans[cur++] = a[i];
}
reverse(ans,ans+cur);
ans[cur] = '\0';
printf("%s\n",ans);
return 0;
}