POJ 3714 最近点对分治

最新推荐文章于 2021-02-25 20:01:33 发布

ljq真的不能再吃了

最新推荐文章于 2021-02-25 20:01:33 发布

阅读量453

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM

本文链接：https://blog.youkuaiyun.com/heucodesong/article/details/83509047

ACM 专栏收录该内容

455 篇文章

订阅专栏

博客介绍了如何运用最经典的分治算法解决POJ 3714题目，即在给定2n个点的情况下，找出前n个点与后n个点中距离最近的点对，关键在于添加id来区分不同组别的点。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

POJ 3714

题意输入n 给你2n个点前n和后n 距离最近的点是？

我们用最经典的最近点分治去做只不过加个id就行了 id区别不同的组别

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
int n;
struct P
{
    int id;
    double x, y;
    bool operator <(const P& B)const { return x < B.x; }
}p[100050*2];
double dis(P A, P B) { if(A.id!=B.id) return (A.x-B.x)*(A.x-B.x) + (A.y-B.y)*(A.y-B.y);else return 3e10+9; }
P Q[100050*2];
double Divide(int l, int r)
{
    if(l == r) return 3e10+9;
    int mid = (l+r)>>1, tx = p[mid].x, tot = 0;
    double d;
    d = min(Divide(l, mid), Divide(mid+1, r));
    for(int i = l, j = mid+1; (i <= mid || j <= r); i++)
    {
        while(j <= r && (p[i].y > p[j].y || i > mid)) Q[tot++] = p[j], j++; //归并按y排序
        if(abs(p[i].x - tx) < d && i <= mid)  //选择中间符合要求的点
        {

            for(int k = j-1; k > mid && j-k < 3; k--) d = min(d, dis(p[i], p[k]));
            for(int k = j; k <= r && k-j < 2; k++)  d = min(d, dis(p[i], p[k]));
        }
        if(i <= mid) Q[tot++] = p[i];
    }
    for(int i = l, j = 0; i <= r; i++, j++) p[i] = Q[j];
    return d;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
    scanf("%d",&n);
    for(int i = 1; i <= n; i++) scanf("%lf%lf",&p[i].x,&p[i].y),p[i].id = 1;
    for(int i = n+1;i<=2*n;++i) scanf("%lf%lf",&p[i].x,&p[i].y),p[i].id = 2;
    sort(p+1, p+1+2*n);
    printf("%.3f\n",sqrt(Divide(1,2*n)));
    }
    return 0;
}