最长上升子序列

题目

给出一个由 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 )