【bzoj 1013】球形空间产生器sphere

传送门~

解题思路

用距离公式以后式子里有$r^2$，所以相邻两个式子相减消去$r^2$，之后高斯消元。
代码

#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
using namespace std;
double c[15][2],f[15][15];
int n;
void gauss(){
    for(int i=1;i<=n;i++){
        int pre=i;
        for(int j=i+1;j<=n;j++)
        if(f[j][i]>f[pre][i])
        pre=j;  
        for(int j=1;j<=n+1;j++) swap(f[i][j],f[pre][j]);
        for(int j=n+1;j>=i;j--) f[i][j]/=f[i][i];
        for(int j=1;j<=n;j++)
        if(i!=j) for(int k=n+1;k>=i;k--) 
        f[j][k]-=f[i][k]*f[j][i];
    }
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){scanf("%lf",&c[i][0]);c[n+1][0]+=c[i][0]*c[i][0];}
    for(int i=1;i<=n;i++){
        c[n+1][1]=0;
        for(int j=1;j<=n;j++){
            scanf("%lf",&c[j][1]);
            c[n+1][1]+=c[j][1]*c[j][1];
            f[i][j]=2*(c[j][1]-c[j][0]);
        }
        f[i][n+1]=c[n+1][1]-c[n+1][0];
        for(int j=1;j<=n+1;j++) c[j][0]=c[j][1];
    }
    gauss();
    for(int i=1;i<n;i++) printf("%.3lf ",f[i][n+1]);
    printf("%.3lf",f[n][n+1]);
    return 0;
}