JZOJ 4933. 【NOIP2017提高组模拟12.24】C

题目

这里写图片描述

题解

这是一道想都能想到方法的码农题
然而,我们要注意几个要点。
在维护值的时候,我们只用打2个LAZY标志。一个更新最新值,另一个判断是修改还是覆盖。
还有,以后在打线段树的时候,不要开二维数组,尽量用一维。这样,各数组的名称可以直观看出维护哪个值。
重点讲一下第7个操作。
我们维护这样几个值:最左边的数ln,最右边的数rn,左边的lwz,右边的rwz,这段区间的答案lg。
这里写图片描述
当合并的时候,比一下两端区间的lg,然后新的区间的ln,ls等于做区间的ln,ls,新的区间的rn,rs等于做区间的rn,rs。如果左区间的rn=右区间的ln,那么尝试用ls_左区间+rs _右区间 更新区间的lg。
在查询的时候,递归要变成函数,比较一下左边lg大还是右边lg大,如果此区间能够合并,那么看看合并后能不能更新lg(注意越界)。
我说的越界如下图:
这里写图片描述
还有其他很多地方要注意的,十分恶心,希望大家打代码时注意一下。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define N 100010
#define LL long long
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
LL a[N];
LL mx[N*15],mn[N*15],sum[N*15],sum1[N*15],la[N*15],rp[N*15];
LL ln[N*15],rn[N*15],ls[N*15],rs[N*15],lg[N*15];
LL i,j,K,L,R,n,q,op,ans,cnt;
inline LL read()
{
    LL p,data=0;
    char ch=0;
    while ((ch!='-') && ch<'0' || ch>'9') ch=getchar();
    if (ch=='-')
    {
        p=-1;
        ch=getchar();
    } else p=1;
    while (ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
    return data*p;
}
void down(LL s,LL l,LL r)
{
    if (rp[s]==1)
    {
        sum[s]+=la[s]*(r-l+1);
        mx[s]+=la[s];
        mn[s]+=la[s];
        ln[s]+=la[s];
        rn[s]+=la[s];
        if (!rp[s*2]) rp[s*2]=1;
        if (!rp[s*2+1]) rp[s*2+1]=1;
        la[s*2]+=la[s];
        la[s*2+1]+=la[s];
        rp[s]=la[s]=0;
    } else
    if (rp[s]==2)
    {
        sum[s]=la[s]*(r-l+1);
        mx[s]=mn[s]=ln[s]=rn[s]=la[s];
        lg[s]=ls[s]=rs[s]=r-l+1;
        rp[s*2]=rp[s*2+1]=2;
        la[s*2]=la[s*2+1]=la[s];
        rp[s]=la[s]=0;
    }
}
void up(LL s,LL l,LL r)
{
    sum[s]=sum[s*2]+sum[s*2+1];
    mx[s]=max(mx[s*2],mx[s*2+1]);
    mn[s]=min(mn[s*2],mn[s*2+1]);
    ls[s]=ls[s*2];
    rs[s]=rs[s*2+1];
    ln[s]=ln[s*2];
    rn[s]=rn[s*2+1];
    lg[s]=max(lg[s*2],lg[s*2+1]);
    if (rn[s*2]==ln[s*2+1])
    {
        lg[s]=max(lg[s],rs[s*2]+ls[s*2+1]);
        if (ls[s*2]==(l+r)/2-l+1) ls[s]+=ls[s*2+1];
        if (rs[s*2+1]==r-(l+r)/2) rs[s]+=rs[s*2];
    }
}
void build(LL s,LL l,LL r)
{
    if (l==r)
    {
        mx[s]=mn[s]=sum[s]=ln[s]=rn[s]=a[l];
        ls[s]=rs[s]=lg[s]=1;
        return;
    }
    LL mid=(l+r)/2;
    build(s*2,l,mid);
    build(s*2+1,mid+1,r);
    up(s,l,r);
}
void inspls(LL s,LL l,LL r,LL x,LL y,LL z)
{
    LL wz=(l+r)/2;
    down(s*2,l,wz);down(s*2+1,wz+1,r);
    if (l==x && r==y)
    {
        la[s]=z;rp[s]=1;
        down(s,l,r);
        return;
    }
    if (y<=wz) inspls(s*2,l,wz,x,y,z);
        else if (x>wz) inspls(s*2+1,wz+1,r,x,y,z);
            else
            {
                inspls(s*2,l,wz,x,wz,z);
                inspls(s*2+1,wz+1,r,wz+1,y,z);
            }
    up(s,l,r);
}
void replace(LL s,LL l,LL r,LL x,LL y,LL z)
{
    LL wz=(l+r)/2;
    down(s*2,l,wz);down(s*2+1,wz+1,r);
    if (l==x && r==y)
    {
        la[s]=z;rp[s]=2;
        down(s,l,r);
        return;
    }
    if (y<=wz) replace(s*2,l,wz,x,y,z);
        else if (x>wz) replace(s*2+1,wz+1,r,x,y,z);
            else
            {
                replace(s*2,l,wz,x,wz,z);
                replace(s*2+1,wz+1,r,wz+1,y,z);
            }
    up(s,l,r);
}
void find(LL s,LL l,LL r,LL x,LL y,LL op)
{
    LL wz=(l+r)/2;
    down(s*2,l,wz);down(s*2+1,wz+1,r);
    if (l==x && r==y)
    {
        if (op==4) ans+=sum[s];
        if (op==5) ans=min(ans,mn[s]);
        if (op==6) ans=max(ans,mx[s]);
        return;
    }
    if (y<=wz) find(s*2,l,wz,x,y,op);
        else if (x>wz) find(s*2+1,wz+1,r,x,y,op);
            else
            {
                find(s*2,l,wz,x,wz,op);
                find(s*2+1,wz+1,r,wz+1,y,op);
            }
    up(s,l,r);
}
LL findlg(LL s,LL l,LL r,LL x,LL y)
{
    LL wz=(l+r)/2;
    down(s*2,l,wz);down(s*2+1,wz+1,r);
    if (l==x && r==y) return lg[s];
    if (y<=wz) return findlg(s*2,l,wz,x,y);
        else if (x>wz) return findlg(s*2+1,wz+1,r,x,y);
            else
            {
                LL ans;
                ans=max(findlg(s*2,l,wz,x,wz),findlg(s*2+1,wz+1,r,wz+1,y));
                if (rn[s*2]==ln[s*2+1]) 
                    ans=max(ans,min(rs[s*2],wz-x+1)+min(ls[s*2+1],y-wz));
                return ans;
            }
    up(s,l,r);
}
int main()
{
    n=read();q=read();
    fo(i,1,n) a[i]=read();
    build(1,1,n);
    fo(i,1,q)
    {
        op=read();
        if (op<4)
        {
            L=read();R=read();K=read();
            if (op==1) inspls(1,1,n,L,R,K);
            if (op==2) inspls(1,1,n,L,R,-K);
            if (op==3) replace(1,1,n,L,R,K);
        } else
        {
            L=read();R=read();
            if (op==4) ans=0;
            if (op==5) ans=2147483647;
            if (op==6) ans=-2147483647;
            if (op==7) ans=findlg(1,1,n,L,R);
                  else find(1,1,n,L,R,op);
            printf("%lld\n",ans);
        }
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值