暑期dp46道（16）hdoj 1025 Constructing Roads In JGShining's Kingdom最长上升子序列（二分）

题目链接：hdu 1025

题意：这题很容易想到最长上升子序列的模型，但是高高兴兴写出来也会TLE。。懵逼

  不过也是，看这个复杂度500000，o(n^2)铁定挂啊，所以就需要优化，二分法就可以。

话不多说，上模板：

void Do()
{
  dp[1]=data[1];
    int r,l,len=1;
    for(int i=2;i<=n;i++)
    {
        l=1;
        r=len;
        while(l<=r)
        {
            int mid=(l+r)/2;
            if(dp[mid]<data[i])
            {
                l=mid+1;
            }
            else
            {
                r=mid-1;
            }
        }
        dp[l]=data[i];
        if(l>len)
        {
            len++;
        }
    }
    printf("Case %d:\nMy king, at most %d %s can be built.\n\n",caseNum,len,(len==1)?"road":"roads");
}

这里dp[i]存储的是当前上升子序列长度为i的序列的最小元素值，所以初始化dp[1]=a[1];没有疑问，接着从1~len(len为当前所求得最大序列长度)二分遍历，从这些最小的元素中找到a[i]最适合插入的位置，更新dp[l];

如果，这时l>len说明最大序列长可以更新，而且l<=len+1；

这样时间复杂度就是O(nlogn)了，然后结束。。

如果还不是很懂，看下这篇解析：二分LIS解析

代码：

#include<cstdio>
#include<cstring>
#include<string>
#include<queue>
#include<algorithm>
using namespace std;
#define debug 0
#define M(a) memset(a,0,sizeof(a))
#define Max(a,b) ((a>b)?a:b)
#define Min(a,b) ((a<b)?a:b)
const int maxn=500000+5;
int data[maxn],dp[maxn],n,caseNum=0;
void Do()
{
	dp[1]=data[1];
	int r,l,len=1;
	for(int i=2;i<=n;i++)
	{
		l=1;
		r=len;//当前最长子序列长度 
		while(l<=r)//找a[i]适合插入的地方 
		{
			int mid=(l+r)/2;
			if(dp[mid]<data[i])
			{
				l=mid+1;
			}
			else
			{
				r=mid-1;
			}
		}
		dp[l]=data[i];//更新 
		if(l>len)
		{
			len++;
		}
	}
	printf("Case %d:\nMy king, at most %d %s can be built.\n\n",caseNum,len,(len==1)?"road":"roads");//注意这里复数。。。。 
}
int main()
{
#if debug 
	freopen("in.txt","r",stdin);
#endif//debug
	int x,y;
	while(~scanf("%d",&n))
	{
		caseNum++;
		M(dp);
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&x,&y);
			data[x]=y;//赋值时候注意下 
		}
		Do();
	}
	return 0;
}