问题标签 [pari]

我有一个函数可以遍历其输入并产生零个或多个结果，每个结果由三个数字组成。我想将这些结果保存在数据结构中（例如，矩阵或向量向量），但我不知道在循环终止之前会有多少条目。我需要能够轻松地提取一列结果（例如每个条目的第一个变量）。

我对数论很感兴趣，想在 pari/gp 中测试我的一些想法，但对这个软件不熟悉。具体来说，我想定义一个 'primeorder' 函数，它将整数 n 映射到我称之为素数的顺序，当且仅当 n 是复合数并且等于最小整数 k 使得第 k 次迭代否则，在 n 处计算的素数计数函数是复合的。

如何在 pari/gp 中定义这样的函数？

我是 PARI/GP 的新手，正在学习不同的数字字段。我正在尝试k = Q(\zeta_23)使用bnfinit(23rd cyclotomic field) 构建该字段来运行以下脚本：

运行到 GP 会出现此错误，

这应该找到素数 p，使得域 K 中有代数整数，范数 p。

请问有什么帮助吗？谢谢。

有没有一种简单方便的方法可以将 Pari/GP 中的输出定向到文件？2^400000-1我的目标是在屏幕上或文本文件中获得完整的十进制扩展？

GP 终端输出给出了这个，这不是目标。基本输出重定向也不起作用。有任何想法吗？谢谢。

网上有一本手册，关于输出没说太多，除了变量TeXstyle。我不确定如何使用它。

PARI/GP 是否具有查找 a 的最小素因数t_INT或以其他方式执行整数的部分因式分解的功能？

例如，如果我有号码：

它需要很长时间才能完成，factor(a)因为a包含两个巨大的主要因素。但是，很容易发现17是 的一个除数a。

当然，在这种情况下，我可以只使用forprime(p=2,,a % p == 0 && return(p))或类似的试验部门来找到因素。但是，如果最小的因素有 20 个小数位数，那将是不切实际的，在这种情况下，我可能想使用复杂的方法factor。

因此，如果我可以factor用某种标志打电话说我会对任何部分分解感到满意，或者说我只关心最小的非平凡除数等，那将是理想的。

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

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

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

我在处理 Pari 中的多项式和在文档中找到正确的命令时遇到了一些问题。

1. 是否可以定义具有多个变量的多项式，例如 f(x,y)=x^2+y^2-1
2. 如何以特定值评估先前定义的多项式（通过使用 Pol() ）？

我有两个（对）要绘制的列表。我知道，我可以使用 plothraw 函数分别绘制每个图。但是我怎样才能将它们绘制在同一张图片中，这样我最终会得到两条不同颜色的曲线呢？

我想计算超过 1-1/p 的乘积，其中 p 在素数上运行高达 10^10

我知道近似值 exp(-gamma)/ln(10^10) ，其中 gamma 是 Euler-Mascheroni 常数， ln 是自然对数，但我想计算确切的乘积以查看近似值有多接近。

问题是 PARI/GP 需要很长时间来计算从大约 4.2 * 10^9 到 10^10 的素数。prodeuler-command 也需要很长时间。

有什么方法可以加快 PARI/GP 的计算速度？

我得到了一个（非常）简单的 DSA 问题，并且已经找到了关键和其他变量。为了验证签名，我需要以某种方式转换等式：

V = [( y^u1*h^u2 )mod p] mod q

进入 BigInteger 操作。这在java上甚至可能吗？到目前为止，我一直在成功使用 modPow，但是到目前为止所有问题都在以下形式：

r.modPow（指数，模数）；

我不知道如何通过 BigInteger 完成上述等式（特别是粗体部分），我想知道它是否可能。有没有人有任何想法？

如果 BigInteger 做不到，我将如何将这个等式通过 Pari？