高等代数(一)-多项式04：从辗转相除法到贝祖等式

1. 从最大公因式到辗转相除法

多项式理论中，最大公因式是个绕不开的核心概念。简单来说，如果多项式φ(x)能同时整除f(x)和g(x)，那它就是这两个多项式的公因式。而最大公因式（GCD）就是所有公因式中次数最高的那个，它不仅要是公因式，还得能被其他所有公因式整除。

我第一次接触这个概念时，总觉得"最大"这个词有点误导——它其实不是指系数大小，而是指次数最高。举个例子，对于f(x)=x²-1和g(x)=x³+x²-x-1，它们的公因式有x+1和x-1，其中x+1就是最大公因式。

求最大公因式最经典的方法就是辗转相除法，这个方法我在实际计算时发现特别像小学数学的除法运算。基本思路就是不断用余式去除前一个除数，直到余数为零。具体来说：

1. 用g(x)除f(x)得余式r₁(x)
2. 用r₁(x)除g(x)得余式r₂(x)
3. 用r₂(x)除r₁(x)得余式r₃(x)
4. 重复这个过程，直到余式为0，此时最后一个非零余式就是最大公因式

这个方法的精妙之处在于，每次除法后，新的两个多项式（除数和余式）和原来的两个多项式有着完全相同的公因式集合。我在推导这个性质时，发现它本质上依赖于多项式除法的一个基本性质：如果d(x)能整除f(x)和g(x)，那么它一定能整除f(x)-q(x)g(x)，也就是余式。

2. 手把手演示辗转相除法的计算过程

让我们通过一个具体例子来感受辗转相除法的威力。考虑f(x)=x⁴+3x³-x²-4x-3和g(x)=3x³+10x²+2x-3。

第一步计算f(x)÷g(x)：
f(x) = (1/3x - 1/9)g(x) + (-5/9x² - 25/9x - 10/3)
这里商是(1/3x - 1/9)，余式是(-5/9x² - 25/9x - 10/3)

第二步用余式去除g(x)：
g(x) = (-27/5x + 9)(-5/9x² - 25/9x - 10/3) + (9x + 27)
这次的余式简化为9x + 27

第三步继续除：
(-5/9x² - 25/9x - 10/3) = (-5/81x - 10/81)(9x + 27) + 0
余式终于归零，过程结束

最终我们得到最大公因式是x + 3（9x + 27的首一化形式）。这个计算过程虽然看起来繁琐，但每一步都是机械化的操作，非常适合编程实现。我在实现这个算法时发现，保持多项式首项系数为1可以大大简化计算。

3. 贝祖等式：最大公因式的线性表示

辗转相除法不仅帮我们找到了最大公因式，还隐含着一个更深刻的结论——贝祖等式（Bézout's Identity）。这个等式告诉我们，任何两个多项式f(x)和g(x)的最大公因式d(x)，都可以表示为它们的线性组合：
d(x) = u(x)f(x) + v(x)g(x)

回到之前的例子，我们已经得到了最大公因式x+3。现在要找到u(x)和v(x)使得这个等式成立。这个过程需要"回代"：

从第二步的等式：
9x + 27 = g(x) - (-27/5x + 9)(-5/9x² - 25/9x - 10/3)

将第一步的余式表达式代入：
= g(x) - (-27/5x + 9)[f(x) - (1/3x - 1/9)g(x)]
= (27/5x - 9)f(x) + [1 - (27/5x - 9)(1/3x - 1/9)]g(x)

展开整理后得到：
x + 3 = (3/5x - 1)f(x) + (-1/5x² + 2/5x)g(x)

这就是我们要的贝祖等式。在实际操作中，我发现系数回代的过程最容易出错，特别是符号和分数运算要格外小心。

4. 互素多项式与贝祖等式的应用

当两个多项式的最大公因式是1时，我们称它们互素。根据贝祖等式，f(x)和g(x)互素当且仅当存在u(x)和v(x)使得：
u(x)f(x) + v(x)g(x) = 1

这个性质在多项式理论中有很多重要应用。比如，如果f(x)与g(x)互素，且f(x)整除g(x)h(x)，那么f(x)必定整除h(x)。这个结论在解多项式方程和分式分解时特别有用。

我在研究多项式同余方程时发现，贝祖等式实际上给出了方程f(x)u(x) ≡ 1 mod g(x)的解存在的充要条件——即f(x)和g(x)互素。这个性质类似于整数中的模逆元概念，是构建多项式域扩展的基础。

另一个有趣的应用是，如果f₁(x)|g(x)，f₂(x)|g(x)，且f₁(x)与f₂(x)互素，那么f₁(x)f₂(x)|g(x)。这个结论可以推广到多个多项式的情况，为处理复杂的分式分解提供了理论依据。