【牛客网】Finding Hotel
【牛客网】Finding Hotel
忘记K远点对的剪枝的我有点自闭
事实上我们只要先建一棵KD树出来,维护一下所在的矩形,和子树里的最小值
每次查询的时候如果最小值比查询的值要大的话就退出
当前的答案构成了一个圆,若圆和矩形没有交就退出(不一定很严格,可以认为是以圆心为中心向上下左右延伸半径长度的一个正方形和矩形有交)
然后看当前点在哪个子树的矩形里,先搜那个子树,如果回来后在另一个子树里可能达到的最小值都没有答案大就不搜索另一棵子树
#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define ba 47
#define MAXN 200005
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
int N,dimension,M,rt,Ncnt;
struct node {
int d[2],c,id;
}p[MAXN];
bool cmp(node a,node b) {
return a.d[dimension] < b.d[dimension];
}
struct KD {
node p;
int r[4],lc,rc,mc;
}tr[MAXN];
#define lc(u) tr[u].lc
#define rc(u) tr[u].rc
void build(int &u,int l,int r,int d) {
u = 0;
if(l > r) return;
u = ++Ncnt;
dimension = d;
int mid = (l + r) >> 1;
nth_element(p + l,p + mid,p + r + 1,cmp);
tr[u].p = p[mid];
for(int i = 0 ; i < 4 ; ++i) tr[u].r[i] = tr[u].p.d[i & 1];
build(lc(u),l,mid - 1,d ^ 1);build(rc(u),mid + 1,r,d ^ 1);
tr[u].mc = min(p[mid].c,min(tr[lc(u)].mc,tr[rc(u)].mc));
for(int i = 0 ; i < 2 ; ++i) tr[u].r[i] = min(tr[u].r[i],min(tr[lc(u)].r[i],tr[rc(u)].r[i]));
for(int i = 2 ; i < 4 ; ++i) tr[u].r[i] = max(tr[u].r[i],max(tr[lc(u)].r[i],tr[rc(u)].r[i]));
}
node ans;
int64 o(int64 x) {return x * x;}
int64 dis(node a,node b) {
return o(a.d[0] - b.d[0]) + o(a.d[1] - b.d[1]);
}
bool checkin(node x,int u) {
for(int i = 0 ; i <= 1 ; ++i) {
if(x.d[i] < tr[u].r[i] || x.d[i] > tr[u].r[i | 2]) return false;
}
return true;
}
int64 min_possible(node x,int u) {
for(int i = 0 ; i <= 1 ; ++i) {
if(x.d[i] < tr[u].r[i] || x.d[i] > tr[u].r[i | 2]) {
return min(o(tr[u].r[i] - x.d[i]),o(tr[u].r[i | 2] - x.d[i]));
}
}
return 1e18;
}
void Query(int u,node pos) {
if(!u) return;
if(tr[u].mc > pos.c) return;
if(tr[u].p.c <= pos.c) {
if(dis(tr[u].p,pos) < dis(ans,pos) || (dis(tr[u].p,pos) == dis(ans,pos) && tr[u].p.id < ans.id)) ans = tr[u].p;
}
int64 d = dis(ans,pos);
for(int i = 0 ; i <= 1 ; ++i) {
if(pos.d[i] < tr[u].r[i] && 1LL * (tr[u].r[i] - pos.d[i]) * (tr[u].r[i] - pos.d[i]) > d) return;
if(pos.d[i] > tr[u].r[i | 2] && 1LL * (tr[u].r[i | 2] - pos.d[i]) * (tr[u].r[i | 2] - pos.d[i]) > d) return;
}
int s = checkin(pos,lc(u)) ? lc(u) : rc(u);
int t = s == lc(u) ? rc(u) : lc(u);
Query(s,pos);
if(min_possible(pos,t) <= dis(pos,ans)) Query(t,pos);
}
void Solve() {
rt = 0;Ncnt = 0;
read(N);read(M);
for(int i = 1 ; i <= N ; ++i) {
read(p[i].d[0]);read(p[i].d[1]);read(p[i].c);
p[i].id = i;
}
tr[0].mc = 1e9;
tr[0].r[0] = tr[0].r[1] = 1e9;
tr[0].r[2] = tr[0].r[3] = -1;
build(rt,1,N,0);
node pos;
for(int i = 1 ; i <= M ; ++i) {
read(pos.d[0]);read(pos.d[1]);read(pos.c);
ans.d[0] = 1e9;ans.d[1] = 1e9;ans.c = 0;ans.id = 1e9;
Query(1,pos);
out(ans.d[0]);space;out(ans.d[1]);space;out(ans.c);enter;
}
}
int main(){
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
int T;
read(T);
while(T--) Solve();
return 0;
}
【牛客网】Finding Hotel的更多相关文章
- 牛客网暑期ACM多校训练营(第三场) J Distance to Work 计算几何求圆与多边形相交面积模板
链接:https://www.nowcoder.com/acm/contest/141/J来源:牛客网 Eddy has graduated from college. Currently, he i ...
- 牛客网 --java问答题
http://www.nowcoder.com/ 主要是自己什么都不怎么会.在这里可以学习很多的! 第一天看题自己回答,第二天看牛客网的答案! 1 什么是Java虚拟机?为什么Java被称作是“平台无 ...
- 牛客网《BAT面试算法精品课》学习笔记
目录 牛客网<BAT面试算法精品课>学习笔记 牛客网<BAT面试算法精品课>笔记一:排序 牛客网<BAT面试算法精品课>笔记二:字符串 牛客网<BAT面试算法 ...
- C++版 - HDUoj 2010 3阶的水仙花数 - 牛客网
版权声明: 本文为博主Bravo Yeung(知乎UserName同名)的原创文章,欲转载请先私信获博主允许,转载时请附上网址 http://blog.csdn.net/lzuacm. C++版 - ...
- 牛客网第9场多校E(思维求期望)
链接:https://www.nowcoder.com/acm/contest/147/E 来源:牛客网 题目描述 Niuniu likes to play OSU! We simplify the ...
- 牛客网暑期ACM多校训练营(第七场)Bit Compression
链接:https://www.nowcoder.com/acm/contest/145/C 来源:牛客网 题目描述 A binary string s of length N = 2n is give ...
- Beautiful Numbers(牛客网)
链接:https://ac.nowcoder.com/acm/problem/17385来源:牛客网 题目描述 NIBGNAUK is an odd boy and his taste is stra ...
- 牛客网华为机试题之Python解法
牛客网华为机试题之Python解法 第1题 字符串最后一个单词的长度 a = input().split(" ") print(len(a[-1])) 第2题 计算字符个数 a = ...
- 牛客网Wannafly挑战赛25A 因子(数论 素因子分解)
链接:https://www.nowcoder.com/acm/contest/197/A来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言5242 ...
随机推荐
- 实现返回顶部-wepy小程序-前端梳理
<script type="text/javascript" src="http://hovertree.com/ziyuan/jquery/jquery-1.11 ...
- xshel链接linuxl安装nginx
原文链接:https://blog.csdn.net/Sweet__dream/article/details/78256952?utm_source=blogxgwz9 这个连接更详细:https: ...
- 控制层@Value注解取不到值
@Value("${enable-upload-image}") private String enable; 如上所示,同样的代码,写在在业务层,运行时能取到正确的值,但在控制层 ...
- @Transactional 注解参数详解
Transactional参数说明 参数名称 功能描述 readOnly 该属性用于设置当前事务是否为只读事务,设置为true表示只读,false则表示可读写,默认值为false.例如:@Transa ...
- mysql转移Data目录
本来很简单的事,为何这么复杂,转目录后经常不成功.mysql安装时就不能指定数据目录吗?老是默认在C盘 测试环境Windows10, mysql5.7,以下是测试记录: 1. 复制Data到D盘下,如 ...
- 2018-2019-2 20165209 《网络对抗技术》Exp9: Web安全基础
2018-2019-2 20165209 <网络对抗技术>Exp9: Web安全基础 1 基础问题回答和实验内容 1.1基础问题回答 (1)SQL注入攻击原理,如何防御? 原理:SQL注入 ...
- Filename too long Resolution
在git bash中,运行下列命令: git config --global core.longpaths true --global是该参数的使用范围,如果只对本版本库设置该参数,只要在上述命令中去 ...
- Flutter移动电商实战 --(47)详细页_Flutter_html插件的使用
详情里面是hemlt和图片组成的,但是flutter是不支持html的所以需要其他插件 flutter webview plugin:这个不太好用 flutter_html:用这个插件 先解决之前一个 ...
- python @abstractmethod
1.写在前面 由于python 没有抽象类.接口的概念,所以要实现这种功能得abc.py 这个类库 2.@abstractmethod特点 @abstractmethod:抽象方法,含abstract ...
- struts2之Action与JSP相互数据传递
package com.loaderman.crm.action; import com.loaderman.crm.entity.User; import com.loaderman.crm.ser ...