login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Primes which can be represented as the sum of a cube and its reverse.
2

%I #12 Jun 26 2021 15:32:31

%S 2,727,79697,85147,100699,3946493,9715169,10301029,11592961,11851481,

%T 13888793,13913093,17746387,125000521,176232571,358030753,417302813,

%U 433748423,726463627,810090007,817807817,832595227,854121557,875444677

%N Primes which can be represented as the sum of a cube and its reverse.

%H Harvey P. Dale, <a href="/A072384/b072384.txt">Table of n, a(n) for n = 1..1000</a>

%e 727 is a term because it is prime and it is the sum of cube 512 and its reverse 215.

%t Select[#+IntegerReverse[#]&/@(Range[1000]^3),PrimeQ]//Union (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Sep 21 2016 *)

%o (Python)

%o from sympy import isprime

%o def rev(n): return int(str(n)[::-1])

%o def aupto(lim):

%o c = [p**3 for p in range(1, int(lim**(1/3))+2)]

%o s = set(ara for ara in (a + rev(a) for a in c) if ara <= lim)

%o return sorted(filter(isprime, s))

%o print(aupto(10**9)) # _Michael S. Branicky_, Jun 26 2021

%Y Cf. A000578, A004086, A004165, A319603.

%K base,nonn

%O 1,1

%A _Shyam Sunder Gupta_, Jul 20 2002