我正在尝试在 c 中实现最长公共子序列算法,矩阵 c[][] 存储最长公共子序列的长度,row[][] 存储 c[][] 矩阵中的父块行col[][] 存储父块列。
对于解决 LCS 的方法非常不方便且效率低下,我深表歉意,但是 什么都没有打印。请帮忙。
#include<stdio.h>
#include<stdlib.h>
void lcs(char *a,char *b,int i,int j)
{
int x,y,z;
int **c=(int **)malloc((j+1)*sizeof(int *));
for(x=0;x<j+1;x++)
{
*(c+x)=(int *)malloc((i+1)*sizeof(int));
*(*(c+x)+0)=0;
}
for(y=0;y<i+1;i++)
c[0][y]=0;
int **row=(int **)malloc((j+1)*sizeof(int *));
for(x=0;x<j+1;x++)
{
*(row+x)=(int *)malloc((i+1)*sizeof(int));
*(*(row+x)+0)=x-1;
}
for(y=0;y<i+1;y++)
row[0][y]=0;
row[0][0]=0;
int **col=(int **)malloc((j+1)*sizeof(int *));
for(x=0;x<j+1;x++)
{
*(col+x)=(int *)malloc((i+1)*sizeof(int));
*(*(col+x)+0)=0;
}
for(y=0;y<i+1;y++)
col[0][y]=y-1;
col[0][0]=0;
for(x=1;x<j+1;x++)
{
for(y=1;y<i+1;y++)
{
if(a[y-1]==b[x-1])
{
c[x][y]=c[x-1][y-1]+1;
row[x][y]=x-1;
col[x][y]=y-1;
}
else if(c[x-1][y]>
c[x][y-1])
{
c[x][y]=c[x-1][y];
row[x][y]=x-1;
col[x][y]=y;
}
else
{
c[x][y]=c[x][y-1];
row[x][y]=x;
col[x][y]=y-1;
}
}
}
for(x=0;x<j+1;x++)
{
for(y=0;y<i+1;y++)
printf("%d ",c[x][y]);
printf("\n");
}
printf("\n\n");
for(x=0;x<j+1;x++)
{
for(y=0;y<i+1;y++)
printf("%d,%d ",row[x][y],col[x][y]);
printf("\n");
}
}
void main()
{
char a[30]="papayaman";
char b[30]="papmn";
lcs(a,b,9,5);
}
答案 0 :(得分:0)
也许可以简化成这样:
迭代版本(自下而上)
int lcs(char * A, char * B)
{
allocate storage for array L;
for (i = m; i >= 0; i--)
for (j = n; j >= 0; j--)
{
if (A[i] == '\0' || B[j] == '\0') L[i,j] = 0;
else if (A[i] == B[j]) L[i,j] = 1 + L[i+1, j+1];
else L[i,j] = max(L[i+1, j], L[i, j+1]);
}
return L[0,0];
}
====================================================