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

bzoj3156防御准备

时间:2016-07-24 17:38:24      阅读:116      评论:0      收藏:0      [点我收藏+]

标签:

 

bzoj3156防御准备

题意:

n个检查点,在第i个检查点放置塔花费a[i],放置木偶花费为该位置右边最近一个塔离它的距离。求最小花费。n≤1000000

题解:

从右往左处理。在第i个点放塔的费用f[i]=min{f[j]+sigma(k,i+1,j-1)k-i}+a[i],用等差数列求和公式化简后作斜率dp,具体看代码。反思:本弱公式总是推错,要稳!

代码:

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #define ll long long
 5 #define maxn 1000010
 6 #define inc(i,j,k) for(int i=j;i<=k;i++)
 7 using namespace std;
 8 
 9 int n,q[maxn],l,r; ll c[maxn],f[maxn];
10 inline int read(){
11     char ch=getchar(); int x=0,f=1;
12     while(ch<0||ch>9){if(ch==-)f=-1; ch=getchar();} while(ch>=0&&ch<=9)x=x*10+ch-0,ch=getchar();
13     return x*f;
14 }
15 double calc(int j,int k){
16     return (double)(2*f[j]-2*f[k]+j-k+(ll)j*j-(ll)k*k)/(double)(j-k);
17 }
18 int main(){
19     n=read(); inc(i,1,n)c[i]=(ll)read(); l=r=0;
20     inc(i,1,n){
21         while(l<r&&calc(q[l],q[l+1])<2*i)l++; f[i]=f[q[l]]+((ll)i-q[l])*(i-q[l]-1)/2+c[i];
22         while(l<r&&calc(q[r-1],q[r])>calc(q[r],i))r--; q[++r]=i;
23     }
24     printf("%lld",f[n]);
25 }

 

20160623

bzoj3156防御准备

标签:

原文地址:http://www.cnblogs.com/YuanZiming/p/5701086.html

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