最长上升子序列 LIS 〇(nlogn)

以前学过一些最长上升子序列的方法，不过时间复杂度都是挺高的，一般的LIS算法的时间复杂度都是〇(n^2)的，但是有些题目对算法的时间要求比较高，数据量比较大，所以

就需要更高效的算法。今天经tengyu大牛指点以后，学到了一种比较好的算法，算法时间复杂度是〇(nlogn)，算法思想如下：

假设要寻找最长上升子序列的序列是a[n]，然后寻找到的递增子序列放入到数组b中。

（1）当遍历到数组a的第一个元素的时候，就将这个元素放入到b数组中，以后遍历到的元素都和已经放入到b数组中的元素进行比较；

（2）如果比b数组中的每个元素都大，则将该元素插入到b数组的最后一个元素，并且b数组的长度要加1；

（3）如果比b数组中最后一个元素小，就要运用二分法进行查找，查找出第一个比该元素大的最小的元素，然后将其替换。

在这个过程中，只重复执行这两步就可以了，最后b数组的长度就是最长的上升子序列长度。

 

代码如下：

[code]

/*
   LIS(最长上升子序列) 〇(nlogn的算法)
*/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int a[1005];
int b[1005];

int main()
{
   int t;
   while(scanf("%d",&t)!=EOF)
   {
      memset(b,0,sizeof(b));
      int k=0;
      for(int i=0;i<t;i++)
      {
         scanf("%d",&a[i]);
         if(i==0){b[k++]=a[i];continue;}
         if(b[k-1]<a[i]){b[k++]=a[i];continue;}
         //if(b[0]>a[i])
           // b[0]=a[i];
         int mid,l=0,h=k-1;
         while(l<h)
         {
            mid=l+h>>1;
            if(b[mid]>a[i])
               h=mid-1;
            else if(b[mid]<a[i]) l=mid+1;
            else break;
         }
         mid=l+h>>1;
         if(l<h) b[mid]=a[i];
         else
         {
            for(int j=mid;j<k;j++)
               if(b[j]>=a[i]){b[j]=a[i];break;}
         }
      }
      printf("%d\n",k);
   }
  
   return 0;
}

 

 

[/code]