矩阵消元(高斯消元)

安利一波高斯消元的博客，内容很详细。

https://www.cnblogs.com/Dumblidor/p/5751579.html

看完这个相信你已经理解了大概，高斯消元求线性方程组，在学习线性代数（大学课程）的时候我们都接触过。原理是先把线程组转换成矩阵，然后把它等价变换成上三角矩阵，这样从下到上依次可以求出解集。

高斯消元模板:   https://www.luogu.org/problemnew/show/P3389 模板题

#include<bits/stdc++.h>
#define N 205
using namespace std; const double eps=1e-8;
int n; double a[N][N],del;
bool gauss(){
    for(int i=1;i<=n;i++){
        int k=i;
        for(int j=i+1;j<=n;j++) if(fabs(a[j][i])>fabs(a[k][i])) k=j;
        //k对应这一列最大值所在行，如果是0表示无法构成上三角矩阵（下面绝对值都是0）
        //所以有多解或无解，返回不存在。
        if(fabs(del=a[k][i])<eps) return 0;
        //下面是把最大值所在的一行弄上来
        for(int j=i;j<=n+1;j++) swap(a[i][j],a[k][j]);
        for(int j=i;j<=n+1;j++) a[i][j]/=del;
        //下面是求解的过程  咀嚼Ing
        for(k=1;k<=n;k++) if(k!=i){
            del=a[k][i];
            for(int j=i;j<=n+1;j++) a[k][j]-=a[i][j]*del;
        }
    }return 1;
}
int main(){
    cin>>n;
    //输入线性方程组
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n+1;j++) cin>>a[i][j];
    }
    bool flag=gauss();
    if(!flag) return puts("NO solution"),0;
    else{
        for(int i=1;i<=n;i++) printf("%.2lf\n",a[i][n+1]);
    }
    return 0;
}