最长上升子序列, N*logN,九度OJ 1533,二分+DP

最长上升子序列, N*logN,九度OJ 1533,二分+DP

给定n个整数，求最长上升子序列的个数

例如 -1,1, 2, 2, 3 的LIS为4 （-1, 1, 2, 3）

下面算法参考了《编程之美》和别人的博客。利用了二分查找，优化了一般的O(N^2)的算法，时间复杂度O(N*logN)

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <vector>
#include <stack>
#include <queue>
#include <deque>
#include <string>
#include <set>
#include <map>
#include <vector>
#include <list>
using namespace std;
#define INF 0x7fffffff
const int maxn = 1e5 + 7;
int aryInt[maxn];
int lis[maxn]; //以ary[i]为最大元素的最长递增子序列的长度
int maxV[maxn]; //长度为lis[i]的递增序列中的最大元素值
int LIS(int ary[], int n) {
	//记录数组中的递增序列信息
	maxV[1] = ary[0]; //数组中的第一值，边界值
	maxV[0] = -INF; //最小值,边界值
	lis[0] = 1;
	int nMaxLis = 1; //数组最长递增子序列的长度
	for (int i = 1; i < n; ++i) {
		lis[i] = 1;
		//遍历历史最长递增序列信息
		//在maxV数组中，从后向前，查找第一小于ary[i]的元素，返回下标
		//lower_bound()是基于二分查找的寻找第一个不小于ary[i]的元素的下标，区别于upper_bound()
		int j = (int)(lower_bound(maxV, maxV + nMaxLis + 1, ary[i]) - maxV) - 1; //j有长度下标的意思
		lis[i] = j + 1;
		//如果当前最长序列大于最长上升序列长度，更新最长信息
		if (lis[i] > nMaxLis) {
			nMaxLis = lis[i];
			maxV[lis[i]] = ary[i];
		}
		else if (maxV[j] < ary[i] && ary[i] < maxV[j + 1]) {
			maxV[j + 1] = ary[i]; //长度为j + 1的递增子序列的最大元素尽可能小
		}
	}
	return nMaxLis;
}

int main()
{
	int n;
	while (~scanf("%d", &n)) {
		for (int i = 0; i < n; ++i)
			scanf("%d", aryInt + i);
		printf("%d\n", LIS(aryInt, n));
	}
	return 0;
}

zzuwenjie 2017-3-29 21:22:42

二次编辑（对空间进行了优化）

#define _CRT_SECURE_NO_WARNINGS  
#include <algorithm>  
#include <climits>
#include <cstdio>  
using namespace std; 
const int maxn = 1e5 + 7;
int myInts[maxn];
int dp[maxn];//dp[i]意为前序列长度为i的序列的最大元素值
int LIS(int aryInt[], int n) {
	if(n <= 0) return 0; dp[0] = INT_MIN;
	dp[1] = aryInt[0]; //初始化
	int nMaxLis = 1;
	for (int i = 1; i < n; ++i) {
		int j = (int)(lower_bound(dp, dp + nMaxLis + 1, aryInt[i]) - dp);
		if (j > nMaxLis) {
			nMaxLis = j;
			dp[j] = aryInt[i];
		}
		else if (dp[j - 1] < aryInt[i] && aryInt[i] < dp[j])
			dp[j] = aryInt[i];
	}
	return nMaxLis;
}
int main()
{
	int n;
	while (~scanf("%d", &n)) {
		for (int i = 0; i < n; ++i)
			scanf("%d", myInts + i);
		printf("%d\n", LIS(myInts, n));
	}
	return 0;
}

2017/8/2对LIS做了小修改，当参数n<= 0时，返回0。若无此语句，n=0时返回1，显然不对。