P4035 [JSOI2008]球形空间产生器

#include<bits/stdc++.h>
using namespace std;
double a[17][17],b[17][17];
int n;
void build(int x,int y){
    for(int i=1;i<=n;i++) a[x][i]=b[y][i]*2-b[x][i]*2,a[x][n+1]+=-(b[x][i]*b[x][i]-b[y][i]*b[y][i]);
    return ;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)
            scanf("%lf",&b[i][j]);
    for(int i=1;i<=n;i++) build(i,i+1);
    int pl;
    for(int i=1;i<=n;i++){
        pl=i;
        while(a[pl][i]==0&&pl<=n) pl++;
        if(pl==n+1){printf("No Solution\n");return 0;} 
        for(int j=1;j<=n+1;j++) swap(a[i][j],a[pl][j]);
        double k=a[i][i];
        for(int j=1;j<=n+1;j++) 
            a[i][j]/=k;
        for(int j=1;j<=n;j++){
            if(j!=i){
                double ki=a[j][i];
                for(int m=1;m<=n+1;m++)
                    a[j][m]=a[j][m]-ki*a[i][m];
            }
        }
    }
    for(int i=1;i<=n;i++) printf("%.3lf ",a[i][n+1]);
    return 0;
}