10 node p[200000];
11 double r;
12 node O;
13 double dist(node a,node b)
14 {
15 return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) );
16 }
17 void calc(double a,double b,double c,double d,double e,double f) //給出兩條直線ax+by+c=0,dx+ey+f=0 求交點(diǎn)
18 { //注意到三角形里兩條中垂線不可能平行,所以不會(huì)產(chǎn)生除0錯(cuò)誤
19 O.y=(c*d-f*a)/(b*d-e*a);
20 O.x=(c*e-f*b)/(a*e-b*d);
21 }
22 int main()
23 {
24 freopen("HYOJ1337.in","r",stdin);
25 freopen("HYOJ1337.out","w",stdout);
26 scanf("%d",&n);
27 for (int i=1;i<=n;++i) scanf("%lf%lf",&p[i].x,&p[i].y);
28 O=p[1];r=0; //初始C1
29
30 for (int i=2;i<=n;++i) //A
31 if (dist(O,p[i])>r+1e-6)
32 {
33 O=p[i];r=0;
34 for (int j=1;j<=i-1;++j) //B
35 if (dist(O,p[j])>r+1e-6)
36 {
37 O.x=(p[i].x+p[j].x)/2;O.y=(p[i].y+p[j].y)/2;r=dist(O,p[j]);
38 for (int k=1;k<=j-1;++k) //C
39 if (dist(O,p[k])>r+1e-6)
40 {
41 calc(p[j].x-p[i].x,p[j].y-p[i].y,(p[j].x*p[j].x+p[j].y*p[j].y-p[i].x*p[i].x-p[i].y*p[i].y)/2,
42 p[k].x-p[i].x,p[k].y-p[i].y,(p[k].x*p[k].x+p[k].y*p[k].y-p[i].x*p[i].x-p[i].y*p[i].y)/2);
43 r=dist(O,p[k]);
44 }
45 }
46 }
47 printf("%.3lf\n",r);
48 return 0;
49 }
相關(guān)推薦:2011計(jì)算機(jī)等級(jí)考試二級(jí)C輔導(dǎo)實(shí)例編程匯總
計(jì)算機(jī)等級(jí)考試二級(jí)C語(yǔ)言歷年真題匯總(2005-2010)
2011年計(jì)算機(jī)等級(jí)考試二級(jí)C語(yǔ)言常見(jiàn)問(wèn)題匯總
計(jì)算機(jī)等級(jí)考試二級(jí)C語(yǔ)言常見(jiàn)知識(shí)點(diǎn)總結(jié)
北京 | 天津 | 上海 | 江蘇 | 山東 |
安徽 | 浙江 | 江西 | 福建 | 深圳 |
廣東 | 河北 | 湖南 | 廣西 | 河南 |
海南 | 湖北 | 四川 | 重慶 | 云南 |
貴州 | 西藏 | 新疆 | 陜西 | 山西 |
寧夏 | 甘肅 | 青海 | 遼寧 | 吉林 |
黑龍江 | 內(nèi)蒙古 |