第四十三章 动态规划——最长单调序列模型
- 一、最长单调序列模型
- 1、模型母题
- 2、思路分析(两种方法:DP,贪心)
- 二、模型的应用
- 1、AcWing 1017. 怪盗基德的滑翔翼
- (1)问题
- (2)分析
- (3)代码
- 2、AcWing 1014. 登山
- (1)问题
- (2)分析
- (3)代码
- 3、AcWing 482. 合唱队形
- (1)问题
- (2)分析
- (3)代码
- 4、AcWing 1012. 友好城市
- (1)问题
- (2)分析
- (3)代码
- 5、AcWing 1016. 最大上升子序列和
- (1)问题
- (2)分析
- (3)代码
- 6、AcWing 1010. 拦截导弹
- (1)问题
- (2)分析
- (3)代码
一、最长单调序列模型
1、模型母题
2、思路分析(两种方法:DP,贪心)
这道母题作者在之前的文章中做过详细的解答,如果大家对这两种方法不太理解的话,作者建议先去看之前的文章。最长上升子序列模型
二、模型的应用
1、AcWing 1017. 怪盗基德的滑翔翼
(1)问题
(2)分析
这道题就是一个很单纯的单调序列的使用,大家只需要分别求出最长上升子序列和最长下降子序列,然后比较出一个最大值输出即可。
(3)代码
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
const int N=110;
int a1[N],a2[N],f1[N],f2[N];
int main()
{
int t;
cin>>t;
while(t--)
{
int n;
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%d",a1+i);
}
for(int i=n;i>0;i--)a2[i]=a1[n-i+1];
for(int i=1;i<=n;i++)f1[i]=1,f2[i]=1;
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(a1[i]<a1[j])
{
f1[i]=max(f1[i],f1[j]+1);
}
}
}
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(a2[i]<a2[j])
{
f2[i]=max(f2[i],f2[j]+1);
}
}
}
int ma=0;
for(int i=1;i<=n;i++)
{
ma=max(f1[i],ma);
}
for(int i=1;i<=n;i++)
{
ma=max(f2[i],ma);
}
cout<<ma<<endl;
}
return 0;
}
2、AcWing 1014. 登山
(1)问题
(2)分析
通过上面的图,我们只需要求出以a[i]结尾的正序的最长上升子序列长度l,以及以a[i]为结尾的逆序最长上升子序列的长度m,然后我们找到所有l+m的最大值。
(3)代码
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=1e3+10;
int a[N],ra[N];
int f1[N],f2[N],f[N];
int n;
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
scanf("%d",a+i);
for(int i=n;i>0;i--)
ra[n-i+1]=a[i];
for(int i=1;i<=n;i++)f1[i]=1,f2[i]=1;
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(a[j]<a[i])
f1[i]=max(f1[i],f1[j]+1);
}
}
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(ra[j]<ra[i])
f2[i]=max(f2[i],f2[j]+1);
}
}
int ma=0;
for(int i=1;i<=n;i++)
{
f[i]=f1[i]+f2[n-i+1]-1;
ma=max(f[i],ma);
}
cout<<ma<<endl;
return 0;
}
3、AcWing 482. 合唱队形
(1)问题
(2)分析
这道题和登山的题目几乎一样。
(3)代码
#include<iostream>
#include<algorithm>
using namespace std;
const int N=110;
int a[N],ra[N];
int f1[N],f2[N],f[N];
int main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++)scanf("%d",a+i);
for(int i=n;i>0;i--)ra[n-i+1]=a[i];
for(int i=1;i<=n;i++)f1[i]=1,f2[i]=1;
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(a[j]<a[i])
f1[i]=max(f1[i],f1[j]+1);
}
}
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(ra[j]<ra[i])
f2[i]=max(f2[i],f2[j]+1);
}
}
int ma=0;
for(int i=1;i<=n;i++)
{
f[i]=f1[i]+f2[n-i+1]-1;
ma=max(ma,f[i]);
}
cout<<n-ma<<endl;
return 0;
}
4、AcWing 1012. 友好城市
(1)问题
(2)分析
这道题也是单调序列的变形,但是这道题比之前的两道题不容易看出背后的模型,这道题作者之前已经做过详细地讲解,建议不会的读者去看一下:
传送门:
AcWing 1012. 友好城市
(3)代码
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef pair<int ,int> pii;
const int N=5100;
pii a[N];
int f[N];
int n;
bool cmp(pii a,pii b)
{
return a.first<b.first;
}
int main()
{
cin>>n;
for(int i=0;i<n;i++)scanf("%d%d",&a[i].first,&a[i].second);
sort(a,a+n,cmp);
for(int i=0;i<n;i++)f[i]=1;
for(int i=0;i<n;i++)
for(int j=0;j<i;j++)
if(a[j].second<a[i].second)
f[i]=max(f[i],f[j]+1);
int ma=0;
for(int i=0;i<n;i++)
{
ma=max(ma,f[i]);
}
cout<<ma<<endl;
return 0;
}
5、AcWing 1016. 最大上升子序列和
(1)问题
(2)分析
这道题就很简单了,只是从长度变成了最大值,只需要稍作调整即可。将转移方程中的f[j]+1变成f[j]+a[i]。然后初始化的时候初始化为数值本身,而不是1。
(3)代码
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e3+10;
int a[N],f[N];
int main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++)scanf("%d",a+i);
for(int i=1;i<=n;i++)f[i]=a[i];
for(int i=1;i<=n;i++)
{
for(int j=i-1;j>0;j--)
{
if(a[i]>a[j])
{
f[i]=max(f[i],f[j]+a[i]);
}
}
}
int ma=0;
for(int i=1;i<=n;i++)
{
ma=max(f[i],ma);
}
cout<<ma<<endl;
return 0;
}
6、AcWing 1010. 拦截导弹
(1)问题
(2)分析
这道题涉及到了最长上升子序列的贪心做法,在之前的文章中曾经做过详细地讲解,建议读者去看一看那一篇讲解。AcWing 1010. 拦截导弹
(3)代码
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e6+10;
int a[N],f[N],g[N];
int n;
int main()
{
while(cin>>a[n])n++;
int len1=0;
for(int i=0;i<n;i++)
{
int l=0,r=len1;
while(l<r)
{
int mid=l+r+1>>1;
if(f[mid]>=a[i])l=mid;
else r=mid-1;
}
len1=max(len1,r+1);
f[r+1]=a[i];
}
cout<<len1<<endl;
int len2=0;
for(int i=0;i<n;i++)
{
int l=0,r=len2;
while(l<r)
{
int mid=l+r+1>>1;
if(g[mid]<a[i])l=mid;
else r=mid-1;
}
len2=max(len2,r+1);
g[r+1]=a[i];
}
cout<<len2<<endl;
return 0;
}