P4035 [JSOI2008]球形空间产生器

之前自己写过的PDF版的题解：

AC Code：

#include<cstdio>
#include<iostream>
#define ri register int
using namespace std;

const int MAXN=15;
int n;
double fst[MAXN][MAXN],k[MAXN],a[MAXN][MAXN],b[MAXN],x[MAXN];

void Gauss()
{
    for(ri i=1;i<=n;i++)
        for(ri j=i+1;j<=n;j++)
        {
            double chu=a[j][i]/a[i][i];
            for(ri k=i;k<=n;k++)	a[j][k]-=a[i][k]*chu;
            b[j]-=chu*b[i];
        }
    for(ri i=n;i>=1;i--)
    {
        double ans=0;
        for(ri j=i+1;j<=n;j++)	 ans+=x[j]*a[i][j];
        x[i]=(b[i]-ans)/a[i][i];
    }
}

int main()
{
    scanf("%d",&n);
    for(ri i=1;i<=n+1;i++)
        for(ri j=1;j<=n;j++)	scanf("%lf",&fst[i][j]);
    for(ri i=1;i<=n+1;i++)
        for(ri j=1;j<=n;j++)	k[i]+=fst[i][j]*fst[i][j];
    for(ri i=1;i<=n;i++)
    {
        for(ri j=1;j<=n;j++)	a[i][j]=(fst[i][j]-fst[i+1][j])*2;
        b[i]=k[i]-k[i+1];
    }
    Gauss();
    for(ri i=1;i<=n;i++)	printf("%.3f ",x[i]);
    return 0;
}