如果单按照距离相等的话既是高次也没有半径，所以因为给了 (n+1) 组点就想到两两做差。
假如一组点是 ({a_i}) 一组是 ({b_i})，我们能轻易地得出
[sum_{i=1}^n(x_i-a_i)^2=sum_{i=1}^n(x_i-b_i)^2 Rightarrow sum_{i=1}^n 2(b_i-a_i)x_i=sum_{i=1}^n(b_i^2-a_i^2)
]
这是一个 (n) 元一次线性方程。我们能搞出 (n) 组这样的方程，高斯消元解之。
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n;
double tmp1[15], tmp2[15], a[15][15];
void gauss(){
for(int i=1; i<=n; i++){
int maxi=i;
for(int j=i+1; j<=n; j++)
if(fabs(a[j][i])>fabs(a[i][i]))
maxi = j;
double now=a[maxi][i];
swap(a[maxi], a[i]);
for(int j=i; j<=n+1; j++)
a[i][j] /= now;
for(int j=i+1; j<=n; j++){
now = a[j][i];
for(int k=i; k<=n+1; k++)
a[j][k] -= a[i][k] * now;
}
}
for(int i=n; i>=1; i--)
for(int j=1; j<=i-1; j++){
a[j][n+1] -= a[j][i] * a[i][n+1];
a[j][i] = 0;
}
}
int main(){
cin>>n;
for(int i=1; i<=n; i++) scanf("%lf", &tmp1[i]);
for(int i=1; i<=n; i++){
for(int j=1; j<=n; j++)
scanf("%lf", &tmp2[j]);
for(int j=1; j<=n; j++){
a[i][j] = 2 * (tmp2[j] - tmp1[j]);
a[i][n+1] += tmp2[j] * tmp2[j] - tmp1[j] * tmp1[j];
}
swap(tmp1, tmp2);
}
gauss();
for(int i=1; i<=n; i++)
printf("%.3f ", a[i][n+1]);
printf("
");
return 0;
}