

A046433


Numbers such that 3 steps are needed to reach a prime under "Sum of digits raised to its digits' powers" procedure.


2



7, 26, 62, 126, 147, 148, 158, 162, 174, 184, 185, 206, 216, 248, 260, 261, 266, 284, 377, 399, 407, 408, 417, 418, 428, 470, 471, 480, 481, 482, 508, 518, 580, 581, 599, 602, 612, 620, 621, 626, 662, 704, 714, 737, 740, 741, 773, 804, 805, 814, 815, 824
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


LINKS



EXAMPLE

26 > 2^2 + 6^6 = 46660 > 4^4 + 6^6 + 6^6 + 6^6 + 0^0 = 140225 > 1^1 + 4^4 + 0^0 + 2^2 + 2^2 + 5^5 = 3391, which is prime, so 26 is a term.


MATHEMATICA

sop[n_] := Total[#^# & /@ (IntegerDigits[n] /. {0 > 1})]; Select[ Range[825], Boole /@ PrimeQ[Rest[NestList[sop, #, 3]]] == {0, 0, 1} &] (* Jayanta Basu, Jun 25 2013 *)


PROG

(Magma) f:=func<n &+[a[i]^a[i]:i in [1..#a]] where a is Intseq(n) >; [k:k in [1..900] not IsPrime(f(k)) and not IsPrime(f(f(k))) and IsPrime(f(f(f(k))))]; // Marius A. Burtea, Oct 27 2019


CROSSREFS



KEYWORD

nonn,base


AUTHOR



STATUS

approved



