ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

bzoj3156 防御准备(斜率优化)

2019-09-10 21:51:51  阅读:233  来源: 互联网

标签:frac int ll bzoj3156 斜率 防御 Input 整数 getchar


Time Limit: 10 Sec  Memory Limit: 512 MB

 

Input

第一行为一个整数N表示战线的总长度。

第二行N个整数,第i个整数表示在位置i放置守卫塔的花费Ai。

Output

共一个整数,表示最小的战线花费值。

Sample Input

10
2 3 1 5 4 5 6 3 1 2

Sample Output

18

HINT

1<=N<=10^6,1<=Ai<=10^9


 

设$f[i]$为已经完成$i+1$到$n$,在$i$放塔的最小代价

$f[i]=a[i]+min(f[j]+ \frac {(j-i-1)(j-i)}{2}),i<j<=n$

移项转化,得

$f[j]+\frac{j(j-1)}{2}=ij+f[i]-\frac{i(i+1)}{2}$

$k=i,x=j$单调

凸包维护单调队列,复杂度$O(n)$

#include<cstdio>
typedef long long ll;
inline ll min(ll A,ll B){return A<B?A:B;}
inline ll max(ll A,ll B){return A>B?A:B;}
int read(){
    char c=getchar(); int x=0;
    while(c<'0'||c>'9') c=getchar();
    while('0'<=c&&c<='9') x=x*10+c-48,c=getchar();
    return x;
}
#define N 1000005
int n,a[N],h[N],L,R; ll f[N],ans=1e16;
inline ll y(ll j){return f[j]+j*(j-1)/2;}
inline int chk(ll A,ll B,ll kx,ll ky){return (y(A)-y(B))*kx>=ky*(A-B);}
int main(){
    R=n=read();
    for(register int i=1;i<=n;++i) a[i]=read();
    f[n]=a[n]; h[L=R=1]=n;
    for(register ll i=n-1;i;--i){
        while(L<R&&chk(h[L],h[L+1],1,i)) ++L;
        f[i]=a[i]+f[h[L]]+(h[L]-i)*(h[L]-i-1)/2;
        while(L<R&&chk(h[R],h[R-1],h[R]-i,y(h[R])-y(i))) --R;
        h[++R]=i;
        ans=min(ans,f[i]+i*(i-1)/2);
    }printf("%lld",ans);
    return 0;
}

 

标签:frac,int,ll,bzoj3156,斜率,防御,Input,整数,getchar
来源: https://www.cnblogs.com/kafuuchino/p/11503195.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有