题目描述 Description

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

编号为i的士兵的初始战斗力为xi ,一支特别行动队的初始战斗力x为队内士兵初始战斗力之和,即x = xi+ xi+1+ … + xi+k。 
通过长期的观察,你总结出一支特别行动队的初始战斗力x将按如下经验公式修正为x':x' = ax2 + 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,没有其它方案能使修正后的战斗力和更大。

输入描述 Input Description

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

输出描述 Output Description

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

样例输入 Sample Input


-1 10 -20 
2 2 3 4

样例输出 Sample Output

9

数据范围及提示 Data Size & Hint

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。

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1000000+5;
int read()
{int x=0,f=1;char ch;ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') f=-1; ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}return x*f;
}
int n,a,b,c;
int w[maxn],q[maxn];
long long s[maxn],f[maxn];
long long g(int i)
{return a*s[i]*s[i]+b*s[i];
}
long long h(int j)
{return f[j]+a*s[j]*s[j]-b*s[j];
}
double slope(int j1,int j2)
{return (h(j1)-h(j2))/(2.0*a*(s[j1]-s[j2]));
}
long long js(int x)
{return a*x*x+b*x+c;
}
int main()
{int i;n=read();a=read();b=read();c=read();for(i=1;i<=n;i++)w[i]=read();for(i=1;i<=n;i++)s[i]=s[i-1]+w[i];int h=1,p=0;q[++p]=0;for(i=1;i<=n;i++){while(p>h&&slope(q[h],q[h+1])<=s[i]) h++;int t=q[h];f[i]=f[t]+js(s[i]-s[t]);while(p>h&&slope(q[p],i)<slope(q[p-1],q[p])) p--;q[++p]=i;}printf("%lld",f[n]);return 0;
}

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

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

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

  2. P3628 [APIO2010]特别行动队(简单斜率优化)

    LINK 题意 nnn名士兵,你需要把士兵分成连续的一段一段特别行动队 设其中一队的战斗力和为xxx,那么修正战斗力是ax2+bx+cax^2+bx+cax2+bx+c(当然a<0a<0a ...

  3. P3628 [APIO2010]特别行动队(斜率优化DP)

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

  4. 洛谷P3628 [APIO2010]特别行动队(斜率优化)

    传送门 先写出转移方程$$dp[i]=max\{dp[j]+a*(sum[i]-sum[j])^2+b*(sum[i]-sum[j])+c\}$$ 假设$j$比$k$更优,则有$$dp[j]+a*(s ...

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

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

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

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

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

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

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

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

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

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

最新文章

  1. hadoop启动_Mac OS X 上搭建 Hadoop 开发环境指南
  2. 异步IO一定更好吗?
  3. AOP代理配置pointcut
  4. linux yun 安装 redis,Linux系统安装NoSQL(MongoDB和Redis)步骤及问题解决办法(总结篇)...
  5. 实现用户名密码登陆,输入错误3次后锁定用户
  6. 【腾许Bugly干货分享】“HTTPS”安全在哪里?
  7. 如何找到需要的rpm包
  8. os.environ['TF_CPP_MIN_LOG_LEVEL'] = '2' 无效
  9. CDA level 1 业务分析师 推荐学习书目
  10. Matlab中median函数的使用
  11. virtualbox下安装雪豹10.6系统(AMD下亲测成功)
  12. 记一次很坑很坑的报错java.lang.Exception: The class is not public.
  13. 【cvpr2020】Referring Image Segmentation via Cross-Modal Progressive Comprehension
  14. ajax前后端通信的头部消息之请求头与响应头
  15. python 读取excel表格数据,然后写入word模板,批量生成
  16. **传统线上支付 区块链**
  17. 计算机图形学常用算法实现11 扫描线z-buffer算法
  18. 安卓android按宽/高缩放图片到指定大小并进行裁剪得到中间部分图片
  19. 狂欢,不过是一群人的孤单--来自人人
  20. 六款窗口比较器电路设计原理

热门文章

  1. 移动支付商如何赚你的钱
  2. 如何快速批量导入非Oracle DB格式的数据--sqlloader
  3. 工业大数据分析算法及工具
  4. 统一用户认证和单点登录(SSO)解决方案
  5. BZOJ 1064 假面舞会(图论-连通分量)
  6. 要不起对三队 - alpha冲刺(第三天)
  7. Linux查看磁盘和系统资源占用
  8. M1卡/CUID卡格式化清卡操作方法
  9. WebService全面详解
  10. 切换tab页,页面局部刷新,地址栏路径修改