平面分治(求最近点对)

本文深入探讨了如何运用分治策略解决平面几何中的最近点对问题。通过对平面上的点集进行划分,递归地查找每个子区域内的最近点对,并结合边界上的点来确定全局最近点对。这种方法有效减少了计算复杂性,提高了求解效率。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include<cstdio>
#include<algorithm>
#include<cmath>
#define MAXN 10100 
#define INF 0x3f3f3f3f
using namespace std;
struct node{
	int x;
	int y;
}point[MAXN];
bool cmpx(node a,node b){
	return a.x < b.x;
}
bool cmpy(node c,node d){
	return c.y < d.y;
}
double min(double w,double v){
	return w <= v ? w : v;
}
double dis(node g,node p){
	return sqrt( (g.x-p.x)*(g.x-p.x) + (g.y-p.y)*(g.y-p.y) );
}
double point_b(int left,int right){
	double d=10000;
	if(left == right)	return d;
	if(left+1 == right)
		return dis(point[left], point[right]);
	double d1 , d2 ;
	int l = ( left + right ) / 2;
	d1 = point_b(left, l );
	d2 = point_b(l+1, right );
	d = min(d1, d2 );
	int i=left , j=right;
	while(i < l && point[l].x - point[i].x >= d)	i++;
	while(j > l && point[j].x - point[l].x >= d)	j--;
	sort(point + i + 1, point + j + 1, cmpy);
	for(int k=i ; k<=j ; k++)
	{
		for(int m=k+1 ; m<=j && point[m].y - point[k].y < d ; m++)
		{
			double d3=dis(point[m], point[k]);
			if(d > d3)
				d=d3;
		}
	}
	return d;
}
int main()
{
	int n;
	while(~scanf("%d",&n) && n)
	{
		for(int i=0 ; i<n ; i++)  
			scanf("%d %d",&point[i].x,&point[i].y);
		sort(point,point+n,cmpx);
		double temp=point_b(0,n-1);
		if(temp >= 10000)
			printf("INFINITY");
		else
			printf("%.4lf\n",temp);
	}
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值