最长公共子序列,英文缩写为LCS(Longest Common Subsequence)。
其定义是,一个序列 S ,如果分别是两个或多个已知序列的子序列,
且是所有符合此条件序列中最长的,则 S 称为已知序列的最长公共子序列。
而最长公共子串(要求连续)和最长公共子序列是不同的
设序列 X 存储在数组 x[m]中,序列 Y 存储在数组 y[n]中, 二维数组 L[m+1][n+1]
存储最长公共子序列的长度的迭代过程, S[ m + 1] [n + 1]存储相应的状态, 最长公共子
序列存储在数组 z[ k]中
int CommonOrder(int m, int n, int x[ ] , int y[ ] , int z[ ])
{
int i , j , k ;
for (j = 0; j <= n; j ++ ) // 初始化第 0 行
L[0][j] = 0;
for (i = 0; j <= m; i++ ) // 初始化第 0 列
L[i] [0] = 0;
for (i = 1; i<= m; i++ )
for (j = 1; j <= n; j++ )
if (x[i] == y[j]) { L[i][j] = L[i - 1] [j - 1] + 1; S[i][j] = 1; }
else if (L[i][j - 1] >= L[i - 1][j] ) { L[i][j] = L[i] [j - 1]; S[i][j] = 2; }
else {L[i][j] = L[i - 1][j]; S[i][j] = 3; }
i = m; j = n; k = L[m] [n];
while(i > 0 && j > 0)
{
if (S[i][j] == 1) { z[k] = x[i] ; k -- ; i -- ; j -- ; }
else if (S[i][j] == 2) j -- ;
else i -- ;
}
return L[m][n] ;
}
int main(void)
{
int i,j;
int x[10] ;
int y[10] ;
int z[10] ;
i = CommonOrder(10,10,x,y,z) ;
printf("%d\n",i) ;
putchar('\n');
for(i = 0 ; i < 10 ; i++){
for(j = 0 ; j < 10 ; j++){
printf("%d ",L[i][j]);
}
putchar('\n');
}
return 0 ;
}
运行结果:
2
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 1
1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1
2 2 2 2 2 2 2 2 2 2