LCS问题一直都知道如何去求长度,但是还真是没真正研究过如何去构造。今天算法课老师一直在讲一些我们不可能现在去研究的东西,我就看了看。大致算法为这样:
void LCSLength(int m,int n,char * x, char * y,int * * c,int * * b)
{
int i,j;
for(i=1;i<=m;i++) c[i][0]=0;
for(i=1;i<=n;i++) c[0][i]=0;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++){
if(x[i]==y[j]) {c[i][j]=c[i-1][j-1]+1;b[i][j]=1;}
else if(c[i-1][j]>=c[i][j-1]) {c[i][j]=c[i-1][j];b[i][j]=2;}
else {c[i][j]=c[i][j-1];b[i][j]=3;}
}
}
void LCS(int i,int j, char * x,int * * b)
{
if(i==0 || j==0) return ;
if(b[i][j]==1) { LCS(i-1,j-1,x,b);cout<<x[i];}
else if(b[i][j]==2) LCS(i-1,j,x,b);
else LCS(i,j-1,x,b);
}
这里LCSLength里,当b[i][j]为1时,证明此点是相同的,2的话代表不相同且取上面的,3代表取前面的。由此可以得出从后往前如果b[i][j]为3,则向前回溯,2则向上回溯,1的话则输出,这样能保证最后最长公共子序列的情况的取法为真。
但是呢,如果想构造,那么就不能用滚动数组了,意思就是空间复杂度必为O(mn),时间复杂度自然也为O(mn)。滚动数组的空间复杂度为O(min{m,n})。
以上。