bzoj3437 小P的牧场

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

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

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

Sample Input
4

2424

3142

Sample Output
9

样例解释

选取牧场1，3，4建立控制站，最小费用为2+(2+1*1)+4=9。

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

HINT

Source
KpmCup#0 By Greens
基本同锯木厂选址那题只不过这题距离变成1 更简单一些的
f[i]=min {f[j]+c[i]-c[j]-sw[i]*(i-j)| 0<=j< i +a[i]}。表示在i建控制站，上一个控制站建在j所需的最小花费。假设k1< k2且k1优于k2，则有(f[k2]-f[k1]+c[k1]-c[k2]+sw[k2]*k2-sw[k1]*k1)/(sw[k2]-sw[k1])>i。因此我们维护一个下凸曲线，保证队列中的斜率都是单增的。
设sw[i]表示1~i的放养量总和 那么我可以o(n)算出cost[i]也就是我把1~i的只在i号位置建立一个牧场 我的总代价是多少 cost[i]=cost[i-1]+sw[i-1]
有了这样两个数组我就可以求出我将j~i全部用i控制的代价了 那就是cost[i]-cost[j]-sw[j]*(i-j) 我首先将他们认为1~j全部由j号点控制 那么我只想要j~i由i控制 那么显然对于1~i全部由i控制我还多了sw[j]由i控制的代价

#include<cstdio>
#define ll long long
#define N 1100000
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x;
}
int n,a[N],q[N];ll f[N],cost[N],sw[N];
inline double slope(int j1,int j2){
    return double(f[j1]-f[j2]+cost[j2]-cost[j1]+sw[j1]*j1-sw[j2]*j2)/(sw[j1]-sw[j2]);
}
int main(){
    freopen("bzoj3437.in","r",stdin);
    n=read();int b;
    for (int i=1;i<=n;++i) a[i]=read();
    for (int i=1;i<=n;++i) b=read(),sw[i]=sw[i-1]+b;
    for (int i=1;i<=n;++i) cost[i]=cost[i-1]+sw[i-1];
    int l=1,r=0;q[++r]=0;
    for (int i=1;i<=n;++i){
        while(l<r&&slope(q[l],q[l+1])<i) ++l;int j=q[l];
        f[i]=f[j]+cost[i]-cost[j]-sw[j]*(i-j)+a[i];
        while(l<r&&slope(q[r-1],q[r])>slope(q[r],i)) --r;q[++r]=i;
    }printf("%lld",f[n]);
    return 0;
}