BZOJ 2648: SJY摆棋子(K-D Tree)
Submit: 6051 Solved: 2113
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
1 1
2 3
2 1 2
1 3 3
2 4 2
Sample Output
1
2
HINT
kdtree可以过
Source
K-D Tree裸题
洛谷上需要拍扁重构才能过
丧心病狂
#include<cstdio> #include<algorithm> #define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 17, stdin), p1 == p2) ? EOF : *p1++) using namespace std; const int MAXN = 6 * 1e5 + 10, INF = 1e9 + 10; const double delat = 0.60; char buf[1 << 17], *p1 = buf, *p2 = buf; inline int read() {char c = getchar(); int x = 0, f = 1;while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();return x * f; } int N, M, WD, root, Ans; inline int abs(int x) {return x < 0 ? -x : x; } #define ls(x) T[x].ls #define rs(x) T[x].rs struct Point {int x[2];bool operator < (const Point rhs) const {return x[WD] < rhs.x[WD];} }p[MAXN]; struct Node {int ls, rs, siz, mi[2], mx[2];Point tp; }T[MAXN]; int rub[MAXN], top, cur; int NewNode() {return top ? rub[top--] : ++cur;//tag } void update(int k) {T[k].siz = T[ls(k)].siz + T[rs(k)].siz + 1;for(int i = 0; i <= 1; i++) {T[k].mi[i] = T[k].mx[i] = T[k].tp.x[i];if(ls(k)) T[k].mi[i] = min(T[k].mi[i], T[ls(k)].mi[i]), T[k].mx[i] = max(T[k].mx[i], T[ls(k)].mx[i]);if(rs(k)) T[k].mi[i] = min(T[k].mi[i], T[rs(k)].mi[i]), T[k].mx[i] = max(T[k].mx[i], T[rs(k)].mx[i]);} } int Build(int l, int r, int wd) {if(l > r) return 0;int k = NewNode(), mid = l + r >> 1;WD = wd, nth_element(p + l, p + mid, p + r + 1);T[k].tp = p[mid];T[k].ls = Build(l, mid - 1, wd ^ 1);T[k].rs = Build(mid + 1, r, wd ^ 1);update(k);return k; } inline void Apart(int k, int num) {if(T[k].ls) Apart(ls(k), num);p[num + T[ls(k)].siz + 1] = T[k].tp, rub[++top] = k;if(T[k].rs) Apart(rs(k), num + T[ls(k)].siz + 1); } inline int check(int &k, int wd) {if(T[k].siz * delat < T[ls(k)].siz || T[k].siz * delat < T[rs(k)].siz) Apart(k, 0), k = Build(1, T[k].siz, wd); } void Insert(Point a, int &k, int wd) {if(k == 0) {k = NewNode(); T[k].tp = a; update(k); return ;}if(a.x[wd] < T[k].tp.x[wd]) Insert(a, ls(k), wd ^ 1);else Insert(a, rs(k), wd ^ 1);update(k); check(k, wd); } inline int dis(Point a, Point b) {return abs(a.x[0] - b.x[0]) + abs(a.x[1] - b.x[1]); } inline int Manha(Point a, int b) {int rt = 0;for(int i = 0; i <= 1; i++) rt += max(0, a.x[i] - T[b].mx[i]) + max(0, T[b].mi[i] - a.x[i]);return rt; } int Query(Point a, int k) {Ans = min(Ans, dis(T[k].tp, a));int disl = INF, disr = INF;if(ls(k)) disl = Manha(a, T[k].ls);if(rs(k)) disr = Manha(a, T[k].rs);if(disl < disr) {if(disl < Ans) Query(a, ls(k));if(disr < Ans) Query(a, rs(k));}else {if(disr < Ans) Query(a, rs(k));if(disl < Ans) Query(a, ls(k));} } int main() {#ifdef WIN32freopen("a.in", "r", stdin);#endifN = read(); M = read();for(int i = 1; i <= N; i++) p[i].x[0] = read(), p[i].x[1] = read();root = Build(1, N, 0);while(M--) {int opt = read(), x = read(), y = read();if(opt == 1) Insert((Point){x, y}, root, 0);else Ans = INF + 1, Query((Point){x, y}, root), printf("%d\n", Ans);}return 0; }
转载于:https://www.cnblogs.com/zwfymqz/p/9062795.html
BZOJ 2648: SJY摆棋子(K-D Tree)相关推荐
- BZOJ 2648 SJY摆棋子(KD-Tree)
题目链接:BZOJ 2648 SJY摆棋子 题意: 在一个棋盘上,有N个黑色棋子.他每次要么放到棋盘上一个黑色棋子,要么放上一个白色棋子,如果是白色棋子,他会找出距离这个白色棋子最近的黑色棋子.此处 ...
- bzoj 2648: SJY摆棋子2716: [Violet 3]天使玩偶 --kdtree
2648: SJY摆棋子&&2716: [Violet 3]天使玩偶 Time Limit: 20 Sec Memory Limit: 128 MB Description 这天,S ...
- bzoj 2648 SJY摆棋子 cdq分治+树状数组
题面 题目传送门 解法 同bzoj2716 自己cdq写的还是丑啊,别人A掉了我T飞了 代码 #include <bits/stdc++.h> #define inf 1 << ...
- 2716: [Violet 3]天使玩偶/2648: SJY摆棋子
题目链接 题目大意:平面动态加点,求与给出的点曼哈顿距离最近点 题解:KD树模板题 我的收获:23333 #include <cstdio> #include <iostream&g ...
- BZOJ2648: SJY摆棋子2716: [Violet 3]天使玩偶
BZOJ2648: SJY摆棋子 BZOJ2716: [Violet 3]天使玩偶 BZOJ氪金无极限... 其实这两道是同一题. 附上2648的题面: Description 这天,SJY显得无聊. ...
- 天使玩偶/SJY摆棋子
P4169 [Violet]天使玩偶/SJY摆棋子 CDQ分治的题目. 我们发现题目要我们求的\(|A_x-B_x|+|A_y-B_y|\)的绝对值号比较恶心. 试想一下怎么去掉 如果所有的点都在我们 ...
- [Luogu4169] [Violet]天使玩偶/SJY摆棋子 [cdq分治/k-d tree]
[Link\frak{Link}Link] 考虑用 cdq 分治解决? 显然依旧分成修改跟查询,考虑一维排序二维cdq?啊但是那显然不太对,cdq能解决三维偏序, 但是她不能解决最近点对啊.平面最近点 ...
- [BZOJ2716/2648][Violet 3]天使玩偶/SJY摆棋子[KDtree]
KDtree干这个复杂度是不对的,重构不一定有作用 解释一下的话,因为复杂度是跟size相关的,所以重构作用不大,KDtree在查询最近点对中的作用仅仅是剪枝,可以构造数据使得他遍历O(n)个节点 h ...
- 【bzoj2648】SJY摆棋子 KD-tree
题目描述 这天,SJY显得无聊.在家自己玩.在一个棋盘上,有N个黑色棋子.他每次要么放到棋盘上一个黑色棋子,要么放上一个白色棋子,如果是白色棋子,他会找出距离这个白色棋子最近的黑色棋子.此处的距离是 ...
- BZOJ2648: SJY摆棋子
Description 这天,SJY显得无聊.在家自己玩.在一个棋盘上,有N个黑色棋子.他每次要么放到棋盘上一个黑色棋子,要么放上一个白色棋子,如果是白色棋子,他会找出距离这个白色棋子最近的黑色棋子. ...
最新文章
- 真正的maven私服搭建器--Artifactory
- 【Flink】Flink 与数据库的集成最佳实践 【视频笔记】
- 计算系统和计算机系统是同一个概念嘛,计算机操作系统的基本概念
- FileOprSer.class.php(文件上传与下载类)
- java redis设置过期时间_Redis的内存回收原理,及内存过期淘汰策略详解
- .net RestSharp使用
- badbody下_badboy下载_badboy测试工具2.0.5官方免费版 - 系统之家
- 【数值优化之线搜索方法】
- MongoDB的一些基本使用
- Altium Designer -- 精心总结
- java自学笔记(4)-Stanford CS106A 弹球动画 20.9.9
- Android获取经纬度以及地址信息
- 项目实战:小米商城官网及秒杀高仿(html+css+js)
- VMware安装苹果系统
- 【未解决】vcs在服务器无法使用——22-03-07工作笔记
- 【LeetCode】934. 最短的桥
- Android 接入支付宝在手机未安装支付宝客户端的情况下掉不起支付宝sdk的h5页面
- 艾默生手操器TREXLFPNAWS1S
- 证券交易的基本流程是怎样的?
- 简明Linux-Linux下的常用文件操作(1)