題目描述

給定平面上n個點，找出其中的一對點的距離，使得在這n個點的所有點對中，該距離為所有點對中最小的

輸入輸出格式

輸入格式：

第一行：n；2≤n≤200000

接下來n行：每行兩個實數：x y，表示一個點的行坐標和列坐標，中間用一個空格隔開。

輸出格式：

僅一行，一個實數，表示最短距離，精確到小數點後面4位。

 1 //分治基礎題
 2 #include<cmath>
 3 #include<vector>
 4 #include<cstdio>
 5 #include<cstring>
 6
 
 #include<algorithm>
 7 using namespace std;
 8 int n,m,q,cnt;
 9 struct node{
10     int x;
11     int y;
12 }poi[200005],tmp[200005];
13 int cmp1(node a,node b){
14     if(a.x==b.x){
15         return a.y<a.y;
16     }
17     return a.x<b.x;
18 }
19 int cmp2(node a,node b){
20     return
 
  a.y<b.y;
21 }
22 double cal(node a,node b){
23     long long x=1ll*abs(a.x-b.x)*abs(a.x-b.x);
24     long long y=1ll*abs(a.y-b.y)*abs(a.y-b.y);
25     return sqrt(x+y);
26 }
27 double solve(int l,int r){
28     double mn=0x3f3f3f3f;
29     if(l==r){
30         return mn;
31     }
32     if
 
(l+1==r){
33         return cal(poi[l],poi[r]);
34     }
35     int mid=(l+r)>>1;
36     int tp=0;
37     double mnl=solve(l,mid);
38     double mnr=solve(mid+1,r);
39     mn=min(mnl,mnr);
40     for(int i=l;i<=r;i++){
41         if(abs(poi[i].x-poi[mid].x)<=mn)tmp[++tp]=poi[i];
42     }
43     sort(tmp+1,tmp+tp+1,cmp2);
44     for(int i=1;i<=tp;i++){
45         for(int j=1;j<=tp;j++){
46             if(i==j)continue;
47             if(abs(tmp[i].y-tmp[j].y)>mn)continue;
48             mn=min(mn,cal(tmp[i],tmp[j]));
49         }
50     }
51     return mn;
52 }
53 int main(){
54     scanf("%d",&n);
55     for(int i=1;i<=n;i++){
56         scanf("%d%d",&poi[i].x,&poi[i].y);
57     }
58     sort(poi+1,poi+n+1,cmp1);
59     double ans=solve(1,n);
60     printf("%.4lf",ans);
61     return 0;
62 }

平面最近點對（分治）