2019icpc南京 D - Country Meow （三分套三分套三分）

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;
}