[bzoj3437][dp+斜率优化]小P的牧场

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
2 4 2 4
3 1 4 2

Sample Output

9

HINT

样例解释

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

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

题解

这题一看其实和3156防御准备差不多。。只不过这里第n个点一定要建站
方程设f[i]为第i个点建站 1~i的最小花费
再开两个数组 sum[i]表示b[i]*i的前缀和 q[i]表示b[i]的前缀和
那么f[i]=min(f[j]+(sum[i-1]-sum[j])* i - j*(q[i-1]-q[j]))
化式子和3156还是差不多的。。
画出来是这个样子(f[k]-f[j]+sum[k]-sum[j])/(q[k]-q[j])＜i
然后开一下longlong就过了

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
LL a[1110000],b[1110000];
LL q[1110000],sum[1110000],n;
LL f[1110000],list[1110000],head,tail;
double slop(int k,int j){return double((LL)(f[k]-f[j]+sum[k]-sum[j]))/double((LL)(q[k]-q[j]));}
int main()
{
    q[0]=sum[0]=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&b[i]);
        q[i]=q[i-1]+b[i];
        sum[i]=sum[i-1]+b[i]*i;
    }
    memset(f,63,sizeof(f));
    f[0]=0;
    head=tail=1;list[1]=0;
    for(int i=1;i<=n;i++)
    {
        while(head<tail && slop(list[head+1],list[head])<double(i))head++;
        f[i]=f[list[head]]+(q[i-1]-q[list[head]])*i-(sum[i-1]-sum[list[head]])+a[i];
        while(head<tail && slop(list[tail],list[tail-1])>slop(i,list[tail]))tail--;
        list[++tail]=i;
    }
    printf("%lld\n",f[n]);
    return 0;
}