不好做的最长上升子序列 O(nlogn)树状数组解法

题目

实际上本题和普通的最长上升子序列的题目一样，只是测试的数据更强
描述
一个数的序列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 <= 300000)。第二行给出序列中的N个整数，这些整数的取值范围都在0到100000000之间。
输出
最长上升子序列的长度。
样例输入
7
1 7 3 5 9 4 8
样例输出
4
来源
test cases by Gong Linyuan

思路与分析

贴上课的PPT



主要的思路和做法上面的截图已经说得很清楚了，但说实话第一眼看到觉得不太好理解这种做法为什么能够得到正确的结果，我是对照着样例的运行过程思考有段时间才搞明白，我把对每个函数的解释贴在程序中的注释里了，可以看着运行过程理解一下

代码

#include<iostream>
#include<algorithm>
using namespace std;
#pragma warning(disable:4996)
class number {
public:
	int v;
	int posi;
};
int lowbit(int x) {
	return x & (-x);
}
int N;
number Num[300010];
int C[300010] = { 0 };
bool great(number l, number r) {
	if (l.v != r.v)return l.v < r.v;
	return l.posi > r.posi;
}/*对于值相同的元素，我们要求posi大的排在num数组的前面，这样先更新后面的值
	 不会对前面值的查询造成影响，在更新时通过max比较能够消去前值对后值的影响*/
int query(int k) {
	int res = 0;
	while (k>0)
	{
		if (C[k] > res)
			res = C[k];
		k -= lowbit(k);
	}
	return res;
}/*由于list数组中list[i]表示终点在a[i-lowbit(i)+1]到a[i]的lis，
	 query要查询从a[1]到a[i]的lis，所以为得到终点在a[1]~a[i]的lis的最大值，
	 需要对list[]中的对应元素做一遍检索，最好的设计就是检索完p后检索p-lowbit(p),
	 这样可以查询最少的元素并且使list[i]表示的lis终点区间满足要求
	*/
void update(int k,int v) {
	while (k<=N)
	{
		if (C[k] < v)C[k] = v;
		k += lowbit(k);
	}
}/*这波是单点更新，已经对值做过排序，所以在num数组中每向后检索一位，都必定会使
 涉及该区间的lis增长一位，即树状数组中的单点更新，通过与之前记录的LIS长度进行比较
 决定是否要更新LIS的值
 */
int main() {

	cin >> N;
	for (int i = 0; i < N; i++)
	{
		scanf("%d", &Num[i].v);
		Num[i].posi = i;
	}
	sort(Num, Num+N, great);
	int ans = 0;
	for (int i = 0; i < N; i++)
	{
		int temp = 0;
		temp = query(Num[i].posi+1) + 1;
		ans = max(temp, ans);
		update(Num[i].posi+1, temp);
	}
		//最后得到的C数组的最大值满足要求
	cout << ans << endl;
	return 0;
}