POJ 3903 Stock Exchange (LIS:最长上升子序列)

标签：acm   algorithm   算法   dp   

POJ 3903 Stock Exchange (LIS:最长上升子序列)

http://poj.org/problem?id=3903

题意:

       给你一个长度为n (n<=100000) 的数字序列, 要你求该序列中的最长(严格)上升子序列的长度.

分析:

       由于n的规模达到10W, 所以只能用O(nlogn)的算法求.

       令g[i]==x表示当前遍历到的长度为i的所有最长上升子序列中的最小序列末尾值为x.(如果到目前为止, 根本不存在长i的上升序列, 那么x==INF无穷大)

       假设当前遍历到了第j个值即a[j], 那么先找到g[n]数组的值a[j]的下确界(即第一个>=a[j]值的g[k]的k值). 那么此时表明存在长度为k-1的最长上升子序列且该序列末尾的位置<j且该序列末尾值<a[j].

       那么我们可以令g[k]=a[j] 且 dp[i]=k (dp含义如解法1).

       (上面一段花时间仔细理解)

       最终我们可以找出下标最大的i使得: g[i]<INF 中i下标最大. 这个i就是LIS的长.

AC代码: O(n*logn)复杂度

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=100000+5;
const int INF=1e8;

int n;
int a[maxn];
int dp[maxn];
int g[maxn];

int main()
{
    while(scanf("%d",&n)==1)
    {
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]);

        for(int i=1;i<=n;i++)
            g[i]=INF;

        int ans=0;
        for(int i=1;i<=n;i++)
        {
            int k=lower_bound(g+1,g+n+1,a[i])-g;
            dp[i]=k;
            g[k]=a[i];
            ans=max(ans,k);
        }

        printf("%d\n",ans);
    }
    return 0;
}

POJ 3903 Stock Exchange (LIS:最长上升子序列)

标签：acm   algorithm   算法   dp   

原文地址：http://blog.csdn.net/u013480600/article/details/40823227