题目
给出两个字符串A B,求A与B的最长公共子序列(子序列不要求是连续的)。
比如两个串为:abcicba
abdkscabab是两个串的子序列,abc也是,abca也是,其中abca是这两个字符串最长的子序列。
Input
第1行:字符串A
第2行:字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列,如果有多个,随意输出1个。
Sample Input
abcicba
abdkscab
Sample Output
abca
题目大意:
两个字符串中最长的公共子字符串,此串不一定连续。
设str1,str2, dp[i][j] 表示str1前i个字符与str2前j个字符的最长公共子字符串。
那么转移的方程是
if (str[i]==str[j]) dp[i][j] = dp[i][j]+1;
else dp[i][j] = max(dp[i-1][j],dp[i][j-1]);
这个题需要打印子字符串,那么在状态转移的时候,记录一下上一个状态。
易错点
第一次记录路径的时候。
for ( int i=1; i<=la; i++ ) {
for( int j=1; j<=lb; 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]);
}
}
if ( cnt+1== dp[i][lb] ) {
ans[cnt] = A[i-1];
cnt++;
}
}
这种记录字符串的方式,错误的以为当前的状态会延续到最后的状态。
例如 str1 = “abcc” str2 = “accb”;
在i==2的时候,str1[1]==str2[3],这种情况下是dp[2][4]的最优解,但并不是全局的最优解。
看了网上的题解之后写出了正确的记录前一个状态的方法。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
const int N = 1000+10;
char A[N];
char B[N];
char ans[N];
int dp[N][N];
int path[N][N][2];
void print( int x,int y ) {
if ( x==0 || y==0 )
return;
print( path[x][y][0],path[x][y][1] ) ;
if ( A[x-1]==B[y-1] ) {
printf("%c",A[x-1]);
}
}
int main()
{
cin>>A>>B;
int la = strlen(A);
int lb = strlen(B);
int cnt = 0;
// dp[i][j] 是 最长的公共字符串的长度
for ( int i=1; i<=la; i++ ) {
for( int j=1; j<=lb; j++ ) {
if ( A[i-1]==B[j-1] ) {
dp[i][j] = dp[i-1][j-1] + 1 ;
path[i][j][0] = i-1;
path[i][j][1] = j-1;
} else {
int index1 = dp[i][j-1];
int index2 = dp[i-1][j];
if (index1>index2) {
dp[i][j] = index1;
path[i][j][0] = i;
path[i][j][1] = j-1;
} else {
dp[i][j] = index2;
path[i][j][0] = i-1;
path[i][j][1] = j;
}
}
}
}
print( la,lb );
return 0 ;
}