平面最近点对——杨子曰题目

最新推荐文章于 2022-06-02 16:18:43 发布

杨子曰

最新推荐文章于 2022-06-02 16:18:43 发布

阅读量292

点赞数 2

分类专栏：恶心的题目

本文链接：https://blog.csdn.net/HenryYang2018/article/details/88549603

版权

恶心的题目专栏收录该内容

23 篇文章 0 订阅

订阅专栏

平面最近点对——杨子曰题目

题目：我是题目

就是说给你一堆平面上的点，让你求出这些点中最近的两个点之间的距离

杨子曰：这道题应该用分治
大佬发问：Why？
“因为oj的分类里面写的是分治”（逃

为什么用分治呢？因为如果你把平面一刀切成两半，计算每一半里面的最小距离是不受另一半影响的，而且合并时我们可以有所取舍，走起：

对于一个平面，找到一个mid（敲黑板：这里的mid是点按横坐标排序后编号的mid，如果是左右距离的mid，那么一堆点全在左边，只有一个孤零零的点在很右边的话，Oh TLE），我们把它一刀切，然后算左边，算右边，取一个min，我们不妨称这个min为d

在这里插入图片描述
但是我们还要考虑合并，也就是我们还要算横跨左右的点之间的距离，怎么办呢？
“全部扫一遍!!”
——这不就是暴力吗？
在我们已经知道左右的最小值d时，所有距离大于d的点对都不可能成为答案，那么那些点对距离肯定大于d呢？
我们作这样两条线：
在这里插入图片描述
有没有发现只有在虚线里面的点才有可能成为答案（显而易见）
“对对对，这个时候就可以暴枚”
——这时候从天而降一张这样的图：

成功退化为O（n^2）,TLE!

我们再来看，对于左边的某个点，是不是右边只有这些区域的点，才可能与它一起成为答案（显而易见）：
在这里插入图片描述
一定有大佬质疑我了，如果有一堆的点都堆在绿色部分那不就凉了吗？
杨子曰：杨子曰的永远是真理

可以证明绿色部分不会超过6个点：
首先你要清楚d已经是左右两边选出来的最小值了，说明绿色部分里的点再怎么小也小不过d
然后你可以去试试看这样一来绿色部分还能摆几个点，你会发现很难摆
如果我们这样把它分成6个完全一样的矩形：
在这里插入图片描述
这样一来每个矩形的对角线就是d5/6<d，嗯，这就说明每个矩形里一旦有一个点了，第二个点就放不下了，也就是说每个矩形里面最多只有一个点，最多6个点，得证！

打代码时注意：

最开始按x排序
归并时的mid的是横坐标排序后的编号mid
合并前对两条线间的点按y排序
合并暴枚时，两重循环可以让i<j,更快
合并时由于已经对纵坐标排序了，所以一单发现纵坐标之差大于了d，就break
OK，完事

c++代码：

#include<bits/stdc++.h>
#define inf 2000000100;

using namespace std;

const int maxn=200005;

struct node{
    double x,y;
}a[maxn*2]; 

int n,q[maxn];

int cmpx(node a,node b){
    return a.x<b.x;
}

int cmpy(const int& n, const int& m)  
{  
    return a[n].y < a[m].y;  
}  

double sqr(double a){
    return a*a;
}

double calc(int i,int j){
    return (sqrt(sqr(a[i].x-a[j].x)+sqr(a[i].y-a[j].y)));
}

double g(int l,int r){
    if (r-l==1) return calc(l,r);
    if (r-l==2) return min(calc(l,r),min(calc(l+1,r),calc(l,l+1)));
    int mid=(l+r)/2;
    double d=min(g(l,mid),g(mid+1,r));
    if (d==0) return 0;
    int num=0;
    for (int i=mid;a[mid].x-a[i].x<d && i>=l;i--){
        q[++num]=i;
    }
    for (int i=mid+1;a[i].x-a[mid].x<d && i<=r;i++){
        q[++num]=i;
    }
    sort(q+1,q+num+1,cmpy);
    for (int i=1;i<num;i++){
        for (int j=i+1;j<=num;j++){
        	if (a[q[j]].y-a[q[i]].y>d) break;
            if (calc(q[i],q[j])<d) d=calc(q[i],q[j]);
        }
    } 
    return d;
}

int main(){
    scanf("%d",&n);
    for (int i=1;i<=n;i++){
        scanf("%lf%lf",&a[i].x,&a[i].y);
    }
    sort(a+1,a+n+1,cmpx);
    printf("%.4f",g(1,n));
    return 0;
}