最长不降子序列LIS.docxVIP

LIS。。Longest Increasing Subsequence ；最长上升/不下降子序列；　　最长上升子序列：　　有两种基本方法：两个时间复杂度分别为O(n^2)和O(nlogn)　　先来看第1种：　　动态规划：　　对于给定数列E，元素个数为n,最长上升子序列Q满足对任意1=ij=n，有Q[i]Q[j]，且E[i]E[j]。　　容易得出O(n^2)的DP递推公式：　　D[i]=max{D[j]}+1;(1=ji且E[j]E[i])　　D[i]为以元素i结尾的最长子序列个数。　　这样经过两重循环一次遍历可以得到最长上升子序列。　　代码如下：　　#includeiostream　　using namespace std;　　int n;　　int a[1001],b[1001];　　int main()　　{ int i,j;　　scanf(%d,n);　　for(i=1;i=n;i++)　　{　　scanf(%d,a[i]);　　b[i]=1;　　}　　for(i=1;i=n;i++)　　{ int maxx=0;　　for(j=1;j=i;j++)　　{?　　if(a[i]a[j]b[j]maxx)　　{ maxx=b[j];　　}　　}　　b[i]=maxx+1;　　}　　int best=0;　　for(i=1;i=n;i++)　　{ if(bestb[i])　　best=b[i];　　}　　printf(%d\n,best);　　system(pause);　　return 0;　　}　　第2种方法时间复杂度为O(nlogn)，用到了二分查找.(又称作CMI算法)　　其操作如下：　　开辟一个栈b，每次取栈顶元素s和读到的元素a做比较，如果as，则置为栈顶；如果as，则二分查找栈中的比a大的第1个数，并替换。最终栈的大小即为最长递增子序列为长度　　考察b栈内每个元素的含义,b[i] 表示所有长度为i的上升子序列中最小的最后一个数.　　?　　举例：原序列为3，4，5，2，4，2　　栈为3，4，5，此时读到2，则用2替换3，得到栈中元素为2，4，5，再读4，用4替换5，得到2，4，4，再读2，得到最终栈为2，2，4,最终得到的解是:?　　长度为1的上升子序列中最小的最后一个数是2 (2)　　长度为2的上升子序列中最小的最后一个数是2 (2,2)长度为3的上升子序列中最小的最后一个数是4 (3,4,4)?　　可知没有长度为4的上升子序列,最长递增子序列长度为3. (3,4,4)　　CMI本质是LIS问题的另一种动态规划思路　　注意:CMI只能求LIS的长度和最后一个数,不能求LIS的序列!　　代码如下：　　#includeiostream　　using namespace std;　　int n;　　int a[1001],b[1001];　　int rear;　　int solve(int t)　　{ int l=1,r=rear;　　while(l=r)　　{ int mid=(l+r)1;　　if(b[mid]=t)//若为非递减序列，则为b[mid]t　　r=mid-1;　　else?　　l=mid+1;　　}　　if(lrear)　　rear=l;　　return l;　　}　　int main()　　{ int i,j;　　scanf(%d,n);　　rear=0;　　for(i=1;i=n;i++)　　{　　scanf(%d,a[i]);　　b[solve(a[i])]=a[i];　　}　　printf(%d\n,rear);　　system(pause);　　return 0;　　}