在应用中,常用诸如点、圆等简单的几何对象代表现实世界中的实体。在涉及这些几何对象的问题中,常需要了解其邻域中其他几何对象的信息。例如,在空中交通控制问题中,若将飞机作为空间中移动的一个点来看待,则具有最大碰撞危险的2架飞机,就是这个空间中最接近的一对点。这类问题是计算几何学中研究的基本问题之一。下面我们着重考虑平面上的最接近点对问题。
最接近点对问题的提法是:给定平面上n个点,找其中的一对点,使得在n个点的所有点对中,该点对的距离最小。
严格地说,最接近点对可能多于1对。为了简单起见,这里只限于找其中的一对。
首先,对于n^2的枚举,100000的数据量肯定不能通过,于是考虑是否能够在O(nlogn)或O(n)时间复杂度内解决
O(n)肯定不能实现,于是继续考虑nlogn
想到是否能通过分治来解决
一维点对最小距离很好想的,但是如何扩展到二维,这里参考下面博客,其对于平面二维点对的最小距离做法做出了比较严谨的推导和证明
参考博客:https://blog.csdn.net/matricer/article/details/53012797.
知道了这些基础知识,再来解决这道题就方便多了,直接上代码
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstdlib>
using namespace std;
const int MAX = 1000010;
const double inf = 1e50;
struct node
{
double x, y;
bool flag;
};
int n;
inline bool cmp1(node A, node B)
{
return A.x < B.x;
}
inline bool cmp2(node A, node B)
{
return A.y < B.y;
}
inline double dist(node A, node B)
{
return sqrt((A.x - B.x) * (A.x - B.x) + (A.y - B.y) * (A.y - B.y));
}
node point[MAX];
node temp[MAX];
double merge(int l, int r)
{
//边界条件
//左右相等,两点同性质返回inf
//两点不同性质返回距离
if (l == r)
return inf;
if (l + 1 == r)
{
if (point[l].flag == point[r].flag)
return inf;
else
return dist(point[l], point[r]);
}
//分治求解
int mid = (l + r) >> 1;
double ll = merge(l, mid);
double rr = merge(mid + 1, r);
//求解当前最小dist
double ans = min(ll, rr);
//合并中间dist的条件一定为距离边界小于dist的值
//k记录当前的小于ans的中间值有多少
int i, j, k = 0;
for (i = l; i <= r; i++)
{
if (fabs(point[i].x - point[mid].x) <= ans)
temp[++k] = point[i];
}
//暴力枚举求解
sort(temp + 1, temp + k + 1, cmp2);
for (i = 1; i <= k; i++)
for (j = i + 1; j <= k; j++)
{
//如果当前排序的距离已经大于了当前最小答案,肯定不符合要求
if (temp[j].y - temp[i].y > ans)
break;
//否者直接更新
if (temp[i].flag != temp[j].flag)
ans = min(ans, dist(temp[i], temp[j]));
}
return ans;
}
int main()
{
int T;
scanf("%d", &T);
while (T--)
{
scanf("%d", &n);
for (int i = 1; i <= n; i++)
{
scanf("%lf%lf", &point[i].x, &point[i].y);
point[i].flag = false;
}
for (int i = 1; i <= n; i++)
{
scanf("%lf%lf", &point[n + i].x, &point[n + i].y);
point[n + i].flag = true;
}
sort(point + 1, point + 2 * n + 1, cmp1);
printf("%.3f\n", merge(1, 2 * n));
}
return 0;
}
重要部分代码中已经给出了注释,应该不难理解