LCS(Longest Common Subsequence),最长公共子序列,典型的dp问题。
问题:设字符串A、B,求A和B的LCS。
二维dp:用dp[i][j]数组表示以A[i]结尾、B[j]结尾的字符串的LCS长度,那么
(1)dp[i][j] = dp[i - 1][j - 1] + 1 if A[i] == B[j]
(2)dp[i][j] = max(dp[i - 1][j], dp[i][j - 1] if A[i] != B[j]
这种方法的时间复杂度是O(nm),空间复杂度也是O(nm),记录路径时,可以记录dp[i][j]是由哪一步转化而来。
#define N 1005
#define DIAGONAL 0
#define LEFT 1
#define UP 2
char str1[N],str2[N];
int dp[N][N],path[N][N];
int dir[3][2] = {{-1,-1},{-1,0},{0,-1}};
void _print(int i,int j)
{
if(i == 1)
{
printf("%c",str1[i]);
return;
}
if(j == 1)
{
printf("%c",str2[j]);
return;
}
int ni,nj;
ni = i + dir[path[i][j]][0];
nj = j + dir[path[i][j]][1];
_print(ni,nj);
if(path[i][j] == DIAGONAL)
printf(" %c",str1[i]);
}
void print()
{
_print(strlen(str1) - 1,strlen(str2) - 1);
puts("");
}
int lcs()
{
memset(dp,0,sizeof(dp));
memset(path,0,sizeof(path));
for(int i = 1;str1[i];i++)
for(int j = 1;str2[j];j++)
if(str1[i] == str2[j])
dp[i][j] = dp[i - 1][j - 1] + 1;
else
{
dp[i][j] = max(dp[i - 1][j],dp[i][j - 1]);
if(dp[i][j] == dp[i - 1][j])
path[i][j] = LEFT;
else
path[i][j] = UP;
}
return dp[strlen(str1) - 1][strlen(str2) - 1];
}