cp's OEIS Frontend

A161354 Cubes whose reversal is a prime number.

%I A161354 #25 Dec 28 2021 04:17:35
%S A161354 125,125000,140608,704969,1643032,1815848,3511808,3723875,3869893,
%T A161354 7301384,9393931,10360232,11543176,14526784,15069223,15252992,
%U A161354 15813251,16777216,19902511,34328125,34645976,35287552,70444997,92345408
%N A161354 Cubes whose reversal is a prime number.
%H A161354 Chai Wah Wu, <a href="/A161354/b161354.txt">Table of n, a(n) for n = 1..4050</a>
%e A161354 1815848 is a term because it is a cube and 8485181 is a prime.
%p A161354 rev := proc (n) local nn: nn := convert(n, base, 10): add(nn[j]*10^(nops(nn)-j), j = 1 .. nops(nn)) end proc: a := proc (n) if isprime(rev(n^3)) = true then n^3 else end if end proc: seq(a(n), n = 1 .. 460); # _Emeric Deutsch_, Jun 27 2009
%t A161354 Select[Range[500]^3,PrimeQ[IntegerReverse[#]]&] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Jul 03 2019 *)
%o A161354 (Python)
%o A161354 from sympy import isprime
%o A161354 A161354_list, i, j = [], 0, 0
%o A161354 while j < 10**15:
%o A161354     p = int(str(j)[::-1])
%o A161354     if isprime(p):
%o A161354         A161354_list.append(j)
%o A161354     j += 3*i*(i+1)+1
%o A161354     i += 1
%o A161354 A161354_list = sorted(A161354_list) # _Chai Wah Wu_, Dec 20 2015
%o A161354 (PARI) lista(nn) = { for(n=1, nn, if(ispseudoprime(eval(concat(Vecrev(Str(n^3))))), print1(n^3, ", ")); ); } \\ _Altug Alkan_, Dec 20 2015
%Y A161354 Cf. A004165, A057699.
%K A161354 nonn,base
%O A161354 1,1
%A A161354 _Claudio Meller_, Jun 07 2009