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:以下程序中打了注释的语句都是错过的。。。

  1 #include<cstdio>
  2 #include<cstring>
  3 #include<algorithm>
  4 #define mid ((l+r)>>1)
  5 #define lc (num<<1)
  6 #define rc (num<<1|1)
  7 using namespace std;
  8 struct Th
  9 {
 10     int a,b,c;
 11     bool d;//是否全部是1
 12     int len;
 13     Th(int aa=0,int bb=0,int cc=0,bool dd=1,int ll=0):a(aa),b(bb),c(cc),d(dd),len(ll){}
 14 };
 15 int num1[270000],maxnum[270000][2][3],settag[270000];
 16 bool revtag[270000];
 17 int n,m;
 18 //定义tag同时存在时settag先进行操作
 19 int a[100010];
 20 int L,R,x;
 21 //maxnum[][0/1][0/1/2]表示最多的连续0/1,接左侧/接右侧/中间
 22 void pd(int l,int r,int num)
 23 {
 24     if(settag[num]!=-1)
 25     {
 26         num1[lc]=settag[num]*(mid-l+1);
 27         maxnum[lc][0][0]=maxnum[lc][0][1]=maxnum[lc][0][2]=(1-settag[num])*(mid-l+1);
 28         maxnum[lc][1][0]=maxnum[lc][1][1]=maxnum[lc][1][2]=settag[num]*(mid-l+1);
 29         num1[rc]=settag[num]*(r-mid);
 30         maxnum[rc][0][0]=maxnum[rc][0][1]=maxnum[rc][0][2]=(1-settag[num])*(r-mid);
 31         maxnum[rc][1][0]=maxnum[rc][1][1]=maxnum[rc][1][2]=settag[num]*(r-mid);
 32         revtag[lc]=revtag[rc]=0;
 33         settag[lc]=settag[rc]=settag[num];
 34         settag[num]=-1;
 35     }
 36     if(revtag[num])
 37     {
 38         num1[lc]=mid-l+1-num1[lc];
 39         num1[rc]=r-mid-num1[rc];
 40         swap(maxnum[lc][0][0],maxnum[lc][1][0]);
 41         swap(maxnum[lc][0][1],maxnum[lc][1][1]);
 42         swap(maxnum[lc][0][2],maxnum[lc][1][2]);
 43         swap(maxnum[rc][0][0],maxnum[rc][1][0]);
 44         swap(maxnum[rc][0][1],maxnum[rc][1][1]);
 45         swap(maxnum[rc][0][2],maxnum[rc][1][2]);
 46         revtag[lc]^=1;revtag[rc]^=1;
 47         revtag[num]=0;
 48     }
 49 }
 50 void update(int l,int r,int num)
 51 {
 52     num1[num]=num1[lc]+num1[rc];
 53     maxnum[num][0][0]=maxnum[lc][0][0];
 54     if(num1[lc]==0)    maxnum[num][0][0]=(mid-l+1)+maxnum[rc][0][0];
 55     maxnum[num][0][1]=maxnum[rc][0][1];
 56     if(num1[rc]==0)    maxnum[num][0][1]=(r-mid)+maxnum[lc][0][1];
 57     maxnum[num][1][0]=maxnum[lc][1][0];
 58     if(num1[lc]==mid-l+1)    maxnum[num][1][0]=(mid-l+1)+maxnum[rc][1][0];
 59     maxnum[num][1][1]=maxnum[rc][1][1];
 60     if(num1[rc]==r-mid)    maxnum[num][1][1]=(r-mid)+maxnum[lc][1][1];
 61     maxnum[num][0][2]=max(max(maxnum[lc][0][2],maxnum[rc][0][2]),maxnum[lc][0][1]+maxnum[rc][0][0]);
 62     maxnum[num][1][2]=max(max(maxnum[lc][1][2],maxnum[rc][1][2]),maxnum[lc][1][1]+maxnum[rc][1][0]);
 63 }
 64 void build(int l,int r,int num)
 65 {
 66     if(l==r)
 67     {
 68         num1[num]=(a[l]==1);
 69         maxnum[num][0][0]=maxnum[num][0][1]=maxnum[num][0][2]=(a[l]==0);
 70         maxnum[num][1][0]=maxnum[num][1][1]=maxnum[num][1][2]=(a[l]==1);
 71         settag[num]=-1;
 72         return;
 73     }
 74     build(l,mid,lc);
 75     build(mid+1,r,rc);
 76     settag[num]=-1;update(l,r,num);
 77 }
 78 void setto(int l,int r,int num)
 79 {
 80     if(L<=l&&r<=R)
 81     {
 82         revtag[num]=0;settag[num]=x;
 83         num1[num]=x*(r-l+1);//mid-l+1
 84         maxnum[num][0][0]=maxnum[num][0][1]=maxnum[num][0][2]=(1-x)*(r-l+1);//mid-l+1
 85         maxnum[num][1][0]=maxnum[num][1][1]=maxnum[num][1][2]=x*(r-l+1);//mid-l+1
 86         return;
 87     }
 88     pd(l,r,num);
 89     if(L<=mid)    setto(l,mid,lc);
 90     if(mid<R)    setto(mid+1,r,rc);
 91     update(l,r,num);
 92 }
 93 void revx(int l,int r,int num)
 94 {
 95     if(L<=l&&r<=R)
 96     {
 97         revtag[num]^=1;
 98         num1[num]=r-l+1-num1[num];
 99         swap(maxnum[num][0][0],maxnum[num][1][0]);
100         swap(maxnum[num][0][1],maxnum[num][1][1]);
101         swap(maxnum[num][0][2],maxnum[num][1][2]);
102         return;
103     }
104     pd(l,r,num);
105     if(L<=mid)    revx(l,mid,lc);
106     if(mid<R)    revx(mid+1,r,rc);
107     update(l,r,num);
108 }
109 int query1(int l,int r,int num)
110 {
111     if(L<=l&&r<=R)    return num1[num];
112     pd(l,r,num);
113     int ans=0;
114     if(L<=mid)    ans+=query1(l,mid,lc);
115     if(mid<R)    ans+=query1(mid+1,r,rc);
116     return ans;
117 }
118 Th query2(int l,int r,int num)
119 {
120     //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);
121     if(L<=l&&r<=R)    return Th(maxnum[num][1][0],maxnum[num][1][1],maxnum[num][1][2],num1[num]==r-l+1,r-l+1);//num1[num]==1
122     pd(l,r,num);
123     Th ans,t1,t2;
124     if(L<=mid)    t1=query2(l,mid,lc);
125     if(mid<R)    t2=query2(mid+1,r,rc);
126     ans.a=t1.a;
127     //if(t1.d||t1.len==0)    ans.a=t1.len+t2.a;
128     if(t1.d)    ans.a=t1.len+t2.a;
129     ans.b=t2.b;
130     //if(t2.d||t2.len==0)    ans.b=t2.len+t1.b;
131     if(t2.d)    ans.b=t2.len+t1.b;
132     ans.c=max(max(t1.c,t2.c),t1.b+t2.a);
133     ans.d=t1.d&&t2.d;
134     ans.len=t1.len+t2.len;
135     return ans;
136 }
137 int main()
138 {
139     int i,idx;Th ttt;
140     scanf("%d%d",&n,&m);
141     for(i=1;i<=n;i++)    scanf("%d",&a[i]);
142     build(1,n,1);
143     for(i=1;i<=m;i++)
144     {
145         scanf("%d%d%d",&idx,&L,&R);L++,R++;
146         if(idx==0)    x=0,setto(1,n,1);
147         else if(idx==1)    x=1,setto(1,n,1);
148         else if(idx==2)    revx(1,n,1);
149         else if(idx==3)    printf("%d
",query1(1,n,1));
150         else if(idx==4)    {ttt=query2(1,n,1);printf("%d
",max(max(ttt.a,ttt.b),ttt.c));}
151     }
152     return 0;
153 }
原文地址:https://www.cnblogs.com/hehe54321/p/8525697.html