高斯消元原理

高消一直是ACM中高层次经常用到的算法，虽然线性代数已经学过，但高消求解的问题模型及高消模板的应用变化是高消的最复杂之处。

先介绍一下高消的基本原理：引入互联网czyuan的帖子：

高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。

高斯消元法的原理是：

若用初等行变换将增广矩阵 化为 ，则AX = B与CX = D是同解方程组。

所以我们可以用初等行变换把增广矩阵转换为行阶梯阵，然后回代求出方程的解。

 

以上是线性代数课的回顾，下面来说说高斯消元法在编程中的应用。

 

首先，先介绍程序中高斯消元法的步骤：

(我们设方程组中方程的个数为equ，变元的个数为var，注意：一般情况下是n个方程，n个变元，但是有些题目就故意让方程数与变元数不同)

 

1. 把方程组转换成增广矩阵。

 

2. 利用初等行变换来把增广矩阵转换成行阶梯阵。

枚举k从0到equ – 1，当前处理的列为col(初始为0) ，每次找第k行以下(包括第k行)，col列中元素绝对值最大的列与第k行交换。如果col列中的元素全为0，那么则处理col + 1列，k不变。

 

3. 转换为行阶梯阵，判断解的情况。

 

① 无解

当方程中出现(0, 0, …, 0, a)的形式，且a != 0时，说明是无解的。

 

② 唯一解

条件是k = equ，即行阶梯阵形成了严格的上三角阵。利用回代逐一求出解集。

 

③ 无穷解。

条件是k < equ，即不能形成严格的上三角形，自由变元的个数即为equ – k，但有些题目要求判断哪些变元是不缺定的。

    这里单独介绍下这种解法：

首先，自由变元有var - k个，即不确定的变元至少有var - k个。我们先把所有的变元视为不确定的。在每个方程中判断不确定变元的个数，如果大于1个，则该方程无法求解。如果只有1个变元，那么该变元即可求出，即为确定变元。

 

以上介绍的是求解整数线性方程组的求法，复杂度是O(n3)。浮点数线性方程组的求法类似，但是要在判断是否为0时，加入EPS，以消除精度问