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

bzoj1013: [JSOI2008]球形空间产生器sphere
想了想既然写了就发上来吧。

---

题解

展开距离公式变成n元一次方程组，直接高斯消元即可，保证有解。

---

奇怪的代码

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n,t;
double e[12][12],ans[11];
void rev(double a[],double b[])
{for(int i=0;i<12;i++)swap(a[i],b[i]);}
void sub(double a[],double b[])
{for(int i=0;i<12;i++)a[i]-=b[i];}
void mul(double a[],double k)
{for(int i=0;i<12;i++)a[i]*=k;}
bool zer(double k)
{return (k+1e-8>0)&&(k-1e-8<0);}
int main()
{
    scanf("%d",&n);
    for(int j=0;j<n;j++)
    scanf("%lf",&e[0][j]),e[0][n]+=e[0][j]*e[0][j],e[0][j]*=2;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<n;j++)
        scanf("%lf",&e[i][j]),e[i][n]+=e[i][j]*e[i][j],e[i][j]*=2;
        for(int j=0;j<=n;j++)
        e[i][j]-=e[0][j];
    }
    for(int i=0;i<n-1;i++)
    {
        t=i;
        while(zer(e[++t][i])&&t<=n);
        if(t>n)continue;
        rev(e[i+1],e[t]);
        for(int j=i+2;j<=n;j++)
        if(!zer(e[j][i]))
        mul(e[j],e[i+1][i]/e[j][i]),sub(e[j],e[i+1]);
    }
    ans[n-1]=e[n][n]/e[n][n-1];
    for(int i=n-1;i>0;i--)
    {
        for(int j=i;j<n;j++)
        e[i][n]-=ans[j]*e[i][j];
        ans[i-1]=e[i][n]/e[i][i-1];
    }
    printf("%.3lf",ans[0]);
    for(int i=1;i<n;i++)
    printf(" %.3lf",ans[i]);
    printf("\n");
}