Description
给出N个点,让你画一个最小的包含所有点的圆。
【题目分析】
最小园覆盖问题,就是选择边界以外的点不断加入,然后把之前的点维护一下就可以了。代码基本上是抄的,尤其是求三个圆的确定的最小覆盖的圆心那一段不是很懂。
【代码】
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const double eps=1e-8;
const int mxn=100000;
int n;
struct point{
double x,y;
friend point operator +(const point a,const point b){return (point){a.x+b.x,a.y+b.y};}
friend point operator -(const point a,const point b){return (point){a.x-b.x,a.y-b.y};}
friend point operator /(const point a,double b){return (point){a.x/b,a.y/b};}
}p[mxn];
inline double dis(point a,point b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
point center(point a,point b,point c){
double a1,a2,b1,b2,c1,c2;
point ans;
a1=2*(b.x-a.x);b1=2*(b.y-a.y);c1=(b.x*b.x)-(a.x*a.x)+(b.y*b.y)-(a.y*a.y);
a2=2*(c.x-a.x);b2=2*(c.y-a.y);c2=(c.x*c.x)-(a.x*a.x)+(c.y*c.y)-(a.y*a.y);
if(fabs(a1)<eps) ans.y=c1/b1,ans.x=(c2-ans.y*b2)/a2;
else if(fabs(b1)<eps) ans.x=c1/a1,ans.y=(c2-ans.x*a2)/b2;
else ans.x=(c2*b1-c1*b2)/(a2*b1-a1*b2),ans.y=(c2*a1-c1*a2)/(b2*a1-b1*a2);
return ans;
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
random_shuffle(p+1,p+n+1);
point t=p[1];
double r=0.0;
for(int i=2;i<=n;i++)
if(dis(t,p[i])>r+eps){
t=(p[i]+p[1])/2,r=dis(p[i],t);
for(int j=2;j<i;j++)
if(dis(t,p[j])>r+eps){
t=(p[i]+p[j])/2,r=dis(t,p[i]);
for(int k=1;k<j;k++)
if(dis(p[k],t)>r+eps)
t=center(p[i],p[j],p[k]),r=dis(p[i],t);
}
}
printf("%.10lf\n%.10lf %.10lf",r,t.x,t.y);
return 0;
}