BZOJ2716 [Violet 3]天使玩偶 【CDQ分治】
题目
输入格式
输出格式
输入样例
//样例太长就不贴了。。。。
输出样例
//见原题
提示
题解
我们将曼哈顿距离式子中的绝对值去掉,每次只考虑x,y比当前点小的更新答案。
为了使所有点都对答案进行更新,将坐标轴旋转三次再算三次
每一次对于点(x,y),找到(x’,y’)【x’<=x,y’<= y且 时间t’ < t】使得x+y−(x′+y′)最小
类似三维偏序的东西,可以用CDQ分治
树状数组维护最大值
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define lbt(x) (x & -x)
using namespace std;
const int maxn = 1000005,maxm = 2010005,INF = 0x7fffffff;
inline int RD(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57) {out = (out << 1) + (out << 3) + c - '0'; c = getchar();}
return out * flag;
}
int ans[maxn],n,m,N = 0,ansi = 0,Qi = 0,S[maxm];
struct Que{int x,y,t,id;}Q[maxn],T[maxn];
inline bool operator <(const Que& a,const Que& b){
if (a.x == b.x && a.y == b.y) return a.t < b.t;
if (a.x == b.x) return a.y < b.y;
return a.x < b.x;
}
void add(int u,int v){while (u <= N) S[u] = max(S[u],v),u += lbt(u);}
int query(int u){int an = 0; while (u) an = max(an,S[u]),u -= lbt(u); return an;}
void cls(int u){while (u <= N) S[u] = 0,u += lbt(u);}
void CDQ(int l,int r){
if (l == r) return;
int mid = l + r >> 1,l1 = l,l2 = mid + 1,temp;
for (int i = l; i <= r; i++){
if (Q[i].t <= mid && !Q[i].id)
add(Q[i].y,Q[i].x + Q[i].y);
else if (Q[i].t > mid && Q[i].id){
temp = query(Q[i].y);
if (temp) ans[Q[i].id] = min(ans[Q[i].id],Q[i].x + Q[i].y - temp);
}
}
for (int i = l; i <= r; i++){
if (Q[i].t <= mid){
T[l1++] = Q[i];
if (!Q[i].id) cls(Q[i].y);
}else T[l2++] = Q[i];
}
for (int i = l; i <= r; i++) Q[i] = T[i];
CDQ(l,mid); CDQ(mid + 1,r);
}
int main(){
//freopen("in.txt","r",stdin);
//freopen("out1.txt","w",stdout);
n = RD();m = RD(); Qi = n + m; int opt;
REP(i,n){
Q[i].x = RD() + 1,Q[i].y = RD() + 1,Q[i].t = i,Q[i].id = 0;
N = max(N,Q[i].x); N = max(N,Q[i].y);
}
for (int i = n + 1; i <= Qi; i++){
opt = RD(); Q[i].x = RD() + 1; Q[i].y = RD() + 1,Q[i].t = i;
N = max(N,Q[i].x); N = max(N,Q[i].y);
if (opt & 1) Q[i].id = 0;
else Q[i].id = ++ansi;
}N++;
fill(ans,ans + maxn,INF);
sort(Q + 1,Q + 1 + Qi); CDQ(1,Qi);
REP(i,Qi) Q[i].x = N - Q[i].x;
sort(Q + 1,Q + 1 + Qi); CDQ(1,Qi);
REP(i,Qi) Q[i].y = N - Q[i].y;
sort(Q + 1,Q + 1 + Qi); CDQ(1,Qi);
REP(i,Qi) Q[i].x = N - Q[i].x;
sort(Q + 1,Q + 1 + Qi); CDQ(1,Qi);
REP(i,ansi) printf("%d\n",ans[i]);
return 0;
}
BZOJ2716 [Violet 3]天使玩偶 【CDQ分治】的更多相关文章
- [BZOJ2716] [Violet 3]天使玩偶(CDQ分治)
[BZOJ2716] [Violet 3]天使玩偶(CDQ分治) 题面 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里, ...
- BZOJ 2716: [Violet 3]天使玩偶( CDQ分治 + 树状数组 )
先cdq分治, 然后要处理点对答案的贡献, 可以以询问点为中心分成4个区域, 然后去掉绝对值(4种情况讨论), 用BIT维护就行了. --------------------------------- ...
- BZOJ 2716 Violet 3 天使玩偶 CDQ分治
题目大意:初始给定平面上的一个点集.提供两种操作: 1.将一个点增加点集 2.查询距离一个点最小的曼哈顿距离 K-D树是啥...不会写... 我仅仅会CDQ分治 对于一个询问,查询的点与这个点的位置关 ...
- BZOJ 2716 [Violet 3]天使玩偶 (CDQ分治、树状数组)
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=2716 怎么KD树跑得都那么快啊..我写的CDQ分治被暴虐 做四遍CDQ分治,每次求一个 ...
- bzoj2648SJY摆棋子&&bzoj2716[Violet 3]天使玩偶*
bzoj2648SJY摆棋子 bzoj2716[Violet 3]天使玩偶 题意: 棋盘上有n个棋子,现在有m个操作,一种是加棋子,一种是查询离某个点最近的棋子.n,m≤500000. 题解: 先将已 ...
- [bzoj] 2716 天使玩偶 || CDQ分治
原题 已知n个点有天使玩偶,有m次操作: 操作1:想起来某个位置有一个天使玩偶 操作2:询问离当前点最近的天使玩偶的曼哈顿距离 显然的CDQ问题,三维分别为时间,x轴,y轴. 但是这道题的问题在于最近 ...
- bzoj2716: [Violet 3]天使玩偶
#include <iostream> #include <cstdio> #include <cstring> #include <cmath> #i ...
- 洛谷P4169 天使玩偶 CDQ分治
还是照着CDQ的思路来. 但是有一些改动: 要求4个方向的,但是可爱的CDQ分治只能求在自己一个角落方向上的.怎么办?旋转!做4次就好了. 统计的不是和,而是——max!理由如下: 设当前点是(x,y ...
- CH 4701 - 天使玩偶 - [CDQ分治]
题目链接:传送门 关于CDQ分治(参考李煜东<算法竞赛进阶指南>): 对于一系列操作,其中的任何一个询问操作,其结果必然等价于:初始值 + 此前所有的修改操作产生的影响. 假设共有 $m$ ...
随机推荐
- 谈谈对bug的一点想法,说说做好日志记录的重要性
说起程序猿,总绕不开的一个话题就是bug,估计每个程序猿听到某某测试跑过来一脸淫笑的告诉你这个功能有个bug的时候,总有种恨不得掐死他的想法.其实程序猿跟bug的关系,感觉有点像父亲和儿子的关系,自己 ...
- pycharm中文乱码问题 总结
前言: 这几天刚刚开始学习python,然后就安装了pycharm,但是那个中文乱码的问题真是让人心烦,在网上找了好久,都写得好乱,今天终于让我解决了,在这里总结一下经验,希望可以帮到你们 问题:如下 ...
- 笔记--tslib 编译
tslib 是qt启动时的一个触屏校正检验程序. 它的配置以及编译比较简单. 第一步, 下载tslib源码包: http://download.csdn.net/detail/MKNDG/329156 ...
- PAT (Basic Level) Practice 1004 成绩排名
个人练习 读入n名学生的姓名.学号.成绩,分别输出成绩最高和成绩最低学生的姓名和学号. 输入格式:每个测试输入包含1个测试用例,格式为\ 第1行:正整数n 第2行:第1个学生的姓名 学号 成绩 第3行 ...
- POJ:2431-Expedition
Expedition Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 20089 Accepted: 5786 Descripti ...
- 适配IE8+等浏览器的适配播放插件
function myBrowser(){ var userAgent = navigator.userAgent; //ȡ���������userAgent�ַ� var isOpera = us ...
- 笔记-pyrhon-lib-requests
笔记-pyrhon-lib-requests 1. 简介 Requests is the only Non-GMO HTTP library for Python, safe for hum ...
- 大话CNN经典模型:VGGNet
2014年,牛津大学计算机视觉组(Visual Geometry Group)和Google DeepMind公司的研究员一起研发出了新的深度卷积神经网络:VGGNet,并取得了ILSVRC20 ...
- Python logging 模块简介
Table of Contents 1. Logging 模块 1.1. 简介 1.2. 简单输出日志 1.3. 输入日志到文件 1.4. 几个基本概念 1.4.1. loggers 1.4.2. h ...
- 《Cracking the Coding Interview》——第6章:智力题——题目5
2014-03-20 01:08 题目:扔鸡蛋问题.有一个鸡蛋,如果从N楼扔下去恰好会摔碎,低于N楼则不碎,可以继续扔.给你两个这样的鸡蛋,要求你一定得求出N,怎么扔才能减少最坏情况下的扔的次数? 解 ...