%我
%S 1,2,3,6,4,5,7,8,10,11,9,13,12,15,16,14,17,18,21,19,20,23,22,25,24,26,
%电话:31,27,28,30,29,32,34,33,35
%N形式为6*k+1的N个连续素数序列的第一次出现的顺序。
%C该序列通过增加幅度排序后给出了指数A055625(n)。
%H J.K.Andersen,<a href=“http://primerrecords.dk/congruent-primes.htm“>连续同余素数。
%H Hugo Pfoertner,<a href=“http://www.randomwalk.de/sequences/a085515.txt“>记录6k+-1形式素数的生成运行</a>
%H Hugo Pfoertner,<a href=“http://www.enginemonitoring.net/math/isprim.txt“>FORTRAN函数ISPRIM检查素性</a>
%e a(3)=3,a(4)=6,a(5)=4,因为形式6*k+1(174117471753175917771783)的6个连续素数的第一次出现发生在形式6*k+1(151157163)的3个连续素物的第一次发生之后,以及4个连续“6k+1”-素数(3049306130673079)的第一次发作之前。
%t(*程序不适合计算大量项。*)
%t pp=表[{p=素数[n],Mod[p,6]},{n,10^7}];
%t sp=拆分[pp,Mod[#1[[2]],6]==Mod[#2[[2]],6]&];
%tb[n_]:=SelectFirst[sp,Length[#]==n&&MatchQ[#,{_Integer,1}..}]&][[1,1]];
%t排序[表[{b[n],n},{n,1,16}][[All,2]](*_Jean-François Alcover_,2018年11月21日*)
%o链接处给出的FORTRAN程序。
%Y参考A055625、A085516(排序出现的第一批“6k-1”素数)。
%K更多,nonn
%O 1,2号机组
%A _胡戈·普福尔特纳,2003年7月31日
%E来自Herman Jamke(hermanjamke(AT)fastmail.fm)的更多术语,2006年10月27日
%E a(32)-a(35)摘自Giovanni Resta,2013年8月4日
|