P2572 [SCOI2010]序列操作

线段树神仙操作==珂朵莉树基本操作???

珂朵莉树是不可能的,这辈子只会码线段树,只有线段树神仙操作才刺激,debug之后AC才最快乐

这题我折腾了半个下午加半个晚上。维护的东西太多了。

如果没有区间反转,这题很简单,但是有反转,所以既要维护1,又要维护0。

tot记录区间中1的个数

len记录区间长度

l记录区间左端点

r记录区间右端点

sum1表示区间中最长连续的1的个数

lmax1表示区间中以左端点为起点从左往右的连续1的个数

rmax1表示区间中以右端点为起点从右往左的连续1的个数

0同理

laz1表示区间覆盖,0表示无操作,1表示0覆盖,2表示1覆盖

laz2表示区间反转,0表示无操作,1表示反转

定义了这么多就够了,然后考虑建树和update维护

update:

1.tot就直接加就行

2.lmax1更新要看一下左区间是否全为1,若全为1,则应加上右区间lmax1

3.rmax1更新要看一下右区间是否全为1,若全为1,则应加上左区间rmax1

4.sum1更新有三个来源,左区间sum1,右区间sum1,左区间rmax1+右区间lmax1

0维护的同理。

懒得写那多if,else,所以这里直接采用了三目运算符

build:

1.边找边求len

2.递归到叶子结点,看原序列中是1还是0,去更新叶子结点的东西

3.三目运算符对1和0的更新相反

4.最后update,由叶子结点回溯,维护整个大区间

然后基本的建树完成了,先不想那些个操作,把main框架写出,

change可以分标记表示修改操作,一个就够了

ask返回的东西不同,要分两个

然后main函数也搞定了

开始本题比较难理解的一块:down和laz标记。

1.如果之前有反转标记,但是又有覆盖标记,之前的反转标记就没用了

,所以当要打覆盖标记时,反转标记是可以清空的

2.区间覆盖的话就是那些东西等于len或0的问题

3.区间反转把0和1的东西交换就行了

4.down完标记记得清空

ask1返回tot就行

ask2分开找区间的话,还有可能是左区间的rmax1+右区间的lmax1,

这个还要取min,因为可能分到两个区间的长度到不了lmax1或rmax1

值得注意的是直接下放标记可能叶子结点也下放,可能会越界,所以 开了八倍空间

