https://www.luogu.org/problemnew/show/P3648

同bzoj3675

这题斜率优化+滚动数组就可以了qwq

因为我是在飞机上瞎bb的式子,所以可能会和别的题解的式子不一样(反正 A 题了对吧)

推的过程什么的都写在最底下的注释里了,大家可以选择先查看最底下的过程(您是神犇可以不用看这篇博客)

但这篇题解不知为什么在bzoj上过不了,请dalao们赐教

#include <bits/stdc++.h>
using namespace std; template <typename T>
inline void read(T &f) {
f = 0; T fu = 1; char c = getchar();
while (c < '0' || c > '9') {if (c == '-') fu = -1; c = getchar();}
while (c >= '0' && c <= '9') {f = (f << 3) + (f << 1) + (c & 15); c = getchar();}
f *= fu;
} typedef long long ll;
const int N = 100000 + 10;
const int M = 200 + 10; ll f[2][N], s[N], a[N];
int fr[M][N], q[N], ans[N], head, tail;
int n, m; bool pd(int a, int b, int c, int d) {
// g(b, c) >= s[n] - s[d]
// (f[b] - f[c]) / (s[b] - s[c]) >= s[n] - s[d]
if(f[a][b] - f[a][c] <= (s[b] - s[c]) * (s[n] - s[d])) return 1;
return 0;
} bool pd2(int a, int b, int c, int d) {
// g(b, c) < g(c, d)
// (f[b] - f[c]) / (s[b] - s[c]) < (f[c] - f[d]) / (s[c] - s[d])
// (f[b] - f[c]) * (s[c] - s[d]) < (f[c] - f[d]) * (s[b] - s[c])
if((f[a][b] - f[a][c]) * (s[c] - s[d]) < (f[a][c] - f[a][d]) * (s[b] - s[c])) return 1;
return 0;
} int main() {
read(n); read(m);
for(int i = 1; i <= n; i++) {
read(a[i]);
s[i] = s[i - 1] + a[i];
}
for(int i = 1; i <= m; i++) {
head = 0, tail = 0; q[0] = i - 1;
for(int j = i; j <= n; j++) {
// 保证队列里有两个元素
while(head + 1 <= tail) {
if(pd((i & 1) ^ 1, q[head], q[head + 1], j)) head++;
else break;
}
f[i & 1][j] = f[(i & 1) ^ 1][q[head]] + (s[n] - s[j]) * (s[j] - s[q[head]]);
fr[i][j] = q[head];
while(head + 1 <= tail) {
if(pd2((i & 1) ^ 1, q[tail - 1], q[tail], j)) tail--;
else break;
}
q[++tail] = j;
}
}
ll maxn = -1; int num;
for(int i = n; i >= m; i--) if(f[m & 1][i] > maxn) maxn = f[m & 1][i], num = i;
printf("%lld\n", maxn);
for(int i = m; i >= 1; i--) {
ans[i] = num;
num = fr[i][num];
}
for(int i = 1; i < m; i++) printf("%d ", ans[i]);
printf("%d\n", ans[m]);
return 0;
} // f[i] 表示取到了第 i 个数获得的最大价值
// f[i] = f[j] + (s[n] - s[i]) * (s[i] - s[j])
// f[i] = f[j] + s[n] * s[i] - s[n] * s[j] - s[i] * s[i] + s[i] * s[j]
// 假设 j 比 k 优, j < k
// f[j] + s[n] * s[i] - s[n] * s[j] - s[i] * s[i] + s[i] * s[j] > f[k] + s[n] * s[i] - s[n] * s[k] - s[i] * s[i] + s[i] * s[k]
// f[j] - s[n] * s[j] + s[i] * s[j] > f[k] - s[n] * s[k] + s[i] * s[k]
// f[j] - s[j] * (s[n] - s[i]) > f[k] - s[k] * (s[n] - s[i])
// f[j] - f[k] > (s[j] - s[k]) * (s[n] - s[i])
// (f[j] - f[k]) / (s[j] - s[k]) < s[n] - s[i]
// 令 g(i, j) 表示 i 和 j 的斜率
// i < j < k
// if g(i, j) < g(j, k)
// if g(i, j) < s[n] - s[x] j 比 i 差
// if g(i, j) >= s[n] - s[x] i 比 j 差, j 比 k 差
// 所以 j 就被扔掉了 QAQ
// 令 a, b 为队首两点
// if g(a, b) >= s[n] - s[i]
// a 比 b 差
// 当 i 上升时, s[n] - s[i] 单调不增
// a 永远比 b 差
// 飞机上没有博客只能自己推式子好痛苦啊 QAQ

