BZOJ1911[APIO2010] 特别行动队

原题链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1911

特别行动队

题目描述

你有一支由 n 名预备役士兵组成的部队,士兵从 1 到 n 编号,要将他们拆分 成若干特别行动队调入战场。出于默契的考虑,同一支特别行动队中队员的编号 应该连续,即为形如 (i,i+1,...,i+k) ( i , i + 1 , . . . , i + k ) 的序列。 编号为 i 的士兵的初始战斗力为 xi ,一支特别行动队的初始战斗力 x 为队内 士兵初始战斗力之和,即 x=xi+xi+1+...+xi+k x = x i + x i + 1 + . . . + x i + k

通过长期的观察,你总结出一支特别行动队的初始战斗力 x 将按如下经验公 式修正为 x:x=ax2+bx+cx x : x ′ = a x 2 + b x + c x ,其中 a, b, c 是已知的系数(a < 0)。 作为部队统帅,现在你要为这支部队进行编队,使得所有特别行动队修正后 战斗力之和最大。试求出这个最大和。

例如,你有 4 名士兵, x1=2,x2=2,x3=3,x4=4 x 1 = 2 , x 2 = 2 , x 3 = 3 , x 4 = 4 。经验公式中的参数为 a = –1, b = 10, c = –20。此时,最佳方案是将士兵组成 3 个特别行动队:第一队包含士兵 1 和士兵 2,第二队包含士兵 3,第三队包含士兵 4。特别行动队的初始战斗力分 别为 4, 3, 4,修正后的战斗力分别为 4, 1, 4。修正后的战斗力和为 9,没有其它 方案能使修正后的战斗力和更大。

输入输出格式
输入格式:

输入由三行组成。第一行包含一个整数 n,表示士兵的总数。第二行包含三 个整数 a, b, c,经验公式中各项的系数。第三行包含 n 个用空格分隔的整数 x1,x2,,xn x 1 , x 2 , … , x n ,分别表示编号为 1, 2, …, n 的士兵的初始战斗力。

输出格式:

输出一个整数,表示所有特别行动队修正后战斗力之和的最大值。

输入输出样例
输入样例#1:

4
-1 10 -20
2 2 3 4

输出样例#1:

9

说明

20%的数据中,n ≤ 1000;

50%的数据中,n ≤ 10,000;

100%的数据中,1 ≤ n ≤ 1,000,000,–5 ≤ a ≤ –1,|b| ≤ 10,000,000,|c| ≤ 10,000,000,1 ≤ xi ≤ 100

题解

一看就是裸的斜率优化,直接推公式啊:

val[i] v a l [ i ] 为战斗力前缀和,设 (k<j) ( k < j ) k k 优于j,易得:  


dp[k]+a(val[i]val[k])2+b(val[i]val[k])+c>dp[j]+a(val[i]val[j])2+b(val[i]val[j])+c d p [ k ] + a ( v a l [ i ] − v a l [ k ] ) 2 + b ( v a l [ i ] − v a l [ k ] ) + c > d p [ j ] + a ( v a l [ i ] − v a l [ j ] ) 2 + b ( v a l [ i ] − v a l [ j ] ) + c
dp[k]+a(val2[i]+val2[k]2val[i]val[k])+bval[i]bval[k]+c>dp[j]+a(val2[i]+val2[j]2val[i]val[j])+bval[i]bval[k]+c d p [ k ] + a ( v a l 2 [ i ] + v a l 2 [ k ] − 2 v a l [ i ] v a l [ k ] ) + b v a l [ i ] − b v a l [ k ] + c > d p [ j ] + a ( v a l 2 [ i ] + v a l 2 [ j ] − 2 v a l [ i ] v a l [ j ] ) + b v a l [ i ] − b v a l [ k ] + c
dp[k]dp[j]+a(val2[k]val2[j])+b(val[j]val[k])>2aval[i](val[k]val[j]) d p [ k ] − d p [ j ] + a ( v a l 2 [ k ] − v a l 2 [ j ] ) + b ( v a l [ j ] − v a l [ k ] ) > 2 a v a l [ i ] ( v a l [ k ] − v a l [ j ] )
dp[k]dp[j]+a(val2[k]val2[j])+b(val[j]val[k])2a(val[k]val[j])>val[i] d p [ k ] − d p [ j ] + a ( v a l 2 [ k ] − v a l 2 [ j ] ) + b ( v a l [ j ] − v a l [ k ] ) 2 a ( v a l [ k ] − v a l [ j ] ) > v a l [ i ]

val[i] v a l [ i ] 单增,完结撒花~~~

代码
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=1e6+5;
int n,a,b,c,que[M];
ll val[M],dp[M];
void in()
{
    scanf("%d%d%d%d",&n,&a,&b,&c);
    for(int i=1;i<=n;++i)
    scanf("%d",&val[i]),val[i]+=val[i-1];
}
ll sqr(ll x){return x*x;}
double slop(int k,int j)
{return 1.0*(dp[k]-dp[j]+a*(sqr(val[k])-sqr(val[j]))+b*(val[j]-val[k]))/(2*a*(val[k]-val[j]));}
void ac()
{
    int le=0,ri=0,k;
    for(int i=1;i<=n;++i)
    {
        while(le<ri&&slop(que[le],que[le+1])<val[i])le++;
        k=que[le];
        dp[i]=dp[k]+a*sqr(val[i]-val[k])+b*(val[i]-val[k])+c;
        while(le<ri&&slop(que[ri],i)<slop(que[ri-1],que[ri]))ri--;
        que[++ri]=i;
    }
    printf("%lld",dp[n]);
}
int main()
{
    in();ac();
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

ShadyPi

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值