[BZOJ4553][HEOI2016/TJOI2016]序列

最新推荐文章于 2021-12-10 13:18:55 发布

denglunza622821951

最新推荐文章于 2021-12-10 13:18:55 发布

阅读量97

点赞数

原文链接：http://www.cnblogs.com/ymzqwq/p/11278980.html

版权

传送门

好像是DP再套个裸的CDQ？
~~树套树是不可能写树套树的，这辈子都不可能写树套树的~~

对于一个 \(i\) ，设它最小为 \(a_i\) ，原数为 \(b_i\) ，最大为 \(c_i\)

\(f_i\) 表示 \(i\) 结尾的最长子序列， \(f_i=f_j+1\) ， \(j\) 要满足

\(j<i\)
\(c_j \leq b_i\)
\(b_j \leq a_i\)

这不就CDQ套个树状数组就完了嘛QAQ

具体的话，把 \([L,mid]\) 按 \(c\) 排序，\([mid+1,r]\) 按 \(b\) 排序，然后像求三维偏序一样套个树状数组就行了。
处理区间的顺序可能要注意下，树状数组的清空也要控制好复杂度。

然后交了一发就过了0.0

#include<bits/stdc++.h>
#define LL long long
#define fr(i,x,y) for(int i=(x);i<=(y);i++)
#define rf(i,x,y) for(int i=(x);i>=(y);i--)
#define frl(i,x,y) for(int i=(x);i<(y);i++)
using namespace std;
const int N=100005;
const int M=N<<1;
int n,m,a[N],b[N],c[N];
int f[N];
int mx=100000;
int num[N],t[N];

void read(int &x){
    char ch=getchar();x=0;
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
}

void chkmin(int &x,int y){ if (y<x) x=y; }
void chkmax(int &x,int y){ if (y>x) x=y; }

inline int cmpc(const int &q,const int &w){ return c[q]<c[w]; }
inline int cmpb(const int &q,const int &w){ return b[q]<b[w]; }

inline void ADD(int x,int v){
    for(;x<=mx;x+=x&(-x)) chkmax(t[x],v);
}

inline int QRY(int x){
    int ans=0;
    for(;x;x^=x&(-x)) chkmax(ans,t[x]);
    return ans;
}

inline void DEL(int x){
    for(;x<=mx;x+=x&(-x)) t[x]=0;
}

void solve(int L,int r){
    if (L==r){
        chkmax(f[L],1);
        return;
    }
    int mid=(L+r)>>1;
    solve(L,mid);
    fr(i,L,r) num[i]=i;
    sort(num+L,num+mid+1,cmpc);
    sort(num+mid+1,num+r+1,cmpb);
    int pp1=L,pp2=mid+1;
    while(pp1<=mid&&pp2<=r){
        int p1=num[pp1],p2=num[pp2];
        if (c[p1]<=b[p2]){
            ADD(b[p1],f[p1]);
            pp1++;
        } else{
            chkmax(f[p2],QRY(a[p2])+1);
            pp2++;
        }
    }
    for(;pp2<=r;pp2++) chkmax(f[num[pp2]],QRY(a[num[pp2]])+1);
    fr(i,L,mid) DEL(b[num[i]]);
    solve(mid+1,r);
}

int main(){
    read(n);read(m);
    fr(i,1,n) read(a[i]),b[i]=c[i]=a[i];
    int x,y;
    fr(i,1,m){
        read(x);read(y);
        chkmin(a[x],y);
        chkmax(c[x],y);
    }
    solve(1,n);
    int ans=0;
    fr(i,1,n) chkmax(ans,f[i]);
    cout<<ans<<endl;
    return 0;
}

转载于:https://www.cnblogs.com/ymzqwq/p/11278980.html