【问题标题】:Prime counting function in Pari GPPari GP 中的素数计数功能
【发布时间】:2017-01-22 23:34:57
【问题描述】:

素数计数函数 pi(x) 计算小于 x 的素数。由于 Pari 非常擅长处理素数,我认为这个函数会实现,但我在文档中没有找到任何内容。

  • pi(x) 是否已在 Pari GP 中实现?
  • 如果没有,是否可以找到素数的索引?我知道命令 prime(n) 返回第 n 个素数? (如果可以快速找到它,那么结合 precprime 就可以了。)

【问题讨论】:

    标签: pari


    【解决方案1】:

    primepi,如:

    primepi(10^6)
    

    (返回78498

    【讨论】:

    • 感谢您的回答。我没有设法在文档中找到它。
    • 你知道为什么这个函数对于较大的 n 值很慢,即使有足够的预先计算的素数超出这个限制?
    • @BeniBogosel 我不知道任何细节。似乎 GP 只会预先计算 10^8 素数(10^8==primepi(2038074743) 对应于 prime(10^8)==2038074743),即使在启动时名为 primelimit 的默认值远高于 10^8。我基于您的观察,当您进行大输入时,函数 primepi 和(右)反函数 prime 都会变慢。
    猜你喜欢
    • 2016-09-05
    • 2018-02-25
    • 2014-08-23
    • 1970-01-01
    • 2013-01-25
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2014-03-10
    相关资源
    最近更新 更多