BZOJ3437 小P的牧场【斜率优化DP】

Time Limit: 10 Sec Memory Limit: 128 MB

Description

小P在MC里有n个牧场，自西向东呈一字形排列（自西向东用1…n编号），于是他就烦恼了：为了控制这n个牧场，他需要在某些牧场上面建立控制站，每个牧场上只能建立一个控制站，每个控制站控制的牧场是它所在的牧场一直到它西边第一个控制站的所有牧场（它西边第一个控制站所在的牧场不被控制）（如果它西边不存在控制站，那么它控制西边所有的牧场），每个牧场被控制都需要一定的花费（毕竟在控制站到牧场间修建道路是需要资源的嘛~），而且该花费等于它到控制它的控制站之间的牧场数目（不包括自身，但包括控制站所在牧场）乘上该牧场的放养量，在第i个牧场建立控制站的花费是ai，每个牧场i的放养量是bi，理所当然，小P需要总花费最小，但是小P的智商有点不够用了，所以这个最小总花费就由你来算出啦。

Input

第一行一个整数 n 表示牧场数目
第二行包括n个整数，第i个整数表示ai
第三行包括n个整数，第i个整数表示bi

Output

只有一行，包括一个整数，表示最小花费

HINT

1<=n<=1000000, 0 < a i ,bi < = 10000

---

题目分析

斜率优化DP–详解

一开始被卡在了${\sum }_{k=1}^{i-j-1}k\ast b_{i-k}$的求解
参考了dalao的博客以后豁然开朗

令$sumB[i]={\sum }_{k=1}^{i}b[k]$
$sumM[i]={\sum }_{k=1}^{i}b[k]\ast k$
那么${\sum }_{k=1}^{i-j-1}k\ast b_{i-k}=i\ast (sumB[i]-sumB[j])-(sumM[i]-sumM[j])$

现在开始考虑dp方程，先确定$O(n^2)$做法
$dp[i]=dp[j]+a[i]+i\ast (sumB[i]-sumB[j])-(sumM[i]-sumM[j])$

然后变形为
$dp[j]+sumM[j]=i\ast sumB[j]-i\ast sumB[i]+sumM[i]-a[i]+dp[i]$
到这里就可以令$dp[j]+sumM[j]$为$y$，$sumB[j]$为$x$，$i$为斜率
单调队列维护凸壳即可

---

#include
#include
#include
#include
#include
#include
using namespace std;
typedef double dd;
typedef long long lt;
    
lt read()
{
    lt f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}
    
const int maxn=2000010;
int n;
lt ai[maxn];
lt sumB[maxn],sumM[maxn];
lt dp[maxn];
int q[maxn],ll,rr;
    
dd calc(int j1,int j2)
{
    lt ty=(dp[j2]+sumM[j2])-(dp[j1]+sumM[j1]);
    lt tx=sumB[j2]-sumB[j1];
    return (dd)ty/(dd)tx;
}
    
int main()
{
    n=read();
    for(int i=1;i<=n;++i) ai[i]=read();
    for(lt i=1;i<=n;++i) 
    {
        lt b=read();
        sumB[i]=sumB[i-1]+b; sumM[i]=sumM[i-1]+i*b;
    }
        
    ll=rr=1;
    for(int i=1;i<=n;++i)
    {
        while( ll<rr && calc(q[ll],q[ll+1])<=i ) ++ll;
        dp[i]=dp[q[ll]]-(sumM[i]-sumM[q[ll]])+i*(sumB[i]-sumB[q[ll]])+ai[i];
        while( ll<rr && calc(q[rr-1],q[rr]) >= calc(q[rr],i) ) --rr;
        q[++rr]=i;
    }
        
    printf("%lld",dp[n]);
    return 0;
}