luogu P4035 [JSOI2008]球形空间产生器

背景：

调参吧。

题意：

在一个$n$为的空间里，我们找到一个球心，使得这$n+1$个点在一个球上（一定在表面）。

思路：

考虑每找到一个点作为圆心时的距离这$n+1$个点的最大距离与最小距离，我们肯定希望这两个的差值尽可能小。那就可以模拟退火了。
$PS$：加油调参，反正我调不出来， 只有$80pts$。

代码：

#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define eps 1e-16
using namespace std;
	int n;
	double now[20],ans[20],a[20][20];
double calc_now()
{
	double sum=0,ma=0,mi=1e60;
	for(int i=1;i<=n+1;i++)
	{
		double dis=0;
		for(int j=1;j<=n;j++)
			dis+=(now[j]-a[i][j])*(now[j]-a[i][j]);
		ma=max(ma,dis);
		mi=min(mi,dis);
	}
	return ma-mi;
}
double calc_ans()
{
	double sum=0,ma=0,mi=1e60;
	for(int i=1;i<=n+1;i++)
	{
		double dis=0;
		for(int j=1;j<=n;j++)
			dis+=(ans[j]-a[i][j])*(ans[j]-a[i][j]);
		ma=max(ma,dis);
		mi=min(mi,dis);
	}
	return ma-mi;
}
void give()
{
	for(int i=1;i<=n;i++)
		ans[i]=now[i];
}
void fire()
{
	double t=5000;
	while(t>eps)
	{
		for(int i=1;i<=n;i++)
			now[i]=ans[i]+(rand()*2-RAND_MAX)*t;
		double delta=calc_now()-calc_ans();
		if(delta<0||exp(-delta/t)*RAND_MAX>rand()) give();
		t*=0.99998;
	}
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n+1;i++)
		for(int j=1;j<=n;j++)
			scanf("%lf",&a[i][j]);
	fire();
	for(int i=1;i<=n;i++)
		printf("%.3lf ",ans[i]);
}