18 南京 D

最新推荐文章于 2024-09-10 22:51:14 发布

weixin_30868855

最新推荐文章于 2024-09-10 22:51:14 发布

阅读量75

点赞数

原文链接：http://www.cnblogs.com/MXang/p/10550957.html

版权

裸的最小球覆盖。

坐标范围大一些所以我们把初始的温度也设置的大一些。

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 typedef double db;
 5 const db INF = 1e100;
 6 const db delta = 0.98;
 7 const db eps = 1e-8;
 8 struct p3{
 9     db x,y,z;
10     p3 operator +(const  p3&k1)const { return p3{x+k1.x,y+k1.y,z+k1.z};}
11     p3 operator -(const p3&k1)const { return p3{x-k1.x,y-k1.y,z-k1.z};}
12     p3 operator *(const db k1)const { return p3{x*k1,y*k1,z*k1};}
13     p3 operator /(const db k1)const { return p3{x/k1,y/k1,z/k1};}
14     db abs(){ return sqrt(x*x+y*y+z*z);}
15     db dis(p3 k1){ return (*this-k1).abs();}
16 };
17 int n;
18 p3 p[105];
19 int main(){
20     scanf("%d",&n);
21     for(int i=1;i<=n;i++){
22         scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z);
23     }
24     db T = 100000;
25     db ans = INF;
26     p3 x = {0,0,0};
27     while (T>eps){
28         db mx =0;int id=-1;
29         for(int i=1;i<=n;i++){
30             if((p[i]-x).abs()>mx){
31                 mx=(p[i]-x).abs();
32                 id=i;
33             }
34         }
35         ans = min(ans,mx);
36         x=x+(p[id]-x)/mx*T;
37         T*=delta;
38     }
39     printf("%.11f\n",ans);
40 }