haskell - 如何在 Haskell 中计算精确（整数）行列式？

Question

我想使用hmatrix但它坚持使用Double，这有时会产生舍入错误。查看源代码，我尝试了

type instance DoubleOf Rational = Rational

但未DoubleOf导出（Numeric.ContainerBool隐藏），这可能是有原因的。

Answer 1

bed-and-breakfast您可以使用package对矩阵执行精确操作。

这是示例 ghci 会话：

>>> import Numeric.Matrix
>>> import Data.Ratio
>>> let m = fromList [[1 % 2, 2 % 3], [3 % 4, 4 % 5]] :: Matrix Rational
>>> det m
(-1) % 10
>>> m * m
 3 % 4  13 % 15
 39 % 40  57 % 50

>>> m^4
 563 % 400  819 % 500
 7371 % 4000  10723 % 5000