Lis（最长上升子序列）nlogn求法

我们发现最原始的求lis的方法是O（n^2）的，本着最优的原则，这里介绍O（nlogn）的求法

Lis的状态转移方程是dp[i]=max(dp[i],dp[j]+1)（0<=j<i<n,a[j]<a[i]）边界条件为dp[0]=0;

这里，我们选择用d[]数组维护最大值，其中d[i]表示长度为i的Lis的最小结尾值。

解释： 1 3 2 4 5  d[2]的Lis最小结尾值应该为2，而不是为3。（因为2<3）

由d数组的定义可知，d数组中的元素一定是递增的。

这里枚举一遍初始数组（从1到n），对于每个a[i]，如果a[i]比d[len]大，则就说明，a[i]可以加到d数组len+1的位置，即此时Lis的最长长度为len+1 。

我们反过来想：我们从d数组中二分查找第一个大于a[i]的数，并且记录该数的位置p，则在p之前的数都比a[i]小，而a[i]是比d[p]小的，也就是说a[i]可以放在p-1的位置，此时的a[i]替换了d[p]的值，维护了d数组的性质。

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
int n;
int a[N],d[N];

int getLis(){
	int n;cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i];
	int len=0;d[++len]=a[1];
	for(int i=2;i<=n;i++){
		int p=lower_bound(d+1,d+1+len,a[i])-d;//在d数组中找到第一个大于a[i]的位置
		if(p==len+1) d[++len]=a[i];
		else d[p]=a[i];
	}
	return len;
}

int main(){
	cout<<getLis();
}