HDU 3007 最小覆盖圆

题意：
一些点，求半径最小的圆能覆盖所有点。

题解：
我觉得随机增量特别扯。。

主要还是我弱。。

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

#define N 2222
#define EPS 1e-7

using namespace std;

struct PO
{
    double x,y;
}p[N];

int n;

inline void read()
{
    for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
}

inline int dc(double x)
{
    if(x>EPS) return 1;
    else if(x<EPS) return -1;
    return 0;
}

inline PO operator -(PO a,PO b)
{
    PO c;
    c.x=a.x-b.x; c.y=a.y-b.y;
    return c;
}

inline PO operator +(PO a,PO b)
{
    PO c;
    c.x=a.x+b.x; c.y=a.y+b.y;
    return c;
}

inline double getdis(PO &a,PO &b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

inline double cross(PO &a,PO &b,PO &c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}

inline PO getpoint(PO &a,PO &b,PO &c,PO &d)//直线交点 
{
    PO ans,tp=b-a;
    double k1=cross(a,d,c);
    double k2=cross(b,c,d);
    ans.x=a.x+tp.x*k1/(k1+k2);
    ans.y=a.y+tp.y*k1/(k1+k2);
    return ans;
}

inline PO getcenter(PO a,PO b,PO c)//三角形外接圆圆心 
{
    PO mid1,mid2,t1,t2;
    mid1=b-a; mid1.x*=0.5; mid1.y*=0.5; mid1=mid1+a;
    mid2=c-a; mid2.x*=0.5; mid2.y*=0.5; mid2=mid2+a;
    t1=mid1; t1.x-=b.y-a.y; t1.y+=b.x-a.x;
    t2=mid2; t2.x-=c.y-a.y; t2.y+=c.x-a.x;
    return getpoint(mid1,t1,mid2,t2);
}

inline void mincover(PO &c,double &r)
{
    random_shuffle(p+1,p+1+n);
    c=p[1]; r=0.0;
    for(int i=2;i<=n;i++)
        if(dc(getdis(p[i],c)-r)>0)
        {
            c=p[i]; r=0.0;
            for(int k=1;k<i;k++)
                if(dc(getdis(p[k],c)-r)>0)
                {
                    c.x=(p[i].x+p[k].x)*0.5;
                    c.y=(p[i].y+p[k].y)*0.5;
                    r=getdis(p[k],c);
                    for(int j=1;j<k;j++)
                        if(dc(getdis(p[j],c)-r)>0)
                        {
                            c=getcenter(p[i],p[k],p[j]);
                            r=getdis(p[i],c);
                        }
                }
        }
}

inline void go()
{
    PO c; double r;
    mincover(c,r);
    printf("%.2lf %.2lf %.2lf\n",c.x,c.y,r);
}

int main()
{
    while(scanf("%d",&n),n) read(),go();
    return 0;
}