问题描述:
最长公共子序列(LCS)是一个在一个序列集合中(通常为两个序列)用
来查找所有序列中最长子序列的问题。这与查找最长公共子串的问题同
的地方是:子序列不需要在原序列中占用连续的位置。而最长公共子串
(要求连续)和最长公共子序列是不同的。
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
int longestSequence[50][50] = {0};//最长子序列表
int main(){
puts("--------------动态规划-----最长公共子序列问题----------------\n");
int *sequence1,*sequence2;
int length1,length2;
puts("输入序列1的长度:");
scanf("%d",&length1);
puts("输入序列1每个数字:");
sequence1 = (int *)malloc(sizeof(int)*(length1+1));
int i = 1;
for(i; i <= length1;i++){
scanf("%d",&sequence1[i]);
}
puts("输入序列2的长度:");
scanf("%d",&length2);
puts("输入序列1每个数字:");
sequence2 = (int *)malloc(sizeof(int)*(length2+1));
i = 1;
for(i; i <= length2;i++){
scanf("%d",&sequence2[i]);
}
i = 0;
dynamicPrograming(sequence1,length1,sequence2,length2);
for(i;i<=length1;i++){
int j=0;
for(j;j<=length2;j++){
printf("f[%d][%d]=%d\n",i,j,longestSequence[i][j]);
}
}
printf("序列1和序列2的最长公共子序列长度为:%d\n",longestSequence[length1][length2]);
return 0;
}
void dynamicPrograming( int *sequence1,int length1,int *sequence2,int length2){
//思想主要是两个序列的最后一个字符相等的话,那么问题规模从f(m,n)降为f(m-1,n-1);f(m,n)=f(m-1,n-1)+1,f(m,n)代表两个子序列的最长长度
//否者不相等那么问题f(m,n) = max{f(m,n-1),f(m-1,n)}
//先处理好所有只有f[m][n]中有一个长度为1的序列的的问题,这时候需要比较长度为1的那个序列是否在长度>1的序列存在
int sequence2Count = 1;
int sequence1Count = 1;
for(sequence2Count;sequence2Count<=length2;sequence2Count++){
int i = 1;
//遍历寻找是sequence1[1]否存在当前的长度sequence2中
for(i;i<=sequence2Count;i++){
if(sequence1[1]==sequence2[i]){
longestSequence[1][sequence2Count] = 1;
break;
}
}
}
for(sequence1Count;sequence1Count<=length1;sequence1Count++){
int i = 1;
//遍历寻找是sequence1[1]否存在当前的长度sequence2中
for(i;i<=sequence1Count;i++){
if(sequence2[1]==sequence1[i]){
longestSequence[sequence1Count][1] = 1;
break;
}
}
}
sequence1Count = 2;
for(sequence1Count;sequence1Count <= length1;sequence1Count++){
sequence2Count = 2;
for(sequence2Count; sequence2Count <= length2;sequence2Count++){
if(sequence1[sequence1Count] == sequence2[sequence2Count]){//如果当前的两个序列最后两个数字相等
longestSequence[sequence1Count][sequence2Count] = longestSequence[sequence1Count-1][sequence2Count-1] + 1;
} else {
//有可能存在f[1][2]但是最后两个数字都不相等那么就会让f[1][2] = f[0][1];这时候就等于0了 这是一个bug
int f1 = longestSequence[sequence1Count-1][sequence2Count];
int f2 = longestSequence[sequence1Count][sequence2Count-1];
longestSequence[sequence1Count][sequence2Count] = f1>f2?f1:f2;
}
}
}
}
实现最长公共子序列问题
最新推荐文章于 2022-09-18 23:43:56 发布