【BZOJ 1013】 [JSOI2008]球形空间产生器sphere 高斯消元

最新推荐文章于 2021-01-03 16:35:55 发布

pbihao

最新推荐文章于 2021-01-03 16:35:55 发布

阅读量155

点赞数

分类专栏： bzoj 数学 ac之路

本文链接：https://blog.csdn.net/pbihao/article/details/53737742

版权

ac之路同时被 3 个专栏收录

303 篇文章 0 订阅

订阅专栏

bzoj

236 篇文章 0 订阅

订阅专栏

数学

30 篇文章 0 订阅

订阅专栏

列方程后高斯消元即可

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define eps 1e-6
using namespace std;
int n;
double a[15][15],b[15][15];

void solve(){
	for(int k,i=1;i<=n;i++){
		for(k=i;k<=n;k++)if(fabs(a[k][i])>eps)break;
		if(k>n)continue;
		if(k^i)for(int j=1;j<=n+1;j++)swap(a[i][j],a[k][j]);
		double z=a[i][i];
		for(int j=1;j<=n+1;j++)a[i][j]/=z;
		for(int j=1;j<=n;j++){
			if(j==i)continue;
			z=a[j][i];
			for(int t=1;t<=n+1;t++){
				a[j][t]-=z*a[i][t];
			}
		}
	}
	for(int i=1;i<=n;i++)
		printf("%.3lf%c",a[i][n+1],i==n?'\n':' ');
}

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n+1;i++){
		for(int j=1;j<=n;j++)
		scanf("%lf",&b[i][j]);
		if(i>1){
			for(int j=1;j<=n;j++){
				a[i-1][j]=2*b[i][j]-2*b[i-1][j];
				a[i-1][n+1]+=b[i][j]*b[i][j]-b[i-1][j]*b[i-1][j];
			}
		}
	}
	solve();
	return 0;
}