BZOJ3437: 小P的牧场

3437: 小P的牧场

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 249  Solved: 153
[Submit][Status]

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

4
2424
3142

Sample Output

9

样例解释

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

数据范围与约定


对于100%的数据，1<=n<=1000000,0<ai,bi<=10000

HINT

 

Source

KpmCup#0 By Greens

题解：

设f[i]表示覆盖1-i的最小费用，则状态转移方程：

f[i]=min(f[j]+sigma((i-k)*b[k]))+a[i]  k=j+1..i-1 

我们令 s1[i]=sigma(b[j]) j<=i  s2[i]=sigma(j*b[j]) j<=i

则上式可以化为 f[i]=min(f[j]+i*(s1[i]-s1[j])-(s2[i]-s2[j]))+a[i]

然后考虑j，k两个决策，j比k优，通过变形得到j，k之间满足的条件。

然后就是喜闻乐见的斜率优化了。

代码：

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 1000000+5
#define maxm 500+100
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    return x*f;
}
int n,q[maxn];
ll a[maxn],b[maxn],c[maxn],f[maxn];
inline double k(int i,int j)
{
    return (double)(f[i]+b[i]-f[j]-b[j])/(double)(a[i]-a[j]);
}
int main()
{
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    n=read();
    for1(i,n)c[i]=read();
    for1(i,n)
    {
        ll x=read();
        a[i]=a[i-1]+x;b[i]=b[i-1]+(ll)i*x;
    }
    int l=1,r=1;
    for1(i,n)
    {
        while(l<r&&k(q[l+1],q[l])<i)l++;
        f[i]=f[q[l]]+i*(a[i]-a[q[l]])-(b[i]-b[q[l]])+c[i];
        while(l<r&&k(q[r],q[r-1])>k(i,q[r]))r--;
        q[++r]=i;
    }
    cout<<f[n]<<endl;
    return 0;
}