【Luogu】 P3928 SAC E#1 - 一道简单题 Sequence2

【题目】洛谷10月月赛R1 提高组

【算法】递推DP+树状数组

【题解】列出DP递推方程，然后用树状数组维护前后缀和。

#include<cstdio>
#include<cstring>
#include<cctype>
#include<cmath>
#include<algorithm>
#define lowbit(x) (x&-x)
#define ll long long
using namespace std;
int read(){
    char c;int s=,t=;
    while(!isdigit(c=getchar()))if(c=='-')t=-;
    do{s=s*+c-'';}while(isdigit(c=getchar()));
    return s*t;
}
int min(int a,int b){return a<b?a:b;}
int max(int a,int b){return a<b?b:a;}
int abs(int x){return x>?x:-x;}
//void insert(int u,int v){tot++;e[tot].v=v;e[tot].from=first[u];first[u]=tot;}
/*------------------------------------------------------------*/
const int inf=0x3f3f3f3f,maxn=;

int n,a[maxn][],f[maxn][],ans,tot,c[][maxn],b[maxn];

void modify(int d,int x,int k){if(d&)x=tot-x+;for(int i=x;i<=tot;i+=lowbit(i))c[d][i]=max(c[d][i],k);}
int ask(int d,int x){if(d&)x=tot-x+;int as=;for(int i=x;i>=;i-=lowbit(i))as=max(as,c[d][i]);return as;}
int main(){
    n=read();tot=;
    for(int i=;i<=n;i++)a[i][]=read(),f[i][]=,b[++tot]=a[i][];
    for(int i=;i<=n;i++)a[i][]=read(),f[i][]=,b[++tot]=a[i][];
    for(int i=;i<=n;i++)a[i][]=read(),f[i][]=f[i][]=,b[++tot]=a[i][];
    sort(b+,b+tot+);
    tot=unique(b+,b+tot+)-b-;
    for(int i=;i<=n;i++)a[i][]=lower_bound(b+,b+tot+,a[i][])-b;
    for(int i=;i<=n;i++)a[i][]=lower_bound(b+,b+tot+,a[i][])-b;
    for(int i=;i<=n;i++)a[i][]=lower_bound(b+,b+tot+,a[i][])-b;
    modify(,a[][],);modify(,a[][],);modify(,a[][],);
    modify(,a[][],);modify(,a[][],);modify(,a[][],);
    modify(,a[][],);modify(,a[][],);modify(,a[][],);
    modify(,a[][],);modify(,a[][],);modify(,a[][],);
    for(int i=;i<=n;i++){
        f[i][]=max(f[i][],ask(,a[i][])+);ans=max(ans,f[i][]);
        f[i][]=max(f[i][],ask(,a[i][])+);ans=max(ans,f[i][]);
        f[i][]=max(f[i][],ask(,a[i][])+);ans=max(ans,f[i][]);
        f[i][]=max(f[i][],ask(,a[i][])+);ans=max(ans,f[i][]);
        modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);modify(,a[i][],max(f[i][],f[i][]));
         modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);modify(,a[i][],max(f[i][],f[i][]));
        modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);
        modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);modify(,a[i][],f[i][]);//
    }
    printf("%d",ans);
    return ;
}