[BZOJ2674]Attack（整体二分+树状数组套线段树）

=== ===

这里放传送门

=== ===

题解

第一次见到都整体二分了还要加上树套树这种鬼畜东西的题目。。
不过写起来真是神清气爽啊蛤蛤蛤（逃
关键就是想到把交换操作看成加入一个点再删除一个点，然后初始的那些点都看成是加入操作，就可以搞整体二分了。
二分一个答案然后把小于这个答案的点都加进去，然后考察每个询问，看这个询问的矩形内点的个数是不是大于k。然后就要求有一个东西能支持插入点然后查询一个矩形内点的个数，那个鬼畜的数据结构就是来搞这个问题的。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 2000000000
using namespace std;
int n,m,X[60010],Y[60010],Z[60010],cnt,val[60010],V,hx[60010],hy[60010],tx,ty,size,s[60010],ans[10010];
struct oper{
    int type,x,y,xx,yy,k,num,id;
    oper(int T=0,int X=0,int Y=0,int XX=0,int YY=0,int K=0,int ID=0){
        type=T;x=X;y=Y;xx=XX;yy=YY;k=K;id=ID;num=0;
    }
}q[200010],tmp[200010];
struct segtree{
    int l,r,val;
    segtree(){l=r=val=0;}
}t[20000010];
void insert(int &i,int l,int r,int x,int v){
    if (i==0){i=++size;t[i]=segtree();}
    t[i].val+=v;
    if (l==r) return;
    int mid=(l+r)>>1;
    if (x<=mid) insert(t[i].l,l,mid,x,v);
    else insert(t[i].r,mid+1,r,x,v);
}
int ask(int i,int l,int r,int left,int right){
    if (i==0) return 0;
    if (left<=l&&right>=r) return t[i].val;
    int mid=(l+r)>>1,ans=0;
    if (left<=mid) ans+=ask(t[i].l,l,mid,left,right);
    if (right>mid) ans+=ask(t[i].r,mid+1,r,left,right);
    return ans;
}
int lowbit(int i){return i&(-i);}
void Change(int i,int y,int d){
    while (i<=tx){
        insert(s[i],1,ty,y,d);
        i+=lowbit(i);
    }
}
int Ask(int i,int y,int yy){
    int ans=0;
    while (i!=0){
        ans+=ask(s[i],1,ty,y,yy);
        i-=lowbit(i);
    }
    return ans;
}
int Query(int x,int y,int xx,int yy){
    return Ask(xx,y,yy)-Ask(x-1,y,yy);
}
void Clear(int i,int y){while (i<=tx){s[i]=0;i+=lowbit(i);}}
void Divide(int l,int r,int x,int y){
    if (x>y) return;
    if (l==r){
        for (int i=x;i<=y;i++)
          if (q[i].type==3) ans[q[i].id]=val[l];
        return;
    }
    int mid=(l+r)>>1,r1=0,r2=0,lim,p1,p2;
    lim=val[mid];
    for (int i=x;i<=y;i++)
      if (q[i].type==3){
          int g=Query(q[i].x,q[i].y,q[i].xx,q[i].yy);
          if (g>=q[i].k){++r1;q[i].num=1;}
          else{++r2;q[i].num=2;q[i].k-=g;}
      }else{
          int g=(q[i].type==1)?1:-1;
          if (q[i].k<=lim){
              Change(q[i].x,q[i].y,g);
              ++r1;q[i].num=1;
          }else{++r2;q[i].num=2;}
      }
    size=0;
    for (int i=x;i<=y;i++)
      if (q[i].type!=3&&q[i].k<=lim)
        Clear(q[i].x,q[i].y);
    p1=x;p2=x+r1;
    for (int i=x;i<=y;i++)
      if (q[i].num==1) tmp[p1++]=q[i];
      else tmp[p2++]=q[i];
    for (int i=x;i<=y;i++) q[i]=tmp[i];
    Divide(l,mid,x,x+r1-1);
    Divide(mid+1,r,x+r1,y);
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++){
        scanf("%d%d%d",&X[i],&Y[i],&Z[i]);
        hx[++tx]=X[i];hy[++ty]=Y[i];val[++V]=Z[i];
    }
    sort(hx+1,hx+tx+1);tx=unique(hx+1,hx+tx+1)-hx-1;
    sort(hy+1,hy+ty+1);ty=unique(hy+1,hy+ty+1)-hy-1;
    sort(val+1,val+V+1);V=unique(val+1,val+V+1)-val-1;
    ++V;val[V]=inf;
    for (int i=1;i<=n;i++){
        X[i]=lower_bound(hx+1,hx+tx+1,X[i])-hx;
        Y[i]=lower_bound(hy+1,hy+ty+1,Y[i])-hy;
        q[++cnt]=oper(1,X[i],Y[i],0,0,Z[i],0);
    }
    for (int i=1;i<=m;i++){
        int x,y,xx,yy,k;
        char s[10];
        scanf("%s",s);
        if (s[0]=='Q'){
            scanf("%d%d%d%d%d",&x,&y,&xx,&yy,&k);
            if (x>xx) swap(x,xx);
            if (y>yy) swap(y,yy);
            x=lower_bound(hx+1,hx+tx+1,x)-hx;
            y=lower_bound(hy+1,hy+ty+1,y)-hy;
            xx=lower_bound(hx+1,hx+tx+1,xx+1)-hx-1;
            yy=lower_bound(hy+1,hy+ty+1,yy+1)-hy-1;
            q[++cnt]=oper(3,x,y,xx,yy,k,i);
        }else{
            scanf("%d%d",&x,&y);++x;++y;
            if (x==y) continue;
            q[++cnt]=oper(2,X[x],Y[x],0,0,Z[x],i);
            q[++cnt]=oper(1,X[x],Y[x],0,0,Z[y],i);
            q[++cnt]=oper(2,X[y],Y[y],0,0,Z[y],i);
            q[++cnt]=oper(1,X[y],Y[y],0,0,Z[x],i);
            swap(Z[x],Z[y]);
        }
    }
    memset(ans,-1,sizeof(ans));
    Divide(1,V,1,cnt);
    for (int i=1;i<=m;i++)
      if (ans[i]!=-1)
        if (ans[i]==inf) printf("It doesn't exist.\n");
        else printf("%d\n",ans[i]);
    return 0;
}