高斯消元模板

详解参照其他博客，简单模板如下

/*处理出的倒三角是这个形状的
x1 x2 x3=..
   x2 x3=..
      x3=..
*/
#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++){//处理第i列，行也处理到i 
        int k=i;//i：当前处理的列 和 行 
        for(int j=i+1;j<=n;j++)if(fabs(a[j][i])>fabs(a[k][i]))k=j;//找到i列最大的系数在的行 
        if(fabs(del=a[k][i])<eps)return 0;//最大数的行都为0，无解或多解     del:i列最大系数 
        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;//最大系数化为 1 
        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;//也可看做a[k][j]=a[k][j]/del-a[i][j]
            //即把a[k][i]也系数化为1，再减去i行，显然k的第i列变成0了，即达到了消元的目的 
        }
    }
    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]);
    bool flag=gauss();
    if(!flag) puts("No Solution");
    else 
    {
        for(int i=1;i<=n;i++)
            printf("%.2lf\n",a[i][n+1]);
    }
    return 0;
}