Longest Increasing Subsequence HDU - 6284

/*
首先预处理好f g数组
 fi :以a[i]为结尾的 最长上升子序列的长度
 gi :以a[i]为开始的 最长上升子序列的长度
 mxx : 最长上升子序列的长度
线段树优化 nlogn
(不包含a[i]==0)

显然把所有0换成x  只可能是mxx变成mxx+1 

 然后我们考虑一对 i j (下标)
 若 f[i]+g[j]==mxx 则 所有a[i]+1~~~a[j]-1之间的x
 他们对用的lis长度为mxx+1
 然后枚举i j凉了
 对于一个i 我们只需要找到 他后面的 一个j 满足  f[i]+g[j]==mxx 并且a[j]最大
 然后维护bg[i] 表示长度为g[j]==i的所有的 a[j]中最大的
 从后往前枚举i 然后维护 bg O(1)转移
 上述过程可能 i和bg维护的j之间 他没有0 那就不能转移
 所以 按0分段  遇到0 就把之前的信息更新bg 

 然后没了 

*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#define lc k*2
#define rc k*2+1
#define mid (l+r)/2
#define maxn 400010
#define ll long long
using namespace std;
ll n,a[maxn],f[maxn],s[maxn],g[maxn],as[maxn],bg[maxn],c[maxn][];
void Insert(ll k,ll l,ll r,ll x,ll y){
    if(x==l&&r==x){
        s[k]=max(s[k],y);return;
    }
    if(x<=mid)Insert(lc,l,mid,x,y);
    else Insert(rc,mid+,r,x,y);
    s[k]=max(s[lc],s[rc]);
}
ll Query(ll k,ll l,ll r,ll x,ll y){
    if(x>y)return ;
    if(x<=l&&y>=r)return s[k];
    ll res=;
    if(x<=mid)res=max(res,Query(lc,l,mid,x,y));
    if(y>mid)res=max(res,Query(rc,mid+,r,x,y));
    return res;
}
int main(){
    while(~scanf("%lld",&n)){
        for(ll i=;i<=n*;i++)
            s[i]=f[i]=g[i]=as[i]=;
        for(ll i=;i<=n;i++){
            scanf("%lld",&a[i]);
            //a[i]=rand();
            f[i]=;g[i]=;
        }
        ll mxx=;
        for(ll i=;i<=n;i++){
            if(a[i]==)continue;
            ll mx=Query(,,n,,a[i]-);
            f[i]=mx+;mxx=max(mxx,f[i]);
            Insert(,,n,a[i],f[i]);
        }
        for(ll i=;i<=n*;i++)s[i]=;
        for(ll i=n;i>=;i--){
            if(a[i]==)continue;
            ll mx=Query(,,n,a[i]+,n);
            g[i]=mx+;Insert(,,n,a[i],g[i]);
        }
        for(ll i=;i<=n*;i++)bg[i]=;
        ll cnt=;a[]=-;
        for(ll i=n;i>=;i--){
            if(a[i]==){
                for(ll j=;j<=cnt;j++)
                    bg[c[j][]]=max(bg[c[j][]],c[j][]);
                cnt=;bg[]=n+;
            }
            else{
                ll mx=bg[mxx-f[i]];
                c[++cnt][]=g[i];c[cnt][]=a[i];
                if(mx-<a[i]+)continue;
                as[a[i]+]++;as[mx]--;
            }
        }
        ll ans=;
        for(ll i=;i<=n;i++)as[i]+=as[i-];
        for(ll i=;i<=n;i++){
            if(as[i]>)ans+=i*(mxx+);
            else ans+=i*mxx;
            //("%lld\n",ans);
        }
        printf("%lld\n",ans);
    }
    return ;
}