动态规划之最长公共子序列

具体分析见算法导论,直接上代码了。

#include <iostream>
using namespace std;


int LCS_Length(string X, string Y, int** C, int XLen, int YLen)
{
	int m = XLen, n = YLen;
	for (int i = 1; i <= m; i++)
		C[i][0] = 0;
	for (int i = 1; i <= n; i++)
		C[0][i] = 0;
	C[0][0] = 0;

	for (int i = 1; i <= m; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			if (X[i-1] == Y[j-1])
			{
				C[i][j] = C[i-1][j-1]+1;
			}
			else
			{
				if (C[i-1][j] >= C[i][j-1])
					C[i][j] = C[i-1][j];
				else
					C[i][j] = C[i][j-1];
			}
		}
	}
/*
	for (int i = 0; i <= m; i++)
	{
		for (int j = 0; j <= n; j++)
		{
			cout<<C[i][j]<<" ";
		}
		cout<<endl;
	}*/
	//调试bug,打印输出,将X[i]==Y[j]改成X[i-1][j-1]。
	return C[m][n];
}

void print_LCS(int** C, int i, int j, string X)
{
	if (i == 0 || j == 0)
		return;

	if (C[i][j] == C[i-1][j])
	{
		print_LCS(C, i-1, j, X);
		
	}
	else if (C[i][j] == C[i][j-1])
	{
		print_LCS(C, i, j-1, X);
	}
	else
	{
		print_LCS(C, i-1, j-1, X);
		cout<<X[i-1];
	}
}

int main()
{
	string X = "ABCBDAB";
	string Y = "BDCABA";
	int XLen = X.length();
	int YLen = Y.length();
	int **C;
	C = new int *[XLen+1];
	for (int i = 0; i <= XLen; i++)
	{
		C[i] = new int[YLen+1];
	}
	cout<<LCS_Length(X,Y,C,XLen,YLen)<<endl;
	print_LCS(C, XLen, YLen, X);
	for (int i = 0; i <= XLen; i++)
	{
		delete[]C[i];
	}
	delete[]C;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值