暑假acwing算法总结23：高斯消元-CSDN博客

本文链接：https://blog.csdn.net/cjp_kangkang_/article/details/118979380

线性代数高斯消元(显式)
算法步骤：
枚举每一列C
1、找当前列绝对值最大的一行
2、将该行换到最上面的未处理行
3、将该行的第一个数变成1
4、将下面所有行的第C列清成0
代码实现

int guass()
{
    int c//枚举列;
    int r//看行;
    for(c=0,r=0;c<n;c++)
    {
        int 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])<esp)continue;
        for(int i=c;i<=n;i++)swap(a[t][i],a[r][i]);
//步骤二：把当前行换到未处理行的最上面位置
        for(int i=n;i>=c;i--)a[r][i]/=a[r][c];
//步骤三：将该行的第一个数变成1，由后向前除保证了正确性
        for(int i=r+1;i<n;i++)
        {
            if(fabs(a[i][c])>esp)
            {
                for(int j=n;j>=c;j--)
                {
                    a[i][j]-=a[r][j]*a[i][c];
                }
            }
        }
        //out();
//步骤四：将下面该列所有值清零
        r++;
//换到下一行重复操作
    }
//如果是完美上三角r=n一定有唯一解
//如果r<n证明有某一列值全为0并且该行一定在最下面，因为可行的行已经被换上去了
//如果这样的行的值==0，则有无穷解因为ar*xr中ar等于0，xr可以任意取，每取一个就是一组解
//出现！=0即出现0==！0的情况就无解
    if(r<n)
    {
        for(int i=r;i<n;i++)
        {
            if(fabs(a[i][n]>esp))
                return 2;//无解
        }
        return 1;//无穷解
    }
    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;//唯一解
}

// void out()
// {
//     for(int i=0;i<n;i++)
//     {
//         for(int j=0;j<n+1;j++)
//             cout<<a[i][j]<<' ';
//     }
//     cout<<endl;
// }观察变化矩阵
int main()
{
    cin>>n;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n+1;j++)
            cin>>a[i][j];
    }
    if(guass()==0)
    {
        for(int i=0;i<n;i++)printf("%.2lf\n",a[i][n]);
    }
    else if(guass()==1)puts("Infinite group solutions");
    else puts("No solution");
    return 0;
}