[JSOI2008]球状空間生成球 - ガウス消去

説明

所与\（N + 1 \）\（\ N-）次元座標、解決\（\ N-） -次元のボールの中心

考え

レビューの質問ガウス消去式から、（X ^ 2、y ^ 2 \）は\ 得るために、線形方程式を排除します

順次最初の選択\（私は\）列の除去、便宜上、上三角行列のバンドに戻し、その後消滅していないが、しかし、完全な排除を取ります

#include <bits/stdc++.h>
using namespace std;

double a[20][20], b[20], c[20][20];
int n;

int main() {
    cin >> n;
    for (int i = 1; i <= n + 1; ++ i)
        for (int j = 1; j <= n; ++ j) scanf("%lf", &a[i][j]);
    for (int i = 1; i <= n; ++ i)
        for (int j = 1; j <= n; ++ j) {
            c[i][j] = 2 * (a[i][j] - a[i + 1][j]);
            b[i] += a[i][j] * a[i][j] - a[i + 1][j] * a[i + 1][j];
        }
    for (int i = 1; i <= n; ++ i) {
        for (int j = i; j <= n; ++ j) {
            if (fabs(c[j][i]) > 1e-8) {
                for (int k = 1; k <= n; ++ k) swap(c[i][k], c[j][k]);
                swap(b[i], b[j]);
                break;
            }
        }
        for (int j = 1; j <= n; ++ j) {
            if (i == j) continue;
            double rate = c[j][i] / c[i][i];
            for (int k = i; k <= n; ++ k) c[j][k] -= c[i][k] * rate;
            b[j] -= b[i] * rate;
        }
    }
    for (int i = 1; i < n; ++ i) printf("%.3f ", b[i] / c[i][i]);
    printf("%.3f\n", b[n] / c[n][n]);
}