最长公共子序列(LCS)

问题:给定两个字符串(或者数字序列)A和B,求一个字符串,使得这个字符串是A和B的最长公共部分(子序列可以不连续)。

比如: s a d s t o r y 和a d m i n s o r r y的最长公共子序列为adsory,长度为6。

思路:

          暴力解法的复杂度较高,这里达不到满意效果。这里令dp[i][j]为字符串A的i号位,和字符串B的j号位之前的LCS长度(下标从1开始),如dp[4][5]表示“sads”和“admin”的LCS长度。那么可以根据A[i]和B[j]的情况,分为两种决策:

  1. 若A[i]!=B[j],则字符串A的第i号位和字符串B的j号位之前的LCS无法延长,因此dp[i][j]将会继承dp[i-1[j]与dp[i][j-1]中的较大值,例如样例中dp[3][3]表示“sad”与“adm”的LCS长度,由于’d‘不等于’m‘,这样dp[3][3]无法在原先的基础上延长,因此继承自”sa”与“adm"的LCS、"sad"与”ad“的LCS较大值,这里为2。
  2. 若A[i]==B[j],则字符串A与字符串B的LCS增加了一位,即有dp[i][j]=dp[i-1][j-1]+1。

这样就得到了状态转移方程:

A[i]!=B[j]:dp[i][j]=max(dp[i-1[j],dp[i][j-1]);

A[i]==B[j]:dp[i][j]=dp[i-1][j-1]+1

代码:

#include <iostream>
#include<cstring>
#include<algorithm>
using namespace std;

const int N = 100;
char A[N], B[N];
int dp[N][N];
int main()
{
	int n;
	gets(A + 1);//从下标1开始读取
	gets(B + 1);
	int lenA = strlen(A + 1);//由于下标从1开始,因此读取长度也从+1开始
	int lenB = strlen(B + 1);//在PAT的编译器里面不能用gets,那么可以直接scanf,然后建立两个新的数组,从新数组的下标1开始分别复制两个字符串。
	//边界
	for (int i = 0; i <= lenA; i++)
	{
		dp[i][0] = 0;
	}
	for (int j = 0; j <= lenB; j++)
	{
		dp[0][j] = 0;
	}
	//状态转移方程
	for (int i = 0; i <= lenA; i++)
	{
		dp[i][0] = 0;
	}
	for (int j = 0; j <= lenB; j++)
	{
		dp[0][j] = 0;
	}
	//状态转移方程
	for (int i = 1; i <= lenA; i++)
	{
		for (int j = 1; j <= lenB; j++)
		{
			if (A[i] == B[i])
			{
				dp[i][j] = dp[i - 1][j - 1] + 1;
			}
			else {
				dp[i][j] = max(dp[i][j - 1], dp[i - 1][j]);
			}
		}
	}
	printf_s("%d\n", dp[lenA][lenB]);



	

}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值