【JSOI 2008】 球形空间产生器

【题目链接】

https://www.lydsy.com/JudgeOnline/problem.php?id=1013

【算法】

高斯消元

【代码】

#include<bits/stdc++.h>
using namespace std;
const double eps = 1e-;

int i,j,k,n;
double a[][],c[][];
double b[];
double rate;

int main()
{

    scanf("%d",&n);
    for (i = ; i <= n + ; i++)
    {
        for (j = ; j <= n; j++)
        {
            scanf("%lf",&a[i][j]);
        }
    }
    for (i = ; i <= n; i++)
    {
        for (j = ; j <= n; j++)
        {
            c[i][j] =  * (a[i+][j] - a[i][j]);
            b[i] += a[i+][j] * a[i+][j] - a[i][j] * a[i][j];
        }
    }
    for (i = ; i <= n; i++)
    {
        for (j = i; j <= n; j++)
        {
            if(fabs(c[j][i]) > eps)
            {
                for (k = ; k <= n; k++) swap(c[i][k],c[j][k]);
                swap(b[i],b[j]);
            }
        }
        for (j = ; j <= n; j++)
        {
            if (i == j) continue;
            rate = c[j][i] / c[i][i];
            for (k = i; k <= n; k++) c[j][k] -= c[i][k] * rate;
            b[j] -= b[i] * rate;
        }
    }
    for (i = ; i < n; i++) printf("%.3lf ",b[i]/c[i][i]);
    printf("%.3lf\n",b[n]/c[n][n]);

    return ;
}