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

时间:2021-06-07 12:54:31

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

错误:

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[],qt1[],qt2[];
LL ans[],qnum;
LL n,m;
namespace SegT
{
#define mid (l+((r-l)>>1))
#define lc (num<<1)
#define rc (num<<1|1)
LL dat[],addv[];
LL L,R,x;
void pd(LL l,LL r,LL num)
{
if(addv[num])
{
dat[lc]+=(mid-l+)*addv[num];addv[lc]+=addv[num];
dat[rc]+=(r-mid)*addv[num];addv[rc]+=addv[num];
addv[num]=;
}
}
void _addx(LL l,LL r,LL num)
{
if(L<=l&&r<=R) {dat[num]+=(r-l+)*x;addv[num]+=x;return;}
pd(l,r,num);
if(L<=mid) _addx(l,mid,lc);
if(mid<R) _addx(mid+,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=;
if(L<=mid) ans+=_query(l,mid,lc);
if(mid<R) ans+=_query(mid+,r,rc);
return ans;
}
void addx(LL l,LL r,LL d)
{
L=l;R=r;x=d;_addx(,n,);
}
LL query(LL l,LL r)
{
L=l;R=r;return _query(,n,);
}
#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==)
ans[q[i].num]=l;
return;
}
LL mid=l+((r-l)>>),tlen1=,tlen2=,t;
for(i=lp;i<=rp;i++)
{
if(q[i].type==)
{
if(q[i].c>mid)
{
SegT::addx(q[i].a,q[i].b,);
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==&&q[i].c>mid)
SegT::addx(q[i].a,q[i].b,-);
memcpy(q+lp,qt1+,sizeof(Q)*tlen1);
memcpy(q+lp+tlen1,qt2+,sizeof(Q)*tlen2);
solve(lp,lp+tlen1-,mid+,r);
solve(lp+tlen1,rp,l,mid);
}
int main()
{
LL i;
scanf("%lld%lld",&n,&m);
for(i=;i<=m;i++)
{
scanf("%lld%lld%lld%lld",&q[i].type,&q[i].a,&q[i].b,&q[i].c);
if(q[i].type==) q[i].num=++qnum;
}
solve(,m,-,);
for(i=;i<=qnum;i++) printf("%lld\n",ans[i]);
return ;
}