This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A072384 #16 Apr 02 2025 15:20:00 %S A072384 2,727,79697,85147,100699,3946493,9715169,10301029,11592961,11851481, %T A072384 13888793,13913093,17746387,125000521,176232571,358030753,417302813, %U A072384 433748423,726463627,810090007,817807817,832595227,854121557,875444677 %N A072384 Primes which can be represented as the sum of a cube and its reverse. %H A072384 Harvey P. Dale, <a href="/A072384/b072384.txt">Table of n, a(n) for n = 1..1000</a> (corrected by Sean A. Irvine) %e A072384 727 is a term because it is prime and it is the sum of cube 512 and its reverse 215. %t A072384 Select[#+IntegerReverse[#]&/@(Range[1000]^3),PrimeQ]//Union (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Sep 21 2016 *) %o A072384 (Python) %o A072384 from sympy import isprime %o A072384 def rev(n): return int(str(n)[::-1]) %o A072384 def aupto(lim): %o A072384 c = [p**3 for p in range(1, int(lim**(1/3))+2)] %o A072384 s = set(ara for ara in (a + rev(a) for a in c) if ara <= lim) %o A072384 return sorted(filter(isprime, s)) %o A072384 print(aupto(10**9)) # _Michael S. Branicky_, Jun 26 2021 %Y A072384 Cf. A000578, A004086, A004165, A319603. %K A072384 base,nonn %O A072384 1,1 %A A072384 _Shyam Sunder Gupta_, Jul 20 2002