最长上升子序列 O(nlogn)

   题目描述

   INPUT    ：  输入一个整数n 表示接下来有 n 个数输入。

   OUTPUT：  输出当前数列最长上升子序列的长度。

   直接上代码 最长下降子序列 最长山峰序列都可以以该问题为母问题 进行延伸。

#include <iostream>
#include <stdio.h>
using namespace std;
#define MAX 1000
int seq[MAX+10];//用于存数
int seqlen[MAX+10];//用于存以第i个数为终点的最长上升子序列的长度
int main()
{
    int n,max,maxlen;
    while(cin>>n)
    {
        maxlen=1;
        for(int i=1;i<=n;i++)
            scanf("%d",&seq[i]);

         for(int i=1;i<=9;i++)
            seqlen[i]=1;

         for(int i=2;i<=n;i++)
         {
             max=0;
             for(int j=1;j<=i-1;j++)
             {
                 if(seq[j]<seq[i]&&max<seqlen[j])
                 max=seqlen[j];
             }//因为为了满足每一步都是最优解 需要进行比较
             seqlen[i]=max+1;
             if(maxlen<seqlen[i])
             maxlen=seqlen[i];
         }
         cout<<maxlen<<endl;
    }

    return 0;
}   

博主在点头网上准备水一发 最长上升子序列的问题的时候用上面的写法发现超时 5000+ms 于是百度了一番 找到了更快的方法 复杂度O(nlogn)

考虑两个dp值相同的数的时候可以得到优化  若 x<y && dp[x]==dp[y] 那么应该选择数值较小的ｘ　因为可能存在z  x<z<y

代码如下

/**
最长递增子序列O(nlogn)算法：
状态转移方程：f[i] = max{f[i],f[j]+1},1<=j<i,a[j]<a[i].
分析：加入x<y,f[x]>=f[y],则x相对于y更有潜力。
首先根据f[]值分类，记录满足f[t]=k的最小的值a[t],记d[k]=min{a[t]},f[t]=k.
    1.发现d[k]在计算过程中单调不上升
    2.d[1]<d[2]<...<d[k] (反证) 1 2 3 8 4 7
解法：
1. 设当前最长递增子序列为len,考虑元素a[i];
2. 若d[len]<a[i],则len++，并将d[len]=a[i];
   否则,在d[0-len]中二分查找,找到第一个比它小的元素d[k],并d[k+1]=a[i].()
*/
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N = 50000;
int a[N];       //a[i] 原始数据
int d[N];       //d[i] 长度为i的递增子序列的最小值

int BinSearch(int key, int* d, int low, int high)
{
    while(low<=high)
    {
        int mid = (low+high)>>1;
        if(key>d[mid] && key<=d[mid+1])
            return mid;
        else if(key>d[mid])
            low = mid+1;
        else
            high = mid-1;
    }
    return 0;
}

int LIS(int* a, int n, int* d)
{
    int i,j;
    d[1] = a[1];
    int len = 1;        //递增子序列长度
    for(i = 2; i <= n; i++)
    {
        if(d[len]<a[i])
            j = ++len;
        else
            j = BinSearch(a[i],d,1,len) + 1;
        d[j] = a[i];
    }
    return len;
}

int main()
{
    int n;
    while(cin>>n)
    {
        for(int i=1;i<=n;i++)
          scanf("%d",&a[i]);
        printf("%d\n",LIS(a,n,d));
    }
    return 0;
}