7-10 哪两个点之间的距离最近[增强版]

最新推荐文章于 2023-12-18 00:00:00 发布

乔冠宇

最新推荐文章于 2023-12-18 00:00:00 发布

阅读量470

点赞数 8

文章标签：算法

本文链接：https://blog.csdn.net/2301_77523019/article/details/135042285

版权

设P={(x1,y1),(x2,y2),⋯,(xn,yn)}是平面上散列的n个点的集合。请编写程序找出集合中距离最近的点对。严格地说，相同距离的最近点对可能不止一对，为了简单期间只找出第一对最近点对即可。题目保证输入的所有数据均按照横坐标值进行升序排列。

输入格式：

输入第一行给出一个正整数n，表示平面上点的数数量。随后n行，每行给出一个实数对，每个实数对表示一个点的坐标值，其中第1数表示横坐标，第2数表示纵坐标。

输出格式：

输出最近点对中两个点的坐标和它们之间的距离。如果 x1+y1<=x2+y2则按
(x1,y1),(x2,y2),miniDist=Distance
输出结果，否则按
(x2,y2),(x1,y1),miniDist=Distance
输出结果。
其中x1,y1,x2,y2是保留两位小数的实数，Distance是保留3位小数的实数

输入样例:

5
-1.00 2.00
0.00 2.00
0.50 0.60
1.00 1.00
2.00 0.00

输出样例:

(0.50,0.60),(1.00,1.00),miniDist=0.640

代码长度限制16 KB

时间限制15 ms

内存限制64 MB

答案：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <algorithm>
#define MAXN 200001

double min(double a, double b) { return a < b ? a : b; }

int n, temp[MAXN];
struct point {
    double x, y;
    friend bool operator < (point p1, point p2) {
        if (p1.x == p2.x) return p1.y < p2.y;
        return p1.x < p2.x;
    }
}p[MAXN];

bool cmp(int a, int b) {
    return p[a].y < p[b].y;
}

double getd(int x, int y) {
    return sqrt((p[x].x - p[y].x) * (p[x].x - p[y].x) + (p[x].y - p[y].y) * (p[x].y - p[y].y));
}

struct closest_pair {
    double dist;
    point p1, p2;
}cp;

closest_pair merge(int l, int r) {
    closest_pair res;
    res.dist = 2147483647.0;
    if (l == r) {
        res.p1 = p[l];
        res.p2 = p[l];
        return res;
    }
    if (l == r - 1) {
        res.dist = getd(l, r);
        if (p[l].x + p[l].y <= p[r].x + p[r].y) {
            res.p1 = p[l];
            res.p2 = p[r];
        }
        else {
            res.p1 = p[r];
            res.p2 = p[l];
        }
        return res;
    }
    int mid = (l + r) >> 1;
    closest_pair d1 = merge(l, mid);
    closest_pair d2 = merge(mid + 1, r);
    if (d1.dist < d2.dist) {
        res = d1;
    }
    else {
        res = d2;
    }
    int cnt = 0;
    for (int i = l; i <= r; ++i) {
        if (fabs(p[mid].x - p[i].x) < res.dist) {
            temp[++cnt] = i;
        }
    }
    std::sort(temp + 1, temp + cnt + 1, cmp);
    for (int i = 1; i <= cnt; ++i) {
        for (int j = i + 1; j <= cnt && p[temp[j]].y - p[temp[i]].y < res.dist; ++j) {
            double d = getd(temp[i], temp[j]);
            if (d < res.dist) {
                res.dist = d;
                if (p[temp[i]].x + p[temp[i]].y <= p[temp[j]].x + p[temp[j]].y) {
                    res.p1 = p[temp[i]];
                    res.p2 = p[temp[j]];
                }
                else {
                    res.p1 = p[temp[j]];
                    res.p2 = p[temp[i]];
                }
            }
        }
    }
    return res;
}

int main() {
    scanf("%d", &n);
    for (int i = 1; i <= n; ++i) {
        scanf("%lf %lf", &p[i].x, &p[i].y);
    }
    std::sort(p + 1, p + n + 1);
    cp = merge(1, n);
    printf("(%.2lf,%.2lf),(%.2lf,%.2lf),miniDist=%.3lf\n", cp.p1.x, cp.p1.y, cp.p2.x, cp.p2.y, cp.dist);
    return 0;
}