编程之美读书笔记2.16—求数列中最长递增序列

解法1:  时间复杂度为 O（N^2）

      

.cpp

#include <iostream>
using namespace std;

int Lis(int A[],int length)
{
	int *L=new int[length];

	for(int i=0;i<length;i++)
	{
		L[i]=1;     //初始化默认的长度
		for(int j=0;j<i;j++)      //找出前面最长的序列
		{
			if(A[i]>A[j]  && L[j]+1>L[i])
				L[i]=L[j]+1;
        }
    }

	int max=1;
	for(int i=0;i<length;i++)     
	{
		if(L[i]>max)      //取L的最大值
			max=L[i];
	}

	return max;
}

int main()
{
	int  A[]={1,-1,2,-3,4,-5,6,-7};
	int length=8;

	cout<<Lis(A,length);

	system("pause");
}

解法2：  时间复杂度为 O（N^2）

#include <iostream>
using namespace std;
int Min(int *data,int length)
{
	int max=data[0];
	for(int i=0;i<length;i++)     
	{
		if(data[i]<max)      //取L的最大值
			max=data[i];
	}
	return max;
}

int Lis(int A[],int length)
{
	int *L=new int[length];
	int *MaxV=new int[length+1];
	MaxV[1]=A[0];
	MaxV[0]=Min(A,length)-1;

	for(int i=0;i<length;i++)
	{
		L[i]=1;     //初始化默认的长度
	}

	int  max=1;
	for(int i=1;i<length;i++)
	{
		int j;
		for(j=max;j>=0;j--)
		{
			if(A[i]>MaxV[j])  //A[i]大于长度为j的序列中最大元素的最小值
			{
				L[i]=j+1;
				break;
			}
		}
		
		//更新递增序列中最大元素的最小值
		if(L[i]>max)
		{
			max=L[i];
			MaxV[L[i]]=A[i];
        }
		else if(MaxV[j]<A[i]  && A[i]<MaxV[j+1])
		{
			MaxV[j+1]=A[i];
		}
	}
	

	return max;
}

int main()
{
	int  A[]={1,-1,2,-3,4,-5,6,-7};
	int length=8;

	cout<<Lis(A,length);

	system("pause");
}

解法3： 时间复杂度为 O（ N* log 2 N）

      统计出的MaxV数组具有很好的递增关系（动态规划都是用来解决最优化问题，我们总能通过优化关系对之前统计出的结果进行排序），即如果i<j，那么有MaxV[i]<MaxV[j]，最后用二分查找法来阶段k的元素在MaxV数组中的位置即可。

证明：反证法，假设当i<j<=k，有MaxV[i]>=MaxV[j]，那么存在一个长度为i的递增序列a1a2...ai，且ai是计算到阶段k时所有长度为i的递增序列中最大元素的最小值，以及长度为j的序列b1b2...bj且ai>=bj，由于i<j长度j的序列b1b2...bj包含一个子序列b1b2...bi，且bi<bj<=ai，即长度为i的递增子序列最大元素的最小值不是ai，矛盾。

 二分搜索  

#include <iostream>
#include <algorithm>

using namespace std;
int Lis(int A[],int length)
{
	int i, j, k;
	
	int *MaxV=new int[length+1];
	MaxV[1] = A[0];
	int nmax = 1;	// 目前找到的最长递增子序列的长度
	// 有n个阶段，每个阶段有1个状态
	for (i=1; i<length; i++)
	{
		// 每个状态有nmaxv种决策，以得出以元素i结尾的最长递归子序列的长度
		int u = 1, v = nmax;
		while (u<=v)    //二分搜索
		{
			int mid = (u+v)>>1;
			if (MaxV[mid] < A[i])
				u = mid+1;
			else
				v = mid-1;
		}

		// 每个元素都会对数组MaxV中的某个元素产生影响
		// 使用二分搜索确定其在第v+1个位置
		nmax = max(nmax, v+1);
		MaxV[v+1] = A[i];
	}

	return nmax;
}

int main()
{
	int length=8;
    int  A[]={1,-1,2,-3,4,-5,6,-7};
   
	cout<<Lis(A,length);
    
	system("pause");
}