我对数论很感兴趣,想在 pari/gp 中测试我的一些想法,但对这个软件不熟悉。具体来说,我想定义一个 'primeorder' 函数,它将整数 n 映射到我称之为素数的顺序,当且仅当 n 是复合数并且等于最小整数 k 使得第 k 次迭代否则,在 n 处计算的素数计数函数是复合的。
如何在 pari/gp 中定义这样的函数?
请检查一下:
primeorder(x) = {
if(!isprime(x), return(0));
my(k=1, p=primepi(x));
while(isprime(p), p=primepi(p); k++);
return(k);
}
我强烈推荐阅读PARI/GP 教程。