bzoj1013: [JSOI2008]球形空间产生器sphere （高斯消元）

题意

已知 $n+1$ 个点的坐标，求球心坐标

---

题解

我们可以通过 $2-n$ 个点，与第一个点列出方程组
即 $(x_{i,1}-a_1)^2+(x_{i,2}-a_2)^2+\cdots+(x_{i,n}-a_n)^2=(x_{1,1}-a_1)^2+(x_{1,2}-a_2)^2+\cdots+(x_{1,n}-a_n)^2$
解出 $2(x_{1,1}-x_{i,1})a_1+2(x_{2,1}-x_{i,2})a_2+\cdots+2(x_{1,n}-x_{i,n})a_n=x_{1,1}^2+x_{1,2}^2+\cdots+x_{1,n}-x_{i,1}^2-x_{i,2}^2-x_{n,2}^2$
这样可以列出 n n <script type="math/tex" id="MathJax-Element-12">n</script> 个方程，高斯消元解方程就可以啦

---

代码

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define eps 1e-6
double a[22][22],f[22];
int n;
void gauss(){
    int now=1;
    for(int i=1;i<=n+1;i++){
        int to=now;for(int j=now+1;j<=n;j++) if(fabs(a[i][to])<fabs(a[i][j])) to=j;
        if(fabs(a[i][to])<eps)continue;
        if(to!=now) for(int j=i;j<=n+1;j++) swap(a[now][j],a[to][j]); 
        for(int j=n+1;j>=i;j--) a[now][j]/=a[now][i];
        for(int j=1;j<=n;j++)
            if(j!=now)
                for(int k=n+1;k>=i;k--) a[j][k]-=a[j][i]*a[now][k];
        now++;
    }
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lf",&f[i]);
    for(int i=1;i<=n;i++){
        a[i][n+1]=0;
        for(int j=1;j<=n;j++){
            double x;scanf("%lf",&x);
            a[i][j]=2*(x-f[j]);a[i][n+1]+=x*x-f[j]*f[j];
        }
    }
    gauss();
    printf("%.3lf",a[1][n+1]);for(int i=2;i<=n;i++) printf(" %.3lf",a[i][n+1]);
    return 0;
}