博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU 2993 MAX Average Problem
阅读量:5316 次
发布时间:2019-06-14

本文共 1630 字,大约阅读时间需要 5 分钟。

题目大意:

给定一个长度为 \(N\) 的序列,求一段长度至少为 \(K\) 的区间,使得这段这段区间的平均值最大。

思路:

首先有 \(ans = max(\frac{sum[i] - sum[j]} {i-j}) (j < i)\)\(O(n^2)\) 暴力。

观察式子,答案就是 \((i, sum[i])\)\((j, sum[j])\) 之间的斜率。

所以,题目求的就是平面上 \(N\) 个点之间最大的斜率。

其次,假设现在有3个点 \(i, j, k\)\((k<j<i)\)

ZszR1O.png

若此时求的点在红色区域中,选 \(k\) 比选 \(j\) 优。

若此时求的点在绿色区域中,选 \(i\) 比选 \(j\) 优。

所以 \(j\) 是无用的点,删除。

故维护的是一个下凸。

若此时求的点与一对点对 \((i,j)\) ,\((i<j)\) 形成上凸,选i更优,

反之形成下凸,选 \(j\) 更优。

查询时是选取两点间的斜率判断,不具有单调性,所以要用单调栈 + 二分查询

代码:

#include 
using namespace std;#define REP(i, a, n) for (int i = a, _n = n; i <= _n; ++i)#define DREP(i, a, n) for (int i = a, _n = n; i >= _n; --i)#define FOR(i, a, n) for (int i = a, _n = n; i < _n; ++i)#define EREP(i, a) for (int i = first[a]; i; i = edge[i].nxt)#define debug(x) cout << #x << " = " << x << endlint n, k;int Q[500005], tail;long long S[500005];int main () { scanf ("%d%d", &n, &k); double ans = -2e9; REP (i, 1, n) scanf ("%lld", S + i), S[i] += S[i - 1]; REP (i, k, n) { int j = i - k; while (tail && (S[Q[tail]] - S[Q[tail - 1]]) * (j - Q[tail]) >= (S[j] - S[Q[tail]]) * (Q[tail] - Q[tail - 1])) --tail; Q[++tail] = j; int l = 1, r = tail - 1, res = tail; while (l <= r) { int mid = l + r >> 1; if ((S[i] - S[Q[mid + 1]]) * (Q[mid + 1] - Q[mid]) <= (S[Q[mid + 1]] - S[Q[mid]]) * (i - Q[mid + 1])) // 是上凸 res = mid, r = mid - 1; else l = mid + 1; } ans = max(ans, 1.0 * (S[i] - S[Q[res]]) / (i - Q[res])); } printf ("%.2f\n", ans); return 0;}

转载于:https://www.cnblogs.com/ZPAYAUR/p/11155778.html

你可能感兴趣的文章
freebsd 实现 tab 命令 补全 命令 提示
查看>>
numpy调试
查看>>
struts1和struts2的区别
查看>>
函数之匿名函数
查看>>
shell习题第16题:查用户
查看>>
梯度下降法与方向导数
查看>>
实验4 [bx]和loop的使用
查看>>
Redis常用命令
查看>>
2018.08.22 NOIP模拟 shop(lower_bound+前缀和预处理)
查看>>
2018.11.06 bzoj1040: [ZJOI2008]骑士(树形dp)
查看>>
2019.02.15 bzoj5210: 最大连通子块和(链分治+ddp)
查看>>
redis cluster 集群资料
查看>>
混合模式程序集是针对“v2.0.50727”版的运行时生成的,在没有配置其他信息的情况下,无法在 4.0 运行时中加载该程序集。...
查看>>
jQuery总结或者锋利的jQuery笔记二
查看>>
微软职位内部推荐-Sr. SE - Office incubation
查看>>
微软职位内部推荐-SOFTWARE ENGINEER II
查看>>
centos系统python2.7更新到3.5
查看>>
【Quartz】常用方法的使用方式(三)
查看>>
MVVM模式下关闭窗口的实现
查看>>
C#区域截图——调用API截图
查看>>