【bzoj 3437】小P的牧场

Description

 背景

    小P是个特么喜欢玩MC的孩纸。。。

 描述

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

Input

    第一行一个整数 n 表示牧场数目

    第二行包括n个整数，第i个整数表示ai

    第三行包括n个整数，第i个整数表示bi

Output

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

Sample Input

Sample Output

斜率优化，具体见代码

#include<cstdio>using namespace std;#define maxn 1000050typedef long long ll;ll n;ll a[maxn],b[maxn];ll c[maxn],d[maxn],f[maxn];int l,r,q[maxn<<1];ll read(){    ll ret=0,flag=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')flag=-1;ch=getchar();}    while(ch>='0'&&ch<='9'){ret=ret*10+ch-'0';ch=getchar();}    return ret*flag;}double calc(int j,int k){return (double)(f[j]-f[k]+d[j]-d[k])/(double)(c[j]-c[k]);}int main(){    n=read();    for(int i=1;i<=n;i++)a[i]=read();    for(int i=1;i<=n;i++)b[i]=read();    for(int i=1;i<=n;i++)c[i]=c[i-1]+b[i];    for(int i=1;i<=n;i++)d[i]=d[i-1]+(ll)i*b[i];    for(int i=1;i<=n;i++){        while(l<r && calc(q[l],q[l+1])<i)l++;        int now = q[l];        f[i]=f[now]+(c[i]-c[now])*i-d[i]+d[now]+a[i];        while(l<r && calc(i,q[r]) < calc(q[r],q[r-1])) r--;        q[++r]=i;    }    PRintf("%lld",f[n]);    return 0;}