球形空间产生器（高斯消元）

#include<bits/stdc++.h>
using namespace std;
const int N=20;
double a[N][N],b[N],c[N][N];
int main(){
	int n;scanf("%d",&n);
	for(int i=1;i<=n+1;i++)
		for(int j=1;j<=n;j++) scanf("%lf",&a[i][j]);
	//构建增广矩阵 
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++){
			c[i][j]=2*(a[i][j]-a[i+1][j]);
			b[i]+=a[i][j]*a[i][j]-a[i+1][j]*a[i+1][j];
		}
	//高斯消元
	for(int i=1;i<=n;i++){
		for(int j=i;j<=n;j++){
			if(fabs(c[j][i])>1e-8){
				for(int k=1;k<=n;k++) swap(c[i][k],c[j][k]);
				swap(b[i],b[j]);
				break;
			}
		}
		for(int j=1;j<=n;j++){
			if(i==j) continue;
			double rate=c[j][i]/c[i][i];
			for(int k=i;k<=n;k++) c[j][k]-=c[i][k]*rate;
			b[j]-=b[i]*rate;
		}
	}
	for(int i=1;i<n;i++) printf("%.3lf ",b[i]/c[i][i]);
	printf("%.3lf\n",b[n]/c[n][n]); 
	return 0;
}