最接近点对问题--分支法

最近接点对

一维

问题描述

现有一个一维坐标，上面有若干个点，要找到其中距离最近的点。

问题分析

在按位置升序排列的区间内，选取一个分界点mid，选取左右距离最小的，再选取mid和mid+1的距离，然后这三个值进行比较。最小的就是这个区间内最小距离点对。

在一个区间只有两个或三个点时：要进行特殊处理：

1. 两点：直接返回距离
2. 三点：返回这三点的最小

代码

#include <iostream>
#include <algorithm>
using namespace std;
const int N=100005;

double p[N];

double cal_dis(double a,double b){
	if(a>=b){
		return a-b;
	}
	else {
		return b-a;
	}
}

typedef struct {
	double p1,p2;
	double dis;
}Pair; 

Pair merge_dis(int l,int r){
	if(r-l==1){
		Pair t;
		t.dis=cal_dis(p[l],p[r]);
		t.p1=p[l],t.p2=p[r];
		return t;
	}
	else if(r-l==2){
		double dis[3],mindis;
		dis[1] = cal_dis(p[l],p[l+1]);
		dis[2] = cal_dis(p[l],p[r]);
		
		mindis=min(dis[1],dis[2]);
		Pair t;
		if(mindis==dis[1]){
			t.p1=p[l];
			t.p2=p[l+1];
			t.dis=dis[1];
		}
		else if(mindis == dis[2]){
			t.p1=p[l];
			t.p2=p[r];
			t.dis=dis[2];
		}
		return t;
	}
	int mid=l+r>>1;
	Pair pa1 = merge_dis(l,mid);
	Pair pa2 = merge_dis(mid+1,r);
	Pair pa3 ;
	pa3.p1=p[mid];
	pa3.p2=p[mid+1];
	pa3.dis = cal_dis(pa3.p1,pa3.p2);
	if(pa1.dis<=pa2.dis&&pa1.dis<=pa3.dis){
		return pa1;
	}
	else if(pa2.dis<pa1.dis&&pa2.dis<=pa3.dis){
		return pa2;
	}
	else if(pa3.dis<pa1.dis&&pa3.dis<pa2.dis){
		return pa3;
	}
}



int main(){
	int n;
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>p[i];
	}
	sort(p+1,p+n+1);
	Pair ans=merge_dis(1,n);
	cout<<ans.p1<<" "<<ans.p2<<" "<<ans.dis<<endl;
	return 0;
}

二维

问题描述

现有一个二维平面，上面有若干个点，要找到其中距离最近的点。

问题分析

与一维问题的大体思路一样，先找两边的最小的最小距离点对，然后再去算中间的点对距离。
不过与二维不一样的是算中间的时候并不是拿两个点去算：
首先两边设从两边找到的最近接点对的距离为d，那么易知，在中线两边的大小为d的区域内属于同一边的点对距离肯定大于d，而且两个点的y坐标差距大于等于d的点对肯定总距离也是大于等于d。
而这其中还有一个隐藏的限制条件：

• 若把一边的点确定到中线上，另一边的点肯定在一个以的d/2为半径的半圆内，不过为了在计算机中更好的处理，我们去一个长宽为d和d/2的长方形内。而且因为同一边的点对距离肯定大于等于d，所以这个矩形内顶多有6个点。 这是一个隐藏的降低复杂度地方。

代码

所用到的数据结构

首先每个点是用结构体存数组储的。
其中存储每个可以满足在中线两边d范围的区域内的点是用一个数组去映射到其在结构体的下标。

具体实现

#include <iostream>
#include <cstdio>
#include <map>
#include <cmath>
#include <algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
const int N = 100005;
int a[N];
typedef struct {
	double x ,y ;
}Point;
Point p[N];

double cal_dis(Point a,Point b){
	return sqrt((a.x - b.x)*(a.x - b.x)+(a.y-b.y)*(a.y-b.y));
}

bool cmpx(Point a,Point b){
	return a.x<=b.x;
}

bool cmpy(int i1,int i2){
	return p[a[i1]].y<=p[a[i2]].y;
}


double f(int l,int r){
	if(r - l == 1){
		return cal_dis(p[l],p[r]);
	}
	else if(r - l == 2){
		return min(cal_dis(p[l],p[r]),min(cal_dis(p[l],p[l+1]),cal_dis(p[l+1],p[r])));
	}
	double d;
	int mid = l + r >> 1;
	d = min(f(l,mid),f(mid+1,r));
	// merge
	int k = 1;
	for(int i=l;i<=r;i++){
		if(p[i].x>=p[mid].x-d&&p[i].x<=p[mid].x+d){
			a[k++] = i;
		}
	}
	
	sort(a+1,a+k,cmpy);
	for(int i=1;i<k;i++){
		for(int j=i+1;j<k;i++){
			if(p[a[i]].y - p[a[j]].y >= d){
				break;
			}
			d = min(d, cal_dis(p[a[i]],p[a[j]]));
		}
	}
	return d;
}

int main(){
	int n;
	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,cmpx);
	double d = f(1,n);
	cout<<d<<endl;
	return 0;
}