NLOG 最长不下降子序列

LIS有N2的算法,还有NLOGN的算法,这里说一下NLOGN的算法。

我们设S[]为原始序列,F[i]为以S[i]结束的LIS的最大长度。在求LIS时,对于S[i],我们找一个最大的可能的F[j],(j<i)使S[i]尽可能大,在寻找一个F[i]时,我们可以用二分查找优化。

我们要通过枚举F[1...i-1]找到可作为F[i]前缀的LIS的最大长度,换一种思路,我们就可以枚举可能的长度——用二分——再检查它是否可以连接在S[i]的前面。

于是我们用C[i]记录长度为i的LIS的末位的最小值(这里有点贪心的想法),然后在处理S[i]时,通过二分,查找末尾数小于S[i]的LIS的最大长度,然后更新C[i]数组,最终通过查询C[i]数组得到结果。

编码的难点是二分查找,需要仔细体会。

对于二分:

A. 如果一个有序表中有元素X, 则通过二分查找必能找到X

B. 在二分中,我用v记录可行的一个LIS长度(即存在长度为v的可以连接在S[i]前面的LIS),然后查找是否有更大的v,于是要让mid=low+1,此时若存在更好的v一定可以找到。

/*
 * TYVJ-1049 最长不下降子序列
 * mike-w
 * 2012-10-16
 */
#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define MAXSIZE 5555
#define INF 0x3ffffff
#define xLOL_DEBUG

int s[MAXSIZE];
int c[MAXSIZE];
int n;

int min(int e1, int e2)
{
	return e1<e2?e1:e2;
}

int main(void)
{
	int i, j, high, low, mid, v;

	scanf("%d", &n);
	for(i=1; i<=n; i++)
		scanf("%d", s+i);
	
	for(i=1; i<=n; i++)
		c[i]=INF;
	c[1]=s[1];
	for(i=2; i<=n; i++)
	{
		high=i-1, low=1, v=-1;
		while(low<=high)
		{
			mid=(low+high)/2;
			if(c[mid]<=s[i])
				low=mid+1, v=mid;
			else
				high=mid-1;
		}

		if(v==-1) 		/* s[i] is too small */
			c[1]=s[i];
		else 			/* update c[] */
			c[v+1]=min(s[i], c[v+1]);

#ifdef LOL_DEBUG
		printf("i=%d v=%d\n", i, v);
		puts("c:");
		for(j=1; j<=n; j++)
			printf("%d ", c[j]);
		putchar('\n');
#endif
	}
	for(i=n; i>=1; i--)
		if(c[i]<INF)
			break;
	printf("%d\n", i);
	return 0;
}



希望上述文字能解决大家的一点问题——貌似我有点言语混乱。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值