最长公共子序列

本文探讨了如何使用C/C++实现 longest common subsequence (LCS) 算法,包括基本思路和递归实现,以及一个应用实例。同时涉及了子序列和问题,通过动态规划找到以特定元素结尾的最长上升子序列的和。适合理解字符串算法和动态规划的开发者。
摘要由CSDN通过智能技术生成

这里搬过来老师的ppt
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

参考程序

#include <bits/stdc++.h>
using namespace std;

/*
x[1..m], y[1..n]

x[m]==y[n] the result is (x[1..m-1], y[1..n-1]) + 1
x[m]!=y[n] the max of the two (x[1..m-1], y[1..n]), (x[1..m],y[1..n-1])

*/

int c[10][10];

int lcs(char x[], int m, char y[], int n){
	int t;
	if(c[m][n]!=-1)
		return c[m][n];	

	if(x[m-1]==y[n-1]){
		t=lcs(x,m-1,y,n-1)+1;	
		c[m][n]=t;
		return t;
	}	
	int t1,t2;
	t1=lcs(x,m,y,n-1);
	t2=lcs(x,m-1,y,n);
	
	t=t1>t2?t1:t2;
	c[m][n]=t;
	return t;
}

int main()
{
	char x[]="ABCBDAB", y[]="BDCABA";
	int i,j;
	for(i=0;i<10;i++)
		for(j=0;j<10;j++)
			c[i][j]=-1;
	
	for(i=0;i<10;i++)	c[i][0]=0;
	for(j=0;j<10;j++)	c[0][j]=0;
	cout<<lcs(x,strlen(x), y, strlen(y));
	
	/*for(i=1;i<=strlen(x);i++){
		for(j=1;j<=strlen(y);j++){
			if(x[i]==y[j])
				c[i][j]=c[i-1][j-1]+1;
			else{
				c[i][j]=c[i-1][j]>c[i][j-1]?c[i-1][j]:c[i][j-1];
			}
		}
	}*/
	//cout<<c[strlen(x)][strlen(y)];
	return 0;
}

例题:
在这里插入图片描述

#include <bits/stdc++.h>
using namespace std;
int c[1001][1001];
char x[100010],y[100010];
int f(char x[], int m, char y[], int n) {
	int t;
	if(c[m][n]!=-1)
		return c[m][n];
	if(x[m-1]==y[n-1]) {
		t=f(x,m-1,y,n-1)+1;
		c[m][n]=t;
		return t;
	}
	int t1,t2;
	t1=f(x,m,y,n-1);
	t2=f(x,m-1,y,n);

	t=t1>t2?t1:t2;
	c[m][n]=t;
	return t;
}

int main() {
	gets(x);
	gets(y);
	int len1=strlen(x);
	int len2=strlen(y);
	int i,j;
	int len=max(len1,len2);
	for(i=0; i<200; i++)
		for(j=0; j<200; j++)
			c[i][j]=-1;

	for(i=0; i<200; i++)	c[i][0]=0;
	for(j=0; j<200; j++)	c[0][j]=0;
	cout<<f(x,strlen(x), y, strlen(y));
	return 0;
}

另一个关于子序列的题,直接在这里补充了
在这里插入图片描述
比较简单的一种算法

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5;
int a[N],sum[N];//sum[i]记录以i结尾的最长上升子序列的和
int main() {
	int n;
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i];
	}
	int ans=0;
	for(int i=1;i<=n;i++){
		sum[i]=a[i];//最短时是他本身
		for(int j=1;j<i;j++){
			if(a[j]<a[i]&&sum[j]+a[i]>sum[i])
				sum[i]=sum[j]+a[i];
		}
		ans=max(ans,sum[i]);
	}
	cout<<ans;
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值