=== ===
这里放传送门
=== ===
题解
第一次见到都整体二分了还要加上树套树这种鬼畜东西的题目。。
不过写起来真是神清气爽啊蛤蛤蛤(逃
关键就是想到把交换操作看成加入一个点再删除一个点,然后初始的那些点都看成是加入操作,就可以搞整体二分了。
二分一个答案然后把小于这个答案的点都加进去,然后考察每个询问,看这个询问的矩形内点的个数是不是大于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;
}