1.问题
最近对问题要求在包含有n个点的集合S中,找出距离最近的两个点。假设 p1(x1,y1),p2(x2,y2),……,pn(xn,yn)是平面的n个点。
2.解析
1.暴力算法,遍历每一个点对的距离求出最小值
2.分治算法,根据中位点二分左右区间得到最短距离,然后对区间[mid - d, mid + d]求最短距离,两值取小,不断递归。
对于最近对问题最容易最直接想到的就是蛮力算法,我们可以得到每两个点的距离求最小值,但这样的时间复杂度较高。我们可以采用分治的想法解决这一问题,将复杂难问题转化为简单问题解决。但是最短点对可能出现在中位点的左右两侧,所以我们还需要求出中间区间的最短点对。
3.设计
double closeset(int low ,int high){
if(low==high)return Max;
if(low+1==high)return dst(p[low],p[high]);
int mid =(low + high )>>1;
double ans=min(closeset(low,mid),closeset(mid+1,high));
int i,j,c=0;
for(i=low ;i<high;++i){
if(p[i].x>=p[mid].x-ans&&p[i].x<=p[mid].x+ans){
a[c++]=i;
}
}
sort(a,a+c,cmpy);
for(i=0;i<c;++i){
int k=i+7>c ? c : i+7;
for(j=i+1;j<k;++j){
ans = min(dst(p[a[i]], p[a[j]]), ans);
}
}
return ans;
}
double bruteForce()
{
double res = Max * 1.0;
for (int i = 1; i <= n; i++)
{
for (int j = i + 1; j <= n; j++)
{
res = min(res, dst(p[i], p[j]));
}
}
return res;
}
4.源码
https://github.com/vup-kano/suanfa/blob/main/closeset.cpp