gauss消元

gauss削消元法♂

就是利用矩阵，进行加减消元

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<iomanip>
using namespace std;
double map[111][111];
double ans[111];
double eps=1e-7;
int main()
{
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n+1;j++)
			scanf("%llf",&map[i][j]);
	for(int i=1;i<=n;i++)
	{
		int r=i;
		for(int j=i+1;j<=n;j++)
			if(fabs(map[r][i])<fabs(map[j][i]))//为了防止精度损失
				r=j;
		if(fabs(map[r][i])<eps)//如果无解
		{
			printf("No Solution");
			return 0;
		}
		if(i!=r)
			for(int j=1;j<=n+1;j++)
				swap(map[i][j],map[r][j]);
		double div=map[i][i];
		for(int j=i;j<=n+1;j++)
			map[i][j]/=div;//同除，使第i行i列的未知数的系数变为1，方便下面加减消元
		for(int j=i+1;j<=n;j++)
		{
			div=map[j][i];
			for(int k=i;k<=n+1;k++)
				map[j][k]-=map[i][k]*div;//开始消元
		}
	}
	ans[n]=map[n][n+1];
	for(int i=n-1;i>=1;i--)//回溯，求解
	{
		ans[i]=map[i][n+1];
		for(int j=i+1;j<=n;j++)
			ans[i]-=(map[i][j]*ans[j]);
	}
	for(int i=1;i<=n;i++)
		printf("%.2lf\n",ans[i]);
}