分治算法求最近点对

该博客介绍了如何使用分治策略解决寻找n个点中距离最近的点对问题。首先对点按x坐标排序,然后分别计算左半部分和右半部分的最近点对距离,再通过合并阶段优化搜索,排除不可能成为最近点对的点,最终找到全局最小距离。文中还提到了通过限制y坐标范围并以7个点为一组进行优化的方法。
摘要由CSDN通过智能技术生成

     http://acm.hdu.edu.cn/showproblem.php?pid=1007

         先说下题意,很简单,给n个点的坐标,求距离最近的一对点之间距离的一半。第一行是一个数n表示有n个点,接下来n行是n个点的x坐标和y坐标,实数。
      这个题目其实就是求最近点对的距离。《算法导论》上有详细讲解,王晓东的书上也有代码。主要思想就是分治。先把n个点按x坐标排序,然后求左边n/2个和右边n/2个的最近距离,最后合并。合并要重点说一下,比较麻烦。
      首先,假设点是n个,编号为1到n。我们要分治求,则找一个中间的编号mid,先求出1到mid点的最近距离设为d1,还有mid+1到n的最近距离设为d2。这里的点需要按x坐标的顺序排好,并且假设这些点中,没有2点在同一个位置。(若有,则直接最小距离为0了)。
      然后,令d为d1, d2中较小的那个点。如果说最近点对中的两点都在1-mid集合中,或者mid+1到n集合中,则d就是最小距离了。但是还有可能的是最近点对中的两点分属这两个集合,所以我们必须先检测一下这种情况是否会存在,若存在,则把这个最近点对的距离记录下来,去更新d。这样我们就可以得道最小的距离d了。
      关键是要去检测最近点对,

以下是分治算法最邻近对的C语言代码: #include <stdio.h> #include <stdlib.h> #include <math.h> #define INF 1e9 //定义无穷大 typedef struct Point { double x, y; } point; int cmpx(const void *a, const void *b) { //按x坐标排序的比较函数 point *p1 = (point *)a; point *p2 = (point *)b; if (p1->x < p2->x) return -1; else if (p1->x > p2->x) return 1; else return 0; } int cmpy(const void *a, const void *b) { //按y坐标排序的比较函数 point *p1 = (point *)a; point *p2 = (point *)b; if (p1->y < p2->y) return -1; else if (p1->y > p2->y) return 1; else return 0; } double dist(point p1, point p2) { //计算两个之间的距离 return sqrt(pow(p1.x - p2.x, 2) + pow(p1.y - p2.y, 2)); } double min(double x, double y) { //返回两个数中较小的数 return x < y ? x : y; } double closestPair(point *p, int n) { double d, d1, d2, d3, dmin; //d为当前最短距离,d1,d2,d3为分治后的三部分最短距离 int i, j, k, m; point *p1, *p2; if (n <= 1) return INF; //若数小于等于1,则返回无穷大 if (n == 2) return dist(p[0], p[1]); //若数等于2,则返回这两个的距离 qsort(p, n, sizeof(point), cmpx); //按x坐标排序 m = n / 2; //取中 d1 = closestPair(p, m); //递归左半部分的最短距离 d2 = closestPair(p + m, n - m); //递归右半部分的最短距离 d = min(d1, d2); //取左右两部分中的最短距离 point *strip = (point *)malloc(n * sizeof(point)); //分配存储跨越中线的的数组 j = 0; for (i = 0; i < n; i++) { //将横坐标距离中线小于d的复制到strip数组中 if (fabs(p[i].x - p[m].x) < d) { strip[j] = p[i]; j++; } } qsort(strip, j, sizeof(point), cmpy); //按y坐标排序 dmin = d; for (i = 0; i < j; i++) { //在strip数组中找距离不超过d的最近对 for (k = i + 1; k < j && strip[k].y - strip[i].y < dmin; k++) { d3 = dist(strip[i], strip[k]); if (d3 < dmin) { dmin = d3; p1 = &strip[i]; p2 = &strip[k]; } } } free(strip); //释放strip数组 if (dmin < d) return dmin; //返回strip数组中的最近对距离 else return d; //返回左右两部分的最近对距离 } int main() { int n, i; printf("Please input the number of points:\n"); scanf("%d", &n); point *p = (point *)malloc(n * sizeof(point)); //分配存储的数组 printf("Please input the coordinates of the points:\n"); for (i = 0; i < n; i++) { scanf("%lf%lf", &p[i].x, &p[i].y); } printf("The closest pair distance is %lf.\n", closestPair(p, n)); //输出最近对距离 free(p); //释放p数组 return 0; }
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值