帮一个MnZn问题

#include<bits/stdc++.h>
#define lo (nw<<1)
#define ro (nw<<1|1)
#define md ((l+r)>>1)
#define int long long
using namespace std;

const int N=100010,INF=0x7fffffffffffffff;

int n,m;
int a[N];

struct node
{
    int tr,sz;
    int sm1,sm2;
    int fr1,fr2;
    int re1,re2;
    int c1,c2;
}t[N << 2];

struct segment_tree
{
    void upd(int nw)
    {
        t[nw].sz=t[lo].sz+t[ro].sz;
        t[nw].tr=t[lo].tr+t[ro].tr;
        t[nw].sm1=max(t[lo].re1+t[ro].fr1,max(t[lo].sm1,t[ro].sm1));
        if(t[lo].fr1==t[lo].sz)
            t[nw].fr1=t[lo].fr1+t[ro].fr1;
        else
            t[nw].fr1=t[lo].fr1;
        if(t[ro].re1==t[ro].sz)
            t[nw].re1=t[ro].re1+t[lo].re1;
        else
            t[nw].re1=t[ro].re1;
        t[nw].sm2=max(t[lo].re2+t[ro].fr2,max(t[lo].sm2,t[ro].sm2));
        if(t[lo].fr2==t[lo].sz)
            t[nw].fr2=t[lo].fr2+t[ro].fr2;
        else
            t[nw].fr2=t[lo].fr2;
        if(t[ro].re2==t[ro].sz)
            t[nw].re2=t[ro].re2+t[lo].re2;
        else
            t[nw].re2=t[ro].re2;
        return ;
    }
    void mkc1(int nw,int z)
    {
        t[nw].c1=z;
        t[nw].tr=z*t[nw].sz;
        t[nw].sm1=z*t[nw].sz;
        t[nw].fr1=z*t[nw].sz;
        t[nw].re1=z*t[nw].sz;
        t[nw].sm2=(z^1)*t[nw].sz;
        t[nw].fr2=(z^1)*t[nw].sz;
        t[nw].re2=(z^1)*t[nw].sz;
        t[nw].c2=0;
        return ;
    }
    void mkc2(int nw)
    {
        t[nw].c2^=1;
        t[nw].tr=t[nw].sz-t[nw].tr;
        swap(t[nw].sm1,t[nw].sm2);
        swap(t[nw].fr1,t[nw].fr2);
        swap(t[nw].re1,t[nw].re2);
        return ;
    }
    void psd(int nw)
    {
        if(t[nw].c1!=-1)
        {
            mkc1(lo,t[nw].c1);
            mkc1(ro,t[nw].c1);
            t[nw].c1=-1;
        }
        if(t[nw].c2!=0)
        {
            mkc2(lo);
            mkc2(ro);
            t[nw].c2=0;
        }
        return ;
    }
    void bld(int nw,int l,int r)
    {
        t[nw].c1=-1;
        t[nw].c2=0;
        if(l==r)
        {
            t[nw].tr=a[l];
            t[nw].sz=1;
            t[nw].sm1=a[l];
            t[nw].fr1=a[l];
            t[nw].re1=a[l];
            t[nw].sm2=a[l]^1;
            t[nw].fr2=a[l]^1;
            t[nw].re2=a[l]^1;
            return ;
        }
        bld(lo,l,md);
        bld(ro,md+1,r);
        upd(nw);
        return ;
    }
    void atr(int nw,int l,int r,int x,int y,int z)
    {
        psd(nw);
        if(x<=l&&r<=y)
        {
            mkc1(nw,z);
            return ;
        }

        if(x<=md)
            atr(lo,l,md,x,y,z);
        if(y>md)
            atr(ro,md+1,r,x,y,z);
        upd(nw);
        return ;
    }
    void atc(int nw,int l,int r,int x,int y)
    {
        psd(nw);
        if(x<=l&&r<=y)
        {
            mkc2(nw);
            return ;
        }

        if(x<=md)
            atc(lo,l,md,x,y);
        if(y>md)
            atc(ro,md+1,r,x,y);
        upd(nw);
        return ;
    }
    int qrt(int nw,int l,int r,int x,int y)
    {
        if(x<=l&&r<=y)
            return t[nw].tr;
        psd(nw);
        int sm=0;
        if(x<=md)
            sm=sm+qrt(lo,l,md,x,y);
        if(y>md)
            sm=sm+qrt(ro,md+1,r,x,y);
        return sm;
    }
    int qry(int nw,int l,int r,int x,int y)
    {
        if(x<=l&&r<=y)
            return t[nw].sm1;
        psd(nw);
        int sm=-INF;
        if(x<=md)
            sm=max(sm, qry(lo,l,md,x,y));
        if(y>md)
            sm=max(sm, qry(ro,md+1,r,x,y));
        return sm;
    }
}ST;

signed main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    ST.bld(1,1,n);
    for(int i=1;i<=m;i++)
    {
        int op,x,y;
        cin>>op>>x>>y;
        ++x; ++y;
        if(op==0||op==1)
        {
            ST.atr(1,1,n,x,y,op);
        }
        if(op==2)
        {
            ST.atc(1,1,n,x,y);
        }
        if(op==3)
        {
            cout<<ST.qrt(1,1,n,x,y)<<'\n';
        }
        if(op==4)
        {
            cout<<ST.qry(1,1,n,x,y)<<'\n';
        }
    }
    return 0;
}