[BZOJ3110][ZJOI2013]K大数查询（整体二分）

题目：

我是超链接

题解：

这道题目的暴力就是线段树套平衡树了，复杂度$O(nlog^3n)$，根本过不去啊
然后就是整体二分的第一题了
话说我一开始写线段树套平衡树的时候，发现线段树是在位置上二分，在查询有多少个比x大的时候是在平衡树上二分，那不如我们合成一个二分吧
有点类似cdq分治的思想
整体二分是指在二分答案的时候把询问和操作也跟着一起分到左边或者右边
只需要做一遍二分就可以求出所有询问的答案
定义过程Solve(l,r,x,y)，表示[x..y]之间的询问可能的答案在[l..r]之间。
终止条件是Solve(r,r,x,y)。这个时候[x..y]之间的询问的答案已经被唯一确定了。
初始调用Solve(1,Max,1,m)。
如何确定一个询问被分到左边还是右边？
在Solve(l,r,x,y)里，二分答案mid，按时间顺序扫描[x..y]之间的所有询问和操作
对于一个询问(L,R,K)，如果在区间[L,R]里加入的小于等于mid的数val大于等于K，说明这个询问的答案小于等于mid，分到左边
否则val < K，说明这个询问的答案大于mid，K->K-val，将这个询问分到右边
对于一个操作(L,R,K)，如果K<=mid分到左边，否则分到右边并且不处理这个操作。
设分到左边的一共有p1个，则接下来递归调用Solve(l,mid,x,x+p1-1)和Solve(mid+1,r,x+p1,y)
当Solve(l,r,x,y)中的[x,y]是空区间的时候退出函数
时间复杂度O(nlog^2n)
常数大大的优于树套树等复杂的数据结构

代码：

#include <cstdio>
#include <algorithm>
#define LL long long
using namespace std;
const int N=50005;
LL sum[N*4],delta[N*4],ans[N];int n,m;bool ql[N*4];
struct hh{int id,l,r,v,num,wz;}q[N];
int cmp(hh a,hh b){return a.num<b.num;}
void updata(int now){sum[now]=sum[now<<1]+sum[now<<1|1];}
void pushdown(int now,int l,int r,int mid)
{
    if (ql[now])
    {
        sum[now<<1]=sum[now<<1|1]=0;
        delta[now<<1]=delta[now<<1|1]=0;
        ql[now<<1]=ql[now<<1|1]=1;
        ql[now]=0;
    }
    if (delta[now])
    {
        sum[now<<1]+=delta[now]*(LL)(mid-l+1);
        sum[now<<1|1]+=delta[now]*(LL)(r-mid);
        delta[now<<1]+=delta[now];
        delta[now<<1|1]+=delta[now];
        delta[now]=0;
    }
}
void change(int now,int l,int r,int lrange,int rrange)
{
    if (lrange<=l && rrange>=r)
    {
        delta[now]++;
        sum[now]+=(LL)(r-l+1);
        return;
    }
    int mid=(l+r)>>1;pushdown(now,l,r,mid);
    if (lrange<=mid) change(now<<1,l,mid,lrange,rrange);
    if (rrange>mid) change(now<<1|1,mid+1,r,lrange,rrange);
    updata(now);
}
LL qurry(int now,int l,int r,int lrange,int rrange)
{
    if (lrange<=l && rrange>=r) return sum[now];
    int mid=(l+r)>>1;LL ans=0;pushdown(now,l,r,mid);
    if (lrange<=mid) ans+=qurry(now<<1,l,mid,lrange,rrange);
    if (rrange>mid) ans+=qurry(now<<1|1,mid+1,r,lrange,rrange);
    return ans;
}
int cmp1(hh a,hh b){return a.wz<b.wz;}
void ef(int l,int r,int a,int b)
{
    if (l==r)
    {
        for (int i=a;i<=b;i++) 
          if (q[i].id==2) ans[q[i].num]=l;
        return;
    }
    sum[1]=delta[1]=0;
    ql[1]=1;
    int mid=(l+r)>>1,pa=0,pb=b-a+1;
    for (int i=a;i<=b;i++)
      if (q[i].id==1)
      {
        if (q[i].v<=mid) q[i].wz=++pa;
        else change(1,1,n,q[i].l,q[i].r),q[i].wz=++pb;
      }
      else
      {
        LL t=qurry(1,1,n,q[i].l,q[i].r);
        if (t>=q[i].v) q[i].wz=++pb;//因为是求第k大 
        else q[i].v-=t,q[i].wz=++pa;
      }
    sort(q+a,q+b+1,cmp1);
    ef(l,mid,a,a+pa-1);
    ef(mid+1,r,a+pa,b);
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&q[i].id,&q[i].l,&q[i].r,&q[i].v);
        q[i].num=i;
    }
    ef(0,n,1,m);
    sort(q+1,q+m+1,cmp);
    for (int i=1;i<=m;i++) 
      if (q[i].id==2) printf("%lld\n",ans[i]);
}