login
9 steps needed to reach a prime under "Sum of digits raised to its digits' powers" procedure.
2

%I #9 Oct 27 2019 00:31:16

%S 133,138,183,303,308,313,318,330,331,380,381,389,398,446,464,644,789,

%T 798,803,813,830,831,839,879,893,897,938,978,983,987,1005,1015,1024,

%U 1042,1046,1050,1051,1064,1105,1115,1124,1142,1146,1150,1151,1164,1204

%N 9 steps needed to reach a prime under "Sum of digits raised to its digits' powers" procedure.

%t sop[n_] := Total[#^# & /@ (IntegerDigits[n] /. {0 -> 1})]; Select[Range[1205], (z = Boole /@ PrimeQ[Rest[NestList[sop, #, 9]]])[[9]] - Total[Drop[z, -1]] == 1 &] (* _Jayanta Basu_, Jun 25 2013 *)

%Y Cf. A046431.

%K nonn,base

%O 0,1

%A _Patrick De Geest_, Jul 15 1998