zoukankan      html  css  js  c++  java
  • hdu5087 Revenge of LIS II (dp)

    只要理解了LIS,这道题稍微搞一下就行了。

    求LIS(最长上升子序列)有两种方法:

    1.O(n^2)的算法:设dp[i]为以a[i]结尾的最长上升子序列的长度。dp[i]最少也得是1,就初始化为1,则dp[i]=max(dp[i],dp[j]+1)(其中j<i且a[j]<a[i])。

    int gao()
    {
        int ans=0;
        for(int i=0;i<n;i++)
        {
            dp[i]=1;
            for(int j=0;j<i;j++)
            {
                if(a[j]<a[i])
                {
                    dp[i]=max(dp[i],dp[j]+1);
                }
            }
            ans=max(ans,dp[i]);
        }
        return ans;
    }
    O(n^2)

    2.O(n*lgn)算法:有这样一个性质:如果上升子序列长度相同,那么最末尾最小的那种在之后的比较中最有优势(有点贪心的味道)。根据这个,设dp[i]为长度为i+1的上升子序列中末尾元素的最小值。开始dp全初始化为INF,每次更新时前面的dp都是排好序的,所以就可以二分来找。

    int gao()
    {
        fill(dp,dp+n,INF);
        for(int i=0;i<n;i++)
        {
            *lower_bound(dp,dp+n,a[i])=a[i];
        }
        ans=lower_bound(dp,dp+n,INF)-dp;
    }
    O(n*lgn)

    对于本题,要求输出第二长子序列的长度。那么就找LIS,如果LIS只有一种情况可以得到,那么输出LIS-1,否则输出LIS。

    用num数组来记录能达到当前长度的方法数。

    #include<iostream>
    #include<cstdio>
    #include<cstdlib>
    #include<cstring>
    #include<string>
    #include<cmath>
    #include<map>
    #include<set>
    #include<vector>
    #include<algorithm>
    #include<stack>
    #include<queue>
    #include<cctype>
    #include<sstream>
    using namespace std;
    #define pii pair<int,int>
    #define LL long long int
    const int eps=1e-8;
    const int INF=1000000000;
    const int maxn=1000+10;
    const LL mod=1000000007;
    int T,n,a[maxn],dp[maxn],num[maxn];
    int main()
    {
        //freopen("in8.txt","r",stdin);
        //freopen("out.txt","w",stdout);
        scanf("%d",&T);
        while(T--)
        {
            scanf("%d",&n);
            int ans=0;
            for(int i=0; i<n; i++)
            {
                scanf("%d",&a[i]);
            }
            for(int i=0; i<n; i++)
            {
                dp[i]=1;
                num[i]=1;
                for(int j=0; j<i; j++)
                {
                    if(a[j]<a[i])
                    {
                        if(dp[j]+1>dp[i])
                        {
                            num[i]=num[j];
                            dp[i]=dp[j]+1;
                        }
                        else if(dp[j]+1==dp[i])
                        {
                            num[i]+=num[j];
                        }
                    }
                }
                if(dp[i]>ans)
                {
                    ans=dp[i];
                }
            }
            int num1=0;
            for(int i=0;i<n;i++)
            {
                if(dp[i]==ans)
                {
                    num1+=num[i];
                }
            }
            if(num1==1)
            {
                printf("%d
    ",ans-1);
            }
            else
            {
                printf("%d
    ",ans);
            }
        }
        //fclose(stdin);
        //fclose(stdout);
        return 0;
    }
  • 相关阅读:
    在emacs上使用博客园的代码功能
    Programming Pearls笔记之一
    Virtualbox中Archlinux联网问题
    微信公众平台消息接口开发集成解决方案
    发送短信
    基于JMS的数据交换既数据互操作平台的解决方案
    Spring MVC基于注解的Junit测试
    获取设置一个字节某一个位的数值
    NotificationManager
    调用系统联系人列表
  • 原文地址:https://www.cnblogs.com/zywscq/p/4113153.html
Copyright © 2011-2022 走看看