佩林伪素数

伊恩·斯图尔特(Ian Stewart)在六月刊上的数学娱乐专栏属于科学美国人讨论Perrin序列A(n)具有:
A(0)=3,A(1)=0,A(2)=2,A(n+1)=A(n-1)+A(n-2)
n个0124567891011121314151617181920
A(n)022557101217222939516890119158209277
A(n)模块n0002050270509810014017

受E.Lucas定理的启发:如果n是素数,它会精确地除a(n),n的素性是否遵循n除以A(n)的问题已制定1899.他们说,到目前为止,还没有人发现一个可以划分a(n)的复合n。这样的数字称为佩林伪素数。文章引用了史蒂文·阿诺的一项实验马里兰州鲍伊的超级计算研究中心最小Perrin伪素数大小的15位数界1991年获得。

1996年7月3日,我能找到两个最小的佩林伪素数:

稍后:筛选Richard Pinch的246683表格卡迈克尔数最多10 ^16,我数过了150佩林伪素数.

关于计算

A(n)可以直接评估,从初始值开始值并保留最后三个结果。在32位计算机上,A(78)=3354494070是序列的最后一个成员在不使用无限精度算法的情况下计算包裹。人工智能语言,如SICStus公司日食Prolog透明地提供了这种便利。A(271441)是一个33150位数,可以在不改变(平凡的)评估算法。
如果计算的目的是检查一致性模n为0的结果,中间结果当然可以简化为模n。将A(n)mod n的线性计算替换为二进制方法的结果是O(log(n))算法。

你可以马上试试。也许你很幸运,发现了另一个佩林伪素数。

模块计算器(Prolog程序)

输入整数n个>0:输入整数>1(如果为空,则默认为n):(大数字可能被切成几行,非数字材料被丢弃)

此页面的URL:http://www.ai.univie.ac.at/perrin.html

请将问题和意见提交给:

克里斯蒂安·霍尔兹鲍尔http://www.ai.univie.ac.网址:/~克里斯蒂亚