hdu 4604 Deque

题目

枚举每一个点，对于a[i]找出n~i中的最长非升非降的子序列，用nlogn的方法可以得到两个数组(b1,b2)，对于b1的元素个数便是最长非降序列的个数，但b1的元素不是最长子序列的元素，他存的是对应LIS长度的最小末尾，他肯定是有序的，所以我们接下来的进deque是对它操作。

在这两个数组中找到a[i]的次数，这两个数组因为都可以看成单调的，所以可以数组的两头开始找a[i]，得到一个起始下标，一个终点下标，那么在此序列中a[i]的次数(t-s+1)，当然你也可以用二分在序列中找a[i]的次数，二分快些，最后取二者的min。

#include<stdio.h>
#include<map>
#include<string.h>
#include<algorithm>
using namespace std;
#define N 100010

int t,n,a[N],b1[N],b2[N],len1,len2,p,c1,c2;
inline int input()
{
    bool IsN=0;
    int ret=0;
    char c;
    c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-') IsN=1;
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        ret=ret*10+c-'0';
        c=getchar();
    }
    return IsN?-ret:ret;
}

//非递减序列的二分
inline int erfen(int s,int t,int v)
{
    int m;
    while(s<=t)
    {
        m=(s+t)>>1;
        if(b1[m]<=v) s=m+1;
        else t=m-1;
    }
    return s;
}
//在二分得到的非递减序列的b1中找a[i]的次数，注意b1不是整个数组的非递减序列
inline int er(int s,int t,int v)
{
    if(b1[s]==v&&b1[t]==v)
    {
        return t-s+1;
    }
    if(s>=t) return 0;
    int m=(s+t)>>1;
    if(b1[m]==v) return er(s,m-1,v)+er(m+1,t,v)+1;
    if(b1[m]>v) return er(s,m-1,v);
    if(b1[m]<v) return er(m+1,t,v);
}
//非递增序列的二分
inline int erfen1(int s,int t,int v)
{
    int m;
    while(s<=t)
    {
        m=(s+t)>>1;
        if(b2[m]>=v) s=m+1;
        else t=m-1;
    }
    return s;
}
//在二分得到的非递增序列的b2中找a[i]的次数
inline int er1(int s,int t,int v)
{


    if(b2[s]==v&&b2[t]==v) return t-s+1;if(s>=t) return 0;
    int m=(s+t)>>1;
    if(b2[m]==v) return er1(s,m-1,v)+er1(m+1,t,v)+1;
    if(b2[m]<v) return er1(s,m-1,v);
    if(b2[m]>v) return er1(m+1,t,v);
}
int main()
{
    t=input();
    while(t--)
    {
        int ans=0;
        len1=len2=0;
        b1[0]=-100000000;
        b2[0]=100000000;
        int b_1,b_2;
        n=input();
        for(int i=1;i<=n;i++)
        {
            a[i]=input();
        }

        for(int i=n;i>=1;i--)
        {
            c1=0,c2=0;
            if(a[i]>=b1[len1])
            {
                len1++;
                p=len1;
            }
            else
            {
                p=erfen(1,len1,a[i]);
            }
            b1[p]=a[i];
            c1=er(1,len1,a[i]);
            b_1=p;
            if(a[i]<=b2[len2])
            {
                len2++;
                p=len2;
            }
            else
            {
                p=erfen1(1,len2,a[i]);
            }
            b2[p]=a[i];
            c2=er1(1,len2,a[i]);
            b_2=p;
            ans=max(ans,b_1+b_2-min(c1,c2));//减去2个序列中重复的a[i]
        }
        printf("%d\n",ans);
    }
    return 0;
}