记录编号	610886	评测结果	AAAAAAAAAA
题目名称	[JSOI 2008]球形空间产生器sphere	最终得分	100
用户昵称	终焉折枝	是否通过	通过
代码语言	C++	运行时间	0.030 s
提交时间	2026-01-24 09:41:18	内存使用	3.85 MiB

#include<iostream>
#include<cmath>
using namespace std;

const int MAXN = 15;
const double eps = 1e-6;

int n;
double a[MAXN][MAXN];
double p[MAXN][MAXN];

bool Jornh(){
    for(int i = 1;i <= n;i ++){
        int piv = i;
        for(int k = i + 1;k <= n;k ++){
            if(fabs(a[k][i]) > fabs(a[piv][i])){
                piv = k;
            }
        }
        for(int j = 1;j <= n + 1;j ++){
            swap(a[i][j], a[piv][j]);
        }
        if(fabs(a[i][i]) < eps){
            return false;
        }
        for(int k = 1;k <= n;k ++){
            if(k == i) continue;
            double fac = a[k][i] / a[i][i];
            for(int j = i;j <= n + 1;j ++){
                a[k][j] -= fac * a[i][j];
            }
        }
    }
    for(int i = 1;i <= n;i ++){
        a[i][n + 1] /= a[i][i];
    }
    return true;
}

int main(){

    cin >> n;
    for(int i = 1;i <= n + 1;i ++){
        for(int j = 1;j <= n;j ++){
            cin >> p[i][j];
        }
    }

    for(int i = 1;i <= n;i ++){
        double res = 0;
        for(int j = 1;j <= n;j ++){
            a[i][j] = 2.0 * (p[i + 1][j] - p[i][j]);
            res += (p[i + 1][j] * p[i + 1][j] - p[i][j] * p[i][j]);
        }
        a[i][n + 1] = res;
    }

    if(Jornh()){
        for(int i = 1;i <= n;i ++){
            printf("%.3lf ", a[i][n + 1]);
        }
    }
    return 0;
}