|
%I #14 Apr 29 2022 19:54:50
%S 1,1,4,27,256,3125,46656,823543,16777216,387420489,2,2,5,28,257,3126,
%T 46657,823544,16777217,387420490,5,5,8,31,260,3129,46660,823547,
%U 16777220,387420493,28,28,31,54,283,3152,46683,823570,16777243,387420516,257,257,260
%N If decimal expansion of n is ab...d, a(n) = a^a + b^b +...+ d^d.
%H Harvey P. Dale, <a href="/A045503/b045503.txt">Table of n, a(n) for n = 0..1000</a>
%p a:= n-> add(i^i, i=convert(n,base,10)):
%p seq(a(n), n=0..42); # _Alois P. Heinz_, Apr 29 2022
%t f[n_]:=Module[{idn=IntegerDigits[n]/.{0->1}},Total[idn^idn]]; Array[f,100,0] (* _Harvey P. Dale_, May 14 2012 *)
%o (PARI) a(n)=if(n, n=digits(n); sum(i=1,#n,n[i]^n[i]), 1) \\ _Charles R Greathouse IV_, Oct 07 2015
%Y Cf. A045512.
%K nonn,base
%O 0,3
%A _Jeff Burch_
%E Checked by Neven Juric (neven.juric(AT)apis-it.hr), Feb 04 2008
|
