算法导论复习（四）| 平面最近点对

brilliantgby

已于 2023-12-30 15:39:15 修改

阅读量397

点赞数 9

分类专栏：算法导论文章标签：算法

于 2023-12-30 15:34:28 首次发布

本文链接：https://blog.csdn.net/brilliantgby_id/article/details/135302498

版权

算法导论专栏收录该内容

11 篇文章 1 订阅

订阅专栏

文章描述了一个使用分治法解决平面上n个点间最近点对查找的问题，通过递归分割点集并暴力求解跨越分割线的点对，时间复杂度为O(nlogn)。给出了C++代码实现，包括点结构定义、排序和分治函数mergenode()的描述。

摘要由CSDN通过智能技术生成

文章目录

问题描述

已知平面上分布着点集P中的n个点p1,p2,…pn，点i的坐标记为(xi,yi)，1≤i≤n。两点之间的距离取其欧式距离。

问题：找出一对距离最近的点。

算法步骤

采用分治法。

记, d_L：P_L中的最近点对距离；d_R：P_R中的最近点对距离；d_C：跨越分割线的最近点对距离

做分割线（垂线）将点集分为P_L和P_R两部分
递归地在P_L中的找具有d_L的点对
递归地在P_R中的找具有d_R的点对
在跨越分割线的点对中找具有d_C的点对
返回min(d_L,d_R,d_C)对应的点对。

主要思想：

在考虑跨越分割线时注意7个点：

对坐标排序导致O(nlog²n)时间复杂度的改进：

算法实现

贴一份O(nlogn)的算法（实在不会写了QAQ）：

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
int n;
struct node{
	double x,y;
}p[200005],q[200005];

bool cmp(node a,node b){
	return a.x<b.x;
}

double mergenode(int l,int r)
{
	if(l==r)
		return 100000000.0;	
	int mid=(l+r)>>1;
	double ans=min(mergenode(l,mid),mergenode(mid+1,r)); //分治
	int cnt=0;
	for(int i=l;i<=r;++i) 
		if(fabs(p[i].x-p[mid].x)<ans)
			q[++cnt]=p[i];
	for(int i=1;i<cnt;++i) //暴力求解跨越分割线的点
		for(int j=i+1;j<=cnt;++j) 
			ans=min(ans,sqrt((q[i].x-q[j].x)*(q[i].x-q[j].x)+(q[i].y-q[j].y)*(q[i].y-q[j].y)));
	return ans;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;++i)
		scanf("%lf%lf",&p[i].x,&p[i].y);
	sort(p+1,p+n+1,cmp); //排序预处理
	printf("%.4lf",mergenode(1,n));
	return 0;
}