AtCoder Beginner Contest 151 F（最小圆覆盖）-CSDN博客

本文链接：https://blog.csdn.net/CHAI_NIAOJINJE/article/details/108625617

本文介绍了一种利用三点定圆原理寻找能够覆盖所有点的最小圆的算法。通过随机选择点并逐步调整圆的大小和位置，确保最终得到的圆是最小且能覆盖所有指定点的。算法中详细解释了如何通过三点确定一个圆心和半径，并通过代码示例展示了整个过程。

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

模拟，就是三点确定一个圆，然后判断剩下的点是不是在圆内，不是的话找三个圆外的点，构成一个新的圆，这个新的圆一定会比原来的圆大，注意找的时候不是一下找三个点，因为我们要求的是最小的圆，找到第一个点的时候，直接把原来的圆的半径扩大到这个点，第二个也是，但是第三个需要特判一下，因为如果直接把半径扩大到第三个点的话，就有可能直接把上面两个点直接包含了，而我们要求的是最小的圆，所以最后一次找到这个三个点时，直接三点定圆，就一定会包含所有的点
三点定圆方法：
三点定圆
AC代码：

#include<iostream>
#include<cstdio>
#include <stdio.h>
#include<algorithm>
#include<cstring>
#include <string>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<map>
#include<vector>
#include<bits/stdc++.h>
#include <set>
#define ll   long long
#define IO ios::sync_with_stdio(false);cin.tie();cout.tie(0)
#define inf 0x3f3f3f3f3f3f
#define pi 3.1415926535898
using namespace std;
const double eps=1e-12;
struct Point
{
    double x,y;
} a[500005];
Point o;
int n;
double ri;
double dis(Point a,Point b)//求距离
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

void tt(Point p1,Point p2,Point p3)//三点定圆
{
    double a,b,c,d,e,f;
    a=p2.y-p1.y;
    b=p3.y-p1.y;
    c=p2.x-p1.x;
    d=p3.x-p1.x;
    f=p3.x*p3.x+p3.y*p3.y-p1.x*p1.x-p1.y*p1.y;
    e=p2.x*p2.x+p2.y*p2.y-p1.x*p1.x-p1.y*p1.y;
    o.x=(a*f-b*e)/(2*a*d-2*b*c);
    o.y=(d*e-c*f)/(2*a*d-2*b*c);
    ri=dis(o,p1);
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf",&a[i].x,&a[i].y);
    }
    random_shuffle(a+1,a+n+1);//打乱顺序玄学加速
    o=a[1];
    ri=0;
    for(int i=1;i<=n;i++)
    {
        if(dis(a[i],o)>ri+eps)
        {
            o=a[i];
            ri=0;
            for(int j=1;j<=i-1;j++)
            {
                if(dis(o,a[j])>ri+eps)
                {
                    o.x=(a[i].x+a[j].x)/2;
                    o.y=(a[i].y+a[j].y)/2;
                    ri=dis(o,a[j]);
                    for(int k=1;k<=j-1;k++)
                    {
                        if(dis(o,a[k])>ri+eps)
                        {
                            tt(a[i],a[j],a[k]);
                        }
                    }
                }
            }
        }
    }
    printf("%.16lf\n",ri);
    return 0;
}