bzoj1858

1858: [Scoi2010]序列操作

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 444  Solved: 240
[Submit][Status][Discuss]

Description

lxhgww最近收到了一个01序列，序列里面包含了n个数，这些数要么是0，要么是1，现在对于这个序列有五种变换操作和询问操作：0 a b 把[a, b]区间内的所有数全变成01 a b 把[a, b]区间内的所有数全变成12 a b 把[a,b]区间内的所有数全部取反，也就是说把所有的0变成1，把所有的1变成03 a b 询问[a, b]区间内总共有多少个14 a b 询问[a, b]区间内最多有多少个连续的1对于每一种询问操作，lxhgww都需要给出回答，聪明的程序员们，你们能帮助他吗？

Input

输入数据第一行包括2个数，n和m，分别表示序列的长度和操作数目 第二行包括n个数，表示序列的初始状态 接下来m行，每行3个数，op, a, b，（0<=op<=4，0<=a<=b

Output

对于每一个询问操作，输出一行，包括1个数，表示其对应的答案

Sample Input

10 10
0 0 0 1 1 0 1 0 1 1
1 0 2
3 0 5
2 2 2
4 0 4
0 3 6
2 3 7
4 2 8
1 0 5
0 5 6
3 3 9

Sample Output

5
2
6
5

HINT

对于30%的数据，1<=n, m<=1000
对于100%的数据，1<=n, m<=100000

线段树的应用

好久不写普通线段树了，写了竟然wa了……

program bzoj1858;type    group=record                max,maxl,maxr:longint;    end;var   a,b,o,tot,n,m,i,j,k:longint;   stop,sum,left,right,maxl1,maxr1,max1,maxl0,maxr0,max0:array [0..200001] of longint;procedure swap (var a,b:longint);inline;begin     if a=b then exit;     a:=a xor b;     b:=a xor b;     a:=a xor b;end;function max (a,b:longint):longint;inline;begin     if a>b then exit(a) else exit(b);end;procedure build (s,e,now:longint);var   mid:longint;begin     stop[now]:=-1;     maxl0[now]:=e-s+1;     maxr0[now]:=e-s+1;     max0[now]:=e-s+1;     if s=e then exit;     mid:=(s+e) div 2;     inc(tot);     left[now]:=tot;     build(s,mid,tot);     inc(tot);     right[now]:=tot;     build(mid+1,e,tot);end;procedure quick (o,l,r,now:longint);begin     case o of          0:          begin               sum[now]:=0;               maxl0[now]:=r-l+1;               maxr0[now]:=r-l+1;               max0[now]:=r-l+1;               maxl1[now]:=0;               maxr1[now]:=0;               max1[now]:=0;               stop[now]:=0;          end;          1:          begin               sum[now]:=r-l+1;               maxl0[now]:=0;               maxr0[now]:=0;               max0[now]:=0;               maxl1[now]:=r-l+1;               maxr1[now]:=r-l+1;               max1[now]:=r-l+1;               stop[now]:=1;          end;          2:          begin               sum[now]:=r-l+1-sum[now];               swap(maxl0[now],maxl1[now]);               swap(maxr0[now],maxr1[now]);               swap(max0[now],max1[now]);               case stop[now] of                    -1:stop[now]:=2;                    0:stop[now]:=1;                    1:stop[now]:=0;                    2:stop[now]:=-1;               end;          end;     end;end;procedure update (l,r,now:longint);var   mid:longint;begin     if stop[now]=-1 then exit;     mid:=(l+r) div 2;     quick(stop[now],l,mid,left[now]);     quick(stop[now],mid+1,r,right[now]);     stop[now]:=-1;end;procedure fill (o,s,e,l,r,now:longint);var   mid:longint;begin     if (s=l)and(e=r) then        begin             quick(o,l,r,now);             exit;        end;     update(l,r,now);     mid:=(l+r) div 2;     if e<=mid then fill(o,s,e,l,mid,left[now])               else     if s>mid then fill(o,s,e,mid+1,r,right[now])              else              begin                   fill(o,s,mid,l,mid,left[now]);                   fill(o,mid+1,e,mid+1,r,right[now]);              end;     sum[now]:=sum[left[now]]+sum[right[now]];     if maxl0[left[now]]=mid-l+1 then maxl0[now]:=maxl0[left[now]]+maxl0[right[now]]                                 else maxl0[now]:=maxl0[left[now]];     if maxr0[right[now]]=r-mid then maxr0[now]:=maxr0[right[now]]+maxr0[left[now]]                                else maxr0[now]:=maxr0[right[now]];     max0[now]:=max0[left[now]];     if max0[now]<max0[right[now]] then max0[now]:=max0[right[now]];     if max0[now]<maxr0[left[now]]+maxl0[right[now]] then max0[now]:=maxr0[left[now]]+maxl0[right[now]];     if maxl1[left[now]]=mid-l+1 then maxl1[now]:=maxl1[left[now]]+maxl1[right[now]]                                 else maxl1[now]:=maxl1[left[now]];     if maxr1[right[now]]=r-mid then maxr1[now]:=maxr1[right[now]]+maxr1[left[now]]                                else maxr1[now]:=maxr1[right[now]];     max1[now]:=max1[left[now]];     if max1[now]<max1[right[now]] then max1[now]:=max1[right[now]];     if max1[now]<maxr1[left[now]]+maxl1[right[now]] then max1[now]:=maxr1[left[now]]+maxl1[right[now]];end;function qsum (s,e,l,r,now:longint):longint;var   mid:longint;begin     if (s=l)and(e=r) then exit(sum[now]);     update(l,r,now);     mid:=(l+r) div 2;     if e<=mid then exit(qsum(s,e,l,mid,left[now]))               else     if s>mid then exit(qsum(s,e,mid+1,r,right[now]))              else exit(qsum(s,mid,l,mid,left[now])+qsum(mid+1,e,mid+1,r,right[now]));end;function find (s,e,l,r,now:longint):group;var   mid:longint;   ans,q,p:group;begin     if (s=l)and(e=r) then        begin             ans.max:=max1[now];             ans.maxl:=maxl1[now];             ans.maxr:=maxr1[now];             exit(ans);        end;     update(l,r,now);     mid:=(l+r) div 2;     if e<=mid then exit(find(s,e,l,mid,left[now]))               else     if s>mid then exit(find(s,e,mid+1,r,right[now]))              else     begin          q:=find(s,mid,l,mid,left[now]);          p:=find(mid+1,e,mid+1,r,right[now]);          ans.max:=max(max(q.max,p.max),q.maxr+p.maxl);          if q.maxl=mid-s+1 then ans.maxl:=q.maxl+p.maxl                            else ans.maxl:=q.maxl;          if p.maxr=e-mid then ans.maxr:=p.maxr+q.maxr                          else ans.maxr:=p.maxr;          exit(ans);     end;end;begin                         read(n,m);     build(1,n,0);     for i:=1 to n do         begin              read(k);              fill(k,i,i,1,n,0);         end;     for i:=1 to m do         begin              read(o,a,b);              inc(a);              inc(b);              case o of                   0:fill(0,a,b,1,n,0);                   1:fill(1,a,b,1,n,0);                   2:fill(2,a,b,1,n,0);                   3:writeln(qsum(a,b,1,n,0));                   4:writeln(find(a,b,1,n,0).max);              end;         end;end.