配对GP多项式运算符

时间:2017-02-19 16:07:47

标签: polynomial-math pari

PARI / GP存在问题。有没有人知道在PARI / GP中操作正确的函数/命令,以细化

的最小多项式
[y = x^2-x+1 (mod x^6+x^5+x^4+x^3+x^2+x+1)]

PARI / GP给出了这个错误:

gp > minpoly(x^6+x^5+x^4+x^3+x^2+x+1,{v=x^2-x+1})
  ***   at top-level: ...(x^6+x^5+x^4+x^3+x^2+x+1,v=x^2-x+1)
  ***                                             ^----------
  ***   incorrect type in evaluator [variable name expected] (t_INT).

感谢您的帮助。

我也试试:

(11:36) gp > elt = Mod(x^2-x+1, x^6+x^5+x^4+x^3+x^2+x+1)
%52 = Mod(43, 39991)
(11:36) gp > poly = minpoly(elt, v='y)
%53 = Mod(1, 39991)*y + Mod(39948, 39991)
(11:36) gp > subst(poly, variable(poly), elt)
%54 = Mod(0, 39991)
(11:36) gp >

这应该是一个剧本吗?

1 个答案:

答案 0 :(得分:2)

实际上,您需要以下呼叫:

elt = Mod('x^2-'x+1, 'x^6+'x^5+'x^4+'x^3+'x^2+'x+1)
poly = minpoly(elt, v='y)
gp > y^6 - 6*y^5 + 15*y^4 - 20*y^3 + 22*y^2 - 6*y + 1

只是为了验证:

subst(poly, variable(poly), elt)
gp > 0

v的参数minpoly仅代表变量名称,而不是模数。