题目描述

  你有一支由n名预备役士兵组成的部队,士兵从1到n编号,要将他们拆分成若干特别行动队调入战场。出于默契考虑,同一支特别行动队中队员的编号应该连续,即为形如(i,i+1,…,i+k)的序列。
  编号为i的士兵的初始战斗力为xi,一支特别运动队的初始战斗力x为队内士兵初始战斗力之和,即x=(xi)+(xi+1)+…+(xi+k)。
  通过长期的观察,你总结出一支特别行动队的初始战斗力x将按如下经验公式修正为x’:x’=ax^2+bx+c,其中a,b,c是已知的系数(a<0)。
  作为部队统帅,现在你要为这支部队进行编队,使得所有特别行动队修正后战斗力之和最大。试求出这个最大和。
  例如,你有4名士兵,x1=2,x2=2,x3=3,x4=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,分别表示编号为1,2,…,n的士兵的初始战斗力。


输出格式

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


样例数据

样例输入

4
-1 10 -20
2 2 3 4

样例输出

9


数据范围

20%的数据中,n<=1000;
50%的数据中,n<=10000;
100%的数据中,1<=n<=1000000,-5<=a<=-1,|b|<=10000000,|c|<=10000000,1<=xi<=100。


题目分析

被这题坑惨了。
容易写出动规方程
fi=max fj+a(sumi−sumj)2+b(sumi−sumj)+c f_i=max f_j+a(sum_i-sum_j)^2+b(sum_i-sum_j)+c
考虑斜率优化
化简为
fi−asum2i−bsumi−c=fj+asum2j−2asumisumj−bsumj f_i-asum_i^2-bsum_i-c=f_j+asum_j^2-2asum_isum_j-bsum_j
令 fj+asum2j−bsumj=y f_j+asum_j^2-bsum_j =y
sumi=k sum_i=k
2asumj=x 2asum_j=x
因为k单调递增且为正,维护下凸包
当然也可以令 2asumi=k 2asum_i=k sumj=x sum_j=x
但因为a为负,所以k单调递减,维护上凸包

一定要看数据范围啊,a为负


