[bzoj5294]二进制

首先可以发现$2^k$模3意义下有循环节，也就是1,-1,1,-1……
考虑对于x个1，y个0，判断是否存在3的倍数
1.x为偶数时一定可以，选择等量的1和-1即可
2.x为奇数，要满足$xge 3$且$yge 2$，这是可以用3个0*(-1)和3个1*1来抵消掉（如果y=2时也可以，因为此时总共有奇数位），同时$x-3$显然为偶数
看上去难以维护，考虑反过来，求不是3的倍数的区间个数，那么即要求$x=1$或x为奇数且$y=0/1$
线段树维护区间，对于每一个区间维护0的数量，1的数量和一个三维数组a[3][4][3]表示在任意/以左端点为开头/右端点为结尾的区间中，x为0/1/非0偶数/非1奇数，y为0/1/(>1)的数量，转移详见代码（比较丑陋）

#include<bits/stdc++.h>
using namespace std;
#define N 100005
#define L (k<<1)
#define R (L+1)
#define mid (l+r>>1)
#define ll long long
struct ji{
    int s[2];
    ll a[3][4][3];
}o,f[N<<2];
int n,m,p,x,y,a[N];
int fx(int x){
    if (x<2)return x;
    return ((x&1)+2);
}
int fy(int y){
    return min(y,2);
}
ji up(ji x,ji y){
    if (x.s[0]<0)return y;
    if (y.s[0]<0)return x;
    for(int i=0;i<2;i++)o.s[i]=x.s[i]+y.s[i];
    memset(o.a[0],0,sizeof(o.a[0]));
    memcpy(o.a[1],x.a[1],sizeof(o.a[1]));
    memcpy(o.a[2],y.a[2],sizeof(o.a[2]));
    for(int i1=0;i1<4;i1++)
        for(int j1=0;j1<3;j1++){
            o.a[0][i1][j1]+=x.a[0][i1][j1]+y.a[0][i1][j1];
            o.a[1][fx(i1+x.s[1])][fy(j1+x.s[0])]+=y.a[1][i1][j1];
            o.a[2][fx(i1+y.s[1])][fy(j1+y.s[0])]+=x.a[2][i1][j1];
            for(int i2=0;i2<4;i2++)
                for(int j2=0;j2<3;j2++)
                    o.a[0][fx(i1+i2)][fy(j1+j2)]+=x.a[2][i1][j1]*y.a[1][i2][j2];
        }
    return o;
}
void update(int k,int l,int r,int x,int y){
    if (l==r){
        f[k].s[y]=1;
        f[k].s[y^1]=0;
        memset(f[k].a,0,sizeof(f[k].a));
        for(int i=0;i<3;i++)f[k].a[i][y][y^1]=1;
        return;
    }
    if (x<=mid)update(L,l,mid,x,y);
    else update(R,mid+1,r,x,y);
    f[k]=up(f[L],f[R]);
}
ji query(int k,int l,int r,int x,int y){
    if ((l>y)||(x>r))return f[0];
    if ((x<=l)&&(r<=y))return f[k];
    return up(query(L,l,mid,x,y),query(R,mid+1,r,x,y));
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)update(1,1,n,i,a[i]);
    f[0].s[0]=-1;
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&p,&x);
        if (p==1)update(1,1,n,x,a[x]^=1);
        else{
            scanf("%d",&y);
            ll ans=(y-x+2LL)*(y-x+1)/2;
            o=query(1,1,n,x,y);
            for(int j=0;j<3;j++)ans-=o.a[0][1][j];
            for(int j=0;j<2;j++)ans-=o.a[0][3][j];
            printf("%lld
",ans);
        }
    }
}