Codeforces Gym 101981D

题目链接

题意：给出n个点的坐标，求一个最小球的半径，使其能覆盖所有的点。

分析：这个题跟poj-2069是一样的，把那道题搞懂这个题也就懂了。这个题用模拟退火做的话唯一麻烦的是调试参数......QAQ

真的贼麻烦....

代码如下：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=200;
struct Node{
	double x,y,z;
}p[maxn];
Node ans;
double avex,avey,avez;
const double eps=1e-8;
int N;
double dis(Node a,Node b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z));
}
void SA(){
	double T=10050;
	double ret=1e30;
	while(T>eps){
		int id=1;
		double max_dis=dis(ans,p[1]);
		for(int i=2;i<=N;i++){
			double tem=dis(ans,p[i]);
			if(tem>max_dis){
				max_dis=tem;
				id=i;
			}
		}
		ret=min(ret,max_dis);
		ans.x+=(p[id].x-ans.x)*T/max_dis;
		ans.y+=(p[id].y-ans.y)/max_dis*T;
		ans.z+=(p[id].z-ans.z)/max_dis*T;
		T*=0.98;
	}
	printf("%.15f\n",ret);
}
int main()
{
	scanf("%d",&N);
	for(int i=1;i<=N;i++){
		scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].z);
		avex+=p[i].x;avey+=p[i].y;avez+=p[i].z;
	}
	avex/=N;avey/=N;avez/=N;
	ans.x=avex;ans.y=avey;ans.z=avez;
	SA();
	return 0;
}

 

<附>我在调试的时候还遇到了一个玄学的bug。。。

如果找最优解取min是这样写的：

ret=min(max_dis,ret);

会神奇的WA3.... QAQ....