用最少圆覆盖所有点问题，给出n个点的坐标和圆的半径大小，求最少多少个圆可以覆盖所有的点，n ＜= 20.

滑稽2号

已于 2023-01-27 18:17:47 修改

阅读量3.2k

点赞数 6

分类专栏：自研算法文章标签：算法贪心算法剪枝近邻算法抽象代数

于 2022-12-16 16:07:29 首次发布

本文链接：https://blog.csdn.net/weixin_66431784/article/details/128342426

版权

本文介绍了二维欧几里得平面上的最少圆覆盖问题，即给定n个点和圆的半径，如何用最少数量的圆覆盖所有点。提出了两种解法：1) 动态精度分配的近似算法，通过遍历点和生成圆心找到最优解；2) 贪心算法，寻找每次能覆盖最多点的圆。针对小规模数据，两种算法能在限制时间内找到接近答案，但因NP-hard性质，可能存在缺陷，需要多次打乱点的顺序并取最小答案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

中文版背景：最少圆覆盖

在二维欧几里得平面上，给出n个点的坐标(x, y)，和一个圆的半径大小r，求至少多少个这样半径的圆可以覆盖所有的点？(n <=20, |x|, |y| <= 180, 给定的每组(x, y)间的距离不会小于等于0.1)

英文版背景：Minimum Geometric Disk Cover Problem (MGDC):

Given a set P of n points in the Euclidean plane, the minimum geometric disk cover (MGDC) problem is to identify a minimally sized set of congruent disks with prescribed radius r that cover all the points in P.

输入：

一共有t组，每组数据中第一行是点的数量n，接下来第2到第n+1行，每行两个数据表示n个点的坐标，第n+2行是圆的半径大小

输出：

输出t行，每行一个数字，表示该组数据需要的最少的圆的数量

以下为一组样例：

input：

2

2
17.2345 20.5678
15.34 28.342
3.5

2
-180 -180
-1 -1
254

output：

2
1

解法1：暴力动态精度分配近似算法（实现简单）

题目显然是np-hard问题，我们使用近似算法来逼近正确答案，由于数据范围很小只有20，因此遍历点比遍历地图要更优的多，我们规定在1s内完成该算法，因此该程序最多进行约1kw次运算，由于浮点运算有精度的要求，我们不能同时兼顾精度和时间，因此我们根据给出的数据范围动态地安排精度, 我们设精度为d_eps, 即dynamic eps

伪代码：

遍历所有的点

该点的坐标为(x, y)

在(x - r, y - r)到(x + r, y + r)的矩形范围内以d_eps为间隔生成 $\frac{4r^2}{d\_eps^2}$ 个圆心

在这些圆中找到一个能覆盖点最多的圆(最优圆)

选择使用这个圆并移除所有被该圆覆盖的点

继续遍历没有被移除的点

时间复杂度

遍历所有的点复杂度为O(n), 生成圆心复杂度为 $O(4r^{2}\frac{1}{d\_eps^2})$ , 找到其中的最优圆复杂度为O(n), 因此总复杂度为 $O(4n^2r^2\frac{1}{d\_eps^2})$ , 我们令上面的多项式值为1kw，即可计算出预设的动态精度的值 $d\_ \,eps = \frac{1}{\sqrt{10000000/20/20/4/r^2}}$ ，该操作使得该算法在极限数据的情况下，理论运算时间恒为1s左右，实际运行中由于剪枝过程使得运算时间大大减少，若要牺牲时间换取更高的精度，则可以将1kw更换成5kw等更大的数