至此，我们已经了解到怎么得到一个RREF，但是你或许会有一个疑惑，为什么得到RREF要分为【向前步骤】和【向后步骤】？
这样做是因为REF拥有一个重要的性质——【它可以判断线性方程组是否有解】。
显然不是所有线性方程组都是有解的，比如：

我们显然不能找到组x1，x2，x3满足“0=15”，所以之后我们也就没有必要继续进行“向后步骤”了。
我们称这种无解的线性方程组为【不相容的】，有唯一解或无穷解的方程组为【相容的】

坏了，被线代统治的记忆回来了

通过以上的介绍，我们已经有能力解决最开始提出的问题了：【我们怎么解掉五阶层方阵的16元1次方程？】
答案是【构建一个16*17的增广矩阵，然后求出它的“RREF”。】
通过以上的介绍你肯定可以直观的察觉到这种矩阵变换是一个典型的【P问题】，即可以通过多项式时间复杂度解决的问题，这个结论其实非常的有意义。
你可能会觉得，这样变换矩阵好像也没有把计算变简单多少啊，矩阵的变换还是很吃力。
虽然矩阵确实难算，但是通过上面的解释你可以看到，这一套规则得价值在于，提供了一套纯“代数”的符号体系来研究这种“线性问题”，后来它被作为独立的数学对象研究，发展出的【矩阵论】在いるいる的领域都有着重要意义，但这就不是这个帖子可以介绍的了，我们还是原意是怎么破解五阶层方阵。
其实还有很多其他规则来帮助进行矩阵的运算（比如著名的Cramer's rule），但是这些规则与我们无用，因为我们将采用一种最先进的，所有人都在用的方法解决这个矩阵，通过计算机计算RREF。
由于这是一个“P问题”，所以我们人类总是不缺算力解决这种问题。

由于前面提到的，五阶层方阵构造的线性方程组前的系数都是1，所以这个增广矩阵中只有1,0还有每一行列的右值（315-xx-xx-xx得到的）
由于楼主支付不起matlab的昂贵费用，于是利用了python提供的免费的numpy+sympy来计算RREF：

只要通过短短两行就能算出RREF，都什么年代了还用matlab？

然后结果就华丽的失败了（没算出来）：

显然五阶层方阵是没有负数的，而且还有3个零行，给出的解的数量也不够。

为什么？这就得提到前面楼主埋的一个坑，16个方程真的就能解出16个未知数吗？
答案是——（不完全正确），这里就必须提到两个重要概念【基本变量】和【自由变量】
形如：

其中x3没有对应的主元列（主元位置的行是零行），其中x1，x2是基本变量，x3是自由变量。

给出的16个线性方程中，有不少是没有意义的“多余方程”，
注意，并不是所有“多余方程”都是形如（3*x1+x2=5， 6*x1+2*x2=10）这样显而易见的形式，他们可能藏得很好，但是经过“高斯消元法”变换后都会变成零行。
这就是16个线性方程解不出五阶层方阵的原因——提供的有意义的线性方程不够。

于是我们补上侧面的6个线性方程，两层平面对角线中的2个方程，以及没写全的垂直方程，共28个线性方程。

对应矩阵：

出去多余的零行，于是得到了对应格子的值：124,84,99，...35。

至此，五阶层方阵就被计算出来了。

🤣我在g吧学线代
大一的时候线代还拿了90+，现在就只知道一些基础的概念了，令人感慨

完事了补个妹妹。

密麻麻的，g吧看线代是吧

不这
需并
要非
立线
志代
解的
出故
方事
程，
的
人