P3389 【模板】高斯消元法

/*
落谷：p3389
给定一个线性方程组，对其求解 
*/ 
#include<bits/stdc++.h>
using namespace std;
const double eps=1e-8;
double a[110][110];
int n;
void Gauss()
{
	for(int i=1;i<=n;i++){//枚举列 
		int maxn=i;
		for(int j=i+1;j<=n;j++){//选择一个当前位置系数绝对值最大的调换过来，防止误差
			if(abs(a[j][i])>abs(a[j][maxn])){
				maxn=j;
			}
		}
		for(int j=1;j<=n+1;j++){//把第i列最大的系数放到第i行 
			swap(a[i][j],a[maxn][j]);
		}
		if(abs(a[i][i])<=eps){//如果该位置系数等于零，则0x=a，一定无解
			printf("No Solution\n");
			exit(0);
		}
		for(int j=1;j<=n;j++){
			if(i!=j){
				double t=a[j][i]/a[i][i];
				for(int k=i;k<=n+1;k++){
					a[j][k]-=a[i][k]*t;//将其他方程用加减法减去系数值
				}
			}
		}
	}
	for(int i=1;i<=n;i++){//把系数消掉 
		a[i][n+1]/=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("%.2f\n",a[i][n+1]);
	}
	return 0;
}