POJ 2533 Longest Ordered Subsequence - from lanshui_Yang

 

        题目大意：求一个数列的最长上升子序列（严格上升）。

        解题思路：

方法一：O(n^2)

dp[i]：表示处理到第i个位置，序列的最长上升子序列末尾为i的长度； a[]数组存储原序列

dp[i] = max{dp[j]+1}，a[i]>a[j]，0≤j≤i

方法二：O(nlogn)

     复杂度降低其实是因为这个算法里面用到了二分搜索。本来有N个数要处理是O(n)，每次计算要查找N次还是O(n)，一共就是O(n^2)；现在搜索换成了O(logn)的二分搜索，总的复杂度就变为O(nlogn)了。

     这个算法的具体操作如下:

     开一个栈，每次取栈顶元素top和读到的元素temp做比较，如果temp > top 则将temp入栈；如果temp < top则二分查找栈中的比temp大的第1个数，并用temp替换它。最长序列长度即为栈的大小top。

    道理：对于数列中的a[i] 和a[j] ，i < j , 假设a[ i ]已在栈stap中，a[ j ] 未在栈中，这时读到元素a[ j ] , 如果a[ j ] < a[ i ](此时a[j]必然小于栈顶元素), 将a[j]与 a[i] 互换，此时这个栈的大小没有变化，但这个栈的“潜力”变大了，因为如果存在a[ j ] < a[ z ] < a[ i ] (i < j < z) ,当a[ i ]为栈顶元素时，a[ j ] 替换 a[ i ]后成为栈顶元素，此后在读到a[ z  ] 时就能把a[z] 压入栈，栈的大小就增加 1 ， 即最长上升子序列长度就增加了1。

    举例：原序列为1，5，8，3，6，7

    栈为1，5，8，此时读到3，用3替换5，得到1，3，8； 再读6，用6替换8，得到1，3，6；再读7，得到最终栈为1，3，6，7。最长递增子序列为长度4。

请看代码：

 

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<algorithm>
#include<cmath>
using namespace std ;
const int MAXN = 1e3 + 5 ;
int s[MAXN] ;
int dp[MAXN] ;
int n ;
int stap[MAXN] ;
int top ;
void init()
{
    int i , j ;
    for(i = 1 ; i <= n ; i ++)
    {
        scanf("%d" , &s[i]) ;
    }
}
void solve()  // O(n^2) 算法
{
    int i , j ;
    dp[1] = 1 ;
    for(i = 2 ; i <= n ; i ++)
    {
        dp[i] = 1 ;
        for(j = 1 ; j < i ; j ++)
        {
            if(s[i] > s[j] && dp[i] < dp[j] + 1)
            {
                dp[i] = dp[j] + 1 ;
            }
        }
    }
    int MAX = -1 ;
    for(i = 1 ; i <= n ; i ++)
    {
        if(MAX < dp[i])
        MAX = dp[i] ;
    }
    printf("%d\n" , MAX) ;
}

void solve2()  // O(n log n) 算法
{
    top = 0 ;
    stap[++ top] = s[1] ;
    int i ;
    for(i = 2 ; i <= n ; i ++)
    {
        if(s[i] > stap[top])
        {
            stap[++ top] = s[i] ;
        }
        else if(s[i] < stap[top])
        {
            int left , right , mid ;
            left = 1 ;
            right = top ;
            while (left < right)
            {
                mid = (left + right) / 2 ;
                if(stap[mid] < s[i])
                {
                    left = mid + 1 ;
                }
                else if(stap[mid] == s[i])
                {
                    break ;
                }
                else
                {
                    right = mid ;  // 注意不是 mid - 1  ！！
                }
            }
            mid = (left + right) / 2 ;
            stap[mid] = s[i] ;
        }
    }
    printf("%d\n" , top) ;
}
int main()
{
    while (scanf("%d" , &n) != EOF)
    {
        init() ;
        //solve() ;
        solve2() ;
    }
    return 0 ;
}