OFFSET
1,1
COMMENTS
Note that all first 4 listed terms are the palindromes. Corresponding palindromic cubes a(n)^3 are listed in A135067 = {8, 343, 1331, 1030301, ...}. PrimePi[ a(n) ] = {1, 4, 5, 26, ...}.
No further terms less than 1.29 * 10^10. - Michael S. Branicky, Feb 07 2021
LINKS
Patrick De Geest, Palindromic Cubes
FORMULA
a(n) = A135067(n)^(1/3).
EXAMPLE
a(3) = 11 because 11^3 = 1331 is a palindrome.
MATHEMATICA
Do[ p = Prime[n]; f = p^3; If[ f == FromDigits[ Reverse[ IntegerDigits[ f ] ] ], Print[ {n, p, f} ]], {n, 1, 200000} ]
PROG
(Python)
from sympy import nextprime
def ispal(n): s = str(n); return s == s[::-1]
p = 2
while True:
if ispal(p**3): print(p)
p = nextprime(p) # Michael S. Branicky, Feb 07 2021
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Alexander Adamchuk, Nov 16 2007
STATUS
approved