如何定义＆＃34;素数顺序＆＃34;在pari / gp？

Question

我对数论非常感兴趣，并希望用pari / gp测试我的一些想法，但我不熟悉这个软件。具体来说，我想定义一个“勋章”＆＃39;将整数n映射到我称之为素数顺序的函数，当且仅当n是复合并且等于最小整数k时才为0，否则在n处计算的素数计数函数的第k次迭代是复合的。

如何在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。