学习笔记--最小圆覆盖

随机化各点后暴力枚举，发现三点定的圆不能包含当前点就更新。。复杂度期望O（n）（不会证。。）

模板题：bzoj1336

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define db double
using namespace std;
const int N=5e5+10;

db eps=1e-10;

int dcmp(db x)
{return x<-eps?-1:x>eps;}

struct Pt{
    db x,y;
    Pt(){};
    Pt(db x,db y):x(x),y(y){};
}pt[N];
int n,H,R,D,cost[N];

bool operator <(const Pt &x,const Pt &y)
{return dcmp(x.x-y.x)==0?dcmp(x.y-y.y)<0:dcmp(x.x-y.x)<0;}
db sqr(db x)
{return x*x;}
db dis(Pt x,Pt y)
{return sqrt(sqr(x.x-y.x)+sqr(x.y-y.y));}
void get2(Pt &O,db &r,int i,int j)
{
    O.x=(pt[i].x+pt[j].x)/2.0;
    O.y=(pt[i].y+pt[j].y)/2.0;
    r=dis(O,pt[i]);
    return;
}
db cj(db x,db y,db xx,db yy)
{return x*yy-y*xx;}
void get3(Pt &O,db &r,int i,int j,int k)
{
    if(cj(pt[i].x-pt[j].x,pt[i].y-pt[j].y,pt[i].x-pt[k].x,pt[i].y-pt[k].y)==0)
    {
        Pt tmp1,tmp2;
        tmp1=Pt(pt[i].x,pt[i].y),tmp2=Pt(pt[i].x,pt[i].y);
        if(pt[j]<tmp1)tmp1=pt[j];
        if(pt[k]<tmp1)tmp1=pt[k];
        if(tmp2<pt[j])tmp2=pt[j];
        if(tmp2<pt[k])tmp2=pt[k];
        O.x=(tmp1.x+tmp2.x)/2.0,O.y=(tmp1.y+tmp2.y)/2.0,r=dis(O,tmp1);
        return;
    }
    db a1,a2,b1,b2,c1,c2,tp;
    a1=2*(pt[i].x-pt[j].x),b1=2*(pt[i].y-pt[j].y),c1=sqr(pt[i].x)-sqr(pt[j].x)+sqr(pt[i].y)-sqr(pt[j].y);
    a2=2*(pt[i].x-pt[k].x),b2=2*(pt[i].y-pt[k].y),c2=sqr(pt[i].x)-sqr(pt[k].x)+sqr(pt[i].y)-sqr(pt[k].y);
    if(a1==0)O.y=c1/b1,O.x=(c2-b2*O.y)/a2;
    else tp=a2/a1,O.y=(c2-c1*tp)/(b2-b1*tp),O.x=(c1-b1*O.y)/a1;
    r=dis(O,pt[i]);
}
int idx[N];

int main()
{   
    srand(517517);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)idx[i]=i;
    random_shuffle(idx+1,idx+n+1);
    for(int i=1;i<=n;i++)
        scanf("%lf%lf",&pt[idx[i]].x,&pt[idx[i]].y);
    db r;Pt O;
    if(n==1){O.x=pt[1].x,O.y=pt[1].y,r=0;goto lb;}
    get2(O,r,1,2);
    for(int i=3;i<=n;i++)
    {
        if(dcmp(dis(O,pt[i])-r)>0)
        {
            O=pt[i],r=0;
            for(int j=1;j<i;j++)
            {
                if(dcmp(dis(O,pt[j])-r)>0)
                {
                    get2(O,r,i,j);
                    for(int k=1;k<j;k++)
                        if(dcmp(dis(O,pt[k])-r)>0)get3(O,r,i,j,k);
                }
            }
        }
    }
    lb:;printf("%.10lf\n%.10lf %.10lf\n",r,O.x,O.y);
    return 0;
}