GMOJ5673 爬山法 题解
Solution
显然先想到处理出每个点能看到的最高的顶点。
然后考虑模拟题目的过程,一段一段走时间复杂度显然不够优秀。
考虑我们要求什么,我们需要求出\(u\)到\(v\)的最近的一个点,使得这个点能看到的点比\(v\)能看到的点更高。
然后这个东西可以直接线段树,当然也可以二分+st表
复杂度\(O(n\log n)\)
Code
#include <cstdio>
#include <iostream>
#define LL long long
#define RE register
#define IN inline
using namespace std;
IN int read() {
int res = 0; char ch = getchar();
for(; !isdigit(ch); ch = getchar());
for(; isdigit(ch); ch = getchar()) res = (res << 1) + (res << 3) + (ch ^ 48);
return res;
}
int n, le[2000010], ri[2000010], nxt[2000010], st[2000010][20], lg[2000010];
LL f[2000010];
struct Point {
int x, y;
}p[2000010];
inline double slope(int x, int y) {return 1.0 * (p[x].y - p[y].y) / (p[x].x - p[y].x);}
inline bool cmp(int x, int y) {
if(x == y) return false;
if(p[x].y > p[y].y) return true;
if(p[x].y == p[y].y) return p[x].x > p[y].x;
return false;
}
int stk[2000010];
inline int pmax(int x, int y) {return cmp(x, y) ? x : y;}
void preWork() {
int top = 0;
for(int i = 1; i <= n; ++i) {
while(top > 1 && slope(stk[top], i) >= slope(stk[top], stk[top - 1])) -- top;
if(top) le[i] = stk[top]; stk[++top] = i;
}
top = 0;
for(int i = n; i; --i) {
while(top > 1 && slope(stk[top], i) <= slope(stk[top], stk[top - 1])) -- top;
if(top) ri[i] = stk[top]; stk[++top] = i;
}
for(int i = 1; i <= n; ++i) {
nxt[i] = pmax(le[i], ri[i]);
if(cmp(i, nxt[i])) nxt[i] = 0;
st[i][0] = nxt[i];
}
for(int i = 1; i <= lg[n]; ++i)
for(int j = 1; j <= n; ++j)
st[j][i] = pmax(st[j][i - 1], st[j + (1 << i - 1)][i - 1]);
return ;
}
inline LL abs(int x) {return x < 0 ? -x : x;}
inline int query(int l, int r) {int len = r - l + 1; return pmax(st[l][lg[len]], st[r - (1 << lg[len]) + 1][lg[len]]);}
LL solve(int k) {
if(f[k] || !nxt[k]) return f[k];
int res = nxt[k];
if(nxt[k] < k) {
for(int l = nxt[k], r = k - 1, mid; l <= r;) {
mid = l + r >> 1;
if(cmp(query(mid, k), nxt[k])) res = mid, l = mid + 1;
else r = mid - 1;
}
}
else {
for(int l = k + 1, r = nxt[k], mid; l <= r;) {
mid = l + r >> 1;
if(cmp(query(k, mid), nxt[k])) res = mid, r = mid - 1;
else l = mid + 1;
}
}
f[k] = solve(res) + abs(res - k);
return f[k];
}
void out(int x) {printf("%d %d\n",p[x].x,p[x].y);}
int main() {
freopen("mountain.in","r",stdin);
freopen("mountain.out","w",stdout);
n = read();
for(int i = 2; i <= n; ++i) lg[i] = lg[i >> 1] + 1;
for(int i = 1; i <= n; ++i) p[i].x = read(), p[i].y = read();
preWork();
for(int i = 1; i <= n; ++i) solve(i);
for(int i = 1; i <= n; ++i) printf("%lld\n",f[i]);
return 0;
}```
GMOJ5673 爬山法 题解的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
随机推荐
- typora的第一天
一级标题 二级标题 三级标题 ..... 表格 java spring mybatis 代码 java代码 public void Hello(){ } 字体 hello word! hello wo ...
- React报错之React Hook 'useEffect' is called in function
正文从这开始~ 总览 为了解决错误"React Hook 'useEffect' is called in function that is neither a React function ...
- flutter系列之:构建Widget的上下文环境BuildContext详解
目录 简介 BuildContext的本质 BuildContext和InheritedWidget BuildContext的层级关系 总结 简介 我们知道Flutter中有两种Widget,分别是 ...
- [HDU3976]Electric resistance(电阻)(信竞&物竞)(高斯消元)
题面 Problem Description Now give you a circuit who has n nodes (marked from 1 to n) , please tell abc ...
- 「题解报告」P7301 【[USACO21JAN] Spaced Out S】
原题传送门 神奇的5分算法:直接输出样例. 20分算法 直接把每个点是否有牛的状态DFS一遍同时判断是否合法,时间复杂度约为\(O(2^{n^2})\)(因为有判断合法的剪枝所以会比这个低).而在前四 ...
- Docker 安装 MySQL、Redis
1 Docker 中安装 Redis 1.1 创建目录 在硬盘上创建 redis 的数据目录: mkdir -p /Users/yygnb/dockerMe/redis/data 为该目录添加权限: ...
- 第六十五篇:Vue的过滤器
好家伙, 过滤器,vue3取消了,只有vue2能用 1.过滤器 过滤器(Filters)是vue为开发者提供的功能,常用于文本的格式化. 过滤器可以用在两个地方:插值表达式和v-bind属性绑定. 过 ...
- ubuntu下安装python
一.安装python3.6 sudo add-apt-repository ppa:jonathonf/python-3.6 如显示不能添加"'ppa:~jonathonf/ubuntu/p ...
- flex常用布局
公共样式: <style> * { margin: 0; padding: 0; } .has-flex { display: flex; } </style> 垂直居中 子元 ...
- 硬核解析MySQL的MVCC实现原理,面试官看了都直呼内行
1. 什么是MVCC MVCC全称是Multi-Version Concurrency Control(多版本并发控制),是一种并发控制的方法,通过维护一个数据的多个版本,减少读写操作的冲突. 如果没 ...