如何在 pari/gp 中定义“素数顺序”？

Question

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

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

Answer 1

请检查一下：

primeorder(x) = {
  if(!isprime(x), return(0));

  my(k=1, p=primepi(x));
  while(isprime(p), p=primepi(p); k++);
  return(k);
}

我强烈推荐阅读PARI/GP tutorial。