Aierbude的博客

词不达意,事与愿违。

利用矩阵求两个序列的最长公共子序列

/**************最长公共子序列*****************/ 
#include<stdio.h>
#define M 7
#define N 6
void INIT_LIST(int L[][N+1],const int arg1[],const int arg2[]){
	for(int i=0;i<=M;i++){
		L[i][0]=0;
	}
	for(int i=0;i<=N;i++){
		L[0][i]=0;
	}
	for(int i=1;i<=M;i++){
		for(int j=1;j<=N;j++){
			if(arg1[i-1]==arg2[j-1]){
				L[i][j]=L[i-1][j-1]+1;
			}else{
				L[i][j]=L[i][j-1]>L[i-1][j]?L[i][j-1]:L[i-1][j];
			}
		}
	}
}

void INIT_ARGS(int arg1[],int arg2[]){
	printf("请输入第一个数组的%d个内容\n",M);
	for(int i=0;i<M;i++){
		scanf("%d",&arg1[i]);
	}
	printf("请输入第二个数组的%d个内容\n",N);
	for(int i=0;i<N;i++){
		scanf("%d",&arg2[i]);
	}
}

void SHOW_LIST(const int L[][N+1]){
	for(int i=0;i<=M;i++){
		for(int j=0;j<=N;j++){
			printf("%d ",L[i][j]);
		}
		printf("\n");
	}
}

void SHOW_ARGS(const int arg1[],const int arg2[]){
	printf("第一个数组的%d个内容是\n",M);
	for(int i=0;i<M;i++){
		printf("%d ",arg1[i]);
	}
	printf("\n");
	printf("第二个数组的%d个内容是\n",N);
	for(int i=0;i<N;i++){
		printf("%d ",arg2[i]);
	}
	printf("\n");
}

void procs(const int L[][N+1],const int arg1[],const int arg2[]){
//	printf("执行");
	int i=M-1;
	int j=N-1;
	int total=0;
	while(i>=0&&j>=0){ 
		//printf("执行i=%d,j=%d\n",i,j);
		if(arg1[i]==arg2[j]){
			printf("%d ",arg1[i]);
			i--;
			j--;
		}else{
			if(L[i][j+1]==L[i+1][j]){
				total++;
				i--;
			}else if(L[i][j+1]<L[i+1][j]){
				j--;
			}else{
				i--;
			}
		}
	}
	int summary=1;
	for(int i=1;i<=total;i++){
		summary*=2;
	} 
	printf("一共有%d种可能",summary);
}
void main(){
	int arg1[M];
	int arg2[N];
	INIT_ARGS(arg1,arg2);
	SHOW_ARGS(arg1,arg2);
	int L[M+1][N+1];
	INIT_LIST(L,arg1,arg2);
	SHOW_LIST(L);
	procs(L,arg1,arg2);
	
}
关于矩阵 其中元素如下图所示



测试数据和运行结果




阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/amanicspater/article/details/49924947
上一篇最小生成树 prim算法实现(利用图的邻接矩阵来存放图)
下一篇最长单调递增子序列O(NlogN)算法
想对作者说点什么? 我来说一句

最长公共子序列

LCS

Code_Mart Code_Mart

2017-10-19 20:30:08

阅读数:525

没有更多推荐了,返回首页

关闭
关闭