luoguP3648 [APIO2014]序列分割的更多相关文章

  1. 【斜率DP】BZOJ 3675:[Apio2014]序列分割

    3675: [Apio2014]序列分割 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 1066  Solved: 427[Submit][Statu ...

  2. BZOJ 3675: [Apio2014]序列分割( dp + 斜率优化 )

    WA了一版... 切点确定的话, 顺序是不会影响结果的..所以可以dp dp(i, k) = max(dp(j, k-1) + (sumn - sumi) * (sumi - sumj)) 然后斜率优 ...

  3. bzoj3675[Apio2014]序列分割 斜率优化dp

    3675: [Apio2014]序列分割 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 3508  Solved: 1402[Submit][Stat ...

  4. BZOJ_3675_[Apio2014]序列分割_斜率优化

    BZOJ_3675_[Apio2014]序列分割_斜率优化 Description 小H最近迷上了一个分隔序列的游戏.在这个游戏里,小H需要将一个长度为n的非负整数序列分割成k+1个非空的子序列.为了 ...

  5. 斜率优化入门学习+总结 Apio2011特别行动队&Apio2014序列分割&HZOI2008玩具装箱&ZJOI2007仓库建设&小P的牧场&防御准备&Sdoi2016征途

    斜率优化: 额...这是篇7个题的题解... 首先说说斜率优化是个啥,额... f[i]=min(f[j]+xxxx(i,j)) ;   1<=j<i (O(n^2)暴力)这样一个式子,首 ...

  6. P3648 [APIO2014]序列分割(斜率优化dp)

    P3648 [APIO2014]序列分割 我们先证明,分块的顺序对结果没有影响. 我们有一个长度为3的序列$abc$ 现在我们将$a,b,c$分开来 随意枚举一种分块方法,如$(ab)(c)$,$(a ...

  7. [luogu P3648] [APIO2014]序列分割

    [luogu P3648] [APIO2014]序列分割 题目描述 小H最近迷上了一个分隔序列的游戏.在这个游戏里,小H需要将一个长度为n的非负整数序列分割成k+1个非空的子序列.为了得到k+1个子序 ...

  8. 洛谷 P3648 [APIO2014]序列分割 解题报告

    P3648 [APIO2014]序列分割 题目描述 你正在玩一个关于长度为\(n\)的非负整数序列的游戏.这个游戏中你需要把序列分成\(k+1\)个非空的块.为了得到\(k+1\)块,你需要重复下面的 ...

  9. [APIO2014]序列分割 --- 斜率优化DP

    [APIO2014]序列分割 题目大意: 你正在玩一个关于长度为\(n\)的非负整数序列的游戏.这个游戏中你需要把序列分成\(k+1\)个非空的块.为了得到\(k+1\)块,你需要重复下面的操作\(k ...

随机推荐

  1. pandas数据对齐

    Pandas的对齐运算 是数据清洗的重要过程,可以按索引对齐进行运算,如果没对齐的位置则补NaN,最后也可以填充NaN Series的对齐运算 1. Series 按行.索引对齐 示例代码: s1 = ...

  2. Python中常见的异常处理

    异常和错误 part1:程序中难免出现错误,而错误分成两种 1. 语法错误(这种错误,根本过不了Python解释器的语法检测,必须在程序执行前就改正) # 语法错误示范一 if # 语法错误示范二 d ...

  3. Java面向对象-代码块

    Java面向对象-代码块 代码块主要就是通过{}花括号 括起来的代码: 主要分为 普通代码块 构造块 静态代码块三类.后面学到线程还有一个同步代码块,到时候再说: 普通代码块:仅仅是花括号括起来的代码 ...

  4. 摄影之HDR

    摄影之HDR 高动态范围图像(High-Dynamic Range,简称HDR),相比普通的图像,可以提供更多的动态范围和图像细节,根据不同的曝光时间的LDR(Low-Dynamic Range)图像 ...

  5. 浏览器get请求到java后台的值是乱码

     get方式提交的参数编码,只支持iso8859-1编码. 因此,如果里面有中文,在后台就需要转换编码,如下 String zhongwen = request.getParameter(" ...

  6. 【bzoj1019】[SHOI2008]汉诺塔

    1019: [SHOI2008]汉诺塔 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 1427  Solved: 872[Submit][Status] ...

  7. 前端性能分析:分析百度和sogou

    先用httpwatch录制这两个网站:www.baidu.com  www.sogou.com 由上图可以看到: 百度用时0.278s 发送7831B 接收36620B 13个请求 搜狗       ...

  8. Mask_RCNN caffe2

    cd DETECTRON/detectron/ python2 tools/infer_simple_ip_camera.py \    --cfg configs/12_2017_baselines ...

  9. IDEA maven打包时跳过测试

    配置这个install -Dmaven.test.skip=true, 可以跳过business项目本地启动自动跑测试用例

  10. hadoop 2.7.3 (hadoop2.x)使用ant制作eclipse插件hadoop-eclipse-plugin-2.7.3.jar

    为了做mapreduce开发,要使用eclipse,并且需要对应的Hadoop插件hadoop-eclipse-plugin-2.7.3.jar,首先说明一下,在hadoop1.x之前官方hadoop ...