【动态规划】 之最长公共子序列LCS

时间:2023-03-09 18:20:53
【动态规划】 之最长公共子序列LCS
int lcs_len(char *a, char *b, int c[][N]){
int aLen=strlen(a),
bLen=strlen(b),
i,j;
for(i=; i<=aLen; i++)
c[i][]=;
for(j=; j<=bLen; j++)
c[][j]=; for(i=;i<=aLen;i++)
for(j=;j<=bLen;j++)
if(a[i-]==b[j-])
c[i][j]=c[i-][j-]+;
else
c[i][j]=MAX( c[i-][j],c[i][j-] );
return c[i][j]; //返回的是长度
}
char *build_LCS(char *s,char *a, char *b){ //重建公共子序列
int k,i=strlen(a),
j=strlen(b), c[N][N]; k=lcs_len(a,b,c);
//s[k]='\0'; while(k>)
if(c[i][j]==c[i-][j])
i--;
else if(c[i][j]==c[i][j-])
j--;
else{
s[--k]=a[i-];
i--;
j--;
}
return s;
}