高斯消元（模板）

传送门：https://www.luogu.org/problem/P3389

其他的博客会讲的很清晰本蒟蒻就退下了

#include<bits/stdc++.h>
using namespace std;
int n;
double b[109][109];
const double esp=1e-18;
double ans[109];
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n+1;j++)
        {
            scanf("%lf",&b[i][j]);
        }
    }
    for(int i=1;i<=n;i++)
    {
        int pos=i;
        for(int j=i;j<=n;j++)
        {
            if(fabs(b[pos][i])-fabs(b[j][i])<esp)
            {
                pos=j;
            }
        }
        if(fabs(b[pos][i])<esp)
        {
            puts("No Solution");
            return 0;
        }
        double div=b[i][i];
        for(int j=i;j<=n+1;j++)
        {
            b[i][j]/=div;//将当前系数化一 
        }
        for(int j=i+1;j<=n;j++)
        {
            div=b[j][i];
            for(int k=i;k<=n+1;k++)
            {
                b[j][k]-=b[i][k]*div;
            }
        }
    }
    /*
    这个时候的b数组是这个样子的
    1 x x x | x
    0 1 x x | x
    0 0 1 x | x
    0 0 0 1 | x 
    */
    ans[n]=b[n][n+1];//所以只能先确定最后一个解
    for(int i=n-1;i>=1;i--)
    {
        ans[i]=b[i][n+1];    
        for(int j=i+1;j<=n;j++)
        {
            ans[i]-=(b[i][j]*ans[j]);//减去前面的解和它相应的系数 
        }
    }
    for(int i=1;i<=n;i++)
    {
        printf("%.2lf
",ans[i]);
    }
}