%I A125724
%S A125724 1,13,25,45,58,88,98
%N A125724 If "sumdigit" denotes the sum of the digits of a number then these are 
               the numbers n such that n=sumdigit(sumdigit(n)^sumdigit(n)).
%e A125724 sumdigit(13)=4; 4^4 = 256; sumdigit(256)=13
%p A125724 P:=proc(n) local i,j,k,w; for i from 1 by 1 to n do w:=0; k:=i; while 
               k>0 do w:=w+k-trunc(k/10)*10; k:=trunc(k/10); od; k:=w^w; w:=0; while 
               k>0 do w:=w+k-trunc(k/10)*10; k:=trunc(k/10); od; if (i=w) then print(i,
               w); fi; od; end: P(100);
%Y A125724 Cf. A125526.
%Y A125724 Sequence in context: A017533 A026058 A032697 this_sequence A163551 A005696 
               A147145
%Y A125724 Adjacent sequences: A125721 A125722 A125723 this_sequence A125725 A125726 
               A125727
%K A125724 easy,fini,nonn,base
%O A125724 0,2
%A A125724 Paolo P. Lava & Giorgio Balzarotti (ppl(AT)spl.at), Feb 02 2007