0

我尝试解决计算第 5 分圆域中理想范数的问题，（表示为 K5） polcyclo(5) 定义该域的多项式在 PARI/GP 中取得了成功。我还能够找到它的主要生成器并进行更多计算（第三行和第四行的命令）。

(09:24) gp > J = idealhnf(bnfinit(polcyclo(5)),11,x-3);
(09:24) gp > idealnorm(bnfinit(polcyclo(5)),J)
%16 = 11
(09:24) gp > Jp = bnfisprincipal(bnfinit(polcyclo(5)),J)
%17 = [[]~, [1, -1, -1, 0]~]
(09:24) gp > u = nfbasistoalg(bnfinit(polcyclo(5)),Jp[2])
%18 = Mod(-x^2 - x + 1, x^4 + x^3 + x^2 + x + 1)
(09:24) gp >

我想为 K179 以相同的速度执行相同的任务，但这不起作用。PARI/GP 有什么问题？我需要另一个程序来执行此操作吗？

(09:25) gp > J = idealhnf(bnfinit(polcyclo(179)),359,x-2);
  ***   at top-level: J=idealhnf(bnfinit(polcyclo(179
  ***                            ^--------------------
  *** bnfinit: the PARI stack overflows !
  current stack size: 16000000 (15.259 Mbytes)
  [hint] you can increase GP stack with allocatemem()

(09:26) gp > idealnorm(bnfinit(polcyclo(179)),J)
  ***   at top-level: idealnorm(bnfinit(polcyclo(179
  ***                           ^--------------------
  *** bnfinit: the PARI stack overflows !
  current stack size: 16000000 (15.259 Mbytes)
  [hint] you can increase GP stack with allocatemem()

(09:27) gp > Jp = bnfisprincipal(bnfinit(polcyclo(179)),J)
  ***   at top-level: Jp=bnfisprincipal(bnfinit(polcyclo(179
  ***                                   ^--------------------
  *** bnfinit: the PARI stack overflows !
  current stack size: 16000000 (15.259 Mbytes)
  [hint] you can increase GP stack with allocatemem()

(09:28) gp > u = nfbasistoalg(bnfinit(polcyclo(179)),Jp[2])
  ***   at top-level: u=nfbasistoalg(bnfinit(polcyclo(179
  ***                                ^--------------------
  *** bnfinit: the PARI stack overflows !
  current stack size: 16000000 (15.259 Mbytes)
  [hint] you can increase GP stack with allocatemem()

(09:29) gp >

我还将堆栈大小增加到 1024000000（最大值），它仍然给出相同的溢出错误。有谁知道我应该怎么做才能正确执行所有四个计算？提前致谢。

4

0 回答 0