BZOJ3156 防御准备(动态规划+斜率优化)

  设f[i]为在i放置守卫塔时1~i的最小花费。那么显然f[i]=min(f[j]+(i-j)*(i-j-1)/2)+a[i]。

  显然这是个斜率优化入门题。将不与i、j同时相关的提出,得f[i]=min(f[j]+j*(j+1)/2-ij)+i*(i-1)/2+a[i]。

  套路地,假设j>k且j转移优于k,则f[j]+j*(j+1)/2-ij<f[k]+k*(k+1)/2-ik,(f[j]+j*(j+1)/2-f[k]-k*(k+1)/2)/(j-k)<i。

  维护下凸壳即可。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
#define N 1000010
#define ll long long
int n,a[N],q[N];
ll f[N];
long double calc(int j,int k)
{
    return (long double)(f[j]+(1ll*j*(j+1)>>1)-f[k]-(1ll*k*(k+1)>>1))/(j-k);
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3156.in","r",stdin);
    freopen("bzoj3156.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=1;i<=n;i++) a[i]=read();
    f[0]=0;
    int head=1,tail=1;q[1]=0;
    for (int i=1;i<=n;i++)
    {
        while (head<tail&&calc(q[head],q[head+1])<i) head++;
        f[i]=f[q[head]]+(1ll*q[head]*(q[head]+1)>>1)-1ll*i*q[head]+(1ll*i*(i-1)>>1)+a[i];
        while (head<tail&&calc(q[tail-1],q[tail])>calc(q[tail],i)) tail--;
        q[++tail]=i;
    }
    cout<<f[n];
    return 0;
}

 

转载于:https://www.cnblogs.com/Gloid/p/9665314.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值