%I#18 2017年12月15日17:35:49
%S 1,2,3,4,6,10,11,12,15,17,29,48,63,77,8818719033811311832,
%电话:2782290731803272539817530
%N对N进行编号,使(素数(N)#+4)/2是素数,其中x#是素数A034386(x)。
%C编号n,使[A002110(n)/2]+2为素数。
%C这些素数是连续奇数素数加上2:2+[3.5.7…..p(n)]的乘积,如果n在这里。
%C a(19)-a(22)是费马和卢卡斯的PRP。(素数(2782)#+4)/2有10865位数字。适用于Windows的PFGW 1.2.0版[FFT v23.8]基本测试(p(2782)#+4)/2[N-1/N+1,Brillhart-Lehmer-Selfriceridge]使用基数5运行N-1测试使用判别式13运行N+1测试,基数1+sqrt(13)(p(2882)#/4)/2是Fermat和Lucas PRP!-_Jason Earls_,2006年12月12日
%C a(28)>25000。-_Robert Price_,2017年9月29日
%t p=1;Do[p=p*Prime[n];如果[PrimeQ[(p+4)/2],打印[n]],{n,1400}]
%t展平[Position[FoldList[Times,Prime[Range[3000]],_?(PrimeQ[(#+4)/2]&)]](*哈维·P·戴尔,2015年5月24日*)
%o(PARI)n=0;pr=1/2;对于素数(p=2,1e4,n++;pr*=p;if(ispseudoprime(pr+2),print1(n“,”))
%Y参考A002110、A067024、A065026。
%K nonn公司
%O 1,2号机组
%A _Labos Elemer,2001年12月29日
%E更多条款,摘自Robert G.Wilson v_,2001年12月30日
%E a(19)-a(22)来自_Jason Earls_,2006年12月12日
%2013年6月16日,雷·钱德勒的E a(23)
%2017年9月29日,Robert Price_的E a(24)-a(27)
|