最长公共子序列

最新推荐文章于 2023-06-01 19:25:01 发布

jinghuaboy

最新推荐文章于 2023-06-01 19:25:01 发布

阅读量312

点赞数

分类专栏： Algorithms Implementation

本文链接：https://blog.csdn.net/jinghuaboy/article/details/39186047

版权

Algorithms Implementation 专栏收录该内容

21 篇文章 0 订阅

订阅专栏

问题描述：

设X=<x₁,x₂,…,x_n>，Y=<y₁,y₂,…,y_m>分别是长度为n和m的序列，X和Y的最长公共子序列是这样一个子序列lcs=<x_k1,x_k2,…,x_kp>，其中1<=k1<k2<…<kp<=n且有x_k1=y_s₁,x_k2=y_s2,…,x_kp=y_sp，1<=s1<s2<...<sp<=m。求p的最大值。

解决方法：

动态规划方法：

f(i,j)=0, i=0 or j=0,

f(i,j)=f(i-1,j-1), xi=yj and i,j>0

f(i,j)=max{f(i-1,j),f(i,j-1)}, xi!=yj and i,j>0

算法实现：

#include <iostream>
#include <cstring>

template<typename T> int LongestCommonSequence(T *, T *, const int size1, const int size2);

using namespace std;
int main(void)
{
	const int SIZE = 30;
	char arr1[SIZE]="asdfsdfasdfsdf";
	char arr2[SIZE]="dfasdfdfasdfdf";
	int len_max = 0;
	
	len_max = LongestCommonSequence(arr1, arr2, strlen(arr1), strlen(arr2));

	cout << len_max << endl;
	return 0;
}

template<typename T> int LongestCommonSequence(T *arr1, T *arr2, const int size1, const int size2)
{
	if (size1 < 0 || size2 < 0)
	{
		return 0;
	}
	int len_max;
	int *len = new int[size1+1];
	for (int i = 0; i <= size1; ++i)
	{
		len[i] = 0;
	}
	for (int i = 0; i < size2; ++i)
	{
		for (int j = 0, tmp1, tmp2 = 0; j < size1 ; ++j)
		{
			tmp1 = tmp2;
			tmp2 = len[j+1];//save f(i-1,j-1) temporarily
			if (arr2[i] == arr1[j])
			{
				len[j+1]=tmp1+1;
			}
			else
			{
				len[j+1]=(len[j+1]>=len[j])?len[j+1]:len[j];
			}
		}
	}
	len_max = len[size1];
	
	delete[] len;
	return len_max;
}