【CF765F】Souvenirs
【CF765F】Souvenirs
题面
题解
我们可以发现,对于某个右端点\(i\),左端点\(j\)在由\(i\rightarrow 1\)的过程中,每一段的答案是单调不增的,由这个性质,我们想办法维护出加入右端点\(i\)后的答案。
我们只考虑形如\(j<i,a_i<a_j\)的答案,因为其他的情况我们只需要将值域翻转即可得到答案。
根据我们上面的说法,你首先需要找到一个\(j\),使得\(a_j>a_i\)且\(j\)最大,这个可以用值域线段树维护。然后我们又要找到一个一个\(j'\),使得\(a_{j'}-a_i<a_j-a_i\),那么我们可以令\(j=j'\),继续更新答案。
如果仅仅是这样的话,我们的复杂度还是不满足要求的,但是限制\(a_{j'}-a_i<a_j-a_i\)可以转化为\(a_{j'}-a_i<a_j-a_{j'}\),因为大于的情况之前肯定以及更新过了,那么我们没必要再次更新,而这个限制又可以变为\(a_{j'}-a_i<\frac 12(a_j-a_i)\),这样子每次值域减半的话复杂度可以保证。
而维护某个左端点的话可以用树状数组维护后缀取\(\min\)即可。
复杂度\(O(n\log ^2 10^9)\)。
代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;
inline int gi() {
register int data = 0, w = 1;
register char ch = 0;
while (!isdigit(ch) && ch != '-') ch = getchar();
if (ch == '-') w = -1, ch = getchar();
while (isdigit(ch)) data = 10 * data + ch - '0', ch = getchar();
return w * data;
}
const int INF = 1e9;
const int MAX_N = 3e5 + 5;
vector<pair<int, int> > vec[MAX_N];
int N, M, a[MAX_N], ans[MAX_N];
int c[MAX_N];
inline int lb(int x) { return x & -x; }
void Add(int x, int v) { while (x) c[x] = min(c[x], v), x -= lb(x); }
int Min(int x) { int res = INF; while (x <= N) res = min(res, c[x]), x += lb(x); return res; }
struct Node { int ls, rs, v; } t[MAX_N << 5];
int rt[MAX_N], tot;
void modify(int &o, int p, int l, int r, int pos, int v) {
o = ++tot, t[o] = t[p], t[o].v = max(t[o].v, v);
if (l == r) return ;
int mid = (l + r) >> 1;
if (pos <= mid) modify(t[o].ls, t[p].ls, l, mid, pos, v);
else modify(t[o].rs, t[p].rs, mid + 1, r, pos, v);
}
int query(int o, int l, int r, int ql, int qr) {
if (!o || l > r) return 0;
if (ql <= l && r <= qr) return t[o].v;
int mid = (l + r) >> 1, res = 0;
if (ql <= mid) res = max(res, query(t[o].ls, l, mid, ql, qr));
if (qr > mid) res = max(res, query(t[o].rs, mid + 1, r, ql, qr));
return res;
}
void solve() {
for (int i = 1; i <= N; i++) c[i] = INF;
for (int i = 1; i <= tot; i++) t[i] = (Node){0, 0, 0};
tot = 0, memset(rt, 0, sizeof(rt));
for (int i = 1; i <= N; i++) modify(rt[i], rt[i - 1], 0, INF, a[i], i);
for (int i = 1; i <= N; i++) {
int j = query(rt[i - 1], 0, INF, a[i], INF);
while (j) {
Add(j, a[j] - a[i]);
j = query(rt[j - 1], 0, INF, a[i], ((a[i] + a[j]) >> 1) - (~(a[i] + a[j]) & 1));
}
for (auto it : vec[i]) ans[it.second] = min(ans[it.second], Min(it.first));
}
}
int main () {
#ifndef ONLINE_JUDGE
freopen("cpp.in", "r", stdin);
#endif
N = gi();
for (int i = 1; i <= N; i++) a[i] = gi();
M = gi();
for (int i = 1; i <= M; i++) {
int l = gi(), r = gi();
vec[r].push_back(make_pair(l, i));
ans[i] = INF;
}
solve();
for (int i = 1; i <= N; i++) a[i] = INF - a[i];
solve();
for (int i = 1; i <= M; i++) printf("%d\n", ans[i]);
return 0;
}
【CF765F】Souvenirs的更多相关文章
- 【CF765F】Souvenirs 主席树
[CF765F]Souvenirs 题意:给你一个长度为n的序列{ai},有m个询问,每次询问给出l,r,问在所有$l\le x < y\le r$中,$|a_x-a_y|$的最小值是多少. $ ...
- Python高手之路【六】python基础之字符串格式化
Python的字符串格式化有两种方式: 百分号方式.format方式 百分号的方式相对来说比较老,而format方式则是比较先进的方式,企图替换古老的方式,目前两者并存.[PEP-3101] This ...
- 【原】谈谈对Objective-C中代理模式的误解
[原]谈谈对Objective-C中代理模式的误解 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 这篇文章主要是对代理模式和委托模式进行了对比,个人认为Objective ...
- 【原】FMDB源码阅读(三)
[原]FMDB源码阅读(三) 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 FMDB比较优秀的地方就在于对多线程的处理.所以这一篇主要是研究FMDB的多线程处理的实现.而 ...
- 【原】Android热更新开源项目Tinker源码解析系列之一:Dex热更新
[原]Android热更新开源项目Tinker源码解析系列之一:Dex热更新 Tinker是微信的第一个开源项目,主要用于安卓应用bug的热修复和功能的迭代. Tinker github地址:http ...
- 【调侃】IOC前世今生
前些天,参与了公司内部小组的一次技术交流,主要是针对<IOC与AOP>,本着学而时习之的态度及积极分享的精神,我就结合一个小故事来初浅地剖析一下我眼中的“IOC前世今生”,以方便初学者能更 ...
- Python高手之路【三】python基础之函数
基本数据类型补充: set 是一个无序且不重复的元素集合 class set(object): """ set() -> new empty set object ...
- Python高手之路【一】初识python
Python简介 1:Python的创始人 Python (英国发音:/ˈpaɪθən/ 美国发音:/ˈpaɪθɑːn/), 是一种解释型.面向对象.动态数据类型的高级程序设计语言,由荷兰人Guido ...
- 【开源】简单4步搞定QQ登录,无需什么代码功底【无语言界限】
说17号发超简单的教程就17号,qq核审通过后就封装了这个,现在放出来~~ 这个是我封装的一个开源项目:https://github.com/dunitian/LoTQQLogin ————————— ...
随机推荐
- Entity Framework Core今日所得:避免 IEnumerable 以及 IQueryable 陷阱
避免 IEnumerable 以及 IQueryable 陷阱: IEnumerable示用Linq会先去数据库查询所有记录,然后再条件查询. IQueryable接口派生自IEnumerable,但 ...
- [大数据学习研究] 错误排查,Hadoop集群部分DataNode不能启动
错误现象 不知道什么原因,今天发现我的hadoop集群启动后datanode只有一台了,我的集群本来有三台的,怎么只剩一台了呢? 用jps命令检查一下,发现果然有两台机器的DataNode没有启动. ...
- 8 Traits of an Experienced Programmer that every beginner programmer should know
Referrence: http://whats-online.info/guides-and-info/36/Traits-of-Experienced-Programmer-that-every- ...
- windows 2003 域控服务器导出全部hash的方法
天下文章一大抄,我也是醉了... 一份“错误”的文章一遍又一遍的被转载,盲目转载,根本不细看.只会误导新手. 谈下windows2003域控下如何导出全部的hash信息. 1. 使用备份还原向导 2. ...
- 单词demantoite翠榴石demantoite英语
一般认为翠榴石demantoite的形成条件是: (1)围岩组成应该是贫铝富铁,且附近有钙质碳酸盐出露地区,即有利于“纯度高”的钙铁榴石结晶环境.否则,若钙铝榴石端员分子比增多,2价Fe必将直接影响晶 ...
- Android-----解析xml文件的三种方式
SAX解析方法介绍: SAX(Simple API for XML)是一个解析速度快并且占用内存少的XML解析器,非常适合用于Android等移动设备.SAX解析XML文件采用的是事件驱动,也就是说, ...
- 以服务方式启动tomcat无法访问NFS共享盘
用startup.bat方式启动tomcat,程序的可以访问NFS共享盘的文件.但用 1).以服务的方式启动tomcat 2).或者用windows的任务计划去执行startup.bat的方式启动to ...
- php workerman定时任务
一.下载workerman https://www.workerman.net/download 二.下载workerman/mysql http://doc3.workerman.net/64020 ...
- 转:win10完美去除快捷方式小箭头
1.去掉小箭头 reg add "HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Windows\CurrentVersion\Explorer\Shell Ic ...
- Django bootstrap按钮点击后激活active
现在有个需求,就是在导航栏上有若干个按钮,我想实现的功能是当点击某个按钮后修改文字颜色,这样网站会更人性化.现总结方法如下: index.html: <ul class="navbar ...