码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj3437 小P的牧场

时间:2016-01-22 17:12:49      阅读:127      评论:0      收藏:0      [点我收藏+]

标签:

题目链接

斜率优化n^2的DP

 1 #include<cstdio>
 2 #include<algorithm>
 3 using namespace std;
 4 int yp()
 5 {
 6     int x=0,f=0;char ch=getchar();
 7     while(ch<0||ch>9){if(ch==-)f=1;ch=getchar();}
 8     while(ch>=0&&ch<=9)x*=10,x+=ch-0,ch=getchar();
 9     return !f?x:-x;
10 }
11 int n,l=1,r,q[1000005];
12 long long a[1000005],b[1000005],tot,ans,h[1000005],sum[1000005];
13 double js(int a,int b)
14 {
15     return (1.*(h[b]-h[a]))/(1.*(a-b));
16 }
17 int main()
18 {
19     n=yp();
20     for(int i=1;i<=n;i++)a[i]=yp();
21     for(int i=1;i<=n;i++)b[i]=yp(),sum[i]=sum[i-1]+b[i];
22     for(int i=1;i<n;i++)tot+=b[i]*(n-i);
23     tot+=a[n],q[++r]=n;
24     for(int i=n-1;i;i--)
25     {
26         while(l<r&&js(q[l],q[l+1])>sum[i])l++;
27         int j=q[l];
28         h[i]=h[j]+sum[i]*(j-i)-a[i];
29         ans=max(ans,h[i]);
30         while(l<r&&js(q[r],i)>js(q[r-1],q[r]))r--;
31         q[++r]=i;
32     }
33     printf("%lld",tot-ans);
34     return 0;
35 }

 

bzoj3437 小P的牧场

标签:

原文地址:http://www.cnblogs.com/HugeGun/p/5151442.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!