//题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3007
//求点集最小覆盖圆
//增量法
#include<iostream>
#include<cstdio>
#include<math.h>
#define eps 1e-8
const int maxn=10000;
struct point {double x,y;}cen,p[maxn],po;
double r;
int n;
double distance(point p1,point p2)
{
return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double dist2(point p1,point p2)
{
return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
//直线a1X+ b1Y+c1=0与a2X+b2Y+c2=0的交点
void calc(double a1,double b1,double c1,double a2,double b2,double c2)
{
cen.y=(c1*a2-c2*a1)/(b1*a2-b2*a1);
cen.x=(c1*b2-c2*b1)/(a1*b2-b1*a2);
}
void cir()
{
cen=p[0];r=0;
int i,j,k;
for(i=1;i<n;i++)
if(distance(cen,p[i])-r>eps)
{
cen=p[i];
r=0;
for(j=0;j<=i-1;j++)
if(distance(cen,p[j])-r>eps)
{
cen.x=(p[i].x+p[j].x)/2.0;
cen.y=(p[i].y+p[j].y)/2.0;
r=distance(cen,p[j]);
for(k=0;k<=j-1;k++)
if(distance(cen,p[k])-r>eps)
{
calc(p[j].x-p[i].x,
p[j].y-p[i].y,
(dist2(p[j],po)-dist2(p[i],po))/2.0,
p[k].x-p[i].x,
p[k].y-p[i].y,
(dist2(p[k],po)-dist2(p[i],po))/2.0
);
r=distance(cen,p[k]);
}
}
}
}
int main()
{
po.x=0;
po.y=0;
while(scanf("%d",&n),n)
{
int i;
for(i=0;i<n;i++)
{
scanf("%lf%lf",&p[i].x,&p[i].y);
}
cir();
printf("%.2lf %.2lf %.2lf\n",cen.x,cen.y,r);
}
return 0;
}
分享到:
相关推荐
hdu2215 最简单的最小圆覆盖
HDU的1250,主要是利用高精度加法,但是代码有点繁琐,效率不是很高
杭电ACMhdu1163
HDU1059的代码
hdu1001解题报告
hdu 1574 passed sorce
HDU的一题........HDU DP动态规
hdu2101AC代码
hdu acm 教案 搜索入门 hdu acm 教案 搜索入门
搜索 dfs 解题代码 hdu1241
hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。
hdu acm 教案 动态规划(1) hdu acm 教案 动态规划(1)
hdu 1166线段树代码
自己做的HDU ACM已经AC的题目
ACM HDU题目分类,我自己总结的大概只有十来个吧
HDU最全ac代码
hdu动态规划算法集锦
hdu题目分类
杭电ACM课件2014版之 (HDUACM201403版_06)并查集(最小生成树)
(HDUACM2010版_06)并查集(最小生成树(HDUACM2010版_06)并查集(最小生成树(HDUACM2010版_06)并查集(最小生成树(HDUACM2010版_06)并查集(最小生成树