A048393 Replacing digits d in decimal expansion of n with d^3 yields a prime.

11, 13, 23, 31, 41, 43, 53, 61, 73, 101, 107, 109, 121, 137, 143, 149, 151, 157, 161, 169, 173, 181, 191, 211, 217, 221, 229, 233, 241, 253, 257, 259, 271, 277, 281, 299, 307, 311, 313, 319, 323, 331, 421, 427, 431, 449, 469, 493, 511, 527, 541, 577, 589

Offset: 1

Author

Patrick De Geest, Mar 15 1999

Keywords

• nonn
• base

Examples

			313 = (3)(1)(3) -> (27)(1)(27) = 27127, which is a prime.
		

Crossrefs

Cf. A048390, A048394.

Programs

Mathematica

Select[Range[600],PrimeQ[FromDigits[Flatten[IntegerDigits/@ (IntegerDigits[ #]^3)]]]&] (* Harvey P. Dale, Jan 27 2015 *)

Python

from sympy import primerange, isprime
for n in primerange(2,600):
    t=int(''.join(str(int(i)**3) for i in str(n)))
    if sympy.isprime(t):
         print(n)
# Abhiram R Devesh, Feb 09 2015

Extensions

Offset corrected by Michel Marcus, Oct 19 2016