%I#30 2023年1月10日20:17:10
%S 4,6,12,18,30,42,60,72,93102108138140150180182190192198,
%电话:2282402422502702822842923123383463433393405410418,
%电话420422430432453462484834955225325485565705786600618
%N与素数邻居相差一个平方的数。
%夹在双素数之间的C数构成子序列。
%这个序列的第一个素数是9551_阿洛伊斯·海因茨,2022年4月22日
%H Alois P.Heinz,n的表格,n=1..10000的a(n)</a>
%93的主要邻居是89和97,他们都是93乘4的平方差。因此,93就是这个序列。
%e 140人的主要邻居是149人和139人。它们与140的差异分别为9和1。两个差值都是平方,因此,140在这个序列中。
%pq:=n->与映射(issqr,[n-前素数(n),下素数(n)-n]):
%p选择(q,[$3.700])[];#_阿洛伊斯·海因茨,2022年4月22日
%t选择[Range[3],2000],IntegerQ[Sqrt[NextPrime[#]-#]]和&IntegerQ[Sqrt[#-Prime[PrimePi[Next Prime[#1]-1]]&]
%o(PARI)isok(k)=(k>1)&&issquare(下一素数(k+1)-k)&&isquare(k-precprime(k-1));\\_米歇尔·马库斯,2022年4月22日
%o(Python)
%o从itertools导入islice,count
%o从sympy导入integer_ntroot、nextprime、prevprime
%o定义A353073_gen(startvalue=3):#术语生成器>=startvalue
%o q=下一个质数(最大值(起始值,3)-1)
%o p,r=前素数(q),下素数(q)
%o为True时:
%o如果integer_ntroot(q-p,2)[1]和integer_Ntroot(r-q,2)[1]:
%o产量q
%o t=q
%o对于计数(1)中的i:
%o t+=2*i-1
%o如果t>=r:
%o中断
%o如果integer_ntroot(r-t,2)[1]:
%o产量t
%o p,q,r=q,r,下一素数(r)
%o A353073_list=list(岛屿(A353073-gen(),30))#_Chai Wah Wu_,2022年4月22日
%Y参考A353072。
%K nonn公司
%O 1,1
%A _Tanya Khovanova,2022年4月21日
|