导弹拦截

转载自洛谷@yhwtd。
题目：导弹拦截

一个手写二分查找求最长不下降子序列和最长上升子序列的方法

要过加强的数据就只能用NlogN的方法，第一个问就是求最长不下降子序列长度，第二个问求最长不下降子序列的个数其实（我也不能证明）就是求最长上升子序列长度，而NlogN求最长上升子序列的方法是用一个数组d[]存放一组“潜力最大的数字”，它存的并不是最长上升子序列，长度一样，并且“潜力”更大，而且每次更新长度时只与最后一个数作比较，所以更新长度时是绝对正确的，以最长上升子序列为例，每次用一个数更新d[]的方法就是查找d[]中第一个比它大的位置并替换掉来提升“潜力“，若比最后一个数还大，则把该数放在最后，更新d[]的长度。全部数都用来更新过后d[]的长度就是最长上升子序列的长度。最长不下降子序列的方法类似，每次更新替换掉第一个比他小的数即可。因为d[]是单调的，所以可以用二分查找的方法。

代码：

#include<cstdio>
#include<cstring>

int a[100005],d[100005],n;

int lis1()
{
    int l,r=n,mid,cnt=0,ans;
    d[++cnt]=a[1];
    for(int i=2;i<=n;i++)
    {
        if(a[i]<=d[cnt]) d[++cnt]=a[i];
        else
        {
            l=1;
            r=cnt;
            while(l<=r)
            {
                mid=(l+r)>>1;
                if(d[mid]<a[i])
                {
                    ans=mid;
                    r=mid-1;
                }
                else l=mid+1;
            }
            d[ans]=a[i];
        }
    }
    return cnt;
}

int lis2()
{
    memset(d,0,sizeof(d));
    int l,r=n,mid,cnt=0,ans;
    d[++cnt]=a[1];
    for(int i=2;i<=n;i++)
    {
        if(a[i]>d[cnt]) d[++cnt]=a[i];
        else
        {
            l=1;
            r=cnt;
            while(l<=r)
            {
                mid=(l+r)>>1;
                if(d[mid]>=a[i])
                {
                    ans=mid;
                    r=mid-1;
                }
                else l=mid+1;
            }
            d[ans]=a[i];
        }
    }
    return cnt;
}

int main()
{
    int ans;
    while(scanf("%d",&a[++n])==1);
    n-=1;
    printf("%d\n%d",lis1(),lis2());
}