%I #20 May 24 2024 14:24:22
%S 1,288,288,288,288,288,288,288,288,1,288,288,288,288,50119,3439,288,
%T 288,288,288,288,288,288,288,288,288,288,288,288,288,288,288,288,288,
%U 288,288,3439,288,288,288,288,288,288,288,288,288,288,288,288,288,50119
%N Smallest limiting value of n under iteration of "Sum of its digits raised to its digits power" (A045512).
%C Numbers that yield a fixed value after adding its digits raised to its digits power are called "Munchhausen Numbers" (A046253); i.e. 3435 = 3^3 + 4^4 + 3^3 + 5^5.
%C Since (9^9)*n < 10^n for n > 12; every initial value should eventually reach a "Munchhausen Number" or a cycle.
%C This sequence assigns to n that Munchhausen number or the lowest member of the cycle.
%C About 80% of the numbers reach the 288 cycle; 15% the 3439 cycle; 5% the 50119 cycle and less than 0.1% the fixed value 438579088; only few sporadic numbers reach 1 or 3435.
%C It has been proved that there are only four Munchhausen numbers (0, 1, 3435 and 438579088) with the convention of 0^0 = 0.
%C Open questions: are there any other cycles than those described here? What is the % of numbers < 10^n reaching a specific limiting value for n = 1, 2, 3, etc...?
%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/MuenchhausenNumber.html">Munchhausen Number</a>.
%t f[n_] := Module[{d = IntegerDigits[n]}, Sum[If[i == 0, 0, i^i], {i, d}]]; Table[s = NestWhileList[f, n, UnsameQ[##] &, All]; Min[Drop[s, Position[s, s[[-1]], 1, 1][[1, 1]]]], {n, 100}] (* _T. D. Noe_, Jul 12 2012 *)
%Y Cf. A046253, A045512.
%K base,nonn
%O 1,2
%A _Sergio Pimentel_, Jul 11 2012