[JSOI2008]球形空间产生器

#include<cstdio>
#include<cstring>
#include<algorithm>
const int maxn = 17;

int n;
double loc[maxn][maxn];
double st[maxn],x[maxn];
void gauss() {
    int t;
    for(int i=1;i<=n;++i) {
        t=i;
        for(int j=i+1;j<=n;++j) if(loc[j][i]>loc[t][i]) t=j;
        if(t!=i) for(int j=i;j<=n+1;++j) std::swap(loc[t][j],loc[i][j]);
        for(int j=i+1;j<=n;++j) {
            double tmp=loc[i][i]/loc[j][i];
            for(int k=i+1;k<=n+1;++k) loc[j][k]=loc[i][k]-tmp*loc[j][k];
        }
    }
    for(int i=n;i>=1;--i) { 
        for(int j=i+1;j<=n;++j) loc[i][n+1]-=x[j]*loc[i][j];
        x[i]=loc[i][n+1]/loc[i][i];
    }
}
int main() {
    scanf("%d",&n);
    for(int i=1;i<=n;++i)scanf("%lf",&st[i]);
    for(int i=1;i<=n;++i) 
        for(int j=1;j<=n;++j) {
            double tmp;
            scanf("%lf",&tmp);
            loc[i][j]=2*(st[j]-tmp);
            loc[i][n+1]+=st[j]*st[j]-tmp*tmp;
        }
    gauss();
    for(int i=1;i<=n;++i) printf("%.3lf ",x[i]);
    return 0;
}