最长公共子序列

来源:互联网 发布:cg建模软件 编辑:程序博客网 时间:2024/05/02 01:58
#include
#include
int m,n;
int a[100],b[100];
int c[100][100];//存储对应的a[i]与b[j]的最大子序列
char d[100][100];//存储该位置到最大子序列走向
void lecLength(){//求出对应的c[][],d[][];
    int i,j;
   for(i=0;i<=m;i++)
c[i][0]=0;
for(j=0;j<=n;j++)
c[0][j]=0;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++){
if(a[i]==b[j]){
c[i][j]=c[i-1][j-1]+1;
d[i][j]='o';//斜向上转
}
else if(c[i-1][j]>=c[i][j-1]){
c[i][j]=c[i-1][j];
d[i][j]='u';//向上转
}
else{
c[i][j]=c[i][j-1];
d[i][j]='l';//向左转
}
}
}
int lookup(int i,int j){做备忘录,功能与lecLength()一样。
if(c[i][j]!=-1)
return c[i][j];
if(i==0||j==0){
c[i][j]=0;
return 0;
    }
if(a[i]==b[j]){
c[i][j]=lookup(i-1,j-1)+1;
d[i][j]='o';
}
else if(lookup(i-1,j)>=lookup(i,j-1)){
  c[i][j]=lookup(i-1,j);
   d[i][j]='u';
}
else{
c[i][j]=lookup(i,j-1);
d[i][j]='l';
}
return c[i][j];
}
void showLec(int i,int j){//根据d[][]输出最长公共子序列
    if(i==0||j==0){
return;
}
if(d[i][j]=='o'){
showLec(i-1,j-1);
printf("%d ",a[i]);
}
else if(d[i][j]=='u')
showLec(i-1,j);
else
showLec(i,j-1);
}
void printLec(int i,int j){//根据c[][]输出最长公共子序列
if(c[i][j]==0)
return;
if(a[i]==b[j]){
printLec(i-1,j-1);
printf("%d ",a[i]);
}
else if(c[i-1][j]>=c[i][j-1]){
printLec(i-1,j);
}
else
printLec(i,j-1);

}
int main(){
int i,j;
scanf("%d",&m);
   for(i=1;i<=m;i++)
scanf("%d",&a[i]);
scanf("%d",&n);
for(i=1;i<=n;i++)
scanf("%d",&b[i]);
   memset(c,-1,sizeof(c));
//lecLength();
    lookup(m,n);
printf("\n");
for(i=0;i<=m;i++){
for(j=0;j<=n;j++)
printf("%d ",c[i][j]);
printf("\n");
}
printf("1.最长公共子序列为:");
showLec(m,n);
printf("\n");
printf("2.最长公共子序列为:");
printLec(m,n);
printf("\n");
return 0;
}
0 0
原创粉丝点击