最长上升子序列模型

Vic.GoodLuck

已于 2023-09-29 18:03:53 修改

阅读量41

点赞数

分类专栏：程序设计算法——动态规划文章标签：算法数据结构

于 2023-01-02 17:28:32 首次发布

本文链接：https://blog.csdn.net/weixin_50533561/article/details/128520154

版权

程序设计算法——动态规划专栏收录该内容

12 篇文章 0 订阅

订阅专栏

怪盗基德的滑翔翼

1017. 怪盗基德的滑翔翼

以某个位置为右端点的最长上升子序列

以某个位置为左端点的最长下降子序列（倒序求最长上升子序列）

两者的最大值

int f1[N];//前边比w[i]小的 
int f2[N];//后边比w[i]小的 
int dp(){
	for(int i=1;i<=n;i++){
		f1[i]=1;
		for(int j=1;j<i;j++){
			if(w[i]>w[j]) f1[i]=max(f1[i],f1[j]+1);
		}
	}
	for(int i=n;i>=1;i--){
		f2[i]=1;
		for(int j=n;j>i;j--){
			if(w[i]>w[j]) f2[i]=max(f2[i],f2[j]+1);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f1[i]);
		maxn=max(maxn,f2[i]);
	}
	return maxn;
}

登山

1014. 登山

以某个位置为右端点的最长上升子序列

以某个位置为左端点的最长下降子序列（倒序求最长上升子序列）

两者之和的最大值-1

int dp(){
	for(int i=1;i<=n;i++){
		f1[i]=1;
		for(int j=1;j<i;j++){
			if(w[i]>w[j]) f1[i]=max(f1[i],f1[j]+1);
		}
	}
	for(int i=n;i>=1;i--){
		f2[i]=1;
		for(int j=n;j>i;j--){
			if(w[i]>w[j]) f2[i]=max(f2[i],f2[j]+1);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f1[i]+f2[i]-1);
	}
	return maxn;
}

合唱队形

482. 合唱队形 - AcWing题库

和上一题相同，不过求的是非最长上升/下降子序列中的元素个数

int dp(){
	for(int i=1;i<=n;i++){
		f1[i]=1;
		for(int j=1;j<i;j++){
			if(w[i]>w[j]) f1[i]=max(f1[i],f1[j]+1);
		}
	}
	for(int i=n;i>=1;i--){
		f2[i]=1;
		for(int j=n;j>i;j--){
			if(w[i]>w[j]) f2[i]=max(f2[i],f2[j]+1);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f1[i]+f2[i]-1);
	}
	return n-maxn;
}

友好城市

1012. 友好城市

按一边的坐标进行排序，求另一边的坐标序列的最长上升子序列

PII w[N];
int f[N];
int dp(){
	sort(w+1,w+n+1);
	for(int i=1;i<=n;i++){
		f[i]=1;
		for(int j=1;j<i;j++){
			if(w[i].second>w[j].second) f[i]=max(f[i],f[j]+1);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f[i]);
	}
	return maxn;
}

最长上升子序列（二分优化）

f[i]表示长度为i的上升子序列中上一次加入的数，对于w[i]，二分找到能继续放w[i]的最后一个子序列的位置l，if(能放) l=mid;else r=mid-1，将其放在l+1；如果找不到则新开一个队，最后的队数即为答案。举例：如果为严格单调上升序列，那么当w[i]>f[mid]时即能放

注意：初始化f[0]要保证第一个数能被放置，上升子序列则初始化为INF，下降子序列则初始化为-INF

int f[N];//长度为i的上升子序列的最后一位 
int dp(){
	int len=0;
	f[0]=-0x3f3f3f3f;
	for(int i=1;i<=n;i++){
		int l=0,r=len;
		while(l<r){
			int mid=l+r+1>>1;
			//找能在队尾放w[i]的最后一个 
			if(f[mid]<w[i]){
				l=mid;
			}else{
				r=mid-1;
			}
		}
		len=max(len,l+1);
		f[l+1]=w[i];
	}
	return len;
}

最大上升子序列和

1016. 最大上升子序列和

f[i]表示以i为右端点的上升子序列的和的最大值

int f[N];//以i为右端点的上升子序列的和 的最大值 
int dp(){
	for(int i=1;i<=n;i++){
		f[i]=w[i];
		for(int j=1;j<i;j++){
			if(w[i]>w[j]) f[i]=max(f[i],f[j]+w[i]);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f[i]);
	}
	return maxn;
}

最大上升子序列和（树状数组优化）

3662. 最大上升子序列和 - AcWing题库

数据范围比较大，朴素方法会超时，可以用树状数组来维护前缀区间（上升序列和）的最大值

首先将原数组W[]离散化并去重，得到位置数组，对于原数组每个数W[i]，在位置数组中找到比它小（此处为严格上升序列）的最大数的位置L，在树状数组查询得到[1~L]的最大值X，则以该数结尾的严格上升子序列的最大和为X+W[i]，并插入到树状数组L+1位置，更新[1,L+1]的最大值

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e5+5;
int n,m;
int w[N];
ll tr[N];
ll res;
int lowbit(int x){
	return x&-x;
}
void modify(int x,ll c){
	for(int i=x;i<=m;i+=lowbit(i)){
		tr[i]=max(tr[i],c);
	}
}
ll query(int x){
	ll res=0;
	for(int i=x;i>0;i-=lowbit(i)){
		res=max(res,tr[i]);
	}
	return res;
}
int main(){
	cin>>n;
	for(int i=1;i<=n;i++) cin>>w[i];
	int q[N];
	memcpy(q,w,sizeof w);
	sort(q+1,q+n+1);
	m=unique(q+1,q+n+1)-(q+1);
	for(int i=1;i<=n;i++){
		int l=1,r=m;
		while(l<r){
			int mid=l+r+1>>1;
			if(q[mid]<=w[i]){
				l=mid;
			}else{
				r=mid-1;
			}
		}
		ll t=query(l)+w[i];
		res=max(res,t);
		modify(l+1,t);
	}
	cout<<res;
}

