高斯消元Gauss

模板题

题目转载———caioj.cn

题目描述

有$n$个未知数$x_1,x_2,x_3\dots \dots x_n$，满足：
$a_1^1\times x_1+a_2^1\times x_2+a_3^1\times x_3\dots \dots +a_n^1\times x_n=b_1$
$a_1^2\times x_1+a_2^2\times x_2+a_3^2\times x_3\dots \dots +a_n^2\times x_n=b_2$
$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$
$a_1^n\times x_1+a_2^n\times x_2+a_3^n\times x_3\dots \dots +a_n^n\times x_n=b_n$
求$x_1,x_2,x_3\dots \dots x_n$的值。（保证有解）

输入格式

第一行给出$n(1\le n\le 100)$
下来$n$行，每行给出$n+1$个实数，分别是$a_1^i{a}_n^i和b_i$。

输出格式

输出$x_1,x_2,x_3\dots \dots x_n$的值（相邻两个用一个空格隔开，每个数保留$3$位小数）

样例输入

3
2.5 5.0 3.0 32.5
1.0 4.5 2.0 22
4.0 3.5 1.5 26.5

样例输出

3.000 2.000 5.000

讲解

————来源百度经验

代码

#include<cstdio>
#include<cstring>
using namespace std;
double a[110][110],f[110];int n;
void gauss()
{
	for(int i=1;i<n;++i)
		for(int i2=i+1;i<=n;++i)
		{
			double bili/*(哲♂学)比例*/=a[i2][i]/a[i][i];
			for(int j=1;i<=n+1;++i)a[i2][j]-=bili*a[i][j];
		}
	for(int i=n;i>=1;--i)
	{
		double t=a[i][n+1];
		for(int j=i+1;j<=n;++j)t-=f[j]*a[i][j];
		f[i]=t/a[i][i];
	}
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;++i)for(int j=1;j<=n+1;++j)scanf("%lf",&a[i][j]);
	gauss();
	for(int i=1;i<n;++i)printf("%.3lf ",f[i]);
	printf("%.3lf\n",f[n]);
	return 0;
}