bzoj 1858: [Scoi2010]序列操作 || 洛谷 P2572

记一下：线段树占空间是$2^{ceil(log2(n))+1}$

这个就是一个线段树区间操作题，各种标记的设置、转移都很明确，只要熟悉这类题应该说是没有什么难度的。

由于对某区间set之后该区间原先待进行的取反操作失效（被覆盖），因此规定tag同时存在时set的标记先进行操作，这样对区间加上set标记时要去掉原有的取反标记。

对于操作4，线段树上每个区间维护6个值，分别表示：该区间最多的连续0/1，从左侧数起/从右侧数起/其中任意位置。题目中不问连续0，为什么连续0也要维护呢？这是为了在进行取反操作时，O(1)完成标记的维护（直接交换对0、1维护的3个值就行了）。

对于操作3，线段树上每个区间维护该区间1的个数。此时不需要额外记录0的个数，0的个数可以由区间总数-1的个数得到。

规定：set的标记为-1时表示不需要set操作，为0/1时表示需要set为该标记的值。

注意一点：对于取反的tag，每次取反操作的时候都是将tag取反，而不是置为1，不然在同一个位置取反两次就不对了。

各个标记含义同普通线段树：当前节点已经操作，其所有（任意级的）子节点待操作。维护的时候转移都是递推，细的就不讲了，有点多，要细心比对。

然而...

感觉自己药丸啊...调了三个小时才调出来，各种各样诡异的错误...如果省选考这个肯定没法做了..

幸好这个样例够强。。有一些奇怪的错误通过样例就查出来了，过了样例就1A了，不然真的不知道要调到什么时候。。。

