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

高斯消元模板

设一个N维球心（a,b,c,d...）

多个方程左右平方后分别减去第一个等式 对得到的式子高斯消元即可(一定有解

#include<iostream>
#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#include<set>
#include<map>
#include<iomanip>
#include<cstdlib>
#include<cctype>
const int NN=100;
const double eps=1e-10;
using namespace std;

int N,M,T;
double A[NN][NN],f[NN];

void Gauss(int N){
    int i,j,k,now=1;double t;
    for(i=1;i<=N;i++){
        for(j=now;j<=N;j++) if(fabs(A[j][i])>eps) break;
        if(j>N) continue;
        if(j!=now) for(k=1;k<=N+1;k++) swap(A[now][k],A[j][k]);
        t=A[now][i];
        for(k=1;k<=N+1;k++) A[now][k]/=t;
        for(j=1;j<=N;j++){
            if(now==j) continue;
            t=A[j][i];
            for(k=1;k<=N+1;k++) A[j][k]-=t*A[now][k];
        }
        now++;
    }
    //for(i=now;i<=N;i++) if(A[i][N+1]!=0) //无解
    //if(now!=N+1)  //无数解
	
	for (int i=1;i<N;i++) printf("%.3f ",A[i][N+1]);
	printf("%.3f\n",A[N][N+1]); 
}


int main()
{
	scanf("%d",&N);
	for(int i=1;i<=N;i++)
	{
		scanf("%lf",&f[i]);		
	}
	for(int i=1;i<=N;i++)
	{
		for(int j=1;j<=N;j++)
		{
			double x;
			scanf("%lf",&x);
			A[i][j]=2.0*(-f[j]+x);
			A[i][N+1]+=x*x-f[j]*f[j];
		}
	}
	Gauss(N);
	
	return 0;
}