|
%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
|
