bzoj 4552: [Tjoi2016&Heoi2016]排序——二分+线段树

Description

在2016年，佳媛姐姐喜欢上了数字序列。因而他经常研究关于序列的一些奇奇怪怪的问题，现在他在研究一个难题

，需要你来帮助他。这个难题是这样子的：给出一个1到n的全排列，现在对这个全排列序列进行m次局部排序，排

序分为两种：1:(0,l,r)表示将区间[l,r]的数字升序排序2:(1,l,r)表示将区间[l,r]的数字降序排序最后询问第q

位置上的数字。

Input

输入数据的第一行为两个整数n和m。n表示序列的长度，m表示局部排序的次数。1 <= n, m <= 10^5第二行为n个整

数，表示1到n的一个全排列。接下来输入m行，每一行有三个整数op, l, r, op为0代表升序排序，op为1代表降序

排序, l, r 表示排序的区间。最后输入一个整数q，q表示排序完之后询问的位置, 1 <= q <= n。1 <= n <= 10^5

，1 <= m <= 10^5

Output

输出数据仅有一行，一个整数，表示按照顺序将全部的部分排序结束后第q位置上的数字。

Sample Input

6 3
1 6 2 5 3 4
0 1 4
1 3 6
0 2 4
3

Sample Output

5

————————————————————————————————

二分p位置的值，把大于mid的数改为1，小于等于mid的数改为0，

变成01串后就可以用线段树实现排序了，像降序升序什么的操作就把1和0各堆到一边就可以辣

排序后如果p的位置上的数为0，说明答案比mid小，如果为1，说明答案比mid大。

至于为什么可以用二分呢

你想如果p位置上是1，说明mid较小，v[p]>mid，所以把v[p]给标记成了1。

如果p位置上是0，就是把v[p]<=mid，所以把v[p]标记成了0，

但是这样还有一些大于v[p]的位置也是0，所以继续往小的地方逼近答案。

满足单调所以就可以这么写辣

#include<cstdio>
#include<cstring>
#include<algorithm>
const int M=5e5+;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
int p,n,m,v[M];
int lx[M],rx[M],op[M];
int L,R,mx;
struct pos{int s,h[];}tr[M];
void up(int x){tr[x].s=tr[x<<].s+tr[x<<^].s;}
void down(int x,int l,int r){
    if(l==r) return ;
    int ls=x<<,rs=x<<^,mid=(l+r)>>;
    if(tr[x].h[]){
        tr[x].h[]=;
        tr[ls].h[]=tr[rs].h[]=;
        tr[ls].s=tr[rs].s=;
        tr[ls].h[]=tr[rs].h[]=;
    }
    else if(tr[x].h[]){
        tr[x].h[]=; tr[ls].h[]=tr[rs].h[]=;
        tr[ls].s=mid-l+; tr[rs].s=r-mid;
        tr[ls].h[]=tr[rs].h[]=;
    }
}
void build(int x,int l,int r){
    tr[x].h[]=tr[x].h[]=;
    if(l==r){
        tr[x].s=(v[l]>mx);
        return ;
    }
    int mid=(l+r)>>;
    build(x<<,l,mid);
    build(x<<^,mid+,r);
    up(x);
}
void modify(int x,int l,int r,int v){
    if(L<=l&&r<=R){
        tr[x].h[v]=;
        tr[x].h[v^]=;
        tr[x].s=(r-l+)*v;
        return ;
    }
    down(x,l,r);
    int mid=(l+r)>>;
    if(L<=mid) modify(x<<,l,mid,v);
    if(R>mid)  modify(x<<^,mid+,r,v);
    up(x);
}
int query(int x,int l,int r){
    if(L<=l&&r<=R) return tr[x].s;
    down(x,l,r);
    int mid=(l+r)>>,sum=;
    if(L<=mid) sum+=query(x<<,l,mid);
    if(R>mid)  sum+=query(x<<^,mid+,r);
    return sum;
}
bool check(int k){
    mx=k; build(,,n);
    for(int i=;i<=m;i++){
        L=lx[i]; R=rx[i];
        int ly=query(,,n);
        if(op[i]){
            if(ly) L=lx[i],R=lx[i]+ly-,modify(,,n,);
            if(lx[i]+ly<=rx[i]) L=lx[i]+ly,R=rx[i],modify(,,n,);
        }
        else{
            if(lx[i]<=rx[i]-ly) L=lx[i],R=rx[i]-ly,modify(,,n,);
            if(ly) L=rx[i]-ly+,R=rx[i],modify(,,n,);
        }
    }
    L=R=p;
    return !query(,,n);
}
int main(){
    n=read(); m=read();
    for(int i=;i<=n;i++) v[i]=read();
    for(int i=;i<=m;i++) op[i]=read(),lx[i]=read(),rx[i]=read();
    p=read();
    int l=,r=n;
    while(l<r){
        int mid=(l+r)>>;
        if(check(mid)) r=mid;
        else l=mid+;
    }printf("%d\n",l);
    return ;
}