这个是最长公共子序列的升级版,不仅需要能算出公共子序列的最大长度
还要能输出公共子序列,并且还要能把两个序列不同的地方分别输出来,这里用了算法导论的LCS算法 ,然后加上后面的输出字符的算法
通过形象的字符表来看到一个一个的字符输出,直到每个字符串的边缘值都被输出来,如此便得到了两个字符的完整合并值
这里为了表示算法的直观方便,我把数组和字符图都输出来了,便于观看效果
#include<stdio.h>
#include<string.h>
#define MAX_LEN 1000
char sz1[MAX_LEN];
char sz2[MAX_LEN];
int aMaxLen[MAX_LEN][MAX_LEN];
char b[MAX_LEN][MAX_LEN];
void printLCS(int i,int j)
{
if(i==0&&j==0)
return;
if(b[i][j]=='\\')
{
printLCS(i-1,j-1);
printf("%c",sz1[i]);
}
else if(b[i][j]=='|')
{
printLCS(i-1,j);
printf("%c",sz1[i]);
}
else
{
printLCS(i,j-1);
printf("%c",sz2[j]);
}
}
int main()
{
while(scanf("%s%s",sz1+1,sz2+1)>0)
{
int nlen1=strlen(sz1+1);
int nlen2=strlen(sz2+1);
int nTmp;
int i,j;
for(i=0;i<=nlen1;i++)
{
aMaxLen[i][0]=0;
b[i][0]='|';
}
for(j=0;j<=nlen2;j++)
{
aMaxLen[0][j]=0;
b[0][j]='-';
}
for(i=1;i<=nlen1;i++)
{
for(j=1;j<=nlen2;j++)
{
if(sz1[i]==sz2[j])
{
aMaxLen[i][j]=aMaxLen[i-1][j-1]+1;
b[i][j]='\\';
}
else if(aMaxLen[i][j-1]>aMaxLen[i-1][j])
{
aMaxLen[i][j]=aMaxLen[i][j-1];
b[i][j]='-';
}
else
{
aMaxLen[i][j]=aMaxLen[i-1][j];
b[i][j]='|';
}
}
}
printf("%d\n",aMaxLen[nlen1][nlen2]);
printLCS(nlen1,nlen2);
printf("\n");
for(i=1;i<=nlen1;i++)
{
for(j=1;j<=nlen2;j++)
printf("%d ",aMaxLen[i][j]);
printf("\n");
}
for(i=0;i<=nlen1;i++)
{
for(j=0;j<=nlen2;j++)
printf("%c ",b[i][j]);
printf("\n");
}
}
return 0;
}
看到上面的字符方向图 ,这里便是字符输出的方向和值,斜杠是公共部分,横杠输出水平字符,竖杠输出垂直字符