%I#28 2021年5月13日01:40:44
%S 3,31,41,5314159141594159,59,29,21592653653,531411592653589,89,
%电话:41592653589789759265358976535897358975897,97,735897958979,79,
%电话:58979392653589793239323、2393238462643462643643、434341592653589793238846264338338462638383、83
%N第一个1、2、3、4、…内出现的未被发现的素数,。。。Pi的数字,A000796(根据最后一个数字的位置排序,然后是初始数字)。
%C考虑第一个,然后前两个,然后是前三个。。。,A000796的术语,即Pi的十进制数字。看看将一定数量的后续数字串联起来是否会产生一个之前没有出现过的素数(因此必然会涉及到所考虑的最后一个数字)。如果是这样,把这个素数加到序列中。
%C包含A005042作为子序列。
%e序列的第一个数字是素数a(1)=3。
%e前两位数字“3.1”表示素数a(2)=31。
%e在“3.14”中没有更多的素数。在“3.141”中,素数a(3)=41。
%e在“3.1415”中有素数a(4)=5。
%e在“3.14159”中,我们有素数314159、14159、4159和59。
%o(PARI){my(PI=数字(PI\.1^30),seen=[]);for(i=1,#PI-1,for(j=1,i,my(p=fromdigits(PI[j.i]));!isprime(p)||setsearch(seen,p)|| print1(p“,”)||seen=setunion(seen),[p])))}\\由_Andrew Howroyd_和_M.F.Hasler_编辑为使用当前的PARI语法,2021年5月10日
%o(PARI){my(a=List());对于(m=0,精度(.)-3,my(pi=pi\.1^m,p);对于
%Y参考A198019(“新”素数根据其大小而非起始位置排序)。
%Y参考A000796、A005042、A047658。
%K nonn、base、dumb
%O 1,1号机组
%A _M.F.Hasler,2011年10月20日