传送门

解题思路

  \(KD-Tree\)模板题,\(KD-Tree\)解决的是多维问题,它是一个可以储存\(K\)维数据的二叉树,每一层都被一维所分割。它的插入删除复杂度为\(log^2 n\),它查询最近点对的复杂度为\(O(n^{\frac{k-1}{k}}\),\(k\)代表维数。用堆维护最近点,查询时就先找到它属于的区域,然后回溯时判断一下它到父节点的距离和堆顶的大小,如果比堆顶还大就不递归它的兄弟节点。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm> using namespace std;
const int N=50005; inline int rd(){
int x=0,f=1; char ch=getchar();
while(!isdigit(ch)) f=ch=='-'?0:1,ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
return f?x:-x;
} inline int pw(int x){
return x*x;
} int n,q,K,t; struct Node{
int a[7];
void init() {
memset(a,0,sizeof(a));
}
friend bool operator<(const Node A,const Node B){
return A.a[t]<B.a[t];
}
}node[N],pt[N<<2],ans[25];
priority_queue<pair<double,Node> > Q; struct KD_Tree{
#define mid ((l+r)>>1)
int end[N<<2];
void build(int x,int l,int r,int dep){
if(l>r) return; t=dep%K;
end[x]=0; end[x<<1]=end[x<<1|1]=1;
nth_element(node+l,node+mid,node+r+1);
pt[x]=node[mid];
build(x<<1,l,mid-1,dep+1); build(x<<1|1,mid+1,r,dep+1);
}
void query(int x,int dep,int lim,Node now){
if(end[x]) return;
pair<double,Node> tmp=make_pair(0,pt[x]);
for(int i=0;i<K;i++) tmp.first+=pw(pt[x].a[i]-now.a[i]);
int ls=x<<1,rs=x<<1|1,t=dep%K,flag=0;
if(now.a[t]>=pt[x].a[t]) swap(ls,rs);
if(!end[ls]) query(ls,dep+1,lim,now);
if(Q.size()<lim) Q.push(tmp),flag=1;
else {
if(Q.top().first>tmp.first) Q.pop(),Q.push(tmp);
if(pw(pt[x].a[t]-now.a[t])<Q.top().first) flag=1;
}
if(!end[rs] && flag) query(rs,dep+1,lim,now);
}
#undef mid
}tree; int main(){
while(~scanf("%d%d",&n,&K)){
for(int i=1;i<=n;i++)
for(int j=0;j<K;j++) node[i].a[j]=rd();
tree.build(1,1,n,0);
for(q=rd();q;q--){
Node now; now.init();
for(int i=0;i<K;i++) now.a[i]=rd();
int t=rd(); tree.query(1,0,t,now);
for(int i=1;!Q.empty();i++)
ans[i]=Q.top().second,Q.pop();
printf("the closest %d points are:\n",t);
for(int i=t;i;i--){
printf("%d",ans[i].a[0]);
for(int j=1;j<K;j++)
printf(" %d",ans[i].a[j]);
putchar('\n');
}
}
}
return 0;
}

hud 4347 The Closest M Points(KD-Tree)的更多相关文章

  1. bzoj 3053 HDU 4347 : The Closest M Points kd树

    bzoj 3053 HDU 4347 : The Closest M Points  kd树 题目大意:求k维空间内某点的前k近的点. 就是一般的kd树,根据实测发现,kd树的两种建树方式,即按照方差 ...

  2. BZOJ 3053: The Closest M Points(K-D Tree)

    Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1235  Solved: 418[Submit][Status][Discuss] Descripti ...

  3. hdu 4347 The Closest M Points (kd树)

    版权声明:本文为博主原创文章,未经博主允许不得转载. hdu 4347 题意: 求k维空间中离所给点最近的m个点,并按顺序输出  . 解法: kd树模板题 . 不懂kd树的可以先看看这个 . 不多说, ...

  4. hdu 4347 The Closest M Points(KD树)

    Problem - 4347 一道KNN的题.直接用kd树加上一个暴力更新就撸过去了.写的时候有一个错误就是搜索一边子树的时候返回有当前层数会被改变了,然后就直接判断搜索另一边子树,搞到wa了半天. ...

  5. 数据结构(KD树):HDU 4347 The Closest M Points

    The Closest M Points Time Limit: 16000/8000 MS (Java/Others)    Memory Limit: 98304/98304 K (Java/Ot ...

  6. 【HDOJ】4347 The Closest M Points

    居然是KD解. /* 4347 */ #include <iostream> #include <sstream> #include <string> #inclu ...

  7. HDU 4347 - The Closest M Points - [KDTree模板题]

    本文参考: https://www.cnblogs.com/GerynOhenz/p/8727415.html kuangbin的ACM模板(新) 题目链接:http://acm.hdu.edu.cn ...

  8. BZOJ3053:The Closest M Points(K-D Teee)

    Description The course of Software Design and Development Practice is objectionable. ZLC is facing a ...

  9. HDU 4347 The Closest M Points (kdTree)

    赤果果的kdTree. 学习传送门:http://www.cnblogs.com/v-July-v/archive/2012/11/20/3125419.html 其实就是二叉树的变形 #includ ...

随机推荐

  1. 升级到mysql5.7无法启动问题解决

    漏洞扫描,老项目升级到5.7位成功,启动发现报错:unknown option log_error 线备份my.cnf配置文件, 猜测应该是写法有问题,先把log_error  改成log #log_ ...

  2. jmeter逻辑控制详解(1)

    逻辑控制器 Jmeter提供了多种逻辑控制器,下面进行讲解说明: 1.Simple Controller 简单控制器是最基本的控制器,对jmeter测试运行没有任何影响,可以将某些请求归集在一个简单控 ...

  3. python每日一练:0001题

    第 0001 题:做为 Apple Store App 独立开发者,你要搞限时促销,为你的应用生成激活码(或者优惠券),使用 Python 如何生成 200 个激活码(或者优惠券)? import o ...

  4. 【报错】springboot thymeleaf超链接跳转 404

    Whitelabel Error Page This application has no explicit mapping for /error, so you are seeing this as ...

  5. JavaScript 开发的 睡眠状况自测(SRSS)

    Javascript 开发睡眠状况自测程序,手记!2019.11.13日... <script>//初始化fbox = new Findpair('fbox','output');fbox ...

  6. [Python3 练习] 009 利用列表隐藏并找到有用的信息

    题目:利用列表隐藏并找到有用的信息 (1) 描述 1) 题源 鱼 C 论坛中"小甲鱼"老师出的题 链接地址:第020讲:函数:内嵌函数和闭包 | 课后测试题及答案 2) 修改 题中 ...

  7. ubuntu终端代理之proxychains

    命令行代理 安装proxychains sudo apt install proxychains 配置proxychains sudo vim /etc/proxychains.conf 在proxy ...

  8. vue.js(11)--案例--关键字搜索列表

    关键字搜索品牌案例 (1)页面布局 <div class="app"> <div class="panel panel-primary"> ...

  9. 使用JSONP,jQuery的ajax跨域获取json数据

    网上找了很多资料,写的不错,推荐下: 1.深入浅出JSONP--解决ajax跨域问题 (http://www.cnblogs.com/chopper/archive/2012/03/24/240394 ...

  10. day02-css

    技术分析 HTML的块标签: div标签: 默认占一行,自动换行 span标签: 内容显示在同一行 CSS概述: Cascading Style Sheets : 层叠样式表 主要用作用: 用来美化我 ...