LCS待完成

最新推荐文章于 2025-05-27 13:41:40 发布

weixin_30809333

最新推荐文章于 2025-05-27 13:41:40 发布

阅读量44

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/foreverlearn/archive/2012/04/10/2723571.html

本文详细介绍了一个经典的计算机科学问题——最长公共子序列（LCS）问题的解决方法，并提供了完整的C++实现代码。通过该算法可以找出两个字符串之间的最长公共子序列，适用于文本比对、生物信息学等领域。

待续：待研究，比较最长子序列和最长子串，研究后缀数组和后缀树

LCS实现代码：

#include <iostream>
using namespace std;

int c[8][7],b[8][7];

void lcs_length(char *X,char *Y)
{
	int m = strlen(X);
	int n = strlen(Y);
	for ( int i = 0; i <=m; i++)c[i][0] = 0;
	for(int j = 0; j <=n; j++)c[0][j] = 0;
	
	for(int i = 1; i <=m; i++)
		for (int j = 1; j <=n; j++)
		{
			int s_i = i - 1;
			int s_j = j - 1;
			if (X[s_i] == Y[s_j])
			{
					c[i][j] = c[i-1][j-1]+1;
					b[i][j] = 2;
			}
			else
			{
				if (c[i-1][j] >= c[i][j-1])
				{
					c[i][j] = c[i-1][j];
					b[i][j] = 3;
				}
				else
				{
					c[i][j] = c[i][j-1];
					b[i][j] = 1;
				}
			}
		}
}

void print_lcs(char *X,int i,int j)
{
	if (i == 0 || j == 0)return;
	if (b[i][j] == 2)
	{
		print_lcs(X,i-1,j-1);
		int s_i = i - 1;
		cout<<X[s_i];
		//cout<<"*";
	}
	else if(b[i][j] == 1) print_lcs(X,i-1,j);
	else print_lcs(X,i,j-1);
}

int main()
{
	char *s1 = "ABCBDAB";
	char *s2 = "BDCABA";
	lcs_length(s1,s2);
	print_lcs(s1,7,6);
	system("pause");
	return 0;
}

转载于:https://www.cnblogs.com/foreverlearn/archive/2012/04/10/2723571.html