最邻近点对

要求用分治算法（O(nlogn)复杂度）实现寻找n个点中最邻近点对,输出最邻近距离的平方
其中0<=x<10000,0<=y<10000，（x,y取整数或者小数点后一位）
点数1<n<=30000.
数据输入:
第一行表示点的数目n；
接下来的n行中，每一行依次表示点的x坐标，y坐标。
结果输出:
输出n个点中最邻近点对距离的平方(小数点后两位)

输入示例：
8
1 1
2 2
4 4
8 8
2 2.8
5 6
7 9
11 11
输出示例：
0.64

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
struct point
{
	float x;
	float y;
	point(float x = 0, float y = 0) :x(x), y(y) {}
	
}p1[100000],p2[100000];
bool cmpx(point a, point b)                                //按x排序
{
	return a.x < b.x;
}
bool cmpy(point a, point b)                                //按y排序
{
	return a.y < b.y;
}
float distance(point a, point b)                           //两点间距离
{
	return (a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y);
}

float clost(int start, int end)
{
	float d = INFINITY;
	if (start == end)return d;                                  //一个点的情况
	if (start + 1 == end)return distance(p1[start], p1[end]);     //两个点直接返回距离
	int mid = (start + end) / 2;
	float d1 = clost(start, mid);                   //左边最近距离
	float d2 = clost(mid, end);                     //右边最近距离
	d = min(d1, d2);                                //d取d1，d2最小值
	int j = 0;
	for (int i = start; i <= end; i++)
	{
		if (abs(p1[mid].x) - p1[i].x<=d)p2[j++] = p1[i];  //获取处于mid-d到mid+d之间的点，按y排序
	}
	sort(p2, p2 + j, cmpy);                                     
	for (int i = 0; i < j; i++)                               //在任意两点x或者y差值均不超过d的区间内算出每两点的距离，取最小值
	{
		for (int k = i + 1; k < j&&p2[k].y - p2[i].y < d; k++)
		{
			float d3 = distance(p2[i], p2[k]);
			d = min(d, d3);
		}
	}
	return d;


}
int main()
{
	int n; cin >> n;
	for (int i = 0; i < n; i++)
	{
		float a, b;
		cin >> a >> b;
		p1[i] = point(a, b);
	}
	sort(p1, p1 + n, cmpx);
	printf("%.2f", clost(0, n - 1));

}