高斯消元解线性方程组

#include<bits/stdc++.h>

using namespace std;
const int N = 110;
typedef double db;
db a[N][N];
const db eps = 1e-8;
int n;

int gauss()
{
        int c,r;
        for(c=0,r=0;c<n;++c)
        {
                int t = r;
                for(int i=r+1;i<n;++i) if(fabs(a[i][c])>fabs(a[r][c])) t = i;
                if(fabs(a[t][c])<eps) continue;

                for(int i=c;i<n+1;++i) swap(a[r][i],a[t][i]);
                for(int i=n;i>=c;--i) a[r][i]/=a[r][c];
                for(int i=r+1;i<n;++i)
                {
                        if(fabs(a[i][c])>eps)
                        {
                                for(int j=n;j>=c;--j)
                                {
                                        a[i][j]-=a[r][j]*a[i][c];
                                }
                        }
                }

                ++r;
        }
        if(r<n)
        {
                for(int i=r;i<n;++i) if(fabs(a[i][n])>eps) return 1;
                return 2;
        }
        for(int i=n-1;i>=0;--i)
        {
                for(int j=i+1;j<n;++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+1;++j)
                {
                        cin>>a[i][j];
                }
        }

        int t = gauss();
        if(t==2)
        {
                cout<<"Infinite group solutions"<<endl;
        }
        else if(t==1) cout<<"No solution"<<endl;
        else
        {
                for(int i=0;i<n;++i)
                {
                        if(fabs(a[i][n])<eps) a[i][n]=0;
                        printf("%.2lf\n",a[i][n]);
                }
        }

        return 0;
}