平面最近点对

题目链接：http://oj.daimayuan.top/course/22/problem/979

 方法：我们先考虑常规做法，常规做法即是把两两点进行互相比对，最终求得最小距离

分治思想：我们把假想一个最小距离d，然后再对分区间，进行最小化d。

然后我们就找到了当前最小距离d，我们再找中间点M的周围的点，再进行最小化即可。

如何找点M周围的点？

我们只需要找M周围所有距离小于d的点，再进行以y点的排序，再暴力更新即可。

解释：因为距离公式：D=,如果y>d,x<d,D>d,则不可能更新d，当y<d,x<d,则D<d，可更新

 代码：

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=2e5+10;

struct Node{
	double x,y;
	bool operator<(const Node&A)const{//定义排序
		return x<A.x;
	}
}a[N],c[N];

double solve(int l,int r){
	if(l==r) return 1e14;//单点不存在最小值，即返回一个非常的大的值
	int m=l+r>>1;
	double d=(solve(l,m),solve(m+1,r));//选取左右区间的最小d
	int cnt=0;//取得点的数量
	for(int i=l;i<=r;i++){//加入M周围的点
		if(abs(a[i].x-a[m].x)<d){
			c[++cnt].x=a[i].y;
			c[cnt].y=a[i].x;
//因为结构体定义的是以x轴排序，所以把y放在x里面，相当于对y进行排序，点位互换，不影响距离
		}
	}
	sort(c+1,c+cnt+1);
	for(int i=1;i<=cnt;i++){//暴力枚举所有点位的距离，更新d的最小值
		for(int j=i+1;j<=cnt && a[j].x-a[i].x<d;j++){//还需要满足两点y的距离<d
			d=min(d,sqrt((c[i].x-c[j].x)*(c[i].x-c[j].x)+(c[i].y-c[j].y)*(c[i].y-c[j].y)));//更新距离
		}
	}
}

signed main(){
	int n;cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i].x>>a[i].y;
	sort(a+1,a+1+n);
	cout<<fixed<<setprecision(6)<<solve(1,n);
}