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

题目背景

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

题目描述

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

输入格式

第一行：$n$ ，保证 $2\le n\le 200000$ 。

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

输出格式

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

样例 #1

样例输入 #1

3
1 1
1 2
2 2

样例输出 #1

1.0000

提示

数据保证 $0\le x,y\le 10^9$

分治法o(nlog n)解法：视频题解

题解：

#include <iostream>
#include <algorithm>
#include <iomanip>
#include <cmath>
using namespace std;
const int N = 200005;
struct p
{
    long long x;
    long long y;
};
int n;
struct p a[N], b[N];
bool cmp1(struct p a, struct p b)
{
    if (a.x == b.x)
    {
        return a.y < b.y;
    }
    return a.x < b.x;
}
bool cmp2(struct p a, struct p b)
{
    return a.y < b.y;
}
double dist(struct p a, struct p b)
{
    double ans = sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
    return ans;
}
double solve(int l, int r)
{
    if (l == r)
    {               // 表示同一个点
        return 2e9; // 返回无穷大
    }
    if (l == r - 1)
    { // 相邻的两个点
        return dist(a[l], a[r]);
    }
    int mid = l + r >> 1;
    double d = min(solve(l, mid), solve(mid + 1, r)); // 递归寻找左右区间的最小值
    int k = 0;
    for (int i = l; i <= r; i++)
    {
        if (abs(a[i].x - a[mid].x) < d) // 只有小于d才有更优的可能
        {
            b[++k] = a[i];
        }
    }
    sort(b + 1, b + k + 1, cmp2);
    for (int i = 1; i < k; i++)
    {
        for (int j = i + 1; j <= k && b[j].y - b[i].y < d; j++)
        {
            d = min(d, dist(b[i], b[j]));
        }
    }
    return d;
}
int main()
{
    cin >> n;
    for (int i = 1; i <= n; i++)
    {
        cin >> a[i].x >> a[i].y;
    }
    sort(a + 1, a + n + 1, cmp1);
    cout << fixed << setprecision(4) << solve(1, n);
    return 0;
}