最长公共子序列（Lcs）51NOD-1006

最新推荐文章于 2020-07-21 23:30:12 发布

原创最新推荐文章于 2020-07-21 23:30:12 发布 · 530 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#51NOD 1006 #最长公共子序列 #Lcs

数据结构同时被 2 个专栏收录

3 篇文章

订阅专栏

51NOD

1 篇文章

订阅专栏

本文介绍了最长公共子序列的概念及其与最长公共子串的区别，并通过一个具体的编程实例，使用动态规划方法解决最长公共子序列问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一、什么是最长公共子序列

什么是最长公共子序列呢？举个简单的例子吧，一个数列S，若分别是两个或多个已知序列的子序列，且是所有符合条件序列中最长的，则S称为已知序列的最长公共子序列。

举例如下，如：有两个随机数列，1 2 3 4 5 6 和 3 4 5 8 9，则它们的最长公共子序列便是：3 4 5。

一直不明白：最长公共子串和最长公共子序列的区别。

最长公共子串（Longest Common Substirng）和最长公共子序列（Longest Common Subsequence，LCS）的区别为：子串是串的一个连续的部分，子序列则是从不改变序列的顺序，而从序列中去掉任意的元素而获得新的序列；也就是说，子串中字符的位置必须是连续的，子序列则可以不必连续。

二、题目分析

Input
第1行：字符串A
第2行：字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列，如果有多个，随意输出1个。
Input示例
abcicba
abdkscab
Output示例
abca

分析：

此题的切入点就是动态规划，通过动归来确定哪些字符是最长公共子序列中的字符。

mat[i][j] 表示第一个序列的前i个字符和第二个序列的前j个字符的公共子序列。动态转移方程为：

mat[i][j] = max(mat[i-1][j], mat[i][j-1], mat[i-1][j-1] + (A[i]==B[j] ? 1 : 0));

代码如下：

#include<stdio.h>
#include<iostream>
using namespace std;
#define MAXN 1002
char a[MAXN] = { 0 };
char b[MAXN] = { 0 };
char r[MAXN] = { 0 };
short mat[MAXN][MAXN] = { 0 };

short Max1(short a,short b,short c)
{
	if (a > b)
	{
		b = a;
	}
	return b > c ? b : c;
}

int main()
{
	int i, j = 0, k;
	scanf("%s %s", a + 1, b + 1); //输入两字符串 从下标1开始存储入a[] b[]
	for (i = 1; a[i]; i++)
	{
		for (j = 1; b[j]; j++)
		{
			mat[i][j] = Max1(mat[i - 1][j], mat[i][j - 1], mat[i-1][j-1] + (a[i] == b[j] ? 1 : 0));
		}
	} //动态规划
	i--;
	j--;
	k = MAXN - 1;
	while (i > 0 && j > 0)
	{
		if (a[i] == b[j])
		{
			r[k--] = a[i];
			i--;
			j--;
		}
		else if (mat[i][j - 1] < mat[i - 1][j])
		{ i--; }
		else 
		{ j--; }
	} //反向查找最长公共子序列
	printf("%s\n",r+k+1);
	
	system("pause");
	return 0;
}