[bzoj 1911--Apio2010]特别行动队

124 篇文章 2 订阅
26 篇文章 0 订阅

这里写图片描述
这道题,一看就知道是dp了,但范围之大提醒着我们这道题并不简单。对,这题要用斜率优化,还比较容易的那种。
先设i>j>k,
f[j]+a*(sum[j]-sum[i])^2+b*(sum[j]-sum[i])+c>=f[k]+a*(sum[k]-sum[i])^2+b*(sum[k]-sum[i])+c
化简得
[(f[j]+a*sum[j]^2-b*sum[j])-(f[k]+a*sum[k]^2-b*sum[k])]/[2*a*(sum[j]-sum[k])]<=sum[i]
设y[i]=f[i]+a*sum[i]^2-b*sum[i],x[i]=2*a*sum[i]
(y[j]-y[k])/(x[j]-x[k])<=sum[i]
斜率形式出来了,这题便解决了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
int v[1100000];
long long a,b,c,f[1100000],sum[1100000];
long long YY(int x){return f[x]+a*sum[x]*sum[x]-b*sum[x];}
long long Y(int x,int y){return YY(y)-YY(x);}
long long XX(int x){return 2*a*sum[x];}
long long X(int x,int y){return XX(y)-XX(x);}
int main()
{
    int n;
    scanf("%d",&n);
    scanf("%lld%lld%lld",&a,&b,&c);
    for(int i=1;i<=n;i++)
    {
        int x;
        scanf("%d",&x);
        sum[i]=sum[i-1]+x;
    }
    int head=1,tail=1;
    for(int i=1;i<=n;i++)
    {
        while(head<tail && Y(v[head],v[head+1])>=sum[i]*X(v[head],v[head+1]))head++;
        int x=v[head];
        f[i]=f[x]+a*(sum[i]-sum[x])*(sum[i]-sum[x])+b*(sum[i]-sum[x])+c;
        while(head<tail && Y(v[tail],i)*X(v[tail-1],v[tail])<Y(v[tail-1],v[tail])*X(v[tail],i))tail--;
        v[++tail]=i;
    }
    printf("%lld\n",f[n]);
    return 0;
}
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值