「营业日志 2020.12.10 Sec2」斜对称矩阵 Pfaffian 的正确求法

我们这里约定

$$\mathrm{pf}\left(\begin{array}{cc}0& a\\ -a& 0\end{array}\right)=a$$

基于 $\det A=(\mathrm{pf}A{)}^2$ 的方法总是在一般交换环上难以计算。但是接下来我们考虑一个 wiki 上的结论：

对于斜对称矩阵 $A$ 和任意矩阵 $B$，有 $$\mathrm{pf}(BA{B}^{\mathsf{T}})=\det B\cdot \mathrm{pf}A$$

那么我们令 $B$ 为初等行变换，就会得到如下消元过程：

若当前矩阵形如

$$\left(\begin{array}{ccc}0& a& \\ -a& 0& ⋮\\ & \cdots & \\ & & \ddots \end{array}\right)$$

我们用 $-a$ 所在的行可以消去其下方，而所对应的转置自然而然地通过 $a$ 消去了右边的列。

如果一开始 $a$ 不是 $0$，那么我们需要从下面找到一行，将第 $i(i>2)$ 行与 $2$ 行交换，转置地，将第 $i$ 列与第 $2$ 列交换。同时我们需要将 $\mathrm{pf}$ 乘以 $-1$。

最后，消元后的矩阵形如

$$\left(\begin{array}{ccccccc}& a_1& & & & & \\ -a_1& & & & & & \\ & & & a_2& & & \\ & & -a_2& & & & \\ & & & & \ddots & & \\ & & & \ddots & & & \\ & & & & & & a_{n/2}\\ & & & & & -a_{n/2}& \end{array}\right)$$

显然其 $\mathrm{pf}$ 为 $\prod a_i$。

而实现上有一个更加简单的方法，注意到我们将行变换全部进行后，列变换实际上已经不影响答案的计算，所以我们几乎可以直接按照正常行列式计算时候的消元顺序，应当形如

$$\left(\begin{array}{ccccc}-a_1& & & & \\ & a_1& & & \\ & & \ddots & & \\ & & & -a_{n/2}& \\ & & & & a_{n/2}\end{array}\right)$$

这就引出了 唐老师文章 结尾的奇妙结果。