洛谷 P3332 [ZJOI2013]K大数查询 || bzoj3110

用树套树就很麻烦，用整体二分就成了裸题。。。。

错误：

1.尝试线段树套平衡树，码农，而且n*log^3(n)慢慢卡反正我觉得卡不过去

2.线段树pushdown写错。。。加法tag对于区间和的更新应该要乘上区间长度的

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
struct Q
{
    LL type,a,b,c,num;
}q[50100],qt1[50100],qt2[50100];
LL ans[50100],qnum;
LL n,m;
namespace SegT
{
#define mid (l+((r-l)>>1))
#define lc (num<<1)
#define rc (num<<1|1)
LL dat[200100],addv[200100];
LL L,R,x;
void pd(LL l,LL r,LL num)
{
    if(addv[num])
    {
        dat[lc]+=(mid-l+1)*addv[num];addv[lc]+=addv[num];
        dat[rc]+=(r-mid)*addv[num];addv[rc]+=addv[num];
        addv[num]=0;
    }
}
void _addx(LL l,LL r,LL num)
{
    if(L<=l&&r<=R)    {dat[num]+=(r-l+1)*x;addv[num]+=x;return;}
    pd(l,r,num);
    if(L<=mid)    _addx(l,mid,lc);
    if(mid<R)    _addx(mid+1,r,rc);
    dat[num]=dat[lc]+dat[rc];
}
LL _query(LL l,LL r,LL num)
{
    if(L<=l&&r<=R)    return dat[num];
    pd(l,r,num);LL ans=0;
    if(L<=mid)    ans+=_query(l,mid,lc);
    if(mid<R)    ans+=_query(mid+1,r,rc);
    return ans;
}
void addx(LL l,LL r,LL d)
{
    L=l;R=r;x=d;_addx(1,n,1);
}
LL query(LL l,LL r)
{
    L=l;R=r;return _query(1,n,1);
}
#undef mid
#undef lc
#undef rc
}
void solve(LL lp,LL rp,LL l,LL r)
{
    if(lp>rp)    return;
    LL i;
    if(l==r)
    {
        for(i=lp;i<=rp;i++)
            if(q[i].type==2)
                ans[q[i].num]=l;
        return;
    }
    LL mid=l+((r-l)>>1),tlen1=0,tlen2=0,t;
    for(i=lp;i<=rp;i++)
    {
        if(q[i].type==1)
        {
            if(q[i].c>mid)
            {
                SegT::addx(q[i].a,q[i].b,1);
                qt1[++tlen1]=q[i];
            }
            else
                qt2[++tlen2]=q[i];
        }
        else
        {
            t=SegT::query(q[i].a,q[i].b);
            if(t>=q[i].c)
                qt1[++tlen1]=q[i];
            else
                qt2[++tlen2]=q[i],qt2[tlen2].c-=t;
        }
    }
    for(i=lp;i<=rp;i++)
        if(q[i].type==1&&q[i].c>mid)
            SegT::addx(q[i].a,q[i].b,-1);
    memcpy(q+lp,qt1+1,sizeof(Q)*tlen1);
    memcpy(q+lp+tlen1,qt2+1,sizeof(Q)*tlen2);
    solve(lp,lp+tlen1-1,mid+1,r);
    solve(lp+tlen1,rp,l,mid);
}
int main()
{
    LL i;
    scanf("%lld%lld",&n,&m);
    for(i=1;i<=m;i++)
    {
        scanf("%lld%lld%lld%lld",&q[i].type,&q[i].a,&q[i].b,&q[i].c);
        if(q[i].type==2)    q[i].num=++qnum;
    }
    solve(1,m,-50000,50000);
    for(i=1;i<=qnum;i++)    printf("%lld
",ans[i]);
    return 0;
}

---