python 有限域函数库_python – Sympy：在有限域中求解矩阵

sympy的Matrix类支持模块化反转.这是模5的示例：

from sympy import Matrix, pprint

A = Matrix([

[5,6],

[7,9]

])

#Find inverse of A modulo 26

A_inv = A.inv_mod(5)

pprint(A_inv)

#Prints the inverse of A modulo 5:

#[3 3]

#[ ]

#[1 0]

用于查找行减少的梯形形式的rref方法支持关键字iszerofunction,其指示矩阵内的哪些条目应被视为零.我相信预期用途是数值稳定性(将小数字视为零),但我也用它来模块化减少.

这是模5的示例：

from sympy import Matrix, Rational, mod_inverse, pprint

B = Matrix([

[2,2,3,2,2],

[2,3,1,1,4],

[0,0,0,1,0],

[4,1,2,2,3]

])

#Find row-reduced echolon form of B modulo 5:

B_rref = B.rref(iszerofunc=lambda x: x % 5==0)

pprint(B_rref)

# Returns row-reduced echelon form of B modulo 5, along with pivot columns:

# ([1 0 7/2 0 -1], [0, 1, 3])

# [ ]

# [0 1 -2 0 2 ]</