题意:给定你一些点,求能够覆盖掉所有点的最小的圆。
方法:最小圆覆盖,随机增量法
代码:
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
#define eps 1e-8
struct Point
{
double x, y;
};
struct Circle
{
Point c;
double r;
};
int dcmp(double a)
{
if (fabs(a) < eps)return 0;
return a>0?1:-1;
}
//两点之间距离
double dist(Point a, Point b)
{
return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) );
}
//外接圆
Circle calc(Point p1, Point p2, Point p3)
{
Circle temp;
double a,b,c,d,e,f;
a=p2.x-p1.x;
b=p2.y-p1.y;
c=(p2.x*p2.x+p2.y*p2.y-p1.x*p1.x-p1.y*p1.y)/2;
d=p3.x-p1.x;
e=p3.y-p1.y;
f=(p3.x*p3.x+p3.y*p3.y-p1.x*p1.x-p1.y*p1.y)/2;
temp.c.y=(c*d-f*a)/(b*d-e*a);
temp.c.x=(c*e-f*b)/(a*e-b*d);
temp.r = dist(temp.c, p1);
return temp;
}
Circle Minc(Point *p, int n)
{
Circle min;
min.c = p[0];//初始化圆心
min.r = 0;//初始化半径
for (int i=1; i<n; i++)
{
if ( dcmp(dist(min.c, p[i])-min.r)<=0)
continue;//新加入的点p[i]在圆内
min.c = p[i];//新加入的点p[i]不在圆内
min.r = 0;//以新加入的为基础,半径为0
for (int j=0; j<i; j++)
{
if ( dcmp(dist(min.c, p[j])-min.r)<=0 )
continue;//在p[0...i-1]个点中,点在新圆中
min.c.x = (p[i].x + p[j].x) / 2;//有点不在圆中,但必须让点在圆中
min.c.y = (p[i].y + p[j].y) / 2;//以p[i]和p[j]中点为圆心,建立新圆
min.r = dist(min.c, p[j]);//求出半径
for (int k=0; k<j; k++)
{
if ( dcmp(dist(min.c, p[k])-min.r)<=0 )
continue;//在p[0...j-1]个点中,点在圆中
min = calc(p[i], p[j], p[k]);//不在圆中的点,以p[i],p[j],p[k]确定一个外接圆
min.r = dist(min.c, p[k]);
}
}
}
return min;
}
Circle minc(Point *p, int n);
int main()
{
Point p[505];
Circle O;
int n;
while (~scanf("%d", &n) && n)
{
for (int i=0; i<n; i++)
scanf("%lf %lf", &p[i].x, &p[i].y);
O = minc(p, n);
printf("%.2f %.2f %.2f\n", O.c.x, O.c.y, O.r);
}
return 0;
}