省选专练之数学[JSOI2008]球形空间产生器sphere

最新推荐文章于 2019-07-04 22:57:47 发布

LauZiyang

最新推荐文章于 2019-07-04 22:57:47 发布

阅读量132

点赞数

分类专栏：其他省选 2008 高斯消元

本文链接：https://blog.csdn.net/fcb_x/article/details/81389387

版权

其他省选同时被 3 个专栏收录

35 篇文章 0 订阅

订阅专栏

2008

12 篇文章 0 订阅

订阅专栏

高斯消元

4 篇文章 0 订阅

订阅专栏

权当是复习高斯消元了吧

对于一个n维的玩意n+1个点只要两两不共线就可以确定一个东西。

怎么找球心？

易知： $\sum_{i=1}^{n}(A_{i}-Ans_{i})=\sum_{i=1}^{n}(B_{i}-Ans_{i})$

AB为读入的点

拆开式子并且整理有 $\sum_{i=1}^{n}[2*(A_{i}-B_{i})]*Ans_{i}=\sum_{i=1}^{n}(A_{i}^{2}-B_{i}^{2})$

右式明显是常数

左式有n个未知量一共会有n个这样的式子

高斯消元

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-9;
double P[201][201]={};
double a[201][201]={};
int n;
void Gauess(){
	for(int i=1;i<=n;i++){
		int Id=i;
		for(int j=i+1;j<=n;j++){
			if(fabs(a[j][i])-fabs(a[Id][i])>eps)Id=j;
		}
		for(int j=1;j<=n+1;j++)swap(a[i][j],a[Id][j]);
		if(fabs(a[i][i])<eps)continue;
		double t=a[i][i];
		for(int j=1;j<=n+1;j++)a[i][j]/=t;
		for(int j=1;j<=n;j++){
			if(i==j)continue;
			double t=a[j][i];
			for(int k=1;k<=n+1;k++){
				a[j][k]-=t*a[i][k];
			}
		}
	}
	for(int i=1;i<=n;i++){
		if(fabs(a[i][n+1])<eps)cout<<0<<" ";
		else cout<<fixed<<setprecision(3)<<a[i][n+1]<<" ";
	}
}
int main(){
//	freopen("3004.in","r",stdin);
	scanf("%d",&n);
	for(int i=1;i<=n+1;i++){
		for(int j=1;j<=n;j++){
			scanf("%lf",&P[i][j]);
		}
	}
	for(int i=1;i<=n;i++){
		int Id=i+1;
		for(int j=1;j<=n;j++){
			a[i][j]=2.0*(P[1][j]-P[Id][j]);
		}
		double b=0.0;
		for(int j=1;j<=n;j++)b+=(P[1][j]*P[1][j]-P[Id][j]*P[Id][j]);
		a[i][n+1]=b;
	}
	/*for(int i=1;i<=n;i++){
		for(int j=1;j<=n+1;j++)cout<<a[i][j]<<" ";
		cout<<'\n';
	}*/
	Gauess();
	return 0;
}