[JSOI2008] bzoj1013 球星空间产生器sphere - 高斯消元

用一点小学数学知识化简一下就是一个高斯消元。
代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#define N 20
#define gabs(x) (x>0?x:-x)
using namespace std;
double a[N][N],b[N],c[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]),
            c[i]+=a[i][j]*a[i][j];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            a[i][j]=2*(a[i+1][j]-a[i][j]);
    for(int i=1;i<=n;i++) b[i]=c[i+1]-c[i];
    for(int i=1;i<=n;i++)
    {
        int t=i;double x,y;
        for(int j=i;j<=n;j++)
            if(gabs(a[i][j])>gabs(a[t][j])) t=i;
        swap(a[i],a[t]),swap(b[i],b[t]),x=a[i][i];
        for(int j=i;j<=n;j++) a[i][j]/=x;b[i]/=x;
        for(int j=1;j<=n;j++)
            if(j^i)
            {
                y=a[j][i],b[j]-=y*b[i];
                for(int k=i;k<=n;k++)
                    a[j][k]-=y*a[i][k];
            }
    }
    for(int i=1;i<=n;i++)
        if(i<n) printf("%.3f ",b[i]);
        else printf("%.3f",b[i]);
    return 0;
}