LCS算法

1. 问题
   [描述算法问题，首选形式化方式（数学语言），其次才是非形式化方式（日常语言）]
   

2. 解析
   [问题的理解和推导，可用电子版直接在此编写，也可用纸笔推导，拍照嵌入本文档]
   （1）如果 xn=ym，即X的最后一个元素与Y的最后一个元素相同，这说明该元素一定位于公共子序列中。因此，现在只需要找：LCS(Xn-1，Ym-1)
   （2）如果xn != ym，这下要麻烦一点，因为它产生了两个子问题：LCS(Xn-1，Ym) 和 LCS(Xn，Ym-1)
   因为序列X和序列Y的最后一个元素不相等，那说明最后一个元素不可能是最长公共子序列中的元素。
   LCS(Xn-1，Ym)表示：最长公共序列可以在(x1,x2,…x(n-1)) 和 (y1,y2,…yn)中找。
   LCS(Xn，Ym-1)表示：最长公共序列可以在(x1,x2,…xn) 和 (y1,y2,…y(n-1))中找。
   求解上面两个子问题，得到的公共子序列谁最长，那谁就是 LCS（X,Y）
   

3. 设计
   [核心伪代码]
   int c[200][200]; //用c[i][j]记录X[i]与Y[j] 的LCS 的长度
   int b[200][200]; //b[i][j]记录c[i][j]是通过哪一个子问题的值求得的，以决定搜索的方向
   char f[200];

/取c[i-1][j]和c[i][j-1]的最大值,并记录c[i][j]是通过哪一个子问题的值求得的，以决定搜索的方向/
int Max(int m, int n, int i, int j)
{
if (m > n)
{
b[i][j] = -1;
return m;
}
else
{
b[i][j] = 1;
return n;
}
}

/递归打印LCS的元素内容/
void print(int i, int j, int s, char x[], char y[])
{
if (b[i][j] == 0)
{
f[s - 1] = x[i - 1];
i–; j–; s–;
print(i, j, s, x, y);
}
else if (b[i][j] == -1)
{
i–;
print(i, j, s, x, y);
}
else if (b[i][j] == 1)
{
j–;
print(i, j, s, x, y);
}
}
int LCS(char x[], char y[])
{
int i, j;
int x_len, y_len;
x_len = strlen(x);
y_len = strlen(y);
printf(" “);
for (i = 0; i < y_len; i++)
{
printf(”%c “, y[i]);
}
printf(”\n");
for (i = 1; i <= x_len; i++)
{
printf("%c “, x[i - 1]);
for (j = 1; j <= y_len; j++)
{
if (x[i - 1] == y[j - 1])
{
c[i][j] = c[i - 1][j - 1] + 1;
b[i][j] = 0;
printf(”%d “, c[i][j]);
}
else
{
c[i][j] = Max(c[i - 1][j], c[i][j - 1], i, j);
printf(”%d “, c[i][j]);
}
}
printf(”\n");
}
printf(“X和Y的LCS是:”);
print(x_len, y_len, c[x_len][y_len], x, y);
printf("%s", f);
printf("\n");
return c[x_len][y_len];
}
4. 分析
[算法复杂度推导]
O（nlogn）
5. 源码
[github源码地址]
https://github.com/zddhandsome/Algorithm/blob/main/LCS%E7%AE%97%E6%B3%95.cpp