input:
3
0 0 0
3 0 0
0 4 0
output:
4
0 0 0
1 0 0
0 1 0
0 0 1
题目大意:给我们n个三维坐标点,要求我们从三维空间中选择一个点,使得这个点到其他所有点的最大距离最小。
解题思路:三分套三分再套个三分,每次三分选择一个坐标,最终选择出一个(x,y,z),用这个坐标和其他所有点求最大距离即可,至于为什么用三分,我们单看一维,在一维坐标系中选择一个点其最大距离满足一个单峰函数,所以可以用三分。
上代码:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int N=110;
struct node{
int x;
int y;
int z;
}p[N];
int n;
double check(double x,double y,double z)
{
double ans=-1;
for(int i=1;i<=n;i++)
{
double temp=sqrt((p[i].x-x)*(p[i].x-x)+(p[i].y-y)*(p[i].y-y)+(p[i].z-z)*(p[i].z-z));
ans=max(ans,temp);
}
return ans;
}
double check2(double x,double y)
{
double l=100000,r=-100000;
while(fabs(r-l)>1e-9)
{
double mid1=(r+l)/2,mid2=(mid1+r)/2;
if(check(x,y,mid1)>check(x,y,mid2))
l=mid1;
else
r=mid2;
}
return check(x,y,l);
}
double check1(double x)
{
double l=100000,r=-100000;
while(fabs(r-l)>1e-9)
{
double mid1=(r+l)/2,mid2=(mid1+r)/2;
if(check2(x,mid1)>check2(x,mid2))
l=mid1;
else
r=mid2;
}
return check2(x,l);
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
cin>>p[i].x>>p[i].y>>p[i].z;
double l=100000,r=-100000;
while(fabs(r-l)>1e-9)
{
double mid1=(l+r)/2,mid2=(mid1+r)/2;
if(check1(mid1)>check1(mid2))
l=mid1;
else
r=mid2;
}
printf("%.15lf\n",check1(l));
return 0;
}