SymPy中旧多项式环子模块未考虑环内逆元的问题

运行代码 #

from sympy import Symbol, GF
x = Symbol("x")
ring = GF(2).old_poly_ring(x, x**-1)
print(ring)
print([0, 1/x] in ring.free_module(2).submodule([0, x]))
print([0, 1] in ring.free_module(2).submodule([0, x]))
print([0, x**3] in ring.free_module(2).submodule([0, x]))

输出结果 #

GF(2)[x,1/x]
False
False
True

问题核心 #

根据环上自由模子模块的定义，由向量$f$生成的子模块应包含所有$r \cdot f$（$r \in$ 环$R$）。此处$R$为Laurent多项式环$GF(2)[x,1/x]$，生成元$f=[0,x]$，按定义：

• $[0,1] = (1/x) \cdot f$
• $[0,1/x] = (1/x)^2 \cdot f$

这两个向量都应属于该子模块，但SymPy返回False，仅$[0,x^3]$返回True。

原因分析 #

SymPy中old_poly_ring创建的Laurent多项式环，其自由模子模块的实现逻辑存在局限性：

• 子模块的成员判断仅允许使用普通多项式环$GF(2)[x]$中的元素作为标量，而非完整的Laurent多项式环$GF(2)[x,1/x]$元素。
• $x^{3$可通过多项式$x}2$与生成元$f$相乘得到（$x^2 \cdot [0,x] = [0,x^{3]$），因此被判定为属于子模块；而$1/x$、$(1/x)}2$不属于普通多项式环，对应的标量乘法结果未被SymPy的子模块实现认可。
• 该问题本质是old_poly_ring的子模块功能未完全适配Laurent环的特性，仍沿用了多项式环子模块的标量限制逻辑。

内容的提问来源于stack exchange，提问作者Joel Klassen