1013: [JSOI2008]球形空间产生器sphere

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=15;

int n;
double a[N][N];

int main()
{
    scanf("%d",&n);
    ++n;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<n;++j)
        {
            scanf("%lf",&a[i][j]);
            a[i][n+1]+=a[i][j]*a[i][j];
            a[i][j]*=2;
        }
        a[i][n]=1;
    }
//    for(int i=1;i<=n;++i)
//    {
//        for(int j=1;j<=n+1;++j)
//            printf("%.2lf ",a[i][j]);
//        puts("");
//    }
    for(int i=1;i<=n;++i)
    {
        int pos=i;
        for(int j=i;j<=n;++j)
        {
            if(fabs(a[pos][i])<fabs(a[j][i]))
                pos=j;
        }
        if(pos!=i)
            swap(a[i],a[pos]);
        for(int j=i+1;j<=n;++j)
        {
            double ratio=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;++k)
            {
                a[j][k]-=a[i][k]*ratio;
            }
        }
    }
    for(int i=n;i;--i)
    {
        for(int j=i+1;j<=n;++j)
        {
            a[i][n+1]-=a[i][j]*a[j][n+1];
        }
        a[i][n+1]/=a[i][i];
    }
    for(int i=1;i<n;++i)
        printf("%.3lf ",a[i][n+1]);
    return 0;
}