bzoj1911 特别行动队 斜率优化

斜率优化最裸的裸题,但是式子推了好久,忘了形式。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
LL dp[1000005];
LL s[1000005];
LL a,b,c;
LL p(LL x)
{   return x*x;}
double get(int j,int k)
{
    double ans=(double)(dp[j]-dp[k])/(double)(s[j]-s[k]);
    ans+=(double)(a*(s[j]+s[k]));
    return ans;
}
int q[1000005];
int H,T;
int main()
{
    int n;
    scanf("%d",&n);
    cin>>a>>b>>c;
    LL x;
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&x);
        s[i]=s[i-1]+x;
    }
    dp[1]=a*p(s[1])+b*s[1]+c;
    H=1;
    q[++T]=0;
    q[++T]=1;
    for(int i=2;i<=n;i++)
    {
        while(H<T&&get(q[H],q[H+1])>(double)(2*a*s[i]+b))H++;
        dp[i]=dp[q[H]]+a*p(s[i]-s[q[H]])+b*(s[i]-s[q[H]])+c;
        while(H<T&&get(q[T-1],q[T])<get(q[T],i))   T--;
        q[++T]=i;
    }
    printf("%lld",dp[n]);
    return 0;
}
  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值