一、问题分析
令P为笛卡尔平面上n>1个点构成的集合,假设这些点都不一样,且按照x轴升序怕排列,还可以按照y轴升序排列并且将列命名为Q.求点对间的最小距离·
二、解析
当n<3时蛮力求解,当n>3时利用点集在x轴上的中位数m,在此处做一条垂线,将点集分成两个P1.,P2;通过P1,P2求最近对问题,但d不一定是最小距离,需要在进行验证,检查。
三代码设计
double closeset(int low, int high) {
if (low == high)
return MAX;
if (low + 1 == high)
return dist(p[low], p[high]);
int mid = (low + high) >> 1;
double ans = min(closeset(low, mid), closeset(mid + 1, high));
int i, j, cnt = 0;
for (i = low; i <= high; i++) {
if (p[mid].x - ans <= p[i].x && p[i].x <= p[mid].x + ans)
a[cnt++] = i;
}
sort(a, a + cnt, cmpy);
for (i = 0; i < cnt; i++) {
int k = i + 7 > cnt ? cnt : i + 7;
for (j = i + 1; j < k; j++) {
if (p[a[j]].y - p[a[i]].y > ans)
break;
ans = min(dist(p[a[i]], p[a[j]]), ans);
}
}
return ans;
}
四、时间复杂度
O(nlogn);
五、源代码
https://github.com/yttb/hello-world