题目描述:
给出一组点的坐标,求出能够覆盖掉他们的最小圆的坐标及半径.
算法分析:
解决这个问题有一种算法,随机化增量算法,在随机数据下可以O(N)的解决这个问题.
具体实现方法:
1.先将所有的点随机化处理
2.按顺序把点一个一个的加入(一步一步的求前i个点的最小覆盖圆),每加入一个点就进入步骤3
3.判断当前点是否在当前的最下覆盖圆内,如果不在进入4,在进入2
4.将圆心设为当前点 i 半径为 0,然后把前 i-1个点加入进这个最小覆盖圆中,每加入一个点进入5
5.如果发现当前j号点在当前的最小圆的外面,那么说明点j也一定在前j个点(包括i)的最小覆盖圆边界上,我们转到6来再进一步确定这个圆,否则前j个点(包括i)的最小覆盖圆与前i-1个点(包括i)的最小覆盖圆是一样的,则不需要更新,返回
6.此时已经确认点i,j一定在前j个点(包括i)的最小覆盖圆的边界上了,那么我们可以把当前圆的圆心设为第i个点与第j的点连线的中点,半径为到这两个点的距离(就是找一个覆盖这两个点的最小圆),然后重新把前j-1个点加入这个圆中(还是类似上面的步骤,只不过这次我们提前确定了两个点在圆上,目的是求出包含点i,j的前k个点的最小覆盖圆),每加入一个点就进入7
7.如果发现当前k号点在当前的最小圆的外面,那么说明点k也一定在前k个点(包括i,j)的最小覆盖圆边界上,我们不需要再进一步确定这个圆了(因为三个点能确定一个圆!),直接求出这三点共圆,否则前k个点(包括i,j)的最小覆盖圆与前k-1个点(包括i,j)的最小覆盖圆是一样的,则不需要更新。
题目链接:
代码如下:
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <cstdlib>
const int maxm=1005;
const double eps=1e-9;
struct Point{
double x,y;
int tim;
};
int n;
Point a[maxm],o;
double r;
double distance(Point p1,Point p2)
{
return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
void makecir(Point p1,Point p2,Point p3)//三点确定一个圆
{
double x1=2*(p2.x-p1.x),y1=2*(p2.y-p1.y);
double c1=p2.x*p2.x+p2.y*p2.y-p1.x*p1.x-p1.y*p1.y;
double x2=2*(p3.x-p1.x),y2=2*(p3.y-p1.y);
double c2=p3.x*p3.x+p3.y*p3.y-p1.x*p1.x-p1.y*p1.y;
o.x=(y1*c2-y2*c1)/(y1*x2-y2*x1);
o.y=(x2*c1-x1*c2)/(y1*x2-y2*x1);
r=distance(o,p1);
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
if(!n) break;
for(int i=1;i<=n;i++)
scanf("%lf%lf",&a[i].x,&a[i].y);
std::random_shuffle(a+1,a+n+1);
o=a[1],r=0;
for(int i=2;i<=n;i++)
if(distance(o,a[i])>r+eps)
{
o=a[i],r=0;
for(int j=1;j<=i-1;j++)
if(distance(o,a[j])>r+eps)
{
o.x=(a[i].x+a[j].x)/2;
o.y=(a[i].y+a[j].y)/2;
r=distance(o,a[j]);
for(int k=1;k<=j-1;k++)
if(distance(o,a[k])>r+eps)
makecir(a[i],a[j],a[k]);
}
}
printf("%.2lf %.2lf %.2lf\n",o.x,o.y,r);
}
return 0;
}