问题标签 [pari]

我想将 Pari/GP p-adic 数的元素提取到向量中。即，如果这是我的 p 进数：

我希望能够提取这样的向量：

有没有直接的方法可以做到这一点，还是我必须编写自己的函数？

尝试编写 Pari 代码来解决上述问题。

我想并行化一个包含来自 PARI/GP 库的一些操作以及其他一些操作的 for 循环。但是，将 OpenMP 与 PARI/GP 一起使用似乎并不是很简单。我发现的唯一类似的例子是this，它实际上似乎是sections从 OpenMP 中使用的。parallel for正如您在下面看到的，我尝试使用子句做一些事情，但是，它会在 PARI/GP 上引发分段错误。

任何想法如何使用GENinsidestruct和 with OpenMP？

我正在尝试以向量形式查看二变量多项式的单项式。因此，例如，如果我输入x^2 + x^3*y + x*y + y^2 + 1，我想将其视为[[2;0], [3;1], [1;1], [0;2], [0;0]]- 由列向量组成的向量。

如果我Vec在双变量多项式上使用，它只是将第二个变量视为一个数字，给Vec( x^2 + x^3*y + x*y + y^2 + 1 ) = [ y, 1, y, y^2 + 1 ]，我认为它不能被扭曲成对我有用的东西。

关于如何做到这一点的任何想法？

我在高斯整数中工作并factor成功使用它们。但是，有时虚部会消失，但我仍然想要对高斯整数进行因式分解。

例如，程序for(j=-3, 3, print(j,": ",factor(17+j*I)))写道：

我希望输出中的中间行有所不同：

是否有一些设置可以做到这一点？还是我必须编写自己的“包装器”来专门处理这种情况？

模乘逆是一个整数“x”，这样。

ax ≡ 1 (mod m)。

那么，有没有办法在 pari\gp 计算器中计算模乘逆？

如果这个问题是微不足道的或重复的，我真的很抱歉。

我大部分时间都使用 pari 进行计算数论研究。我输入了以下内容：for(x=1,100,print1(eulerphi(n)-1))
它给出了整数 1 到 100 的 φ(n)-1 的值。这里 φ(n) 是小于 n 且与​​ n 互质的整数个数。它给出了以下输出：（
0 0 1 1 3 1 5 3 5 3 ... 65 31 43 23 69 2 43 59 45 71 31 95 41 59 39
我排除了一些整数以使问题简短）我想检查这些数字中是否有数字。我该怎么做？

我非常使用 pari 进行数学研究。但我不知道如何找出一个数字是否可以表示为某个数字 x 和函数 f 的 f(x)。例如，欧拉的总函数 φ(x) 被定义为小于 x 且与 x 互质的自然数的数量。现在我想弄清楚，给定一个数字 n，对于某个正整数 n，它是否可以表示为 φ(x)-1？我无法测试所有正整数，那我该怎么做呢？

请看这个例子

我这样做是为了将结果存储在 e 中，但是当我将函数调用为 f(27,t) 时，程序不会将“3”存储在任何帮助中

我有大量向量，每个形式都是[a,b,c,d].

对于每个向量，我想返回结果[a,b,c,d,a+b^2+c^3+d^4]。这样做的最佳方法是什么？

例如，假设我的向量是V = [ [1,2,3,4], [5,6,7,8], [9, 10, 11, 12] ]. 最好的方法是定义一个函数，例如：

test( W ) = for( i=1, #W, print( [ [W[i][1],W[i][2],W[i][3],W[i][4],W[i][1]+W[i][2]^2+W[i][3]^3+W[i][4]^4] ] ) )?

（请记住，在实践中，我还会收集更多的向量）

这感觉像是一种天真的方法，那么有更好的方法吗？特别是，这是最快的方法吗？有没有更好的方法来存储我的向量集合？