傳送門

Analysis

第一次寫高斯消元，式子倒是很快就推出來了
然而，高斯消元板子寫掛了……注意最後一項啊！！！！

隨便推一下：設球心的座標為 $(x_1,x_2$

Code

#include<bits/stdc++.h>
using namespace std;
int n;
double g[20][20],ans[20],f[20][20];
void gauss(){
	int i,j,k;
	for(i=1;i<=n;++i){
		k=i;
		for(j=i+1;j<=n;++j) if(fabs(g[k][i])<fabs(g[j][i])) k=j;
		for(j=1;j<=n+1;++j) swap(g[k][j],g[i][j]);//n+1
		for(j=i+1;j<=n;++j)
			for(k=i+1;k<=n+1;++k)//n+1
				g[j][k]-=g[j][i]*g[i][k]/g[i][i];
	}
	for(i=n;i>=1;--i){
		for(j=i+1;j<=n;++j)
			g[i][n+1]-=ans[j]*g[i][j];
		ans[i]=g[i][n+1]/g[i][i];
	}
}
int main(){
	scanf("%d",&n);
	int i,j,k;
	for(i=1;i<=n;++i) scanf("%lf",&f[0][i]);
	for(i=1;i<=n;++i){
		for(j=1;j<=n;++j)
		{
			scanf("%lf",&f[i][j]);
			g[i][n+1]+=(f[i][j]*f[i][j])-(f[i-1][j]*f[i-1][j]);
			g[i][j]=2.0*(f[i][j]-f[i-1][j]);
		}
	}
	gauss();
	for(i=1;i<=n;++i) printf("%.3lf ",ans[i]);
	return 0;
}