#include<iostream>
#include<cstdio>
#define R register
#define lson k<<1,l,mid   //减少码量或者为了好看
#define rson k<<1|1,mid+1,r
#define ls k<<1
#define rs k<<1|1
#define mid ((l+r)>>1)
#define sum1(rt) tr[rt].sum1
#define sum0(rt) tr[rt].sum0
#define lmax1(rt) tr[rt].lmax1
#define lmax0(rt) tr[rt].lmax0
#define rmax1(rt) tr[rt].rmax1
#define rmax0(rt) tr[rt].rmax0
#define len(rt) tr[rt].len
#define laz1(rt) tr[rt].laz1
#define laz2(rt) tr[rt].laz2
#define tot(rt) tr[rt].tot
using namespace std;
const int maxn=200005;
int n,m,a[maxn];
struct node{
    int sum1,lmax1,rmax1;
    int tot,len,laz1,laz2;
    int sum0,lmax0,rmax0;
}tr[maxn<<2];
inline int read(){
    R int s=0,w=1;
    R char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){s=s*10+ch-'0';ch=getchar();}
    return s*w;
}
inline void update(R int k){
    lmax1(k)=(sum1(ls)==len(ls))?len(ls)+lmax1(rs):lmax1(ls);
    rmax1(k)=(sum1(rs)==len(rs))?len(rs)+rmax1(ls):rmax1(rs);
    sum1(k)=max(max(sum1(ls),sum1(rs)),rmax1(ls)+lmax1(rs));
    lmax0(k)=(sum0(ls)==len(ls))?len(ls)+lmax0(rs):lmax0(ls);
    rmax0(k)=(sum0(rs)==len(rs))?len(rs)+rmax0(ls):rmax0(rs);
    sum0(k)=max(max(sum0(ls),sum0(rs)),rmax0(ls)+lmax0(rs));
    tot(k)=tot(ls)+tot(rs);
}
void build(R int k,R int l,R int r){
    len(k)=r-l+1;
  if(l==r){ tot(k)=a[l]; lmax1(k)=rmax1(k)=sum1(k)=((a[l]==1)?1:0); lmax0(k)=rmax0(k)=sum0(k)=((a[l]==1)?0:1); return ; } build(lson);build(rson); update(k); } inline void down(R int k){ if(laz1(k)==1){ sum0(ls)=lmax0(ls)=rmax0(ls)=len(ls); tot(ls)=sum1(ls)=lmax1(ls)=rmax1(ls)=0; laz1(ls)=1;laz2(ls)=0; sum0(rs)=lmax0(rs)=rmax0(rs)=len(rs); tot(rs)=sum1(rs)=lmax1(rs)=rmax1(rs)=0; laz1(rs)=1;laz2(rs)=0; laz1(k)=0;laz2(k)=0; } if(laz1(k)==2){ tot(ls)=sum1(ls)=lmax1(ls)=rmax1(ls)=len(ls); sum0(ls)=lmax0(ls)=rmax0(ls)=0; laz1(ls)=2;laz2(ls)=0; tot(rs)=sum1(rs)=lmax1(rs)=rmax1(rs)=len(rs); sum0(rs)=lmax0(rs)=rmax0(rs)=0; laz1(rs)=2;laz2(rs)=0; laz1(k)=0;laz2(k)=0; } if(laz2(k)){ tot(ls)=len(ls)-tot(ls); swap(sum0(ls),sum1(ls)); swap(lmax0(ls),lmax1(ls)); swap(rmax0(ls),rmax1(ls)); if(laz1(ls)==1)laz1(ls)=2; else if(laz1(ls)==2)laz1(ls)=1; else laz2(ls)^=1; tot(rs)=len(rs)-tot(rs); swap(sum0(rs),sum1(rs)); swap(lmax0(rs),lmax1(rs)); swap(rmax0(rs),rmax1(rs)); if(laz1(rs)==1)laz1(rs)=2; else if(laz1(rs)==2)laz1(rs)=1; else laz2(rs)^=1; laz2(k)=0; } } void change(R int k,R int l,R int r,R int x,R int y,R int z){ // printf("?%d %d? ",l,r); down(k); if(l==x&&y==r){ if(z==1){ sum0(k)=lmax0(k)=rmax0(k)=len(k); tot(k)=sum1(k)=lmax1(k)=rmax1(k)=0; laz1(k)=1; laz2(k)=0; } else if(z==2){ tot(k)=sum1(k)=lmax1(k)=rmax1(k)=len(k); sum0(k)=lmax0(k)=rmax0(k)=0; laz1(k)=2; laz2(k)=0; } else { tot(k)=len(k)-tot(k); swap(sum0(k),sum1(k)); swap(lmax0(k),lmax1(k)); swap(rmax0(k),rmax1(k)); laz2(k)^=1; } return ; } if(y<=mid)change(lson,x,y,z); else if(x>mid)change(rson,x,y,z); else change(lson,x,mid,z),change(rson,mid+1,y,z); update(k); } int ask1(R int k,R int l,R int r,R int x,R int y){ down(k); if(l==x&&y==r){ return tot(k); } if(y<=mid)return ask1(lson,x,y); else if(x>mid)return ask1(rson,x,y); else return ask1(lson,x,mid)+ask1(rson,mid+1,y); } int ask2(R int k,R int l,R int r,R int x,R int y){ down(k); if(l==x&&y==r){ return sum1(k); } if(y<=mid)return ask2(lson,x,y); else if(x>mid)return ask2(rson,x,y); else return max(max(ask2(lson,x,mid),ask2(rson,mid+1,y)),min(lmax1(rs),y-mid)+min(rmax1(ls),mid-x+1)); } int main(){ n=read();m=read(); for(R int i=1;i<=n;++i) a[i]=read(); build(1,1,n); int opt,a,b; while(m--){ opt=read();a=read()+1;b=read()+1; if(opt==0)change(1,1,n,a,b,1); else if(opt==1)change(1,1,n,a,b,2); else if(opt==2)change(1,1,n,a,b,3); else if(opt==3)printf("%d ",ask1(1,1,n,a,b)); else if(opt==4)printf("%d ",ask2(1,1,n,a,b)); } return 0; }
原文地址:https://www.cnblogs.com/sanjinliushi/p/11537162.html