【模板】高斯消元法.cpp

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define dd double
using namespace std;
const int maxn=110;
const dd ex=1e-8;//皮皮的精度 
int n;dd f[maxn][maxn];
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)for(int j=1;j<=n+1;j++)scanf("%lf",&f[i][j]);//将万恶方程读入 
	for(int i=1;i<=n;i++)
	{
		dd a=f[i][i];
		//伪装成fabs的abs(绝对值) 
		if(fabs(a)<=ex) 
		{printf("No Solution\n");return 0;}
		//如果当前a为0,则方程无解 至于原因:当a=0时 f[i][j]/a= 1.#J(正无穷) 
		for(int j=i;j<=n+1;j++)f[i][j]/=a;
		for(int j=1;j<=n;j++)if(j!=i)
		{
			dd b=f[j][i];
			for(int k=i;k<=n+1;k++)f[j][k]-=f[i][k]*b;
		}
	}
	for(int i=1;i<=n;i++) printf("%.2lf\n",f[i][n+1]);//保留两位小数输出 
	return 0;
}