[JZOJ5618]【NOI2018模拟3.31】华胥梦天

Description

N个数的序列A
要求支持区间对一个数取min，单点做减法且对0取max，区间求和
$N\leq500000,A_i\leq10^9$

Solution

经典套路题啊（参见吉老师的SegmentTree-Beats）

之前CTY讲课的时候还跟他谈笑风生来着说这不是大暴力么
结果过没多久就忘了qwq

直接用线段树维护区间最大值、最大值有多少个和区间次大值
假设我们要取min的数为c
若c大于等于最大值，那么显然没有影响
若c小于最大值但大于次大值，那么我们只需要将最大值改为c，直接用个数算贡献，然后在这个区间打上取min标记，等待以后下传

若c小于等于次大值，那么线段树继续递归下去
果然暴力。。。

复杂度证明参考吉老师的集训队论文
最后的复杂度是$O(m\log n)$的

Code

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 500005
#define M 2000005
#define INF 1802201963
#define LL long long
using namespace std;
int n,m,n1,mx[2*M],md[2*M],ct[2*M],lz[2*M];
LL sm[2*M],a[N];
LL read(LL &n)
{
    char ch=' ';LL q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
struct node
{
    int l,r;
}tr[2*M];
void update(int k,int c)
{
    if(c>=mx[k]) return;
    if(md[k]<c) 
    {
        sm[k]-=(LL)(mx[k]-c)*(LL)ct[k],mx[k]=c; 
        lz[k]=c;
        return;
    }
}
void down(int k)
{
    if(lz[k]!=INF)
    {
        update(tr[k].l,lz[k]),update(tr[k].r,lz[k]);
        lz[k]=INF;
    }
}
void up(int k)
{
    int l=tr[k].l,r=tr[k].r;
    if(mx[l]>mx[r]) ct[k]=ct[l],mx[k]=mx[l],md[k]=max(md[l],mx[r]);
    else if(mx[l]<mx[r]) ct[k]=ct[r],mx[k]=mx[r],md[k]=max(md[r],mx[l]);
    else ct[k]=ct[l]+ct[r],mx[k]=mx[r],md[k]=max(md[l],md[r]);
    sm[k]=sm[l]+sm[r];
}
int nwp(int k)
{
    return (k)?k:++n1;
}
void ins(int k,int l,int r,int w,int v)
{
    if(l==r&&l==w)
    {
        int q=max(mx[k]+v,0);
        sm[k]+=q-mx[k],mx[k]=q,ct[k]=1;
        return;
    }
    int mid=(l+r)/2;
    tr[k].l=nwp(tr[k].l),tr[k].r=nwp(tr[k].r);
    down(k);
    if(w<=mid) ins(tr[k].l,l,mid,w,v);
    else ins(tr[k].r,mid+1,r,w,v);
    up(k);
}
LL find(int k,int l,int r,int x,int y)
{
    x=max(x,l),y=min(y,r);
    if(x>y||!k||sm[k]==0) return 0;
    if(l==x&&r==y) return sm[k];
    int mid=(l+r)/2;
    down(k);
    return find(tr[k].l,l,mid,x,y)+find(tr[k].r,mid+1,r,x,y);
}
void make(int k,int l,int r,int x,int y,int c)
{
    x=max(x,l),y=min(y,r);
    if(x>y||!k) return;
    if(l==x&&r==y)
    {
        if(c>=mx[k]) return;
        if(md[k]<c) 
        {
            sm[k]-=(LL)(mx[k]-c)*(LL)ct[k],mx[k]=c; 
            lz[k]=c;
            return;
        }
    }
    int mid=(l+r)/2;
    down(k);
    make(tr[k].l,l,mid,x,y,c),make(tr[k].r,mid+1,r,x,y,c);
    up(k);
}
int main()
{
    cin>>n>>m;
    n1=1;
    fo(i,1,n) 
        read(a[i]),ins(1,1,n,i,a[i]);
    LL lst=0;
    memset(lz,107,sizeof(lz));
    fo(i,1,m)
    {
        LL p,x,y,c;
        read(p),read(x),read(y);
        x^=lst,y^=lst;
        if(p==1) ins(1,1,n,x,-y);
        else if(p==2)
        {
            read(c);
            c^=lst;
            make(1,1,n,x,y,c);
        }
        else
        {
            lst=find(1,1,n,x,y);
            printf("%lld\n",lst);
        }
    }
}