BZOJ1336&1337：[Balkan2002]Alien最小圆覆盖

commonc

于 2016-08-23 16:07:29 发布

阅读量472

点赞数

分类专栏： BZOJ 最小圆覆盖数学几何 Balkan 随机化文章标签： BZOJ 数学几何 balkan 最小圆覆盖

本文链接：https://blog.youkuaiyun.com/commonc/article/details/52291842

版权

BZOJ 同时被 3 个专栏收录

142 篇文章

订阅专栏

数学

14 篇文章

订阅专栏

几何

4 篇文章

订阅专栏

本文介绍了一种解决最小圆覆盖问题的算法实现。通过输入一系列坐标点，该算法能够找到一个圆，使得所有给定的点都位于这个圆内或圆上，并且该圆的半径尽可能小。代码中使用了随机洗牌的方法来提高搜索效率。

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

题目大意：给定n个点，求最小圆覆盖

最小圆覆盖裸题....算法讲解戳→这里

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100010
using namespace std;
struct node{double x,y;}b[N];
node O;
double R;
double sqr(double x){return x*x;}
double dis(node x,node y)
{
    return sqrt(sqr(x.x-y.x)+sqr(x.y-y.y));
}
bool incircle(node x)
{
    if(dis(O,x)<=R) return true;
    return false;
}
node solve(double a,double b,double c,double d,double e,double f)
{
    double y=(f*a-c*d)/(b*d-e*a);
    double x=(f*b-c*e)/(a*e-b*d);
    return (node){x,y};
}
int main()
{
    int n;
    scanf("%d",&n);
    int i,j,k;
    for(i=1;i<=n;i++)
    scanf("%lf%lf",&b[i].x,&b[i].y);
    random_shuffle(b+1,b+n+1);
    R=0;
    for(i=1;i<=n;i++)
    if(!incircle(b[i]))
    {
        O.x=b[i].x;O.y=b[i].y;R=0;
        for(j=1;j<i;j++)
        if(!incircle(b[j]))
        {
            O.x=(b[i].x+b[j].x)/2;
            O.y=(b[i].y+b[j].y)/2;
            R=dis(O,b[i]);
            for(k=1;k<j;k++)
            if(!incircle(b[k]))
            {
                O=solve(
                b[i].x-b[j].x,b[i].y-b[j].y,(sqr(b[j].x)+sqr(b[j].y)-sqr(b[i].x)-sqr(b[i].y))/2,
                b[i].x-b[k].x,b[i].y-b[k].y,(sqr(b[k].x)+sqr(b[k].y)-sqr(b[i].x)-sqr(b[i].y))/2 
                );
                R=dis(b[i],O);
            }
        }
    }
    printf("%.10lf\n%.10lf %.10lf",R,O.x,O.y);
}