note：以下程序中打了注释的语句都是错过的。。。

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define mid ((l+r)>>1)
 #define lc (num<<1)
 #define rc (num<<1|1)
 using namespace std;
 struct Th
 {
     int a,b,c;
     bool d;//是否全部是1
     int len;
     Th(int aa=,int bb=,int cc=,bool dd=,int ll=):a(aa),b(bb),c(cc),d(dd),len(ll){}
 };
 int num1[],maxnum[][][],settag[];
 bool revtag[];
 int n,m;
 //定义tag同时存在时settag先进行操作
 int a[];
 int L,R,x;
 //maxnum[][0/1][0/1/2]表示最多的连续0/1,接左侧/接右侧/中间
 void pd(int l,int r,int num)
 {
     if(settag[num]!=-)
     {
         num1[lc]=settag[num]*(mid-l+);
         maxnum[lc][][]=maxnum[lc][][]=maxnum[lc][][]=(-settag[num])*(mid-l+);
         maxnum[lc][][]=maxnum[lc][][]=maxnum[lc][][]=settag[num]*(mid-l+);
         num1[rc]=settag[num]*(r-mid);
         maxnum[rc][][]=maxnum[rc][][]=maxnum[rc][][]=(-settag[num])*(r-mid);
         maxnum[rc][][]=maxnum[rc][][]=maxnum[rc][][]=settag[num]*(r-mid);
         revtag[lc]=revtag[rc]=;
         settag[lc]=settag[rc]=settag[num];
         settag[num]=-;
     }
     if(revtag[num])
     {
         num1[lc]=mid-l+-num1[lc];
         num1[rc]=r-mid-num1[rc];
         swap(maxnum[lc][][],maxnum[lc][][]);
         swap(maxnum[lc][][],maxnum[lc][][]);
         swap(maxnum[lc][][],maxnum[lc][][]);
         swap(maxnum[rc][][],maxnum[rc][][]);
         swap(maxnum[rc][][],maxnum[rc][][]);
         swap(maxnum[rc][][],maxnum[rc][][]);
         revtag[lc]^=;revtag[rc]^=;
         revtag[num]=;
     }
 }
 void update(int l,int r,int num)
 {
     num1[num]=num1[lc]+num1[rc];
     maxnum[num][][]=maxnum[lc][][];
     if(num1[lc]==)    maxnum[num][][]=(mid-l+)+maxnum[rc][][];
     maxnum[num][][]=maxnum[rc][][];
     if(num1[rc]==)    maxnum[num][][]=(r-mid)+maxnum[lc][][];
     maxnum[num][][]=maxnum[lc][][];
     if(num1[lc]==mid-l+)    maxnum[num][][]=(mid-l+)+maxnum[rc][][];
     maxnum[num][][]=maxnum[rc][][];
     if(num1[rc]==r-mid)    maxnum[num][][]=(r-mid)+maxnum[lc][][];
     maxnum[num][][]=max(max(maxnum[lc][][],maxnum[rc][][]),maxnum[lc][][]+maxnum[rc][][]);
     maxnum[num][][]=max(max(maxnum[lc][][],maxnum[rc][][]),maxnum[lc][][]+maxnum[rc][][]);
 }
 void build(int l,int r,int num)
 {
     if(l==r)
     {
         num1[num]=(a[l]==);
         maxnum[num][][]=maxnum[num][][]=maxnum[num][][]=(a[l]==);
         maxnum[num][][]=maxnum[num][][]=maxnum[num][][]=(a[l]==);
         settag[num]=-;
         return;
     }
     build(l,mid,lc);
     build(mid+,r,rc);
     settag[num]=-;update(l,r,num);
 }
 void setto(int l,int r,int num)
 {
     if(L<=l&&r<=R)
     {
         revtag[num]=;settag[num]=x;
         num1[num]=x*(r-l+);//mid-l+1
         maxnum[num][][]=maxnum[num][][]=maxnum[num][][]=(-x)*(r-l+);//mid-l+1
         maxnum[num][][]=maxnum[num][][]=maxnum[num][][]=x*(r-l+);//mid-l+1
         return;
     }
     pd(l,r,num);
     if(L<=mid)    setto(l,mid,lc);
     if(mid<R)    setto(mid+,r,rc);
     update(l,r,num);
 }
 void revx(int l,int r,int num)
 {
     if(L<=l&&r<=R)
     {
         revtag[num]^=;
         num1[num]=r-l+-num1[num];
         swap(maxnum[num][][],maxnum[num][][]);
         swap(maxnum[num][][],maxnum[num][][]);
         swap(maxnum[num][][],maxnum[num][][]);
         return;
     }
     pd(l,r,num);
     if(L<=mid)    revx(l,mid,lc);
     if(mid<R)    revx(mid+,r,rc);
     update(l,r,num);
 }
 int query1(int l,int r,int num)
 {
     if(L<=l&&r<=R)    return num1[num];
     pd(l,r,num);
     int ans=;
     if(L<=mid)    ans+=query1(l,mid,lc);
     if(mid<R)    ans+=query1(mid+,r,rc);
     return ans;
 }
 Th query2(int l,int r,int num)
 {
     //if(L<=l&&r<=R)    return Th(maxnum[num][1][0],maxnum[num][1][1],maxnum[num][1][2],num1[num]==r-l+1,1);
     if(L<=l&&r<=R)    return Th(maxnum[num][][],maxnum[num][][],maxnum[num][][],num1[num]==r-l+,r-l+);//num1[num]==1
     pd(l,r,num);
     Th ans,t1,t2;
     if(L<=mid)    t1=query2(l,mid,lc);
     if(mid<R)    t2=query2(mid+,r,rc);
     ans.a=t1.a;
     //if(t1.d||t1.len==0)    ans.a=t1.len+t2.a;
     if(t1.d)    ans.a=t1.len+t2.a;
     ans.b=t2.b;
     //if(t2.d||t2.len==0)    ans.b=t2.len+t1.b;
     if(t2.d)    ans.b=t2.len+t1.b;
     ans.c=max(max(t1.c,t2.c),t1.b+t2.a);
     ans.d=t1.d&&t2.d;
     ans.len=t1.len+t2.len;
     return ans;
 }
 int main()
 {
     int i,idx;Th ttt;
     scanf("%d%d",&n,&m);
     for(i=;i<=n;i++)    scanf("%d",&a[i]);
     build(,n,);
     for(i=;i<=m;i++)
     {
         scanf("%d%d%d",&idx,&L,&R);L++,R++;
         if(idx==)    x=,setto(,n,);
         else if(idx==)    x=,setto(,n,);
         else if(idx==)    revx(,n,);
         else if(idx==)    printf("%d\n",query1(,n,));
         else if(idx==)    {ttt=query2(,n,);printf("%d\n",max(max(ttt.a,ttt.b),ttt.c));}
     }
     return ;
 }