【爬山】P4035 [JSOI2008]球形空间产生器

题目描述：

给出n维空间中的n个点，已知它们在同一个n维球面上，求出球心。

操作过程

1. 初始化球心为各个给定点的重心
2. 求出每个已知点到这个球心距离的平均值。
3. 遍历所有已知点。记录一个改变值cans（分开每一维度记录）对于每一个点的距离，如果大于平均值，就是要较劲这个点，相当于把改变值加上差值，否则减去
4. 将我们记录的cans乘上醉酒成度，更新球心，反复重复步骤 2，直到醉酒成度小于某个给定值的时候结束。

 

 

代码

#include<bits/stdc++.h>
using namespace std;
struct pointt
{
	double num[12];
}point[12],ans,cans;
int n;
double dis[12];
double calc_dis(pointt x,pointt y)
{
	double ans=0;
	for(int i=1;i<=n;i++)
		ans+=(x.num[i]-y.num[i])*(x.num[i]-y.num[i]);
	return sqrt(ans);
}
int main()
{
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n+1;i++)
		for(int j=1;j<=n;j++)
		{

			scanf("%lf",&point[i].num[j]);
			ans.num[j]+=point[i].num[j]/(n+1);
		}
	for(double t=10000;t>=0.0001;t*=0.9999)
	{
		double aver=0;
		for(int i=1;i<=n+1;i++)
			dis[i]=calc_dis(ans,point[i]),aver+=dis[i];
		aver/=n+1;
		for(int i=1;i<=n;i++) cans.num[i]=0;
		for(int i=1;i<=n+1;i++)
			for(int j=1;j<=n;j++)
				cans.num[j]+=(dis[i]-aver)*(point[i].num[j]-ans.num[j])/aver;
		for(int i=1;i<=n;i++)
			ans.num[i]+=cans.num[i]*t;
	}
	for(int i=1;i<=n;i++)
		printf("%0.3f ",ans.num[i]);
	return 0;
}