acwing883高斯消元

#include<bits/stdc++.h>
using namespace std;
 
const double eps=1e-6;
const int N=105;
 
int c,r,n,t,flag;
double a[N][N];
 
int guass(){
    for(c=0, r=0;c<n;c++){//c列 r行
        t=r;
        for(int i=r;i<n;i++)
            if(fabs(a[i][c])>fabs(a[t][c])) t=i;//找到绝对值最大的行
         
        if(fabs(a[t][c])<eps) continue;//最大是0
         
        for(int i=c;i<=n;i++) swap(a[r][i],a[t][i]);//将该行换到最上面
         
        for(int i=n;i>=c;i--) a[r][i]/=a[r][c];//将该行第一个数变成1
         
        for(int i=r+1;i<n;i++)//将下面所以行的第c列变成0
            if(fabs(a[i][c])>eps)
            for(int j=n;j>=r;j--)
                a[i][j]-=a[r][j]*a[i][c];
         
        r++;
    }
    if(r<n){//如果非完美阶梯型，即正方形矩阵的行还没遍历完
        for(int i=r;i<n;i++)//从r那行开始，即从为0=？的那行开始
            if(a[i][n]>eps) return -1;//0x=1，无解
        return 1;//0x=0，无穷多组解
    }
    else{//完美阶梯型
        for(int i=r-1;i>=0;i--)//转变成单位矩阵
            for(int j=i+1;j<r;j++)
                a[i][n]-=a[i][j]*a[j][n];
        return 0;
    }
     
}
int main(){
    cin>>n;
    for(int i=0;i<n;i++)
        for(int j=0;j<=n;j++)
            scanf("%lf",&a[i][j]);
    flag=guass();
    if(flag==0){
        for(int i=0;i<n;i++) printf("%.2lf\n",a[i][n]);
    }else if(flag==-1) puts("No solution");
    else puts("Infinite group solutions");
    return 0;
}