算法基础课——第四章数学知识（三）

最新推荐文章于 2024-05-22 09:22:48 发布

VIP文章华北理工大学ACM协会

最新推荐文章于 2024-05-22 09:22:48 发布

阅读量977

点赞数 1

分类专栏：算法竞赛——算法基础课文章标签： c++ 数学

本文链接：https://blog.csdn.net/m0_62021646/article/details/122784913

版权

第四章数学知识（三）

如无特殊说明，所有数均为正整数.

高斯消元是用来解方程的，可以在 $O(n^3)$ 时间复杂度内求解一个 $n$ 个方程 $n$ 个未知数的多元线性方程组.
方程组的解有三种情况： $①$ 无解、 $②$ 无穷多组解、 $③$ 唯一解.
样例方程组为 $\begin{cases}x_1+2x_2-x_3=-6\\2x_1+x_2-3x_3=-9\\-x_1-x_2+2x_3=7\end{cases}$ ，答案为 $\begin{cases}x_1=1\\x_2=-2\\x_3=3\end{cases}$
高斯消元法应用了线性代数中矩阵的初等行变换.

这里一些比较复杂的线性代数知识就不涉及了，只介绍高斯消元的做法和原理；
虽然有线性代数的术语，但是所运用的原理还是初中和小学的知识.

将系数抽出来，得到一个 $n\times (n+1)$ 的系数矩阵，对系数矩阵可以做一些初等行变换的操作；
样例方程组的系数矩阵为 $\begin{pmatrix}1 & 2 & -1 & -6\\2 & 1 & -3 & -9\\-1 & -1 & 2 & 7\end{pmatrix}$
初等行变换一共有三种，这三种行列变换可以保证不影响整个方程组的解：
1. 某一行乘一个非零的数（等价于方程两边乘一个相同的数）
2. 交换某两行（等价于把方程组中某两个数交换一下位置）
3. 把某行的若干倍加到另一行上去（等价于某个方程的若干倍加到另一个方程上去）
通过这三种基本操作，把矩阵化简为一个上三角的形式 $\left\{\begin{aligned}a_{11}x_1+a_{12}x_2+\cdots+a_{1n}x_n=b_1 \\ a_{22}x_2+\cdots+a_{2n}x_n=b_2 \\ \cdots\\ a_{nn}x_n=b_n\end{aligned}\right.$

变成上三角有什么好处呢？
首先 $x_n$ 可以直接求出；
$x_{n-1}$ 的方程为 $a_{n-1,n}x_{n-1}+a_{nn}x_n=b_{n-1}$ ，通过 $x_n$ 就可以推出 $x_{n-1}$ 的值；
以此类推，就可求出所有未知数的值了.

什么情况会没有解？
- 如果通过化简，矩阵恰好是一个完美的阶梯型（即 $a_{ii}$ 都存在），完美阶梯型就会有唯一解；
- 如果是不完美的阶梯型，而且出现类似 $0=b_i$ 的方程，则无解；
- 如果是不完美的阶梯型，而且出现类似 $0 = 0$ 的方程，则有无穷多组解.

具体算法

高斯消元是一个迭代的过程，需要枚举 $1\sim n-1$ 列，假设当前枚举到第 $c$ 列：
1. 找到第 $c$ 列中绝对值最大的一行；
2. 将该行换到第 $c$ 行上去；
3. 将该行第一个数变成 $1$ ；
4. 将下面所有行的第 $c$ 列都消成 $0$ .
遍历到第 $1$ 列：
- $\begin{pmatrix}1 & 2 & -1 & -6\\2 & 1 & -3 & -9\\-1 & -1 & 2 & 7\end{pmatrix} \overset{1.2.}{\to} \begin{pmatrix}2 & 1 & -3 & -9\\1 & 2 & -1 & -6\\-1 & -1 & 2 & 7\end{pmatrix} \overset{3.}{\to} \begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\1 & 2 & -1 & -6\\-1 & -1 & 2 & 7\end{pmatrix} \overset{4.}{\to} \begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1.5 & 0.5 & -1.5\\0 & -0.5 & 0.5 & 2.5\end{pmatrix}$
遍历到第 $2$ 列：
- $\begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1.5 & 0.5 & -1.5\\0 & -0.5 & 0.5 & 2.5\end{pmatrix} \overset{3.}{\to} \begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1 & \frac{1}{3} & -1\\0 & -0.5 & 0.5 & 2.5\end{pmatrix} \overset{4.}{\to} \begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1 & \frac{1}{3} & -1\\0 & 0 & \frac{3}{2} & 2\end{pmatrix}$
遍历到第 $3$ 列：
- $\begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1 & \frac{1}{3} & -1\\0 & 0 & \frac{3}{2} & 2\end{pmatrix}\overset{3.}{\to} \begin{pmatrix}1 & 0.5 & -1.5 & -4.5\\0 & 1 & \frac{1}{3} & -1\\0 & 0 & 1 & 3\end{pmatrix}$

关注