[ZJOI2013]K大数查询——整体二分

题目描述

有N个位置，M个操作。操作有两种，每次操作如果是：

• 1 a b c：表示在第a个位置到第b个位置，每个位置加上一个数c
• 2 a b c：表示询问从第a个位置到第b个位置，第C大的数是多少。

思路

　　比较基础的整体二分。我们二分出$mid，对于值域[l,r]对应的操作[L,R]$，若为操作1，则考虑把$val>mid$的插入线段树中，表示比$mid$大的值的个数，若为操作2，先询问$[q[i].l,q[i].r]$中比$mid$大的值的个数，然后把当前询问填到左右区间再处理。讲的很简单，调过来整体二分原理的一些东西，，，毕竟这题还是比较板子的。

code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define I inline
#define smid (l+r>>1)
#define lch (now<<1)
#define rch (now<<1|1)
using namespace std;
const int N=;
typedef long long LL;
int n,m,tot;
LL ans[N];
struct segment
{
    LL sum,pls;
}sgt[N<<];
struct node
{
    int l,r,k,op,id;
}q[N<<],q1[N<<],q2[N<<];

I int read()
{
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}

I LL readll()
{
    LL x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}

I void pushup(int now)
{
    sgt[now].sum=sgt[lch].sum+sgt[rch].sum;
}

I void pushdown(int now,int l,int r)
{
    if(!sgt[now].pls)return;
    sgt[lch].pls+=sgt[now].pls;
    sgt[rch].pls+=sgt[now].pls;
    sgt[lch].sum+=(smid-l+)*sgt[now].pls;
    sgt[rch].sum+=(r-smid)*sgt[now].pls;
    sgt[now].pls=;
}

I void modify(int now,int l,int r,int x,int y,int val)
{
    if(x<=l&&r<=y)
    {
        sgt[now].pls+=val;
        sgt[now].sum+=(r-l+)*val;
        return;
    }
    pushdown(now,l,r);
    if(x<=smid)modify(lch,l,smid,x,y,val);
    if(smid<y)modify(rch,smid+,r,x,y,val);
    pushup(now);
}

I LL query(int now,int l,int r,int x,int y)
{
    if(x<=l&&r<=y)return sgt[now].sum;
    pushdown(now,l,r);
    LL res=;
    if(x<=smid)res+=query(lch,l,smid,x,y);
    if(smid<y)res+=query(rch,smid+,r,x,y);
    pushup(now);
    return res;
}

I void solve(int l,int r,int L,int R)
{
    if(L>R)return;
    if(l==r)
    {
        for(int i=L;i<=R;i++)if(q[i].op==)ans[q[i].id]=l;
        return;
    }
    int mid=(l+r>>),cnt1=,cnt2=;
    for(int i=L;i<=R;i++)
    {
        if(q[i].op==)
        {
            if(q[i].k>mid)
            modify(,,n,q[i].l,q[i].r,),q2[++cnt2]=q[i];
            else q1[++cnt1]=q[i];
        }
        else
        {
            LL tmp=query(,,n,q[i].l,q[i].r);
            if(q[i].k>tmp)q[i].k-=tmp,q1[++cnt1]=q[i];
            else q2[++cnt2]=q[i];
        }
    }
    for(int i=;i<=cnt2;i++)if(q2[i].op==)modify(,,n,q2[i].l,q2[i].r,-);
    for(int i=;i<=cnt1;i++)q[L+i-]=q1[i];
    for(int i=;i<=cnt2;i++)q[L+cnt1+i-]=q2[i];
    solve(l,mid,L,L+cnt1-);solve(mid+,r,L+cnt1,R);
}

int main()
{
    n=read();m=read();
    for(int i=;i<=m;i++)
    {
        q[i].op=read();q[i].l=read();q[i].r=read();q[i].k=readll();
        if(q[i].op==)q[i].id=++tot;
    }
    solve(-n,n,,m);
    for(int i=;i<=tot;i++)printf("%lld\n",ans[i]);
}