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;
}

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