[bzoj1336] [Balkan2002]Alien最小圆覆盖

　　最小圆覆盖。。三个for是O(n)的QAQ。。因为随机化后新的点不在当前圆内的几率不大。。

　　学习了下求中垂线的姿势...

#include<cstdio>
#include<cmath>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#define d double
using namespace std;
const int maxn=1e5+233;
const d eps=1e-8;
struct poi{d x,y;}a[maxn];
struct cir{poi O;d r;}now;
poi operator +(poi a,poi b){return (poi){a.x+b.x,a.y+b.y};}
poi operator -(poi a,poi b){return (poi){a.x-b.x,a.y-b.y};}
poi operator /(poi a,d b){return (poi){a.x/b,a.y/b};}
int i,j,k,n,m;
 
inline d dis(poi a,poi b){
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
inline poi rev(poi a){return (poi){a.y,-a.x};}
inline poi center(poi a,poi b,poi c){
    if(fabs((b.y-a.y)*(c.x-a.x)-(c.y-a.y)*(b.x-a.x))<eps){
        if(a.x>b.x)swap(a,b);
        if(a.x>c.x)swap(a,c);
        if(b.x>c.x)swap(b,c);
        return (a+c)/2.0;
    }
    poi x1=(a+b)/2.0,x2=x1+rev(a-b),y1=(a+c)/2.0,y2=y1+rev(a-c);
    if(fabs(y1.x-y2.x)<eps)swap(x1,y1),swap(x2,y2);
    d k2=(y1.y-y2.y)/(y1.x-y2.x),b2=y2.y-y2.x*k2;
    if(fabs(x1.x-x2.x)<eps)return (poi){x1.x,k2*x1.x+b2};
    d k1=(x1.y-x2.y)/(x1.x-x2.x),b1=x2.y-x2.x*k1;
    d x=(b2-b1)/(k1-k2);
    return (poi){x,k1*x+b1};
}
inline bool oncir(poi a){
    return dis(a,now.O)<now.r+eps;
}
 
int main(){
    scanf("%d",&n);srand(n+233);
    for(i=1;i<=n;i++)scanf("%lf%lf",&a[i].x,&a[i].y);
    random_shuffle(a+1,a+n+1);
    now.O=(a[1]+a[2])/2.0,now.r=dis(now.O,a[1]);
    for(i=3;i<=n;i++)if(!oncir(a[i])){
        now.O=(a[1]+a[i])/2.0,now.r=dis(a[i],a[1])/2.0;
        for(j=2;j<i;j++)if(!oncir(a[j])){
            now.O=(a[i]+a[j])/2.0,now.r=dis(a[i],a[j])/2.0;
            for(k=1;k<j;k++)if(!oncir(a[k])){
                now.O=center(a[i],a[j],a[k]);
                now.r=dis(now.O,a[k]);
            }
        }
    }
    printf("%.6lf
",now.r);
    printf("%.6lf %.6lf
",now.O.x,now.O.y);
    return 0;
}