OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
727 is a term because it is prime and it is the sum of cube 512 and its reverse 215.
MATHEMATICA
Select[#+IntegerReverse[#]&/@(Range[1000]^3), PrimeQ]//Union (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Sep 21 2016 *)
PROG
(Python)
from sympy import isprime
def rev(n): return int(str(n)[::-1])
def aupto(lim):
c = [p**3 for p in range(1, int(lim**(1/3))+2)]
s = set(ara for ara in (a + rev(a) for a in c) if ara <= lim)
return sorted(filter(isprime, s))
print(aupto(10**9)) # Michael S. Branicky, Jun 26 2021
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Shyam Sunder Gupta, Jul 20 2002
STATUS
approved