最长公共子序列

3 篇文章 0 订阅
#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
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值