动态规划思想如下:
1.刻画最长公共子序列的结构特征
2.写出递归表达式
设c[i,j]为LCS的长度,i表示序列X的长度,j表示序列Y的长度
那么c[i,j]= 0 , i=0|| j=0
= c[i-1,j-1]+1; Xi=Yj
=max(c[i-1,j],c[i,j-1]), Xi !=Yi
好了说了一大堆套话,上代码。
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define MAX(A,B) ((A)>(B)? (A):(B))
char * insert(char *str,char c)
{
int n=strlen(str);
char * in_str=malloc(sizeof(char)*(n+1));
in_str[0]=c;
for(int i=1;i<n+1;i++)
in_str[i]=str[i-1];
return in_str;
}
int Longest(char * str1,char *str2)
{
str1=insert(str1,'#');
str2=insert(str2,'$');
printf("%s\n",str1);
printf("%s\n",str2);
int l1=strlen(str1);
int l2=strlen(str2);
int c[l1][l2];
memset(c,0,sizeof(c));
for(int i=1;i<l1;i++)
for(int j=1;j<l2;j++)
{
if(i==0||j==0)
{
c[i][j]=0;
}
else{
if(str1[i]==str2[j])
c[i][j]=c[i-1][j-1]+1;
else
{
c[i][j]=MAX(c[i-1][j],c[i][j-1]);
}
}
}
for(int i=0;i<l1;i++)
{
for(int j=0;j<l2;j++)
{
printf("%d ",c[i][j]);
}
printf("\n");
}
return c[l1-1][l2-1];
}
int main()
{
printf("%d\n",Longest("aaccde","aced"));
}
思路比较简单,就不解释了。