zoj 1450 Minimal Circle

题意：给定你一些点，求能够覆盖掉所有点的最小的圆。

方法：最小圆覆盖，随机增量法

代码：

#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;
}