高斯消元法解线性方程组

高斯消元法

这种方法，可以在接近O(n³)的复杂度下求解线性方程组，忧郁克拉默法则的O(n*n!)

对于一组线性方程组，枚举每一列进行如下步骤：
1、找到绝对值最大的一行
2、将这一行交换到第一行
3、将这一行的第一个数变成1，对当前这一行进行操作，不涉及矩阵的初等变换
4、将下面所有行的当前列全部消成0，利用矩阵的初等变换

对于原线性方程组进行变换后
如果得到的矩阵是一个完美的上三角矩阵，则说明方程组有唯一解
否则会出现两种情况：
1、推出矛盾，即：零==非零，这种情况下表示无解
2、否则，即原线性方程组中有多余的方程，此时有无穷多解

高斯消元解线性方程组

#include <bits\stdc++.h>
using namespace std;
typedef long long ll;
const int N = 110;
const double eps = 1e-6; //解决c++中存储浮点数存在的精度问题
int n;
double a[N][N];
int gauss()
{
    int c, r; //c表示枚举的哪一列，r表示枚举的哪一行
    for (c = 0, r = 0; c < n; c++)
    {
        int t = r; //记录当前列
        for (int i = r; i < n; i++)
        { //abs返回整数，fabs返回浮点数
            if (fabs(a[i][c]) > fabs(a[t][c]))
                t = i;
        }
        if (fabs(a[t][c]) < eps) 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];
        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 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;
}
int main()
{
    cin >> n;
    for (int i = 0; i < n; i++)
        for (int j = 0; j <= n; j++)
            cin >> a[i][j];
    int t = gauss();
    if (t == 0) //表示只有一组解
    {
        for (int i = 0; i < n; i++)
            printf("%.2lf\n", a[i][n]);
    }
    else if (t == 1)
        cout << "Infinite group solutions" << endl;
    else
        cout << "No solution" << endl;
    return 0;
}