球形空间产生器[高斯消元]

传送门

每一行与第一行联立 , 然后高斯消元

---

#include<bits/stdc++.h>
#define N 15
using namespace std;
int n;
double c[N][N],f[N];
double gauss(){
	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+1;k++) swap(c[i][k],c[j][k]);
			}
		}
		for(int j=1;j<=n;j++){
			if(i==j) continue;
			double now = c[j][i] / c[i][i];
			for(int k=i;k<=n+1;k++) c[j][k] -= c[i][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++){
		for(int j=1;j<=n;j++){
			double x;scanf("%lf",&x);
			c[i][j] = 2 * (x-f[j]);
			c[i][n+1] += x * x - (f[j] * f[j]);
		}
	}
	gauss();
	for(int i=1;i<=n;i++) printf("%0.3lf ",c[i][n+1]/c[i][i]);
	return 0;
}