点击这里查看原题
随机增量法的思路:
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;
}