15.7 计算几何（最小圆覆盖模板代码）——【用圆覆盖点】

题目描述

平面直角坐标系中求一个最小圆，使这个圆能包括给定的 n 个点 。

输入描述

输出描述

第一行输出圆的坐标

第二行输出圆的半径。

输出保留6位小数。

输入输出样例

输入：

6
8.0 9.0
4.0 7.5
1.0 2.0
5.1 8.7
9.0 2.0
4.5 1.0

输出：

5.000000 5.000000
5.000000

---

最终代码c/c++

#include <bits/stdc++.h>
using namespace std;
#define eps 1e-8
const int N = 1e5+2;
int sgn(double x)
{
    if(fabs(x) < eps)  return 0;
    else return x<0?-1:1;
}
struct Point{ double x, y; };
Point p[N];      //输入点
//hypot()函数是cmath标头的库函数，用于查找给定数字的斜边，接受两个数字并返回斜边的计算结果，即sqrt(x * x + y * y) 。
double Distance(Point A, Point B){return hypot(A.x-B.x,A.y-B.y);}

//求三角形abc的外接圆的圆心：
Point circle_center(const Point a, const Point b, const Point c)
{
    Point center;
    double x1=a.x,y1=a.y;
    double x2=b.x,y2=b.y;
    double x3=c.x,y3=c.y;
    
    
    double A1 = 2 * (x2 - x1),
       B1 = 2 * (y2 - y1),
       C1 = x2 * x2 - x1 * x1 - y1 * y1 + y2 *y2,
       A2 = 2 * (x3 - x2),
       B2 = 2 * (y3 - y2),
       C2 = x3 * x3 - x2 * x2 - y2 * y2  + y3 * y3;
    double x,y;
    center.x=((C1*B2)-(C2*B1))/((A1*B2)-(A2*B1));
    center.y=((A1*C2)-(A2*C1))/((A1*B2)-(A2*B1));
   return center;
    
    
}
//点数组，点数组大小，返回圆心c，半径r：
void min_cover_circle(Point *p, int n, Point &c, double &r)
{
    //随机函数，打乱所有点。这一步很重要
    random_shuffle(p, p + n);  
    
    c=p[0]; r=0;                  //从第1个点p0开始。圆心为p0，半径为0
    for(int i=1;i<n;i++)         //扩展所有点
        //点pi在圆外部
        if(sgn(Distance(p[i],c)-r)>0)
        { 
            c=p[i]; r=0;                   //重新设置圆心为pi，半径为0
            for(int j=0;j<i;j++)         //重新检查前面所有的点。
                if(sgn(Distance(p[j],c)-r)>0)
                {   //两点定圆
                    c.x=(p[i].x + p[j].x)/2;
                    c.y=(p[i].y + p[j].y)/2;
                    r=Distance(p[j],c);
                    for(int k=0;k<j;k++)
                        if(sgn(Distance(p[k],c)-r)>0){//两点不能定圆，就三点定圆
                            c=circle_center(p[i],p[j],p[k]);
                            r=Distance(p[i], c);
                        }
                }
        }
}
int main(){
    int n;                //点的个数
    Point c; double r;  //最小覆盖圆的圆心和半径
    scanf("%d",&n);
    for(int i=0;i<n;i++)  scanf("%lf%lf",&p[i].x,&p[i].y);
    
    min_cover_circle(p,n,c,r);
    printf("%.6f %.6f\n",c.x,c.y);
    printf("%.6f\n",r);
    return 0;
}

---

补充知识：外心、已知三点求外心的公式

外心是：所形成的三角形的外接圆的圆心，是三条直线的垂直平分线

---

最小圆覆盖思想

最小圆覆盖问题即给定 n 个点的平面坐标，求一个半径最小的圆，把 n 个点全部包围，部分点在圆上。

分析一下：
首先这个最小圆可以由n个点中的两个点或三个点确定。
两点定圆时，圆心是线段AB的中点，半径是AB长度的一半，其它点都在这个圆内。
如果两点不足以包围所有点，就需要三点定圆，此时圆心是A、B、C这3个点组成的三角形的外心，如下图：

所以最小覆盖圆的获得，就是寻找能两点定圆或三点定圆的那几个点？

大佬：是的。一般我们用增量法求最小圆覆盖。算法从一个点开始，每次加入一个新的点，更新最小圆，直到扩展到全部n个点。设前i个点的最小覆盖圆是Ci，过程如下：

1. 加第1个点p1。C1的圆心就是p1，半径为0。

2. 加第2个点p2。新的C2的圆心是线段p1，p2的中心，半径为两点距离的一半。这一步操作是两点定圆。

3. 加第3个点p3。有两种情况：
   p3在C2的内部或圆周上，不影响原来的最小圆，忽略p3
   p3在C2的外部，此时C2已不能覆盖所有3个点，需要更新。而当p3在C2的外部时，因为P3一定在新的C3上，所以现在的任务转换为：在p1、p2中找1个点或2个点，与p3一起两点定圆或三点定圆。重新定圆的过程，相当于回到第1步，把p3作为第1个点加入，然后再加入p1、P2。

4. 加第4个点p4。分析和步骤3类似，为加强理解，这里重复说明一次。
   如果p4在C3内部或圆周上，忽略它。
   如果在C3外部，那么需要求新的最小圆，此时p4肯定在新的C4的圆周上。任务转换为在p1、p2、p3中找1个点或2个点，与p4一起构成最小圆。我们先检查能不能找到一个点，用两点定圆；如果两点不够，就找到第三个点，用三点定圆。重新定圆的过，程和前3个步骤类似，即把p4作为第1个点加入，然后加入p1、p2、p3。

5. 持续进行下去，直到加完所有点。

概括一下算法思路：假设已经求得前 i-1个点的C_i-1，现在加入第i个点，有两种情况：

1. i 在C_i-1的内部或圆周上，忽略i。
2. i 在C_i-1的外部，需要求新的C_i。首先，i肯定在C_i上，然后重新把前面的i-1个点依次加入，根据两点定圆或者三点定圆，重新构造最小圆。

---

模板代码的复杂度