首页 > 代码库 > 【bzoj 3437】小P的牧场

【bzoj 3437】小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
2424
3142

Sample Output

9

HINT

选取牧场1,3,4建立控制站,最小费用为2+(2+1*1)+4=9。
1<=n<=1000000, 0 < a i ,bi < = 10000

 

 

技术分享
 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 const int N=1000050;
 6 long long n,a[N],b[N],sum[N],num[N],f[N],q[N],l=0,r=0;
 7 int read()
 8 {
 9     int x=0,f=1;char c=getchar();
10     while(c<0||c>9){if(c==-)f=-1;c=getchar();}
11     while(c>=0&&c<=9){x=x*10+c-0;c=getchar();}
12     return x*f;
13 }
14 double slope(long long k,long long j)
15 {
16     return (double)(f[j]-sum[j+1]+num[j+1]*n-f[k]+sum[k+1]-num[k+1]*n)/(double)(num[j+1]-num[k+1]);
17 }
18 int main()
19 {
20     n=read();
21     for(int i=1;i<=n;i++)a[i]=read();
22     for(int i=1;i<=n;i++)b[i]=read();
23     for(long long i=1;i<=n;i++)
24     {
25         sum[i]=sum[i-1]+(long long)b[i-1]*(n-i+1);
26         num[i]=num[i-1]+b[i-1];
27     }
28     for(long long i=1;i<=n;i++)
29     {
30         while(l<r&&slope(q[l],q[l+1])<i)l++;
31         long long t=q[l];
32         f[i]=f[t]+sum[i]-sum[t+1]-(num[i]-num[t+1])*(n-i)+a[i];
33         while(l<r&&slope(q[r],i)<slope(q[r-1],q[r]))r--;
34         q[++r]=i;
35     }
36     printf("%lld\n",f[n]);
37     return 0;
38 }
View Code

 

 

 

【bzoj 3437】小P的牧场