最接近点对问题（分治法）

最接近点对问题

  来 自

 http://www.ahhf45.com/info/Data_Structures_and_Algorithms/problems/problem_set/ndp/problem.htm

问题描述

    在应用中，常用诸如点、圆等简单的几何对象代表现实世界中的实体。在涉及这些几何对象的问题中，常需要了解其邻域中其他几何对象的信息。例如，在空中交通控制问题中，若将飞机作为空间中移动的一个点来看待，则具有最大碰撞危险的2架飞机，就是这个空间中最接近的一对点。这类问题是计算几何学中研究的基本问题之一。下面我们着重考虑平面上的最接近点对问题。

最接近点对问题的提法是:给定平面上n个点，找其中的一对点，使得在n个点的所有点对中，该点对的距离最小。

严格地说，最接近点对可能多于1对。为了简单起见，这里只限于找其中的一对。

代码如下：

#include <iostream>
#include <fstream>
#include <math.h>
using namespace std;

ifstream fin("in.txt");
struct point
{
	int x;
	int y;
}*p;


void sort(point *p,int n)
{
	point t;
	for(int i=1;i<n;i++)
	{
		t = p[i];
		for(int j=i-1;j>=0;j--)
		{
			if(t.x < p[j].x) p[j+1]=p[j];
			else break;
		}
		p[j+1]=t;
	}
}

double Distance(int i,int j)
{
	return sqrt((p[j].x-p[i].x)*(p[j].x-p[i].x)+(p[j].y-p[i].y)*(p[j].y-p[i].y));
}

double Closest(int begin,int end,int &one,int &two)
{	
	if(end-begin==1)    //两个点的情况
	{
		one = begin; two = begin+1;
		return Distance(begin,end);
	}
	double min,dist; 
	if(end-begin==2)      //三个点的情况
	{
		min = Distance(begin,begin+1);
		one = begin;	two = begin+1;
		dist = Distance(begin,end);
		if(min>dist){min=dist;two = end;}
		dist = Distance(begin+1,end);
		if(min>dist){min=dist;one = begin+1;}
		return min;
	}

	int mid = (begin+end)/2;               //多余三个点 分治法求解
	dist = Closest(begin,mid,one,two);
	int three=0,four=0;
	double dist2 =  Closest(mid,end,three,four); 
	if(dist<dist2)                           //对结果合并处理
	{min = dist;}
	else{
		min = dist2;
		one = three; two = four;
	}
	int i,j,num;
	double xia,shang;                  //关键部分！
	xia = p[mid].x-min;
	for(i=mid;i>=0;i--)             
	{
		if(p[i].x<xia)break;
		num = mid+6 < end ? mid+6:end;
		shang = p[i].y+min;
		for(j=mid+1;j<num;j++)
		{
			if(p[j].x-p[i].x > min)break;
			if(p[j].y > shang)continue;
			dist = Distance(i,j);
			if(min>dist)
			{
				min=dist;
				one = i;
				two = j;
			
			}
		}
	}
	return min;
}

int main()
{
	int n;
	fin>>n;
	p= (point *)malloc(sizeof(point)*n);
	int i;
	for(i=0;i<n;i++)
	{
		fin>>p[i].x>>p[i].y;
	}
	sort(p,n);
	for(i=0;i<n;i++)
		cout<<p[i].x<<" ";
	cout<<endl;
	int one=0,two=0;
	cout<<"The closest distance is "<<Closest(0,n,one,two)<<endl;
	cout<<"<"<<p[one].x<<","<<p[one].y<<">"<<endl;
	cout<<"<"<<p[two].x<<","<<p[two].y<<">"<<endl;
	return 0;
}

输入文件 in.txt：

10

43 67
99 35
81 36
64 78
45 65
71 94
24 61
21 34
5 29
31 51

结果：

5 21 24 31 43 45 64 71 81 99   //按x轴从小到大排序
The closest distance is 2.82843
<43,67>
<45,65>
Press any key to continue