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

题目描述
给定平面上n个点，找出其中的一对点的距离，使得在这n个点的所有点对中，该距离为所有点对中最小的

输入格式
第一行：n；2≤n≤200000

接下来n行：每行两个实数：x y，表示一个点的行坐标和列坐标，中间用一个空格隔开。

输出格式
仅一行，一个实数，表示最短距离，精确到小数点后面4位。

输入输出样例
输入 #1复制
3
1 1
1 2
2 2
输出 #1复制
1.0000
说明/提示
0<=x,y<=10^9

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=1001000;
const int inf=2 << 20;
int n,temp[maxn];
struct point{
	double x,y;
}s[maxn];
bool cmp(const point &a,const point &b){
	if(a.x==b.x)
		return a.y<b.y;
	else
		return a.x<b.x;
}
bool cmps(const int &a,const int &b) { return s[a].y < s[b].y; }

double min(double a, double b) { return a < b ? a : b; }
double dist(int i, int j) 
{
    double x = (s[i].x - s[j].x) * (s[i].x - s[j].x);
    double y = (s[i].y - s[j].y) * (s[i].y - s[j].y);
    return sqrt(x + y);
}
double merge(int left, int right) 
{
	double d=inf;
	if(left==right)
		return d;
	if(left+1==right)
        return dist(left,right);
    int mid=(left+right)>>1;
	double d1=merge(left,mid);
    double d2=merge(mid+1,right);   
	d=min(d1,d2);
	int k=0;
	for(int i=left;i<=right;i++){
		if(fabs(s[mid].x-s[i].x)<d)
			temp[k++]=i;
	}
	sort(temp,temp+k,cmps);
	for (int i = 0; i < k; i++)
        for (int j = i + 1; j < k && s[temp[j]].y - s[temp[i]].y < d; j++) 
		{
            double d3 = dist(temp[i], temp[j]);
            if (d > d3)
                d = d3;
        }
    return d;
}
int main() 
{
    scanf("%d", &n);
    for (int i = 0; i < n; i++) scanf("%lf%lf", &s[i].x, &s[i].y);
    sort(s, s + n, cmp);
    printf("%.4lf\n", merge(0, n - 1));
    return 0;
} 

学习自洛谷 syksykCCC

merge(left,right) 是返回由编号 left和right 的点构成的最近点对的距离

Temp 包含了离直线距离不超过 d 的所有点

如果离mid位置的竖线大于d则跨越竖线不可能小于d了

&& s[temp[j]].y - s[temp[i]].y < d是剪枝