拦截导弹

1010. 拦截导弹

第一问就是求最长不上升子序列

第二问求所有的数最少可以组成多少个最长不上升子序列，通过贪心求解：

对于第i个数，把它加入前 i - 1 个数构成的不上升子序列组中，结尾元素>=第i个数的数中最小的那个数（即能容纳第i个数的第一个序列，while(不能放) k++ ）

同时对于序列结尾数组的性质，当序列是不上升序列时，序列结尾数组是单调不减的（即序列结尾数组与原序列的单调性相反）

注意：

如果读入数组时不知道个数，如果使用while(cin>>w[++n])来读入，最后n等于个数+1，而不是个数

所以推荐使用while(cin>>x) w[++n]=x;的写法来读入

int dp1(){
	for(int i=1;i<=n;i++){
		f[i]=1;
		for(int j=1;j<i;j++){
			if(w[i]<=w[j]) f[i]=max(f[i],f[j]+1);
		}
	}
	int maxn=-1;
	for(int i=1;i<=n;i++){
		maxn=max(maxn,f[i]);
	}
	return maxn;
}
int dp2(){
	int cnt=0;
	for(int i=1;i<=n;i++){
		int k=0;
		//满足p[k]>=w[i]的第一个数，即能容纳w[i]的第一组数 
		while(k<cnt&&p[k]<w[i]) k++;
		if(k>=cnt) cnt++;
		p[k]=w[i];
	}
	return cnt;
}

导弹防御系统

187. 导弹防御系统 - AcWing题库

题意：最少可以组成多少个上升/下降子序列

分别开两个数组p1[ ] p2[ ]，表示所有上升/下降子序列的结尾元素

对于每个数，分别选择将它放到上升/下降子序列中，进行DFS

注意：①改变p1[ ] p2[ ]前需要保护现场，对应分支搜索完后进行恢复

②当目前解已经劣于当前最优解时，进行剪枝

int p1[N],p2[N];
//dfs(1,0,0)
void dfs(int u,int cnt1,int cnt2){
	if(cnt1+cnt2>=res) return;
	if(u==n+1){
		res=cnt1+cnt2;
		return;
	}
	//上升
	int k=0;
	while(k<cnt1&&w[u]>=p1[k]) k++;
	int t=p1[k];
	p1[k]=w[u];
	if(k>=cnt1) dfs(u+1,cnt1+1,cnt2);
	else dfs(u+1,cnt1,cnt2);
	p1[k]=t;

	//下降 
	k=0;
	while(k<cnt2&&w[u]<=p2[k]) k++;
	t=p2[k];
	p2[k]=w[u];
	if(k>=cnt2) dfs(u+1,cnt1,cnt2+1);
	else dfs(u+1,cnt1,cnt2);
	p2[k]=t;
}

最长公共上升子序列

272. 最长公共上升子序列 - AcWing题库

f[i][j] 表示a[1~i] b[1~j]中，以b[j]结尾的公共子序列的最大长度

根据子序列中是否包含a[i]来划分集合：

①包含a[i]，即 $a[i]=b[j]$ 时， $f[i][j]=max(f[i-1][k])+1$ ， $b[j]> b[k](k\in [1,j-1])$

②不包含a[i]，即 $a[i]\neq b[j]$ 时，

//f[i][j] a[1~i] b[1~j]中，以b[j]结尾的公共子序列的最大长度
//不包含a[i]，即a[i]!=b[j]时，f[i][j]=f[i-1][j]
//包含a[i]时，f[i][j]=max(f[i-1][1~j-1])+1 
int dp(){
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(a[i]!=b[j]){
				f[i][j]=f[i-1][j];
			}else{
				int maxn=0;
				for(int k=1;k<j;k++){
					if(b[j]>b[k]){
						maxn=max(maxn,f[i-1][k]);
					}
				}
				f[i][j]=max(f[i][j],maxn+1);
			}
		}
	}
	int res=-1;
	for(int i=1;i<=n;i++){
		res=max(res,f[n][i]);
	}
	return res;
}

注意：

在转移时需要反复求f[i-1][1~j-1]的最大值，直接在循环时对能够接在a[i]前面的状态的最大值进行维护即可

int dp1(){
	for(int i=1;i<=n;i++){
		int maxn=0;
		for(int j=1;j<=n;j++){
			if(a[i]!=b[j]){
				f[i][j]=f[i-1][j];
			}else{
				f[i][j]=max(f[i][j],maxn+1);
			}
			if(a[i]>b[j]){
				maxn=max(maxn,f[i-1][j]);
			}
		}
	}
	int res=-1;
	for(int i=1;i<=n;i++){
		res=max(res,f[n][i]);
	}
	return res;
}