bzoj4552 [Tjoi2016&Heoi2016]排序 线段树+二分答案

Description

---

在2016年，佳媛姐姐喜欢上了数字序列。因而他经常研究关于序列的一些奇奇怪怪的问题，现在他在研究一个难题
，需要你来帮助他。这个难题是这样子的：给出一个1到n的全排列，现在对这个全排列序列进行m次局部排序，排
序分为两种：1:(0,l,r)表示将区间[l,r]的数字升序排序2:(1,l,r)表示将区间[l,r]的数字降序排序最后询问第q
位置上的数字。
1 <= n, m <= 10^5

Solution

---

可以二分一个答案，然后把数字分成大于mid和小于mid的两部分，这样就能用线段树区间维护排序操作了
感谢loj的数据

Code

---

#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))
const int N=200005;
int lazy[N<<2],c[N<<2];
int opt[N][3],a[N];
int n,m,pos;
int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}
void push_down(int now,int tl,int tr) {
    if (lazy[now]==-1) return ;
    int mid=(tl+tr)>>1;
    lazy[now<<1]=lazy[now<<1|1]=lazy[now];
    c[now<<1]=(mid-tl+1)*lazy[now];
    c[now<<1|1]=(tr-mid)*lazy[now];
    lazy[now]=-1;
}
void modify(int now,int tl,int tr,int l,int r,int v) {
    if (l>r) return ;
    if (tl==l&&tr==r) {
        lazy[now]=v;
        c[now]=(r-l+1)*v;
        return ;
    }
    push_down(now,tl,tr);
    int mid=(tl+tr)>>1;
    if (r<=mid) modify(now<<1,tl,mid,l,r,v);
    else if (l>mid) modify(now<<1|1,mid+1,tr,l,r,v);
    else {
        modify(now<<1,tl,mid,l,mid,v);
        modify(now<<1|1,mid+1,tr,mid+1,r,v);
    }
    c[now]=c[now<<1]+c[now<<1|1];
}
int query(int now,int tl,int tr,int l,int r) {
    if (tl==l&&tr==r) return c[now];
    push_down(now,tl,tr);
    int mid=(tl+tr)>>1;
    if (r<=mid) return query(now<<1,tl,mid,l,r);
    else if (l>mid) return query(now<<1|1,mid+1,tr,l,r);
    else return query(now<<1,tl,mid,l,mid)+query(now<<1|1,mid+1,tr,mid+1,r);
}
int solve(int now) {
    fill(c,0);
    fill(lazy,-1);
    rep(i,1,n) modify(1,1,n,i,i,a[i]>=now);
    rep(i,1,m) {
        int ret=query(1,1,n,opt[i][1],opt[i][2]);
        if (opt[i][0]==0) {
            modify(1,1,n,opt[i][1],opt[i][2]-ret,0);
            modify(1,1,n,opt[i][2]-ret+1,opt[i][2],1);
//          printf("%d %d %d %d\n",opt[i][1],opt[i][2]-ret,opt[i][2]-ret+1,opt[i][2]);
        } else {
            modify(1,1,n,opt[i][1],opt[i][1]+ret-1,1);
            modify(1,1,n,opt[i][1]+ret,opt[i][2],0);
//          printf("%d %d %d %d\n",opt[i][1],opt[i][1]+ret-1,opt[i][1]+ret,opt[i][2]);
        }
    }
    return query(1,1,n,pos,pos);
}
int main(void) {
    n=read(); m=read();
    rep(i,1,n) a[i]=read();
    rep(i,1,m) opt[i][0]=read(),opt[i][1]=read(),opt[i][2]=read();
    pos=read();
    int ans=0;
    int l=1,r=n;
    while (l<=r) {
        int mid=(l+r)>>1;
        if (solve(mid)) l=mid+1,ans=mid;
        else r=mid-1;
    }
    printf("%d\n", ans);
    return 0;
}