HDU 3397 Sequence operation

题目意思很简单，四种操作：全部变成1，全部变成0，或者区间异或，查询区间1的个数，和区间连续的1的最大长度。理解起来没有问题，可是写起代码来发现非常繁琐，共用了9个数组，差不多200行了，第一次写这么长得代码，WA一次，后来发现某个地方rt<<1,写成了rt<<1|1,还有PushUp向上更新时忘记了cover[rt]也要根据sum[rt]的值进行更新。
  1 #include<cstdio>
  2 #include<cstring>
  3 #include<algorithm>
  4 #include<cstdlib>
  5 #define N 111111
  6 #define lson l,m,rt<<1
  7 #define rson m+1,r,rt<<1|1
  8 using namespace std;
  9 
 10 int sum[N<<2],msum[N<<2],lsum[N<<2],rsum[N<<2],XOR[N<<2],cover[N<<2],msum0[N<<2],lsum0[N<<2],rsum0[N<<2];
 11 
 12 
 13 //sum表示区间1的个数，msum表示区间连续1的最大长度，lsum表示从左端开始的1的最大长度，
 14 //XOR为异或标记，cover为覆盖标记，1,0表示区间分别被1,0完全覆盖，-1不完全覆盖
 15 //msum0,lsum0,rsum0表示对应的0的个数的统计。
 16 
 17 
 18 void PushUp(int rt,int m)
 19 {
 20     sum[rt]=sum[rt<<1]+sum[rt<<1|1];
 21     cover[rt]=((sum[rt]==m)?1:((sum[rt]==0)?0:-1));
 22     lsum[rt]=lsum[rt<<1];//我擦
 23     rsum[rt]=rsum[rt<<1|1];
 24     lsum0[rt]=lsum0[rt<<1];
 25     rsum0[rt]=rsum0[rt<<1|1];
 26     if(lsum[rt]==(m-(m>>1)))
 27         lsum[rt]+=lsum[rt<<1|1];
 28     else if(lsum0[rt]==(m-(m>>1)))
 29         lsum0[rt]+=lsum0[rt<<1|1];
 30     if(rsum[rt]==(m>>1))
 31         rsum[rt]+=rsum[rt<<1];
 32     else if(rsum0[rt]==(m>>1))
 33         rsum0[rt]+=rsum0[rt<<1];
 34     msum0[rt]=max(max(msum0[rt<<1],msum0[rt<<1|1]),rsum0[rt<<1]+lsum0[rt<<1|1]);
 35     msum[rt]=max(max(msum[rt<<1],msum[rt<<1|1]),rsum[rt<<1]+lsum[rt<<1|1]);
 36 }
 37 
 38 void FXOR(int rt,int m)
 39 {
 40     if(cover[rt]!=-1)
 41     {
 42         cover[rt]^=1;
 43         sum[rt]=msum[rt]=lsum[rt]=rsum[rt]=m*cover[rt];
 44         msum0[rt]=lsum0[rt]=rsum0[rt]=m*(1-cover[rt]);
 45     }
 46     else
 47     {
 48         XOR[rt]^=1;
 49         sum[rt]=(m-sum[rt]);
 50         int a=lsum[rt],b=rsum[rt];
 51         lsum[rt]=lsum0[rt];
 52         rsum[rt]=rsum0[rt];
 53         lsum0[rt]=a;
 54         rsum0[rt]=b;
 55         a=msum[rt];
 56         msum[rt]=msum0[rt];
 57         msum0[rt]=a;
 58     }
 59 }
 60 
 61 void PushDown(int rt,int m)
 62 {
 63     if(cover[rt]!=-1)
 64     {
 65         lsum[rt<<1]=rsum[rt<<1]=msum[rt<<1]=sum[rt<<1]=(m-(m>>1))*cover[rt];
 66         lsum[rt<<1|1]=rsum[rt<<1|1]=msum[rt<<1|1]=sum[rt<<1|1]=(m>>1)*cover[rt];
 67         lsum0[rt<<1]=rsum0[rt<<1]=msum0[rt<<1]=(m-(m>>1))*(1-cover[rt]);
 68         lsum0[rt<<1|1]=rsum0[rt<<1|1]=msum0[rt<<1|1]=(m>>1)*(1-cover[rt]);
 69         cover[rt<<1]=cover[rt<<1|1]=cover[rt];
 70         XOR[rt<<1]=XOR[rt<<1|1]=0;
 71         cover[rt]=-1;
 72     }
 73     if(XOR[rt])
 74     {
 75         FXOR(rt<<1,(m-(m>>1)));
 76         FXOR(rt<<1|1,(m>>1));
 77         XOR[rt]=0;
 78     }
 79 
 80 }
 81 
 82 void build(int l,int r,int rt)
 83 {
 84     XOR[rt]=0;
 85     if(l==r)
 86     {
 87         scanf("%d",&sum[rt]);
 88         lsum[rt]=rsum[rt]=msum[rt]=cover[rt]=sum[rt];
 89         lsum0[rt]=rsum0[rt]=msum0[rt]=1-sum[rt];
 90         return;
 91     }
 92     int m=(l+r)>>1;
 93     build(lson);
 94     build(rson);
 95     PushUp(rt,r-l+1);
 96 }
 97 
 98 void update(int L,int R,int c,int l,int r,int rt)
 99 {
100     if(L<=l&&R>=r)
101     {
102         if(c==0)
103             cover[rt]=0,lsum[rt]=rsum[rt]=sum[rt]=msum[rt]=0,lsum0[rt]=rsum0[rt]=msum0[rt]=r-l+1,XOR[rt]=0;
104         else if(c==1)
105             cover[rt]=1,lsum[rt]=rsum[rt]=sum[rt]=msum[rt]=r-l+1,lsum0[rt]=rsum0[rt]=msum0[rt]=0,XOR[rt]=0;
106         else FXOR(rt,r-l+1);
107         return;
108     }
109     int m=(l+r)>>1;
110     PushDown(rt,r-l+1);
111     if(L<=m)
112         update(L,R,c,lson);
113     if(R>m)
114         update(L,R,c,rson);
115     PushUp(rt,r-l+1);
116 }
117 
118 int query1(int L,int R,int l,int r,int rt)
119 {
120     if(L<=l&&R>=r)
121         return sum[rt];
122     int m=(l+r)>>1;
123     PushDown(rt,r-l+1);
124     int res=0;
125     if(L<=m)
126         res+=query1(L,R,lson);
127     if(R>m)
128         res+=query1(L,R,rson);
129     return res;
130 }
131 
132 int query2(int L,int R,int l,int r,int rt)
133 {
134     if(L<=l&&R>=r)
135         return msum[rt];
136     int m=(l+r)>>1;
137     PushDown(rt,r-l+1);
138     int ans=0;
139     if(R<=m)
140         ans=max(ans,query2(L,R,lson));
141     else if(L>m)
142         ans=max(ans,query2(L,R,rson));
143     else
144     {
145         ans=max(ans,query2(L,R,lson));
146         ans=max(ans,query2(L,R,rson));
147         int ll,rr;
148         ll=min(m-L+1,rsum[rt<<1]);
149         rr=min(R-m,lsum[rt<<1|1]);
150         ans=max(ans,ll+rr);
151     }
152     return ans;
153 }
154 
155 int main(void)
156 {
157     int tc,n,m;
158     int op,a,b,ans;
159     scanf("%d",&tc);
160     while(tc--)
161     {
162         scanf("%d%d",&n,&m);
163         build(1,n,1);
164         while(m--)
165         {
166             scanf("%d%d%d",&op,&a,&b);
167             switch(op)
168             {
169             case 0:
170             case 1:
171             case 2:
172                 update(a+1,b+1,op,1,n,1);
173                 break;
174             case 3:
175                 ans=query1(a+1,b+1,1,n,1);
176                 printf("%d
",ans);
177                 break;
178             case 4:
179                 ans=query2(a+1,b+1,1,n,1);
180                 printf("%d
",ans);
181                 break;
182             }
183         }
184     }
185     return 0;
186 }