之前自己写过的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;
}