我试图在一个二值有限域{0,1}上构造多项式,并且希望它们使用此设置中存在的某些标识自动简化。
我尝试了以下方法:
from sympy import *
from sympy.polys.domains.finitefield import FiniteField
x, y, z, t = symbols('x y z t')
k = Poly(x+y * z*z + (x + y) + y + 1, domain=FiniteField(2))
这已经简化为:
Poly(y*z**2 + 1, x, y, z, modulus=2)
但是,z ** 2实际上与我要使用的字段中的z相同。它似乎确实会自动识别y + y =0。我如何实现另一个标识z * z = z(幂等)?
最佳答案
想要的东西似乎并没有实现多边形,但是也许可以模拟效果:
In [54]: normalise = lambda p: Poly(ratsimpmodprime(p, [n**2-n for n in p.free_symbols]), modulus=2)
In [55]: e = x+y * z*z + (x + y) + y + 1
In [56]: normalise(e)
Out[56]: Poly(y*z + 1, x, y, z, modulus=2)
关于python - 如何使用Sympy自动简化二值有限域上的表达式?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/56706342/