平面最近对(p1429)

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

题目背景

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

题目描述

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

输入格式

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

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

输出格式

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

输入输出样例 #1

输入 #1

3
1 1
1 2
2 2

输出 #1

1.0000

说明/提示

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

AC代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1e6;

struct Node{
    int x,y;
} s[maxn];

bool cmp(const Node &a,const Node &b){
    if(a.x==b.x) return a.y<b.y;
    return a.x<b.x;
}

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

ll dis(int a,int b){
    return (ll)(s[a].x-s[b].x)*(s[a].x-s[b].x)+(ll)(s[a].y-s[b].y)*(s[a].y-s[b].y);
}

int tmp[maxn];

ll memrge(int l,int r){
    if(l>=r) return LLONG_MAX;
    if(l+1==r) return dis(l,r);
    int mid=(l+r)>>1;
    ll d1=memrge(l,mid);
    ll d2=memrge(mid+1,r);
    ll d=min(d1,d2);
    int t=0;
    for (int i = l; i <= r; i++)
    {
        if ((s[i].x - s[mid].x) * (s[i].x - s[mid].x) < d)
        {
            tmp[t++]=i;
        }
    }
    sort(tmp,tmp+t,cmp1);
    for(int i=0;i<t;i++){
        for (int j = i + 1; j < t && (s[tmp[j]].y - s[tmp[i]].y) * (s[tmp[j]].y - s[tmp[i]].y) < d; j++)
        {
            d=min(d,dis(tmp[i],tmp[j]));
        }
    }
    return d;
}

int main(){
    int n;
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>s[i].x>>s[i].y;
    }
    double ans;
    sort(s+1,s+n+1,cmp);
    ans=memrge(1,n);
    ans=sqrt(ans);
    printf("%.4lf\n",ans);
    return 0;
}