[BZOJ 1336/1337][Balkan2002]Alien最小圆覆盖：随机增量法

点击这里查看原题

随机增量法的思路：
A：如果遇到某个点i在圆外，那么扩圆后这个点一定在圆上。
B：进行扩圆，如果有之前的点j在当前的圆外，那么i,j一定都在圆上，以i,j两点的连线作为圆的直径。
C：继续检查，如果有j之前的点k在圆外，那么i,j,k三个点一定都在圆上，于是做出三角形ijk的最小覆盖圆

/*
User:Small
Language:C++
Problem No.:1336
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 999999999
#define pii pair<int,int>
#define mp make_pair
using namespace std;
const int M=1e5+5;
const double eps=1e-6;
int n;
double r;
struct no{
    double x,y;
}p[M],o;
double sqr(double x){
    return x*x;
}
double dis(no a,no b){
    return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
}
void cal(double a,double b,double c,double d,double e,double f){//计算ax+by+c=0和dx+ey+f=0的交点 
    o.x=(c*e-b*f)/(b*d-a*e);
    o.y=(a*f-c*d)/(b*d-a*e);
}
int main(){
    freopen("data.in","r",stdin);//
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
    for(int i=1;i<=n;i++)
        swap(p[rand()%n+1],p[rand()%n+1]);
    o=p[1];
    r=0;
    for(int i=2;i<=n;i++){
        if(dis(o,p[i])>r+eps){//A
            o=p[i];
            r=0;
            for(int j=1;j<i;j++){
                if(dis(o,p[j])>r+eps){//B
                    o.x=(p[i].x+p[j].x)/2;
                    o.y=(p[i].y+p[j].y)/2;
                    r=dis(o,p[j]);
                    for(int k=1;k<j;k++){
                        if(dis(o,p[k])>r+eps){//C，其中cal传入的是ij连线的中垂线方程和ik连线的中垂线方程，这两条直线的交点即为圆心 
                            cal(p[j].x-p[i].x,p[j].y-p[i].y,-(sqr(p[j].y)+sqr(p[j].x)-sqr(p[i].x)-sqr(p[i].y))/2,
                                p[k].x-p[i].x,p[k].y-p[i].y,-(sqr(p[k].y)+sqr(p[k].x)-sqr(p[i].x)-sqr(p[i].y))/2);
                            r=dis(o,p[k]);
                        }
                    }
                }
            }
        }
    }
    printf("%.6lf\n%.6lf %.6lf\n",r,o.x,o.y);
    return 0;
}