╮(╯▽╰)╭,小球球啊,你的心究竟在哪里
很明显,这是一道高斯消元的题目,虽然我们不知道4维、5维……n维的小球究竟长什么样,但是我们还是可以找到他的心的
(^o^)/YES!
他的心对每个点是平等的,到他们的距离都相等,于是就有了方程啦
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>//double gauss
using namespace std;
const int maxn=10+5;
int n;
double coefficient[maxn][maxn];
double spot[maxn];
void gauss(){
for(int i=1;i<=n;i++){
int k;double maxs=-1.0;
/*
for(k=i;k<=n;k++)
if(fabs(coefficient[k][i])>1e-6)//精度问题
break;
*/
for(int j=1;j<=n;j++)
if(fabs(coefficient[j][i])>maxs)
maxs=coefficient[j][i],k=j;//寻找最大的主元
if(i!=k)
for(int s=1;s<=n+1;s++)
swap(coefficient[i][s],coefficient[k][s]);
double t=coefficient[i][i];
for(int s=i+1;s<=n+1;s++)
coefficient[i][s]/=t;
for(k=1;k<=n;k++)
if(k!=i){
t=coefficient[k][i];
for(int s=1;s<=n+1;s++)
coefficient[k][s]-=coefficient[i][s]*t;
}
}
}
int main(){
scanf("%d",&n);
//cin>>n;
for(int i=1;i<=n;i++)
//cin>>spot[i];
scanf("%lf",&spot[i]);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
double x;scanf("%lf",&x);//cin>>x;
coefficient[i][j]=2*(x-spot[j]);
coefficient[i][n+1]+=(x*x)-(spot[j]*spot[j]);
}
gauss();
for(int i=1;i<n;i++)
printf("%.3f ",coefficient[i][n+1]);
printf("%.3f\n",coefficient[n][n+1]);
return 0;
}
关于为什么寻找主元时要找最大主元这个问题是借鉴的一位大神的代码
~(≧▽≦)/~(感谢大神)
理论上只要不为零都可以
取最大的话后面除的时候,结果就比较小,小数点后面的精确位数就更多,精度误差更小一些
by >o< neighthorn