HDU-1025,最长上升子序列(nlogn)算法+详解

题目连接

分析：

此题就是一个最长上升子序列的一个变形题，没什么难度。关于LIS问题，有两种算法，一种是dp，时间复杂度为n^2，另一个是二分+Dp，之所以能用二分查找，关键在于它用辅助数组保存了具有递增的数据，正好适宜二分查找。而查找的整个序列D中，有错落的递增式，也就必定存在最长的一个，它们的元素就构成了一个递增有序的集合M，第一个必定是在M中最小，而且第一个最小的元素，一定是在D中靠前的。

分析一下nlogn的方法：设辅助数组为dp[]。

可以知道最短len=1；此时dp[1]=a[1]，若第二个元素要大于第一个，把它加入dp[2]，len++,否则就更新dp[1]，len不加，分析第三元素就可以推广到n个了，这样就降低分析问题的

难度了。如果第三个大于地二个加入dp[3]中,如果小于第二个，直接更新dp[2]。要判断j>len（j代表dp数组中更新元素的位置）这样2==2，不需要加。如果小于第一个，更新dp[1],

len=2>1，不需要加。通过进一步分析：我们发现，dp数组中是递增的。它含有的元素的个数是当前上升子序列的个数（除0），当有数据不是加入dp数组中，要更新前面的（用二分法找）.而它不会影响结果，因为dp数组中值（除0以外）的个数并没增加，要使len 增加，有两种可能：

                              1:后面有元素大于dp数组的最后一个元素就有(j>len)。

                              2:被另一个元素连续更新,且更新到j>len,意思就是说：选出另一元素它能比当前元素构造的上升子序列要多。这样就能选出最长上升子序列了。

这就是所谓的nlogn算法。也不过如此吧！

代码如下：(注意格式 别被“road”和“roads”坑了）

#include<stdio.h>

int dp[500005],road[500005];

long search(long n,long m)
{
	int L=1,R=n,mid;
	while(L<=R)
	{
	  mid=(L+R)>>1;
	  if(m>dp[mid]) L=mid+1;
	  else R=mid-1;
	}
	return L;
}

int main()
{
	int x,y;
	long t,i,j,len,ans=0;
	while(scanf("%ld",&t)!=EOF)
	{
	 for(i=0;i<t;i++)
	 {
	  scanf("%d%d",&x,&y);
	  road[x]=y;
	 }
	len=1;
	dp[1]=road[1];
	for(i=2;i<=t;i++)
	{
	    j=search(len,road[i]);
		dp[j]=road[i];
		if(j>len) len++;
	}
    printf("Case %d:\nMy king, at most %d road",++ans,len);  
    if(len!=1) printf("s");  
    printf(" can be built.\n\n");  
    }
    return 0;
}