如何在 pari/gp 中定义 "primality order"？

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。