%我#11 2021年2月7日14:58:02
%序号2,7,11101
%N素数p,使得p^3是回文。
%注意,所有列出的前4个术语都是回文。对应的回文立方体a(n)^3列在A135067={8,343,1331,1030301,…}中。素数Pi[a(n)]={1,4,5,26,…}。
%C没有低于1.29*10^10的其他条款。-_Michael S.Branicky_,2021年2月7日
%H Patrick De Geest,<a href=“http://www.worldofnumbers.com/cube.htm“>回文立方体</a>
%F a(n)=A135067(n)^(1/3)。
%e a(3)=11,因为11^3=1331是一个回文。
%t Do[p=素数[n];f=p^3;如果[f==FromDigits[Reverse[IntegerDigits[f]],打印[{n,p,f}]],{n,1,200000}]
%o(Python)
%o来自sympy import nextprime
%o定义ispal(n):s=str(n);返回s==s[::-1]
%o p=2
%o为True时:
%o如果ispal(p**3):打印(p)
%o p=下一个素数(p)#_Michael S.Branicky_,2021年2月7日
%Y参考A002780(立方体是一个回文)、A069748(n和n^3都是回文),A002781(回文立方体),A135067(素数的回文立方体)。
%K nonn,基础,更多
%O 1,1号机组
%A_Alexander Adamchuk,2007年11月16日
|