最近对问题

1.问题
[描述算法问题，首选形式化方式（数学语言），其次才是非形式化方式（日常语言）]
在二维平面上的n个点中，如何快速的找出最近的一对点。
2.解析
我们先考虑一维情形，即所有的点都在x轴上，x轴上某个点m将集合划分为s1和s2，并且两个集合的点数都是相同的，递归的求解两个子集最近对的点，(p1,p2)和(q1,q2)。如果，两个集合的最近点对都在子集s1和s2中，则 d = min{ |p1-p2| , |q1-q2| }

如果，最近点对分别在s1和s2中，则一定是（p3,q3)。|p3-q3|<d，则p3和q3两者与m的距离都不超过d，且在区间(m-d,d]和(d,m+d]各有且仅有一个点。这样，就可以在线性时间内实现合并。
推广一维，来看二维，我们仿照一维的情况先把所有点按照x(横坐标)从左到右升序排列。以X横坐标中间的点作为分界线.将平面的点分成左边和右边。找到左边的点中最近点对的距离d1,和右边最近点对的距离d2。 当前的d1,d2都是两边各自的点的最近距离,但是没有考虑到两边的点相互配对的情况,即是点对一个在左边区域,一个在右边区域。可以用上面找到的d来限制配对.即是明显超过d的两点不需要配对，如下
以中间线为中心，dmin 为半径划分一个长带，最小点对还有可能存在于s1和s2的交界处，，p点和q点分别位于S1和S1的虚线范围内，只有在这个范围内，p点和q点之间的距离才会小于dmin。

不选取2dmin 带中所有的点进行计算，而是对于S1虚框范围内的p点，只选取S2虚框范围内长2dmin，宽为dmin的中的点进行计算，由dmin的意义可知S2中任何2个S中的点的距离都不小于d。由此可以推出矩形R中最多只有6个S2中的点。

[问题的理解和推导，可用电子版直接在此编写，也可用纸笔推导，拍照嵌入本文档]
3.设计
分解
对所有的点按照x坐标（或者y）从小到大排序（排序方法时间复杂度O ( n l o g n ) O(nlogn)O(nlogn)）。
根据下标进行分割，使得点集分为两个集合。
解决
递归的寻找两个集合中的最近点对。
取两个集合最近点对中的最小值m i n ( d i s l e f t ， d i s r i g h t )
合并
最近距离不一定存在于两个集合中，可能一个点在集合A，一个点在集合B，而这两点间距离小于dis。

#include<iostream>
#include<math.h>  //计算距离
#include<time.h>  //计算运行时间
#include<fstream>  //文件操作
#include<algorithm>  //排序
 
using namespace std;
#define MAX 0x3f3f3f3f  //定义无穷大
#define M 99999
 
struct point {
  double x, y;
}p[M];
 
int a[M];// 保存排序的索引
 
int cmpx(const point& a, const point& b)  //排序升序
{  
  return a.x < b.x;
}
 
int cmpy(int &a, int &b)   //排序升序
{
  return p[a].y < p[b].y;
}
 
inline double min(double a, double b)   //返回两个值中较小的
{
	return a < b ? a : b;
}
inline double dist(const point& a, const point& b)
{
  return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}
double closeset(int low, int high)
{
	if (low == high)   
		return MAX;
	if (low + 1 == high)  //即n=2，返回两点之间的距离
		return dist(p[low], p[high]);
	int mid = (low + high)>>1;  //右移一位，相当于除以2，但右移的运算速度更快，若使用(low+high)/2求中间位置容易溢出
	double ans = min(closeset(low, mid), closeset(mid+1, high));  //递归求出两边最小距离
	int i, j, c = 0;
	for (i = low; i <= high; i++) //统计那些点位于两虚线内，并记录
	{
		if (p[mid].x - ans <= p[i].x && p[i].x <= p[mid].x + ans)
			a[c++] = i;
	}
	sort(a, a + c, cmpy);
	for(i = 0; i < c; i++)//比较s1中虚线内的点和s2中虚线内的点的距离是否存在有小于两侧最小对的
	{
		int k = i+7 > c ? c : i+7;  
		for (j = i+1; j < k; j++)
		{
			if (p[a[j]].y - p[a[i]].y > ans)  //如果位于中位线两侧的点的距离大于anx则跳出第一个循环
				break;
			ans = min(dist(p[a[i]], p[a[j]]), ans);   //如果有两个点小于两侧的最近对，则选出最小的点
		}
	}
	return ans;
}
 
int main()
{
	clock_t start;
	double totaltime;
	start=clock();
	int n;  //一共多少个点
	double dmin;
	ifstream read_in;
	read_in.open("close.txt");
	read_in >> n;
	cout<<"读入点数为："<<n<<endl;
	cout<<"点的坐标为："<<endl;
	for(int i=0;i<n;i++)  //循环读入文件
	{
		cout<<"p"<<i+1<<":";
		read_in>>p[i].x>>p[i].y;
		cout<<p[i].x<<" "<<p[i].y<<endl;
	}
	sort(p,p+n,cmpx); //按照x轴排序
	dmin=closeset(0, n-1);
	cout<<"最近的距离是："<<dmin<<endl;
	clock_t end=clock();
	totaltime=(double)(end-start)/CLOCKS_PER_SEC;
	cout<<"程序运行时间是："<<totaltime<<endl;
	return 0;
}