%I#15 2017年5月29日23:59:20
%S 127142720032713763976519817031396631664317377349589,
%电话:3710276533576960837520337936998835499386179747131150733,
%电话:1176983120722313107791675577170257718803632715169
%N具有p,q,r,s连续素数的p*q+q*r+r*s形式的素数。
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%e a(1)=127=3*5+5*7+7*11=
%e A000040(2)*A000040=
%e A006094(2)+A006093(3)+A006 094(4)。
%p N:=100:#得到a(1)-a(N)
%p p:=2:q:=3:r:=5:s:=7:
%p计数:=0:
%p而计数<N do
%p p:=q;q: =r;r: =秒;s: =下一素数;
%pn:=p*q+q*r+r*s;
%p如果isprime(n),则计数:=计数+1;A[count]:=n fi
%日期:
%p序列(A[i],i=1..N);#_Robert Israel_,2017年5月29日
%o(PARI){p=2;q=3;r=5;s=7;对于(k=11000,if(i素数(a=p*q+q*r+r*s),
%o打印1(“,”);p=q;q=r;r=秒;s=下一素数(1+s)}
%o(Python)
%o从sympy导入nextprime,isprime
%o A287653_列表,pq,qr,rs,s=[],6,15,35,7
%o当s<=10**6时:
%o n=pq+qr+rs
%o如果是质数(n):
%o A287653_列表.附加(n)
%o t=下一素数
%o pq,qr,rs,s=qr,r,s*t,t#_Chai Wah Wu_,2017年5月29日
%Y参考A000040(素数),A006094(两个连续素数的乘积)。
%K nonn公司
%O 1,1号机组
%A _Zak Seidov,2017年5月29日
|