数组中最长递增子序列

写一个时间复杂度尽可能低的程序，求一个一维数组（N个元素）中的最长递增子序列的长度。

例如：在序列1，-1，2，-3，4，-5，6，-7中，其最长的递增子序列为1，2，4，6。

分析与解法

根据题目的要求，求一维数组中的最长递增子序列，也就是找一个标号的序列b[0],b[1],…,b[m](0 <= b[0] < b[1] < … < b[m] < N)，使得array[b[0]]<array[b[1]]<…<array[b[m]]。

解法一

根据无后效性的定义我们知道，将各阶段按照一定的次序排列好之后，对于某个给定的阶段状态来说，它以前各阶段的状态无法直接影响它未来的决策，而只能间接地通过当前的这个状态来影响。换句话说，每个状态都是历史的一个完整总结。

同样的，仍以序列1，-1，2，-3，4，-5，6，-7为例，我们在找到4之后，并不关心4之前的两个值具体是怎样，因为它对找到6没有直接影响。因此，这个问题满足无后效性，可以通过使用动态规划来解决。

可以通过数字的规律来分析目标串：1，-1，2，-3，4，-5，6，-7。

使用i来表示当前遍历的位置

当i=1时，显然，最长的递增序列为（1），序列长度为1.

当i=2是，由于-1<1。因此，必须丢弃第一个值后重新建立串。当前的递增序列为（-1），长度为1。

当i=3时，由于2>1,2>-1。因此，最长的递增序列为（1,2）,(-1,2),长度为2。在这里，2前面是1还是-1对求出后面的递增序列没有直接影响。（但是在其它情况下可能有影响）

依此类推之后，我们得出如下的结论。

假设在目标数组array[]的前i个元素中，最长递增子序列的长度为LIS[i]。那么，

LIS[i+1]=max{1,LIS[k]+1},array[i+1]>array[k],for any k <= i

即如果array[i+1]大于array[k],那么第i+1个元素可以接在LIS[k]长的子序列后面构成一个更长的子序列。于此同时array[i+1]本身至少可以构成一个长度为1的子序列。

根据上面的分析，就可以得到代码清单

C++代码：

int Max(int *a, int n)

{

int max = a[0];

for(int i = 1; i < n; i++)

if(max < a[i])

max = a[i];

return max;

}

int LIS(vector<int> &array)

{

int *a = new int[array.size()];

for(int i = 0; i < array.size(); i++)

{

a[i] = 1; //初始化默认的长度

for(int j = 0; j < i; j++) //前面最长的序列

{

if(array [i] > array [j] && a[j] + 1 > a[i])

{

a[i] = a[j] + 1;

}

}

}

return Max(a, array.size());

}

这种方法的时间复杂度为O(N2 + N) = O(N2)

o(n*logn)算法：

保留递增长度为m的最小值n：b[m]=n。

有b[m] >= b[m-1]：

假设当前最大递增长度为k，当前比较元素为a[j]：

If a[j] > b[k]，k++,b[k]=a[j]；

Else 将a[j]在b[1]到b[k]进行二分查找：If a[j] <= b[1]，b[1] = a[j]

Else a[j] >b[i] && a[j]<b[i+1]，b[i+1] = a[j]

初始化：b[1] = a[0];

代码：

#include<stdio.h>

#include<stdlib.h>

int binary_search(int num, int b[], int k){

int left=1,mid,right=k;

if(num <= b[1]){

return 1;

}

mid = (left+mid)/2;

while(left <= right){

if(num == b[mid]){

return mid;

}else{

if(num > b[mid]){

left = mid + 1;

}else{

right = mid - 1;

}

}

}

return left;

}

int main(){

int i,k,n,ret;

int a[100],b[100];

scanf("%d",&n);

if(n == 0){

printf("no data err\n");

return main();

}

for(i=1; i<=n; i++)

scanf("%d",&a[i]);

k = 1;

b[k] = a[1];

for(i=2; i<=n; i++){

if(a[i] > b[k]){

b[++k] = a[i];

}else{

ret = binary_search(a[i],b,k);

b[ret] = a[i];

}

}

printf("%d",k);

system("pause");

return main();

}