tjut 3007

#include<iostream>  
#include<cmath>  
#include<cstdio>  
#include<algorithm>   
using namespace std;  
const double eps=1e-8;  
struct Point{  
    double x,y;   
}p[505];  
double dis(const Point &a,const Point &b)  
{  
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));   
}   
Point circumcenter(const Point &a,const Point &b,const Point &c)  
{ //返回三角形的外心   
    Point ret;   
    double a1=b.x-a.x,b1=b.y-a.y,c1=(a1*a1+b1*b1)/2;  
    double a2=c.x-a.x,b2=c.y-a.y,c2=(a2*a2+b2*b2)/2;  
    double d=a1*b2-a2*b1;  
    ret.x=a.x+(c1*b2-c2*b1)/d;  
    ret.y=a.y+(a1*c2-a2*c1)/d;  
    return ret;   
}   
void min_cover_circle(Point *p,int n,Point &c,double &r){ //c为圆心,r为半径   
    random_shuffle(p,p+n); //   
    c=p[0]; r=0;  
    for(int i=1;i<n;i++)  
    {  
        if(dis(p[i],c)>r+eps)  //第一个点  
        {   
            c=p[i]; r=0;  
            for(int j=0;j<i;j++)  
                if(dis(p[j],c)>r+eps) //第二个点  
                {  
                    c.x=(p[i].x+p[j].x)/2;  
                    c.y=(p[i].y+p[j].y)/2;  
                    r=dis(p[j],c);  
                    for(int k=0;k<j;k++)  
                        if(dis(p[k],c)>r+eps) //第三个点  
                        {//求外接圆圆心,三点必不共线   
                            c=circumcenter(p[i],p[j],p[k]);   
                            r=dis(p[i],c);   
                        }   
                }    
        }      
    }   
}   
int main(){  
    int n;  
    Point c;  
    double r;   
    while(scanf("%d",&n)==1 && n){  
        for(int i=0;i<n;i++)  
            scanf("%lf%lf",&p[i].x,&p[i].y);  
        min_cover_circle(p,n,c,r);                      
        printf("%.2lf %.2lf %.2lf\n",c.x,c.y,r);                     
    }   
    return 0;  
}  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值