问题标签 [pari-gp]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
63 浏览

pari - How to modify a Pari code to test a conjecture

This is a Pari code.

Now how can I modify this program to see if it is true that if numbers generated with that Pari code are multiple of 43, then they are congruent either to 0 or 344 mod (559)?

0 投票
2 回答
615 浏览

pari - 在 Pari/GP 中将 p 进数转换为向量

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

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

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

0 投票
1 回答
89 浏览

vector - 在 PARI/GP 中将二变量多项式的返回幂作为向量?

我正在尝试以向量形式查看二变量多项式的单项式。因此,例如,如果我输入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 ],我认为它不能被扭曲成对我有用的东西。

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

0 投票
1 回答
159 浏览

pari-gp - 在 Pari/GP 中计算 Goldbach 分区的最快方法

我正在尝试计算哥德巴赫分区的数量,并在 n 很大时碰壁。关于如何尽可能快地编写此代码的任何建议都会有所帮助。这是迄今为止我能做到的最好的:

0 投票
1 回答
40 浏览

complex-numbers - 请求高斯分解的简单方法

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

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

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

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

0 投票
1 回答
280 浏览

pari - 如何计算pari中的模乘法逆?

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

ax ≡ 1 (mod m)。

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

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

0 投票
1 回答
70 浏览

pari-gp - 如何获得特定碱基中分数的表示(在 GP shell 中)?

如何获得特定基数中任意分数的表示?例如,以 2 为底的 1/3 的表示应该是

但是当我输入digits(1/3, 2)GP shell 并按 Enter 时,我看到的只是以下错误:

如何解决问题?

0 投票
1 回答
51 浏览

pari - pari gp函数编程中的指针可选参数

请看这个例子

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

0 投票
1 回答
105 浏览

for-loop - 循环 PARI/GP 中的特定值

我有大量向量,每个形式都是[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] ] ) )?

(请记住,在实践中,我还会收集更多的向量)

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

0 投票
1 回答
48 浏览

pari-gp - 打印 Pari/GP 正在处理的并行计算编号

我正在使用 Pari/GP 来测试数字序列的(伪)素性f(n),其中f(n)有一些函数。我已将代码并行化为

这很好用,但我也想知道正在测试哪个号码。我当然可以print(n)if. 问题是,在非常快的同时ispseudoprime(a(n))可能需要相当长的时间。由于计算是并行的,输出很快就会变得一团糟:我会看到已测试的数字列表,但我不知道在给定时刻正在测试什么数字(想法是,如果程序是花费大量时间在一个数字上它越来越可能是素数,所以我想知道它)。nn+1

有没有一种简单的方法可以做到这一点?谢谢!