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

很直观的一个gauss题；

用的是以前用过的一个模板；

#include<cstdio>

#include<algorithm>

#include<cmath>

#include<cstring>

#define maxn 12

#define eps 0.00001

using namespace std;

double matrix[][];

double ans[];

void exchange_col(int p1,int p2,int n)

{

    double t;

    int i;

    for(int i=; i<=n; i++)

        swap(matrix[p1][i],matrix[p2][i]);

}

bool gauss(int n)

{

    int i,j,k;

    int p;

    double r;

    for(i=; i<n-; i++)

    {

        p=i;

        for(j=i+; j<n; j++)

            if(fabs(matrix[j][i])>fabs(matrix[p][i]))

                p=j;

        if(p!=i)

            exchange_col(i,p,n);

        if(matrix[i][i]==) return false;

        for(j=i+; j<n; j++)

        {

            r=matrix[j][i]/matrix[i][i];

            for(k=i; k<=n; k++)

                matrix[j][k]-=r*matrix[i][k];

        }

    }

    for(i=n-; i>=; i--)

    {

        ans[i]=matrix[i][n];

        for(j=n-; j>i; j--)

            ans[i]-=matrix[i][j]*ans[j];

        if(fabs(matrix[i][i])<eps)return false;

        ans[i]/=matrix[i][i];

    }

    return true;

}

double map[][];

double fang[][];

int main()

{

    int n;

    scanf("%d",&n);

    for(int i=;i<=n;i++)

    {

        for(int j=;j<n;j++)

        {

            scanf("%lf",&map[i][j]);

            if(i>)

            {

                matrix[i-][j]=*(map[i][j]-map[i-][j]);

                fang[i-][j]=map[i][j]*map[i][j]-map[i-][j]*map[i-][j];

            }

        }

    }

    for(int i=;i<n;i++)

    {

        double ret=;

        for(int j=;j<n;j++)

            ret+=fang[i][j];

        matrix[i][n]=ret;

    }

    gauss(n);

    for(int i=;i<n-;i++)

        printf("%.3lf ",ans[i]);

    printf("%.3lf",ans[n-]);

    return ;

}