[BZOJ1911][Apio2010]特别行动队

31 篇文章 0 订阅
8 篇文章 0 订阅

原题地址

斜率优化…

#include <cstdio>
#include <queue>
using namespace std;
typedef long long ll;
const ll N=1000010;
ll n,a,b,c;
ll x[N],s[N],f[N];

struct Vec{
    ll x,y;

    Vec() {}

    Vec(ll x,ll y):x(x),y(y) {}

    friend Vec operator-(Vec a,Vec b){
        return Vec(a.x-b.x,a.y-b.y);
    }

    friend ll operator*(Vec a,Vec b){
        return a.x*b.y-a.y*b.x;
    }
};
deque<Vec> Q;

int main(){
    scanf("%lld%lld%lld%lld",&n,&a,&b,&c);
    for(ll i=1;i<=n;i++) scanf("%lld",&x[i]);
    for(ll i=1;i<=n;i++) s[i]=s[i-1]+x[i];
    for(ll i=1;i<=n;i++){
        ll  A=-s[i];
        Vec tmp=Vec(((a*s[i-1])<<1),f[i-1]+a*s[i-1]*s[i-1]-b*s[i-1]);
        while(Q.size()>1){
            Vec a,b;
            a=Q.front();Q.pop_front();b=Q.front();
            if((tmp-a)*(a-b)<0){
                Q.push_front(a);
                break;
            }
        }
        Q.push_front(tmp);
        while(Q.size()>1){
            Vec a,b;
            a=Q.back();Q.pop_back();b=Q.back();
            if(A*a.x+a.y>A*b.x+b.y){
                Q.push_back(a);
                f[i]=A*a.x+a.y;
                break;
            }
            f[i]=A*b.x+b.y;
        }
        f[i]+=a*s[i]*s[i]+b*s[i]+c;
    }
    printf("%lld\n",f[n]);

    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值