bzoj 4592: [Shoi2015]脑洞治疗仪

线段树维护最长连续段的长度,然后暴力修改,暴力查询,这个复杂度应该是nlogn的。然而我调了好久。。。
   
   
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
 
#define ll long long
#define inf 1e9
#define eps 1e-8
#define md
#define N 200010
using namespace std;
struct Tr { int L,R,sum,ans,sz;} tr[4*N];
int tag[4*N];
 
void tagit(int i,int d)
{
if (d==0) tr[i].L=tr[i].R=tr[i].ans=tr[i].sz,tr[i].sum=0;
else tr[i].L=tr[i].R=tr[i].ans=0,tr[i].sum=tr[i].sz;
tag[i]=d;
}
 
void release(int i)
{
if (tag[i]!=-1)
{
int d=tag[i]; tag[i]=-1;
tagit(i<<1,d); tagit(i<<1|1,d);
}
}
 
Tr update(Tr l,Tr r)
{
Tr ans;
ans.L=l.L+(l.sum==0?r.L:0);
ans.R=r.R+(r.sum==0?l.R:0);
ans.sz=l.sz+r.sz;
ans.sum=l.sum+r.sum;
ans.ans=max(l.R+r.L,max(l.ans,r.ans));
return ans;
}
 
void build(int i,int l,int r)
{
tr[i].L=tr[i].R=tr[i].ans=0; tr[i].sz=tr[i].sum=r-l+1; tag[i]=-1;
if (l==r) return;
int mid=(l+r)>>1;
build(i<<1,l,mid); build(i<<1|1,mid+1,r);
}
 
void modify_0(int i,int l,int r,int ql,int qr)
{
if (ql<=l&&r<=qr) { tagit(i,0); return;}
release(i);
int mid=(l+r)>>1;
if (ql<=mid) modify_0(i<<1,l,mid,ql,qr);
if (mid+1<=qr) modify_0(i<<1|1,mid+1,r,ql,qr);
tr[i]=update(tr[i<<1],tr[i<<1|1]);
}
 
int modify_1(int i,int l,int r,int ql,int qr,int k)
{
//printf("modify_1 [%d %d] now %d %d %d\n",ql,qr,l,r,k);
if (ql<=l&&r<=qr)
{
int t=tr[i].sz-tr[i].sum;
if (t<=k) { tagit(i,1); return t;}
release(i);
int mid=(l+r)>>1;
int p=modify_1(i<<1,l,mid,ql,qr,k);
if (p<k) p+=modify_1(i<<1|1,mid+1,r,ql,qr,k-p);
tr[i]=update(tr[i<<1],tr[i<<1|1]);
return p;
}
release(i);
int mid=(l+r)>>1,t=0;
if (qr<=mid) t=modify_1(i<<1,l,mid,ql,qr,k);
else if (mid+1<=ql) t=modify_1(i<<1|1,mid+1,r,ql,qr,k);
else
{
t=modify_1(i<<1,l,mid,ql,qr,k);
if (t<k) t+=modify_1(i<<1|1,mid+1,r,ql,qr,k-t);
}
tr[i]=update(tr[i<<1],tr[i<<1|1]);
return t;
}
 
Tr query(int i,int l,int r,int ql,int qr)
{
if (ql<=l&&r<=qr) return tr[i];
release(i);
int mid=(l+r)>>1;
if (qr<=mid) return query(i<<1,l,mid,ql,qr);
if (mid+1<=ql) return query(i<<1|1,mid+1,r,ql,qr);
return update(query(i<<1,l,mid,ql,qr),query(i<<1|1,mid+1,r,ql,qr));
}
 
int main()
{
int n,m;
scanf("%d%d",&n,&m);
build(1,1,n);
for (int i=1;i<=m;i++)
{
int ty;
scanf("%d",&ty);
if (ty==0)
{
int l,r;
scanf("%d%d",&l,&r);
modify_0(1,1,n,l,r);
}
else if (ty==1)
{
int l0,r0,l1,r1;
scanf("%d%d%d%d",&l0,&r0,&l1,&r1);
Tr x=query(1,1,n,l0,r0);
modify_0(1,1,n,l0,r0);
if (x.sum) modify_1(1,1,n,l1,r1,x.sum);
}
else
{
int l,r;
scanf("%d%d",&l,&r);
Tr x=query(1,1,n,l,r);
printf("%d\n",x.ans);
}
}
return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值