BZOJ1013: [JSOI2008]球形空间产生器sphere 高斯消元

有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在，你被困在了这个n维球体中，你只知道球面上n+1个点的坐标，你需要以最快的速度确定这个n维球体的球心坐标，以便于摧毁这个球形空间产生器。
n<=10
列出n+1个点到圆心的距离表达式，其值都相等。于是用其中n个分别减去最后一个，由平方差公式即可得到n个线性方程，高斯消元求解即可。
高斯消元的做法：
1.枚举消到第几项
2.枚举第一个要消的项不为0的方程，将其提到当前位置
3.用这个方程将它下面的方程这一项全部消掉
4.从下到上求解每一个方程，将其解代入它上面的方程来消掉它对应的项。

#include<cstdio>
#include<algorithm>
using namespace std;
int n;
double a[11][12],p[11][11];
inline void cosmos()
{
    for(int i=1;i<=n;i++)
    {
        int j=i;
        while(a[j][i]==0.0) ++j;
        if(i!=j) swap(a[i],a[j]);
        for(j=i+1;j<=n;++j)
        {
            double r=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;k++)
            a[j][k]-=r*a[i][k];
        }
    }
    for(int i=n;i>=1;i--)
    {
        a[i][n+1]/=a[i][i];
        for(int j=i-1;j>=1;j--)
        a[j][n+1]-=a[j][i]*a[i][n+1];
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=0;i<=n;i++)
    for(int j=1;j<=n;j++)
    scanf("%lf",&p[i][j]);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    {
        double x1=p[0][j],x2=p[i][j];
        a[i][j]=2*(x2-x1);
        a[i][n+1]-=(x1+x2)*(x1-x2);
    }
    cosmos();
    for(int i=1;i<n;i++)
    printf("%.3lf ",a[i][n+1]);
    printf("%.3lf\n",a[n][n+1]);
}