OFFSET
1,1
COMMENTS
Note that 11 is the only two-digit number in the sequence.
a(n) ~ n. For 414 < n < 10000, 6.38*n - 528 provides an estimate of a(n) to within 6%.
LINKS
Christian N. K. Anderson, Table of n, a(n) for n = 1..10000
EXAMPLE
139 is in the sequence because 1^3 + 3^3 + 9^3 = 757, which is prime.
MATHEMATICA
Select[Range[350], PrimeQ[Total[IntegerDigits[#]^3]]&] (* Harvey P. Dale, Mar 16 2016 *)
PROG
(R)
digcubesum<-function(x) sum(as.numeric(strsplit(as.character(x), split="")[[1]])^3); library(gmp);
which(sapply(1:1000, function(x) isprime(digcubesum(x))>0))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Kevin L. Schwartz and Christian N. K. Anderson, May 10 2013
STATUS
approved