最长上升子序列(LIS)-O(nlogn)算法总结.

最近做了poj-1836这道题,用O(n^2)的LIS算法TLE了,百度了很久大概的学会了LIS的O(nlogn)算法;



最长上升子序列,简称LIS,O(n^2)的算法这里我就不写了,百度一下很容易理解;

1.首先定义一个数组num[n]:即为题目给出的序列;

2.其次定义一个数组LIS[len]:LIS[len]的值表示长度为len的子序列的最小末元素;

3. LIS[len]中的元素是单调递增的;

4.开始时len=1;LIS[1]=num[1];然后对num[i]:若num[i]>LIS[len],那么len++,LIS[len] = num[i];

5.否则,我们要从LIS[1]到LIS[len-1]中找到一个j,满足LIS[j-1]<num[i]<LIS[j],则根据LIS的定义,我们需要更新长度为j的上升子序列的最末元素(使之为最小的)即 LIS[j] = num[i];

6.最终答案就是len,利用LIS的单调性,在查找j的时候可以二分查找,从而时间复杂度为nlogn。


=================================================


下面举个实例:

假设存在一个序列num[9] = 2 6 8 3 5 4 10 9 3;

1.当i=1时,LIS[1]=num[1]=2,及长度为1的上升子序列的最小末元素为2,len=1;

2.当i=2时,此时先判断num[i]与LIS[len]的关系(num[2]与LIS[1]的关系),由于num[2]>LIS[1],所以len++;LIS[len]=num[i](LIS[2]=num[2]=6);即长度为2的上升子序列的最小末元素为6,len=2;LIS[1...2]=2,6;

3.当i=3时,num[3]>LIS[2],所以len++;LIS[len]=num[i](LIS[3]=num[3]=8);即长度为3的上升子序列的最小的末元素为8 ,len=3;LIS[1...3]=2,6,8;

4.当i=4时,num[4]<LIS[3],所以在LIS[1..3]=2,6,8中找到比num[4]=3大的最小的数然后将它替换成num[4]这里就是将LIS[2]=6替换成3,即LIS[2]=3,即长度为2的最小末元素为3,此时len=3,LIS[1...3]=2,3,8;

5.当i=5时,num[5]<LIS[3],在LIS[1..3]=2,3,8中找到比num[5]=5大的最小的数然后将它替换成num[5]这里就是将LIS[3]=8替换成5,即LIS[3]=5,即长度为3的最小末元素为3,此时len=3,LIS[1...3]=2,3,5;

6.当i=6时,num[6]<LIS[3],在LIS[1..3]=2,3,5中找到比num[6]=4大的最小的数然后将它替换成num[6]这里就是将LIS[3]=5替换成4,即LIS[3]=4,即长度为3的最小末元素为4,此时len=3,LIS[1...3]=2,3,4;

7.当i=7时,num[7]>LIS[3],所以len++;LIS[len]=num[7]=10;即长度为4的子序列的最小末元素为10,len=4;LIS[1..4]=2,3,4,10;

8.当i=8时,num[8]<LIS[4],在LIS[1..4]=2,3,4,10中找到....即将LIS[4]替换成8,此时len=4;LIS[1..4]=2,3,4,8;

9.当i=9时,num[9]<LIS[4],在LIS[1..4]=2,3,4,8中找到比num[9]=3大的最小的末元素为然后将它替换成num[9]这里就是将LIS[3]=4替换成3,即LIS=[3],即此时长度为3的子序列的最小末元素更新为3;此时len=4;LIS[1..4]=2,3,3,8;

此时说明这个数列num的最长上升子序列的长度为4


注意LIS表示的并不是那个最长的子序列,它只是存储长度为i的子序列的最小末元素,最后一步看似没有用处,实际上假如给num添加两个元素3,5那么可以将d[4]更新为3 d[5]更新为5这样长度就变成了5;

在B中插入数据是有序的,而且是进行替换而不需要挪动——也就是说,我们可以使用二分查找,将每一个数字的插入时间优化到O(logN)


代码如下:

#include<cstdio>
#include<cstring>
#define MAXN 40005

int arr[MAXN],ans[MAXN],len;

int binary_search(int i){
	int left,right,mid;
	left=0,right=len;
	while(left<right){
		mid = left+(right-left)/2;
		if(ans[mid]>=arr[i]) right=mid;
		else left=mid+1;
	}
	return left;
}

int main()
{
	int T,p,i,j,k;
	scanf("%d",&T);
	while(T--)
        {
		scanf("%d",&p);
		for(i=1; i<=p; ++i)
			scanf("%d",&arr[i]);

		ans[1] = arr[1];
		len=1;
		for(i=2; i<=p; ++i){
			if(arr[i]>ans[len])
				ans[++len]=arr[i];
			else{
				int pos=binary_search(i);
				ans[pos] = arr[i];
                               }
		}
		printf("%d\n",len);
	}
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值