源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
typedef long long LL;
inline const LL Get_Int() {LL num=0,bj=1;char x=getchar();while(x<'0'||x>'9') {if(x=='-')bj=-1;x=getchar();}while(x>='0'&&x<='9') {num=num*10+x-'0';x=getchar();}return num*bj;
}
LL n,f[1000005],Cost[1000005],sum[1000005],Q[1000005],ans=0,Cut=0;
double Slope(int j,int k) {return (double)(f[j]-f[k])/(k-j);
}
int main() {cin>>n;for(int i=1; i<=n; i++)Cost[i]=Get_Int();for(int i=1; i<=n; i++)sum[i]=sum[i-1]+Get_Int();for(int i=1; i<n; i++)ans+=(sum[i]-sum[i-1])*(n-i);ans+=Cost[n];int Left=1,Right=1;Q[1]=n;for(int i=n-1; i>=1; i--) {while(Left<Right&&Slope(Q[Left],Q[Left+1])>=sum[i])Left++; //维护队首(删除非最优决策)int Front=Q[Left];f[i]=f[Front]+sum[i]*(Front-i)-Cost[i];Cut=max(Cut,f[i]);while(Left<Right&&Slope(Q[Right-1],Q[Right])<=Slope(Q[Right],i))Right--; //维护队尾(维护上凸包性质)Q[++Right]=i;}printf("%lld\n",ans-Cut);return 0;
}

[APIO2010] 特别行动队相关推荐

  1. bzoj 1911: [Apio2010]特别行动队 -- 斜率优化

    1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Description Input Output Sample Input 4 ...

  2. bzoj 1911: [Apio2010]特别行动队 2011-12-26

    1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Submit: 892  Solved: 359 [Submit][Statu ...

  3. BZOJ 1911: [Apio2010]特别行动队 [斜率优化DP]

    1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Submit: 4142  Solved: 1964 [Submit][Sta ...

  4. 1911: [Apio2010]特别行动队

    1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Submit: 4061  Solved: 1922 [Submit][Sta ...

  5. [BZOJ] 1911 [Apio2010]特别行动队

    题目描述 你有一支由 n 名预备役士兵组成的部队,士兵从 1 到 n 编号,要将他们拆分 成若干特别行动队调入战场.出于默契的考虑,同一支特别行动队中队员的编号 应该连续,即为形如 (i, i + 1 ...

  6. [APIO2010]特别行动队

    题面在这里 题意 你有一支由\(n\)名预备役士兵组成的部队,士兵从\(1\)到\(n\)编号,要将他们拆分 成若干特别行动队调入战场. 出于默契的考虑,同一支特别行动队中队员的编号应该连续,即为形如 ...

  7. 【斜率优化】[APIO2010]特别行动队

    题目 题目描述 你有一支由 nn 名预备役士兵组成的部队,士兵从 11 到 nn 编号,你要将他们拆分成若干特别行动队调入战场.出于默契的考虑,同一支特别行动队中队员的编号应该连续,即为形如 (i, ...

  8. 洛谷 P3628 [APIO2010]特别行动队

    题目描述 你有一支由 n 名预备役士兵组成的部队,士兵从 1 到 n 编号,要将他们拆分 成若干特别行动队调入战场.出于默契的考虑,同一支特别行动队中队员的编号 应该连续,即为形如的序列. 编号为 i ...

  9. APIO2010 特别行动队

    https://www.luogu.org/problem/show?pid=3628 特别行动队 可以说是斜率优化dp的模板题,50分的方程很容易写出,先推导函数然后用单调队列维护上凸壳. 代码: ...

最新文章

  1. Qt中ui文件的使用
  2. 图论 Make Unique:有向图和无向图的一些算法
  3. windows 执行bat脚本
  4. java xsd 解析 xml文件_Java针对XSD文件验证XML文件的最佳方法是什么?
  5. leetcode-sort-colors
  6. char* 赋值 const char* 释放_Effective C++读书笔记之条款3:尽可能使用const
  7. MailBee.NET Objects接收电子邮件(POP3)教程一:接收简单的电子邮件
  8. VS Code Java 开发指南!
  9. bes配置oracle数据源,安装 ZFS 根文件系统(Oracle Solaris 初始安装)
  10. 调研报告|在线语音识别改进之 RNN-T 训练
  11. 二元函数偏导数公式_高等数学六:(3)复合函数与微分性质
  12. WinNT无盘Win95维护经验
  13. 操作系统-------------------内存空间的分配方式(连续分配和非连续分配和虚拟存储技术)
  14. 在夕阳再晨的日子里(二)----掌管市场部的岁月之合作团队与社区的拓展
  15. DTOJ3704 威士忌(whiskey)
  16. 【PR 基础】轨道遮罩键、交叉溶解的简单使用
  17. 使用余弦相似度算法计算文本相似度
  18. hdfs文件系统和linux文件系统,分布式文件系统HADOOP HDFS与传统文件系统LINUX FS的比较与分析...
  19. 微信ibeacon使用场景步骤总结
  20. 银行加速布局大数据蓝海 激战线上供应链金融

热门文章

  1. Docker中定制镜像的两种方法
  2. 高德地图 SDK 的应用 01:绘制多边形区域图
  3. Linux系统使用--Ubuntu 16.04 安装为知笔记
  4. 教你如何求一个集合的所有非空子集的方差和
  5. 我的世界制作服务器地图,我的世界1.7地图制作必看心得 教你怎么做地图
  6. 水箱建模最小二乘法_大气VOCs在线监测系统评估工作指南(二)
  7. (1)从1开始写一个操作系统
  8. Android 设置黑白滤镜
  9. Linux查看DNS地址以及端口,nslookup命令以及更强大的dig命令
  10. 案例分享:Qt高频fpga采集数据压力位移速度加速度分析系统(通道配置、电压转换、采样频率、通道补偿、定时采集、距离采集,导出excel、自动XY轴、隐藏XY轴、隐藏显示通道,文件回放等等)