gauss elimination

gausscutEradication method♂

It is to use the matrix to add and subtract elements

#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]);
}