login
A126143
Primes whose sum of the decimal digits of their cubes is also a prime.
1
13, 17, 23, 37, 47, 67, 107, 139, 181, 193, 211, 223, 229, 241, 283, 307, 311, 313, 331, 337, 349, 353, 367, 373, 389, 397, 421, 433, 439, 487, 499, 523, 547, 571, 577, 587, 607, 613, 631, 709, 719, 727, 739, 751, 773, 787, 839, 853, 877, 881, 907, 983, 991
OFFSET
1,1
LINKS
EXAMPLE
13^3=2197 2+1+9+7=19
17^3=4913 4+9+1+3=17
23^3=12167 1+2+1+6+7=17
MAPLE
P:=proc(n) local i, k, w; for i from 1 by 1 to n do w:=0; k:=ithprime(i)^3; while k>0 do w:=w+(k-trunc(k/10)*10); k:=trunc(k/10); od; if isprime(w) then print(ithprime(i)); fi; od; end: P(500);
MATHEMATICA
Select[Prime[Range[200]], PrimeQ[Total[IntegerDigits[#^3]]]&] (* Harvey P. Dale, Oct 21 2018 *)
CROSSREFS
Sequence in context: A277689 A240894 A244078 * A145483 A349978 A125262
KEYWORD
easy,nonn,base
AUTHOR
STATUS
approved