最长上升子序列（优化）

描述

一个数的序列bi，当b1 < b2 < … < bS的时候，我们称这个序列是上升的。对于给定的一个序列(a1, a2, …, aN)，我们可以得到一些上升的子序列(ai1, ai2, …, aiK)，这里1 <= i1 < i2 < … < iK <= N。比如，对于序列(1, 7, 3, 5, 9, 4, 8)，有它的一些上升子序列，如(1, 7), (3, 4, 8)等等。这些子序列中最长的长度是4，比如子序列(1, 3, 5, 8).

你的任务，就是对于给定的序列，求出最长上升子序列的长度。

输入

输入的第一行是序列的长度N (1 <= N <= 1000)。第二行给出序列中的N个整数，这些整数的取值范围都在0到10000。

输出

最长上升子序列的长度。

样例输入

7
1 7 3 5 9 4 8

样例输出

4

先写一下dp普通算法

思路

f[i]以第i个数结尾的最长上升子序列的长度
每一个上升子序列包括自己，将每一个值初始化为1，之后遍历之前的值，如果该值大于之前的值，则比较是现在的值大还是之前某个值加1大。最后答案不是为f数组的最后一个数，因为最长上升子序列不一定以最后一个数为结尾。
转移方程
if(w[i]>w[j])
f[i]=max(f[j]+1,f[i])

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;

int w[1001];
int f[1001];

int main()
{
	int n,ans=0;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%d",&w[i]);
	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[j]+1,f[i]);
		}
		ans=max(ans,f[i]);
	}
	cout<<ans;
	return 0;
}

优化

普通算法的时间复杂度为O(n^2)，如果数据容量过大，肯定会超时的。
定义一个辅助数组，用于存储最长上升子序列，且其中的元素均达到能达到的最小状态。先把f[1]设为w[1]，之后遍历w数组的每个元素，判断它是否大于f数组的尾元素，如果大于，则把它也加入到f数组之中；如果不大于，则从头遍历f数组，如果该数小于其中一个数，将其替换，因为f数组中的数越小越好，其越小，后面的数才容易进入。为什么可以替换呢？比如f数组原来是1 5 7 8，需要判断的数为4,4大于1，不能替换，继续，4小于5，可以替换。需要替换的数肯定小于被替换的数，大于被替换的数的前一个数，替换后保证原数组还是上升的。这样就能做到每一个数都尽可能得小，是一种贪心的算法。最后的答案就为f数组的长度
注意
替换的判断条件是小于等于，如果小于的话，碰见一个相等的数，在此处没有替换，它肯定会比下一个数小，然后把下一个数替换掉，从而导致了序列中有了两个相等的数，影响后面所有的判断。
另外，替换掉一个记得退出循环

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;

int w[100001];
int f[100001];

int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%d",&w[i]);
	f[1]=w[1];
	int len=1;							//f数组的长度 
	for(int i=2;i<=n;i++)
	{ 
		if(w[i]>f[len])					//大于末尾的数，则加入 
 			f[++len]=w[i];						
		else
			for(int j=1;j<=len;j++)
				if(w[i]<=f[j])			//等号注意！！！一样也要替换掉，如果不换，它就会换掉后面一个，从而使得这个数与后面那个相等
				{
					f[j]=w[i];			//小于前面的数，则替换 
					break;				//！！！替换完，别忘退出循环 
				}
	} 
	cout<<len;
	return 0;
}

这么算在查找第一个小于的数时，可能还会消耗大量时间，可以使用二分查找优化。这样优化后时间复杂度为O(nlogn)

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;

int w[100001];
int f[100001];

int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%d",&w[i]);
	f[1]=w[1];
	int len=1;							//f数组的长度 
	for(int i=2;i<=n;i++)
	{ 
		if(w[i]>f[len])					//大于末尾的数，则加入 
 			f[++len]=w[i];						
		else
		{
			int l=1,r=len,mid,t=len;
			while(l<=r)
			{
				mid=l+(r-l)/2;
				if(w[i]<=f[mid])		//注意等号
				{
					t=mid;				//二分替换 
					r=mid-1;
				}
				else
					l=mid+1;
			}
			f[t]=w[i];
		}
	} 
	cout<<len;
	return 0;
}

另外，二分查找可以不用写，直接调用lower_bound()函数也可
点击查看 lower_bound( )和upper_bound( )的用法

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;

int w[100001];
int f[100001];

int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		scanf("%d",&w[i]);
	f[1]=w[1];
	int len=1;							//f数组的长度 
	for(int i=2;i<=n;i++)
	{ 
		if(w[i]>f[len])					//大于末尾的数，则加入 
 			f[++len]=w[i];						
		else
		{
			int p=lower_bound(f+1,f+len+1,w[i])-f;	//注意此处是lower  不是upper
			f[p]=w[i];					//替换 
		}
	} 
	cout<<len;
}