editing
proposed
editing
proposed
0, 1, 6, 24, 80, 240, 672, 1792, 4608, 11520, 28160, 67584, 159744, 372736, 860160, 1966080, 4456448, 10027008, 22413312, 49807360, 110100480, 242221056, 530579456, 1157627904, 2516582400, 5452595200, 11777605632, 25367150592, 54492397568, 116769423360, 249644974080, 532575944704
(PARI) A001788_upto(n)=Vec(x/(1-2*x)^3+O(x^n), -n) \\ for illustration. - M. F. Hasler, Oct 05 2024
approved
editing
editing
proposed
Cf. A045512 (excludes zero digits).
proposed
editing
editing
proposed
A variant of Munchausen Münchausen numbers. Cf, cf. A166623.
The sequence is finite, because the sum can't exceed 9^9*L < 10^9*L, where L is the number of digits, and for L > 10 this is less than the number N >= 10^(L-1). - M. F. Hasler, Oct 01 2024
Geoff Bailey, <a href="https://homepage.kranzky.com/puzzles/power_ultra.c">C program for the sequence</a> (cf. Hutchens link for more info), Aug. 1998
Jason Hutchens, <a href="https://homepage.kranzky.com/puzzles/Power.html">power summation</a> (originally at ciips.ee.uwa.edu.au/~hutch), 1997
(C) see Bailey and Hutchens links
(PARI) select( {is_A046253(n)=n==A045512(n)}, [0..10^4]) \\ To find the 4th solution, multiply the set by 51817. - M. F. Hasler, Oct 01 2024
approved
editing
editing
proposed
(PARI) apply( {A045512(n)=vecsum([d^d|d<-digits(n), d])}, [0..44]) \\ M. F. Hasler, Oct 01 2024
See A046253 for fixed points.
approved
editing
editing
proposed