bzoj2388 旅行规划（分块+凸包+三分）

题目链接

分析：
路径权值实际上就是序列的前缀和
我一开始觉得可以直接用树状数组维护前缀和，每一个询问直接在序列上三分即可
但是正解远远没有这么简单

这是一道比较难的题：
路径权值虽然是前缀和的形式，但是需要区间修改
（好像树状数组有区间修改区间查询的操作，但是不会啊。。。）
我们不好直接用简答数据结构维护，我们可以考虑将序列分块（有点像块状链表）
每一块内维护$delta,k,s[x]$

$delta$：整个块中每个位置都需要增加的权值

$k$：每个位置都要增加$(x-l+1)*k$（$l$表示该块第一个元素的位置）

举个例子：

分块：[1,4],[5,8],[9,12]
区间[3,10]+z
对于块2,每一个结点都需要增加(5-3)*z,除此之外,块2中的结点x还要增加(x-5+1)*z
所以delta=(5-3)*z,k=z

$sum[x]$：$x$位置的前缀和

每个块内的元素一定是连续的一段
维护每个块控制的左端点和右端点

对于一个块内的最大值，我们相当于查询$delta+max(sum[x]+k*(x-l+1))$
设$y_i=sum[i],x_i=i$，上式可以化为$z=y_i+k*(x_i-l+1)$
可以看做是一条过$(x_i,y_i)$斜率为$-k$的直线与$y$轴截距的最大值
取最大值的点一定在该区域点的上凸壳上，所以我们对于每个块维护凸壳，查询的时候在凸包上三分

修改一个区间时，暴力修改不完整的块，重建凸壳；对于完整的块，我们直接修改$delta,k$
查询一个区间时，暴力得到不完整块的最大值；对于完整的块，在凸壳上三分得到最大值

tip

一定要头脑清晰。。。代码难写
三分不要写错了
调不出来了，就差抄代码了。。。精神AC
AC代码

普通数据结构难以下手的时候，可以考虑分块

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define ll long long

using namespace std;

const ll INF=1e12;
const int N=100005;
struct po{
    int size,a[350];
};
po c[350];
ll sum[N],delta[N],k[N];
int block,cnt,n,m,belong[N],l[N],r[N],S[400];

double F(int x,int y)
{
    return (double)(sum[x]-sum[y])/(x-y);
}

void build(int x)
{
    int top=0;
    S[++top]=l[x];
    for (int i=l[x]+1;i<=r[x];i++)
    {
        while (top>1&&(F(S[top-1],S[top])<F(S[top-1],i))) top--;
        S[++top]=i;
    }
    S[0]=0; S[top+1]=n+1;
    c[x].size=top;
    for (int i=0;i<=top+1;i++) c[x].a[i]=S[i];
}

void push(int x)
{
    for (int i=l[x];i<=r[x];i++)
        sum[i]=sum[i]+delta[x]+(i-l[x]+1)*k[x];
    delta[x]=k[x]=0;
}

ll cal(int num,int x)
{
    return (ll)delta[num]+sum[x]+k[num]*(x-l[num]+1);
}

ll sanfen(int num)
{
    int L=1,R=c[num].size;
    while (L<=R)
    {
        int M=(L+R)>>1;
        ll t1=cal(num,M-1);
        ll t2=cal(num,M);
        ll t3=cal(num,M+1);
        if (t1<t2&&t2<t3) L=M+1;
        else if (t1>t2&&t2>t3) R=M-1;
        else return t2;
    }
}

ll ask(int x,int y)
{
    ll ans=-INF;
    if (belong[x]==belong[y])    //暴力计算 
    {
        for (int i=x;i<=y;i++) ans=max(ans,cal(belong[i],i));
        return ans;
    }
    for (int i=x;i<=r[belong[x]];i++) ans=max(ans,cal(belong[i],i));
    for (int i=belong[x]+1;i<=belong[y]-1;i++) ans=max(ans,sanfen(i));
    for (int i=l[belong[y]];i<=y;i++) ans=max(ans,cal(belong[i],i));
    return ans;
}

void change(int x,int y,ll z)
{
    if (belong[x]==belong[y])
    {
        for (int i=x;i<=y;i++) sum[i]+=(ll)(i-x+1)*z;
        for (int i=y+1;i<=r[belong[y]];i++) sum[i]+=(ll)(y-x+1)*z;
        for (int i=belong[y]+1;i<=cnt;i++) delta[i]+=(ll)(y-x+1)*z;
        push(belong[x]); build(belong[x]);
    }
    for (int i=belong[x]+1;i<=belong[y]-1;i++) 
        delta[i]+=(l[i]-x)*z,k[i]+=z;    //
    for (int i=x;i<=r[belong[x]];i++) sum[i]+=(ll)(i-x+1)*z;
    push(belong[x]); build(belong[x]);
    for (int i=l[belong[y]];i<=y;i++) sum[i]+=(ll)(i-x+1)*z;
    for (int i=y+1;i<=r[belong[y]];i++) sum[i]+=(ll)(y-x+1)*z;
    for (int i=belong[y]+1;i<=cnt;i++) delta[i]+=(ll)(y-x+1)*z;
    push(belong[y]); build(belong[y]);
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) 
    {
        scanf("%lld",&sum[i]);
        sum[i]+=sum[i-1];
    }
    sum[0]=sum[n+1]=-INF;

    block=sqrt(n);                //分块 
    cnt=ceil((double)n/block);
    for (int i=1;i<=cnt;i++) l[i]=n,r[i]=0;
    for (int i=1;i<=n;i++)
    {
        belong[i]=(i-1)/block+1;
        l[belong[i]]=min(l[belong[i]],i);   //每个块管辖的区间 
        r[belong[i]]=max(r[belong[i]],i);
    }
    for (int i=1;i<=cnt;i++) build(i);

    scanf("%d",&m);
    int opt,x,y,z;
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&opt,&x,&y);
        if (opt==1)
            printf("%lld\n",ask(x,y));
        else
        {
            scanf("%d",&z);
            change(x,y,(ll)z);
        }
    }

    return 0;
}