【bzoj1013】球形空间产生器

高斯消元……

看完线代那一节之后感觉真的是……naive！

线代大法好。

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-;
int n;
double a[][],f[];
inline double sqr(double x){return x*x;}
void init(){
    scanf("%d",&n);
    for(int i=;i<=n;i++)scanf("%lf",f+i);
    for(int i=;i<=n;i++)
    for(int j=;j<=n;j++){
        double t;scanf("%lf",&t);
        a[i][j]=*(t-f[j]);
        a[i][n+]+=sqr(t)-sqr(f[j]);
    }
}
bool gauss(){
    int now=,to;double t;
    for(int i=;i<=n;i++){
        for(to=now;to<=n;to++)if(fabs(a[to][i])>eps)break;
        if(to>n)continue;
        if(to!=now)for(int j=;j<=n+;j++)swap(a[to][j],a[now][j]);
        t=a[now][i];
        for(int j=;j<=n+;j++)a[now][j]/=t;
        for(int j=;j<=n;j++)if(j!=now){
            t=a[j][i];
            for(int k=;k<=n+;k++)a[j][k]-=t*a[now][k];
        }
        ++now;
    }
    for(int i=now;i<=n;i++)if(fabs(a[i][n+])>eps)return ;
    return ;
}
int main(){
    init();gauss();
    for(int i=;i<n;i++)printf("%.3lf ",a[i][n+]);
    printf("%.3lf\n",a[n][n+]);
    return ;
}