[模板] 高斯消元法

题目描述：

求线性方程组的解

题目分析：

高斯消元法，用矩阵来存储线性方程…
利用加减消元法把矩阵消成上三角矩阵，然后我们就可以从后往前递推回带就可以求出每个解了.
回带完毕之后这个上三角矩阵就变成了只有主对角线为1的矩阵，那么n+1的值就是解惹

题目链接：

貌似很水的模板

Ac 代码：

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <cstring>
const double eps=1e-8;
const int maxm=201; 
double a[maxm][maxm],del;
int n;
inline bool Gauss_Eli()
{
    for(int i=1;i<=n;i++)
    {
        int k=i;
        for(int j=i+1;j<=n;j++) if(std::fabs(a[j][i])>std::fabs(a[k][i])) k=j;
        if(fabs(del=a[k][i])<eps) return 0;
        for(int j=i;j<=n+1;j++) std::swap(a[i][j],a[k][j]);
        for(int j=i;j<=n+1;j++) a[i][j]/=del;
        for(int l=1;l<=n;l++)
        if(l!=i)
        {
            del=a[l][i];
            for(int j=i;j<=n+1;j++) a[l][j]-=a[i][j]*del;
        } 
    }
    return 1;
}
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]);
    if(Gauss_Eli())
    {
        for(int i=1;i<=n;i++)
         printf("%.2lf\n",a[i][n+1]);
    }
    else puts("No Solution");
    return 0; 
}