This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A125724 #9 Oct 22 2013 10:10:15
%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)).
%C A125724 Sequence is complete. - _Giovanni Resta_, Mar 21 2013
%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);
%t A125724 sdQ[n_]:=Module[{sd=Total[IntegerDigits[n]]},n==Total[ IntegerDigits[ sd^sd]]]; Select[Range[100],sdQ] (* _Harvey P. Dale_, Oct 22 2013 *)
%Y A125724 Cf. A125526.
%K A125724 easy,fini,full,nonn,base
%O A125724 1,2
%A A125724 _Paolo P. Lava_ and _Giorgio Balzarotti_, Feb 02 2007