题目
给出一个由 n 个不超过 10^6 的正整数组成的序列。请输出这个序列的最长上升子序列的长度。
最长上升子序列是指,从原序列中按顺序取出一些数字排在一起,这些数字是逐渐增大的。
题目链接https://www.luogu.com.cn/problem/B3637
方法一:dp
dp基本思路
定义以第 i 个数结尾的最长上升子序列为dp[ i ]
当第 j 个数小于第 i 个数时有状态转移方程式如下
dp[ i ]=max(dp[ i ],dp[ j ]+1)
即当第 i 个数大于第 j 个数时,以第 i 个数结尾的最长上升子序列可以为以第 j 个数结尾的最长上升子序列加上一
举个例子
a[ i ] = { 1,2,4,1,3,4}
算dp[ 1 ]时,只有一个数,所以dp[ 1 ]=1
算dp[ 2 ]时,因为a[ 2 ]>a[ 1 ],所以dp[ 2 ]=dp[ 1 ]+1=2
算dp[ 3 ]时,因为a[ 3 ]>a[ 1 ]且a[ 3 ]>a[ 2 ],所以dp[ 3 ]=max(dp[ 1 ],dp[ 2 ])+1=3
算dp[ 4 ]时,因为没有比它小的数,所以dp[ 4 ]=1
算dp[ 5 ]时,因为a[ 5 ]>a[ 1 ]且a[ 5 ]>a[ 2 ],所以dp[ 5 ]=max(dp[ 1 ],dp[ 2 ])+1=3
算dp[ 6 ]时,因为a[ 6 ]>a[ 1 ]且a[ 6 ]>a[ 2 ]且a[ 6 ]>a[ 5 ],所以dp[ 6 ]=max(dp[ 1 ],dp[ 2 ], dp[ 5 ])+1=4
最终答案就是上面的数的最大值,也就是4
dp具体代码
#include <iostream>
using namespace std;
int dp[6100],map[6100];
int main()
{
int n,maxn=1;
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>map[i];
dp[i]=1;
}
dp[1]=1;
for(int i=2;i<=n;i++)
{
int f=1;
for(int j=1;j<i;j++)
{
if(map[i]>map[j])
{
dp[i]=max(dp[i],dp[j]+1);//找以第i个数结尾的最长上升子序列
f=0;
}
}
if(f) dp[i]=1;
if(dp[i]>maxn)
{
maxn=dp[i];//找最长上升子序列
}
}
cout<<maxn;
return 0;
}
此算法的时间复杂度为 O ( n^2 ) ,在条件苛刻时容易出现TLE,不是最优解
方法二:贪心+二分
贪心+二分基本思路
对于一个上升子序列,显然其结尾元素越小,越有利于在后面接入其他的元素,也就越可能变得更长。
设 a[ i ] 为第 i 个数,f[ i ] 长 i 的上升子序列长度的最小值
对于每个 a[ i ] ,如果 f[ 当前最长上升子序列长度 ] < a[ i ] ,就将 a[ i ] 接到 f[ 当前最长上升子序列长度 ] 后;否则用 a[ i ] 替换 f 中的第一个大于等于a [ i ]的元素
因为 a[ i ] 加入 f 时,a[ i ]一定大于上一个值,所以 f 一定单调递增
再举个例子
a[ i ] = { 1,2,4,1,3,4}
a[ 1 ]直接放,f[ 1 ]=a[ 1 ]=1
注:红色点为当前点,黑色点为之前的点
a[ 2 ]>f[ 1 ],f[ 2 ]=a[ 2 ]=2
a[ 3 ]>f[ 2 ],f[ 3 ]=a[ 3 ]=4
a[ 4 ]<f[ 2 ],f[ 2 ]=a[ 4 ]=1
a[ 5 ]<f[ 3 ],f[ 3 ]=a[ 5 ]=3
a[ 6 ]>f[ 3 ],f[ 4 ]=a[ 6 ]=4
所以答案就是 f 的最大下标,也就是4
时间复杂度改进
考虑到时间复杂度,在找最小元素时需要使用二分(不然改不改没什么区别还容易错)
二分代码:
int l=0,r=len,mid;
while(l<r){
mid=(l+r+1)/2;
if(a[ i ]>f[ mid ]) l=mid;
else r=mid-1;
}
len=max(len,r+1);
f[ r+1 ]=a[ i ];
贪心+二分具体代码
#include <bits/stdc++.h>
using namespace std;
const int N=5e5+5;
int f[N],a[N];
int main()
{
int n,len=0;
cin>>n;
for(int i=1;i<=n;i++) cin>>a[i];
memset(f,0x3f,sizeof(f));
for(int i=1;i<=n;i++)
{
int l=0,r=len,mid;
while(l<r){//用二分找最小值
mid=(l+r+1)/2;
if(a[i]>f[mid]) l=mid;
else r=mid-1;
}
len=max(len,r+1);//当r+1比len大时,替换len
f[r+1]=a[i];
}
cout<<len;
return 0;
}
时间复杂度 O ( nlogn )