hdu 1025 二分法求LIS

      题目给出两个序列，两个序列间有边相连，我们要选出不交叉的边，使得所选的边最多

     我们发现，如果我们对其中一个序列从小到大排序的话，这个问题就变成了求另一个序列的最长上升子序列的问题。

     由于这个题目的数据比较大，一般的求最长上升子序列的方法是n方的，显然不能满足题目的要求，我们建立一个数组g，g[i]表示长度为i的上升子序列最后一个数的最小值，那么显然g数组是严格单调的，确切的说是严格单调上升的，那么，我们想到利用二分的方法对LIS进行优化，对于一个数x，找到数组g中第一个大于x的数的位置，如果找不到，就将序列的长度增加1，g[len]=x，否则，由于x小于目前找到的这个位置上的数，我们用x更新g数组，即g[第一个大于x的位置]=x

    这样我们就得到了一个新的求LIS的算法，并且这个算法的时间复杂度是O(nlogn)，这样的话就可以很好的解决这个问题了。

    另外，这个题还有一个陷阱：就是当找到的序列长度如果大于1输出的是roads否则是road，我就因为没有注意到这个wa了好几次TAT

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;

struct node
{
    int x,y;
}a[500005];

int n,k,q=0;
int g[500005];

void lis(int x)
{
    int l=0,r=k-1;
    while (l<r)
    {
        int mid=(l+r)>>1;
        if (g[mid]>x) r=mid;
        else l=mid+1;
    }
    if ((l==k-1)&&(g[k-1]<x)) {g[k]=x; k++; return;}
    if (g[l]>x) g[l]=x;
}

bool cmp(node a,node b)
{
    return a.x<b.x;
}

int main()
{
    while (scanf("%d",&n)!=EOF)
    {
        q++;
        for (int i=0;i<n;i++)
        {
            int j,k;
            scanf("%d%d",&j,&k);
            a[i].x=j; a[i].y=k;
        }
        sort(a,a+n,cmp);
        k=1;
        memset(g,0,sizeof(g));
        g[0]=a[0].y;
        for (int i=1;i<n;i++) lis(a[i].y);
        if (k>1) printf("Case %d:\nMy king, at most %d roads can be built.\n\n",q,k);
        else printf("Case %d:\nMy king, at most %d road can be built.\n\n",q,k);
    }
    return 0;
}