【计算方法】高斯消去法（线性方程组的求解）

在本章，我们需要了解一些基本的线性代数基础知识。如非奇异矩阵，矩阵的初等变换等。

高斯消元法

我们以一个例子来讲解高斯消元法

那么，我们可以将它抽象为一个增广矩阵，
该矩阵由系数矩阵和y矩阵组成。我们采用初等变换的方法将系数矩阵转换为上三角矩阵，该做法即为高斯消去法。其中要注意，在做初等交换时，若akk=0，则为需要将k行去与其下面的行进行交换，寻求非零元，若不存在非零元，则该系数矩阵时退化矩阵，该线性方程组也寻求不到唯一解。若化简成功则只需要应用上次讲的上三角线性方程组求解的解法即可解的答案。

有回代的高斯消去法做法

而这里又涉及到了一个选取主元的问题，通常情况下我们有两种选取主元的策略

平凡选主元策略
当且仅当app=0时，往下面去寻找第一个含有非零元的行，然后交换两行
但由于计算机使用固定精度计算，这样每次算术计算中可能引入微小误差，而平凡选主元策略不能很好解决这一点。为了解决这一点我们提出了
偏序选主元策略
为了减少传播误差，我们在选取app的第p行及其下面行中对应第p个元素的绝对值最大的那行。

例题

【问题描述】为求解一个线性方程组，首先构造增广矩阵[A|B]，采用偏序选主元策略的高斯消去法变换成上三角矩阵，再执行回代过程得到解。
【输入形式】在屏幕上依次输入方阵阶数n，系数矩阵A和常数矩阵B。
【输出形式】首先输出上三角矩阵（变换后的增广矩阵），然后每一行输出一个根
【样例1输入】
4
1 2 1 4
2 0 4 3
4 2 2 1
-3 1 3 2
13
28
20
6
【样例1输出】
4 2 2 1 20
0 2.5 4.5 2.75 21
0 0 4.8 3.6 26.4
0 0 0 3.75 7.5
3
-1
4
2
【样例1说明】输入：第1行为方阵阶数4，第2行至5行为系数矩阵A，第6行至9行为常数矩阵B。
输出：首先输出上三角矩阵（变换后的增广矩阵），然后每行依次输出方程解：x1, x2, x3, x4。

ACcode：

/*
 * @Author: csc
 * @Date: 2021-04-02 19:59:11
 * @LastEditTime: 2021-04-24 15:13:28
 * @LastEditors: Please set LastEditors
 * @Description: In User Settings Edit
 * @FilePath: \code_formal\course\cal\gause.cpp
 */
#include <iostream>
#include <vector>
#include <iomanip>
#include <algorithm>
#define pr printf
#define sc scanf
#define sf(n) scanf("%d", &n)
#define sff(n1, n2) scanf("%d %d", &n1, &n2)
#define sfff(n1, n2, n3) scanf("%d %d %d", &n1, &n2, &n3)
#define sl(n) scanf("%lld", &n)
#define sll(n1, n2) scanf("%lld %lld", &n1, &n2)
#define slll(n1, n2, n3) scanf("%lld %lld %lld", &n1, &n2, &n3)
#define for0(i, n) for (i = 0; i < n; i++)
#define