BZOJ5089 最大连续子段和（分块）

　　假设所有操作都是对整个序列的。考虑每个子区间，区间和与其被加的值构成一次函数关系。最大子段和相当于多个子区间取最大值，答案显然就在这些一次函数构成的下凸壳上。如果预处理出凸壳，只要在凸壳上暴力跳就可以回答询问了，因为加的都是正数，并且斜率不同的一次函数数量是O(n)的。暴力建凸壳的复杂度是O(n²)的。

　　那么考虑分块。每个块预处理出凸壳。区间加时，对于整块在凸壳上暴跳，零散部分暴力重构。查询时合并区间，类似于单点加的线段树做法，维护块内最大前缀和、最大后缀和。块大小取n^1/3时最优。

　　造凸壳写挂调了1h，退役了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 50010
#define NUM 2000
#define BLOCK 100
#define inf 10000000000000000ll
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,m,pos[N],num,block;
ll a[N],tmp[N],PRE[N],SUF[N];
struct line
{
    ll k,b;
    ll f(ll x){return k*x+b;}
};
ll cross(line x,line y){return (x.b-y.b-)/(y.k-x.k)+;}
struct hull
{
    line a[BLOCK];int cnt,cur;
    void ins(line x){while (cnt&&x.b>=a[cnt].b||cnt>&&x.f(cross(a[cnt],a[cnt-]))>a[cnt].f(cross(a[cnt],a[cnt-]))) cnt--;a[++cnt]=x;}
    void clear(){cnt=,cur=,ins((line){,});}
    void jump(ll x){while (cur<cnt&&a[cur].f(x)<a[cur+].f(x)) cur++;}
    ll get(ll x){return a[cur].f(x);}
};
struct data
{
    hull pre,suf,seq;int L,R;ll lazy,sum;
    void build()
    {
        for (int i=L;i<=R;i++) a[i]+=lazy;lazy=;
        sum=;for (int i=L;i<=R;i++) sum+=a[i];
        ll x=;pre.clear();
        for (int i=L;i<=R;i++) pre.ins((line){i-L+,x+=a[i]});
        x=;suf.clear();
        for (int i=R;i>=L;i--) suf.ins((line){R-i+,x+=a[i]});
        seq.clear();
        for (int k=;k<=R-L+;k++)
        {
            ll s=-inf;x=;for (int i=L;i<=L+k-;i++) x+=a[i];
            for (int i=L+k-;i<=R;i++) s=max(s,x),x+=a[i+]-a[i-k+];
            seq.ins((line){k,s});
        }
    }
    void add(int x){sum+=1ll*(R-L+)*x,lazy+=x,pre.jump(lazy),suf.jump(lazy),seq.jump(lazy);}
    ll maxpre(){return pre.get(lazy);}
    ll maxsuf(){return suf.get(lazy);}
    ll maxseq(){return seq.get(lazy);}
}f[NUM];
ll getseq(int l,int r)
{
    ll s=,ans=;
    for (int i=l;i<=r;i++)
    {
        s+=tmp[i];
        if (s<) s=;
        ans=max(ans,s);
    }
    return ans;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj5089.in","r",stdin);
    freopen("bzoj5089.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();block=*pow(n,1.0/)+;num=(n-)/block+;
    for (int i=;i<=n;i++) a[i]=read();
    for (int i=;i<=num;i++)
    {
        f[i].L=f[i-].R+,f[i].R=min(n,f[i].L+block-);
        for (int j=f[i].L;j<=f[i].R;j++) pos[j]=i;
        f[i].build();
    }
    while (m--)
    {
        char c=getc();
        if (c=='A')
        {
            int l=read(),r=read(),x=read();
            if (pos[l]==pos[r])
            {
                for (int i=l;i<=r;i++) a[i]+=x;
                f[pos[l]].build();
            }
            else
            {
                for (int i=pos[l]+;i<pos[r];i++) f[i].add(x);
                for (int i=l;i<=f[pos[l]].R;i++) a[i]+=x;f[pos[l]].build();
                for (int i=f[pos[r]].L;i<=r;i++) a[i]+=x;f[pos[r]].build();
            }
        }
        else
        {
            int l=read(),r=read();
            if (pos[l]==pos[r])
            {
                for (int i=l;i<=r;i++) tmp[i]=a[i]+f[pos[l]].lazy;
                printf(LL,getseq(l,r));
            }
            else
            {
                for (int i=l;i<=f[pos[l]].R;i++) tmp[i]=a[i]+f[pos[l]].lazy;
                for (int i=f[pos[r]].L;i<=r;i++) tmp[i]=a[i]+f[pos[r]].lazy;
                ll ans=max(getseq(l,f[pos[l]].R),getseq(f[pos[r]].L,r)),s;
                for (int i=pos[l]+;i<pos[r];i++) ans=max(ans,f[i].maxseq());
                for (int i=pos[l];i<=pos[r];i++) PRE[i]=SUF[i]=;
                s=;for (int i=f[pos[l]].R;i>=l;i--) s+=tmp[i],SUF[pos[l]]=max(SUF[pos[l]],s);
                s=;for (int i=f[pos[r]].L;i<=r;i++) s+=tmp[i],PRE[pos[r]]=max(PRE[pos[r]],s);
                for (int i=pos[l]+;i<pos[r];i++) SUF[i]=max(f[i].maxsuf(),SUF[i-]+f[i].sum);
                for (int i=pos[r]-;i>pos[l];i--) PRE[i]=max(f[i].maxpre(),PRE[i+]+f[i].sum);
                for (int i=pos[l];i<pos[r];i++) ans=max(ans,SUF[i]+PRE[i+]);
                printf(LL,ans);
            }
        }
    }
    return ;
}