luoguP1429 平面最近点对（加强版）

这个东西比较强。。。。。

当然，可以交智商税（开发KD—ree）

然后我并不会。。
于是乎只能看题解。。。。

就是把平面一直分割

然后每次根据左右得出的结论进行剪枝。。。

然后就可以优化过掉这题

大概是$O(NlogN)$上界挺松的

#include<bits/stdc++.h>
#define MAXN 3000001
using namespace std;

int n,q[MAXN],js;
struct node{
	double x,y;
}t[MAXN];

bool cmp(node a , node b){
	if(a.x != b.x)return a.x < b.x;
	else return a.y < b.y;
}

bool cmp2(const int &a , const int &b){
	return t[a].y < t[a].y;
}

double dis(int a , int b){
	return sqrt((t[a].x - t[b].x) * (t[a].x - t[b].x) + 
	(t[a].y - t[b].y) * (t[a].y - t[b].y));
}

void init(){
	cin>>n;
	for(int i = 1 ; i <= n ; i++)cin>>t[i].x>>t[i].y;
	sort(t + 1 , t + 1 + n , cmp);
}

double cal(int l , int r){
	if(l == r)return 1 << 20;
	if(l + 1 == r)return dis(l , r);
	int mid = (l + r) >> 1;
	double zz = min(cal(l , mid) , cal(mid + 1 , r));
	js = 0;
	for(int i = l ; i <= r ; i++)
	if(fabs(t[mid].x - t[i].x) <= zz)
	q[++js] = i;
	sort(q + 1 , q + 1 + js , cmp2);
	for(int i = 1 ; i <= js ; i++){
		for(int j = i + 1 ; j <= js && t[q[j]].y - t[q[i]].y < zz ; j++){
			zz = min(zz , dis(q[i] , q[j]));
		}
	}
	return zz;
}

int main(){
	init();
	double jj = cal(1 , n);
	printf("%.4lf\n" , jj);
}

启发的话。。。
想到暴力怎么做。。有些眉目后，要大胆去做。。。。
一些很像nlogn的数据，大多数都是分治结构。。。

考虑优化暴力。。