1.问题
[描述算法问题,首选形式化方式(数学语言),其次才是非形式化方式(日常语言)]
在二维平面上的n个点中,如何快速的找出最近的一对点。
2.解析
我们先考虑一维情形,即所有的点都在x轴上,x轴上某个点m将集合划分为s1和s2,并且两个集合的点数都是相同的,递归的求解两个子集最近对的点,(p1,p2)和(q1,q2)。如果,两个集合的最近点对都在子集s1和s2中,则 d = min{ |p1-p2| , |q1-q2| }
如果,最近点对分别在s1和s2中,则一定是(p3,q3)。|p3-q3|<d,则p3和q3两者与m的距离都不超过d,且在区间(m-d,d]和(d,m+d]各有且仅有一个点。这样,就可以在线性时间内实现合并。
推广一维,来看二维,我们仿照一维的情况先把所有点按照x(横坐标)从左到右升序排列。以X横坐标中间的点作为分界线.将平面的点分成左边和右边。找到左边的点中最近点对的距离d1,和右边最近点对的距离d2。 当前的d1,d2都是两边各自的点的最近距离,但是没有考虑到两边的点相互配对的情况,即是点对一个在左边区域,一个在右边区域。可以用上面找到的d来限制配对.即是明显超过d的两点不需要配对,如下
以中间线为中心,dmin 为半径划分一个长带,最小点对还有可能存在于s1和s2的交界处,,p点和q点分别位于S1和S1的虚线范围内,只有在这个范围内,p点和q点之间的距离才会小于dmin。
不选取2dmin 带中所有的点进行计算,而是对于S1虚框范围内的p点,只选取S2虚框范围内长2dmin,宽为dmin的中的点进行计算,由dmin的意义可知S2中任何2个S中的点的距离都不小于d。由此可以推出矩形R中最多只有6个S2中的点。
[问题的理解和推导,可用电子版直接在此编写,也可用纸笔推导,拍照嵌入本文档]
3.设计
分解
对所有的点按照x坐标(或者y)从小到大排序(排序方法时间复杂度O ( n l o g n ) O(nlogn)O(nlogn))。
根据下标进行分割,使得点集分为两个集合。
解决
递归的寻找两个集合中的最近点对。
取两个集合最近点对中的最小值m i n ( d i s l e f t , d i s r i g h t )
合并
最近距离不一定存在于两个集合中,可能一个点在集合A,一个点在集合B,而这两点间距离小于dis。
#include<iostream>
#include<math.h> //计算距离
#include<time.h> //计算运行时间
#include<fstream> //文件操作
#include<algorithm> //排序
using namespace std;
#define MAX 0x3f3f3f3f //定义无穷大
#define M 99999
struct point {
double x, y;
}p[M];
int a[M];// 保存排序的索引
int cmpx(const point& a, const point& b) //排序升序
{
return a.x < b.x;
}
int cmpy(int &a, int &b) //排序升序
{
return p[a].y < p[b].y;
}
inline double min(double a, double b) //返回两个值中较小的
{
return a < b ? a : b;
}
inline double dist(const point& a, const point& b)
{
return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}
double closeset(int low, int high)
{
if (low == high)
return MAX;
if (low + 1 == high) //即n=2,返回两点之间的距离
return dist(p[low], p[high]);
int mid = (low + high)>>1; //右移一位,相当于除以2,但右移的运算速度更快,若使用(low+high)/2求中间位置容易溢出
double ans = min(closeset(low, mid), closeset(mid+1, high)); //递归求出两边最小距离
int i, j, c = 0;
for (i = low; i <= high; i++) //统计那些点位于两虚线内,并记录
{
if (p[mid].x - ans <= p[i].x && p[i].x <= p[mid].x + ans)
a[c++] = i;
}
sort(a, a + c, cmpy);
for(i = 0; i < c; i++)//比较s1中虚线内的点和s2中虚线内的点的距离是否存在有小于两侧最小对的
{
int k = i+7 > c ? c : i+7;
for (j = i+1; j < k; j++)
{
if (p[a[j]].y - p[a[i]].y > ans) //如果位于中位线两侧的点的距离大于anx则跳出第一个循环
break;
ans = min(dist(p[a[i]], p[a[j]]), ans); //如果有两个点小于两侧的最近对,则选出最小的点
}
}
return ans;
}
int main()
{
clock_t start;
double totaltime;
start=clock();
int n; //一共多少个点
double dmin;
ifstream read_in;
read_in.open("close.txt");
read_in >> n;
cout<<"读入点数为:"<<n<<endl;
cout<<"点的坐标为:"<<endl;
for(int i=0;i<n;i++) //循环读入文件
{
cout<<"p"<<i+1<<":";
read_in>>p[i].x>>p[i].y;
cout<<p[i].x<<" "<<p[i].y<<endl;
}
sort(p,p+n,cmpx); //按照x轴排序
dmin=closeset(0, n-1);
cout<<"最近的距离是:"<<dmin<<endl;
clock_t end=clock();
totaltime=(double)(end-start)/CLOCKS_PER_SEC;
cout<<"程序运行时间是:"<<totaltime<<endl;
return 0;
}