【NOIP2016模拟7.11】排序

Description

Input

Output

拍完序后第k位上的数字

Sample Input

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

Sample Output

5

Data Constraint

30%：n,m<=300
100% n,m<=100000

Solution

对于30分，直接按照操作暴力即可，然而数据太水，可以水出60分。
对于100%显然是不能直接暴力的。可以想到二分答案。
分析复杂度，发现对于每次操作，必须不大于$log_2(n)$才能，那么想到了什么？线段树！
对于二分的答案，将大于这个的数变为1，小于的变为0。对于一个区间，从小到大排序就是把0堆在前面，1放在后面，从大到小排序则相反。这可以用线段树区间修改来做到。
最后再按照第k位看看是0还是1来修改区间。
注意线段树的归零

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1010000
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
int n,m,a[N],c[N],t[N*2][2],lazy[N*2];
struct noo{
    int o,x,y;
};
noo q[N];
void down(int v,int i,int j)
{
    int m=(i+j)/2;
    t[v*2][lazy[v]]=(m-i+1);t[v*2][1-lazy[v]]=0;
    t[v*2+1][lazy[v]]=(j-m);t[v*2+1][1-lazy[v]]=0;
    lazy[v*2]=lazy[v*2+1]=lazy[v];lazy[v]=-1;
}
void build(int v,int i,int j)
{
    lazy[v]=-1;
    if (i==j) {t[v][c[i]]=1;t[v][1-c[i]]=0;return;}
    build(v*2,i,(i+j)/2);build(v*2+1,(i+j)/2+1,j);
    t[v][0]=t[v*2][0]+t[v*2+1][0];
    t[v][1]=t[v*2][1]+t[v*2+1][1];
}
int find(int v,int i,int j,int x,int y)
{
    if(i==x&&j==y){return t[v][0];}
    int m=(i+j)/2;if(lazy[v]!=-1)down(v,i,j);
    if(y<=m) return find(v*2,i,m,x,y);
    else if(x>m) return find(v*2+1,m+1,j,x,y);
         else return find(v*2,i,m,x,m)+find(v*2+1,m+1,j,m+1,y);
}
void change(int v,int i,int j,int x,int y,int z)
{
    if(i==x&&j==y) {lazy[v]=z;t[v][z]=y-x+1;t[v][1-z]=0;return;}
    int m=(i+j)/2;if (lazy[v]!=-1) down(v,i,j);
    if(y<=m) change(v*2,i,m,x,y,z);
    else if(x>m) change(v*2+1,m+1,j,x,y,z);
         else change(v*2,i,m,x,m,z),change(v*2+1,m+1,j,m+1,y,z);
    t[v][0]=t[v*2][0]+t[v*2+1][0];
    t[v][1]=t[v*2][1]+t[v*2+1][1];
}
int main()
{
    freopen("4605.in","r",stdin);
//  freopen("4605.out","w",stdout);
    scanf("%d%d",&n,&m);
    fo(i,1,n) scanf("%d",&a[i]);
    fo(i,1,m)
    {
        int op,x,y;scanf("%d%d%d",&q[i].o,&q[i].x,&q[i].y);
    }
    int qq;scanf("%d",&qq);
    memset(lazy,-1,sizeof(lazy));
    int l=1,r=n;
    for(;l+1<r;)
    {
        int mid=(l+r)/2;
        fo(i,1,n) c[i]=a[i]>=mid?1:0;
        build(1,1,n);
        fo(i,1,m){
            int k=find(1,1,n,q[i].x,q[i].y);
            int jy=q[i].o;if(jy==1) k=(q[i].y-q[i].x+1)-k;
            if (k!=q[i].y-q[i].x+1&&k!=0) change(1,1,n,q[i].x,q[i].x+k-1,jy),change(1,1,n,q[i].x+k,q[i].y,1-jy);
        }
        int k=1-find(1,1,n,qq,qq);
        if(k==1) l=mid;else r=mid;
    }
    int mid=l;
    fo(i,1,n) c[i]=a[i]>=mid?1:0;
    build(1,1,n);
    fo(i,1,m){
        int k=find(1,1,n,q[i].x,q[i].y);
        int jy=q[i].o;if(jy==1) k=(q[i].y-q[i].x+1)-k;
        if (k!=q[i].y-q[i].x+1&&k!=0) change(1,1,n,q[i].x,q[i].x+k-1,jy),change(1,1,n,q[i].x+k,q[i].y,1-jy);
    }
    int k=1-find(1,1,n,qq,qq);  
    if(k==1) r=l;
    printf("%d",r);
}