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

题目：

我是超链接

题解：

Emmm。整体二分+树套树。
我们先把坐标离散，权值嘛，离不离散的无所谓，因为log1e9也是个常数
这里的交换操作相当于两个修改操作，每个修改操作也要拆成两个：把以前的减去，把新来的加上
那么现在的操作就是修改，询问就是查询第k小
我们依然按照权值二分，每次把小于等于mid值得操作执行，可以完成任务的询问（小于的数字比k多）放在左边
那么为什么不能用二维bit呢？空间啊空间啊，只能写动态加点的线段树了啊
注意一个坑就是不保证给的矩形的两个点的大小，要判断+交换
还有几个容易让你WA/T的问题。比如说你在记录交换操作的时候要记得把值也交换，因为ta可能还被换走
clear的操作比你把操作全都倒过来再加一遍相反数要快得多

代码：

#include <map>
#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=60005;
int sz,sj,ls[N*100],rs[N*100],s[N*10],sum[N*100],root[N*10],tmp[N*10],ans[N],val[N];
struct hh{int x,y,x1,y1,k,f,ans,id,wz;}cz[N*10],tmpl[N*10],tmpr[N*10];
map<int,int>hash;bool flag[N];
int read()
{
    int x=0;char ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch<='9' && ch>='0') x=x*10+ch-'0',ch=getchar();
    return x;
}
void change(int &now,int l,int r,int x,int v)
{
    if (!now) now=++sz;
    sum[now]+=v;
    if (l==r) return;
    int mid=(l+r)>>1;
    if (x<=mid) change(ls[now],l,mid,x,v);
    else change(rs[now],mid+1,r,x,v);
}
int qurry(int now,int l,int r,int lrange,int rrange)
{
    if (lrange<=l && rrange>=r) return sum[now];
    int mid=(l+r)>>1,ans=0;
    if (ls[now] && lrange<=mid) ans+=qurry(ls[now],l,mid,lrange,rrange);
    if (rs[now] && rrange>mid) ans+=qurry(rs[now],mid+1,r,lrange,rrange);
    return ans;
}
void add(int x,int y,int v)
{
    for (int i=x;i<=sj;i+=i&(-i)) change(root[i],1,sj,y,v);
}
int what(int x,int y,int z)
{
    int ans=0;
    for (int i=x;i>=1;i-=i&(-i))
      ans+=qurry(root[i],1,sj,y,z);
    return ans;
}
void clear(int x)
{
    if (ls[x]) clear(ls[x]);
    if (rs[x]) clear(rs[x]);
    sum[x]=ls[x]=rs[x]=0;
}
void qk(int x)
{
    for (int i=x;i<=sj;i+=i&(-i)) 
      if (root[i]) clear(root[i]),root[i]=0;
}
void ef(int l,int r,int a,int b)
{
    if (a>b) return;
    int mid=(l+r)>>1;
    int pa=0,pb=0,pz=a-1,sb=0;sz=0;
    for (int i=a;i<=b;i++)
      if (cz[i].id==1)
      {
        if (cz[i].k<=mid)
        {
            add(cz[i].x,cz[i].y,cz[i].f);
            tmpl[++pa]=cz[i];tmp[++sb]=cz[i].x;
        }else tmpr[++pb]=cz[i];
      }else
      {
        int t1=what(cz[i].x-1,cz[i].y,cz[i].y1);
        int t2=what(cz[i].x1,cz[i].y,cz[i].y1);
        int t=t2-t1;cz[i].ans=t;
        if (t>=cz[i].k) tmpl[++pa]=cz[i];
        else cz[i].k-=t,tmpr[++pb]=cz[i],cz[i].ans=0;
      }
    for (int i=1;i<=sb;i++) qk(tmp[i]);
    if (l==r)
    {
        for (int i=a;i<=b;i++)
          if (cz[i].id==2)
          {
            if (cz[i].ans>=cz[i].k) ans[cz[i].wz]=l; 
            else ans[cz[i].wz]=-1; 
          }
        return;
    }
    for (int i=1;i<=pa;i++) cz[++pz]=tmpl[i];
    for (int i=1;i<=pb;i++) cz[++pz]=tmpr[i];
    //省去sort 
    ef(l,mid,a,a+pa-1);
    ef(mid+1,r,a+pa,b);
}
int cmp(hh a,hh b){return a.wz<b.wz;}
int main()
{
    int n,m,x,y;n=read();m=read();int num=0,maxx=0;
    for (int i=1;i<=n;i++) 
    {
        ++num;cz[num].id=1;cz[num].f=1;
        cz[num].x=read()+1;cz[num].y=read()+1;cz[num].k=read();maxx=max(maxx,cz[num].k);
        s[++sj]=cz[num].x; s[++sj]=cz[num].y; val[i]=cz[num].k;
    }
    for (int i=1;i<=m;i++)
    {
        char st[10];scanf("%s",st);
        if (st[0]=='Q')
        {
            ++num;cz[num].id=2;cz[num].wz=i;flag[i]=1;
            cz[num].x=read()+1;cz[num].y=read()+1;cz[num].x1=read()+1;cz[num].y1=read()+1;cz[num].k=read();
            maxx=max(maxx,cz[num].k);
            if (cz[num].x>cz[num].x1) swap(cz[num].x,cz[num].x1);
            if (cz[num].y>cz[num].y1) swap(cz[num].y,cz[num].y1);
            s[++sj]=cz[num].x; s[++sj]=cz[num].y; s[++sj]=cz[num].x1; s[++sj]=cz[num].y1;
        }
        else
        {
            x=read()+1; y=read()+1;
            if (x==y||val[x]==val[y]) continue;
            ++num; cz[num].x=cz[x].x; cz[num].y=cz[x].y; cz[num].k=val[x]; cz[num].id=1; cz[num].f=-1;
            ++num; cz[num].x=cz[x].x; cz[num].y=cz[x].y; cz[num].k=val[y]; cz[num].id=1; cz[num].f=1;
            ++num; cz[num].x=cz[y].x; cz[num].y=cz[y].y; cz[num].k=val[y]; cz[num].id=1; cz[num].f=-1;
            ++num; cz[num].x=cz[y].x; cz[num].y=cz[y].y; cz[num].k=val[x]; cz[num].id=1; cz[num].f=1;
            swap(val[x],val[y]);
        }
    }
    sort(s+1,s+sj+1);
    int ts=unique(s+1,s+sj+1)-s-1; sj=ts;
    for (int i=1;i<=sj;i++) hash[s[i]]=i;
    for (int i=1;i<=num;i++)
    {
        if (cz[i].id==2) cz[i].x1=hash[cz[i].x1],cz[i].y1=hash[cz[i].y1];
        cz[i].x=hash[cz[i].x],cz[i].y=hash[cz[i].y];
    }
    ef(0,maxx,1,num);
    for (int i=1;i<=m;i++)
      if (flag[i])
      {
        if (ans[i]==-1) printf("It doesn't exist.\n");
        else printf("%d\n",ans[i]);
      }
}