【算法基础课】线性DP-CSDN博客

本文链接：https://blog.csdn.net/weixin_51948235/article/details/119411078

本文介绍了动态规划在解决最长上升子序列、最长公共子序列和字符串编辑距离问题上的应用。通过ACW895和ACW899题目，详细阐述了线性DP的解题思路，包括如何定义状态转移方程，并给出了C++实现代码。此外，还探讨了最大子序列和问题，如最大子序列和与最大子矩阵和的计算方法。

摘要由CSDN通过智能技术生成

文章目录

动态规划 —— 线性 DP

最长上升子序列

Acw895-最长上升子序列

$d p [i]$ 定义为以 $a [i]$ 为结尾的最长上升子序列的长度/总和

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int N=1e3+10;
int a[N],dp[N];
int main()
{
	int n,res=0; cin>>n;
	for(int i=1;i<=n;i++)
	{
	    cin>>a[i]; dp[i]=1;
	    for(int j=1;j<i;j++)
	        if(a[j]<a[i]) dp[i]=max(dp[i],dp[j]+1);//max(dp[i],dp[j]+a[i])
	    res=max(res,dp[i]);
	}
	cout<<res;
	return 0;
}

最长公共子序列

Acw899-最长公共子序列

定义 $d p [i, j]$ 为 $a [i]$ , $b [j]$ 之前的最长公共子序列长度。
转移时，对于 $a [i]$ , $b [j]$ ，选择情况为 00，01，10，11（1代表选其为结尾）

当 $a [i] = b [j]$ 时，11才不为空集。
由定义可知，00空间是01空间的子集，00空间也是10空间的子集.即01，10空间交集为00空间。
求 $m a x$ 的话，00就可以不考虑，只考虑01，10，11.

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int N=1010;
int dp[N][N];
char a[N],b[N];
int main()
{
	int n,m;
	cin>>n>>m;
	cin>>a+1>>b+1;
	
	for(int i=1;i<=n;i++)
	    for(int j=1;j<=m;j++)
	    {
	        dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
	        if(a[i]==b[j]) dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
	    }
	   
	cout<<dp[n][m];
	return 0;
}

字符串距离

动态规划 —— 线性 DP —— 字符串编辑距离

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=2100;

char a[N],b[N];
int dp[N][N];

int main()
{
	cin>>a+1>>b+1;
	
	for(int i=1;i<=strlen(a+1);i++) dp[i][0]=i;
	for(int j=1;j<=strlen(b+1);j++) dp[0][j]=j;
	
	for(int i=1;i<=strlen(a+1);i++)
		for(int j=1;j<=strlen(b+1);j++)
		{
			if(a[i]==b[j]) dp[i][j]=dp[i-1][j-1];
			else dp[i][j]=min(dp[i-1][j],min(dp[i][j-1],dp[i-1][j-1]))+1;
		}
		
	cout<<dp[strlen(a+1)][strlen(b+1)];
	
	return 0;
}