Search: id:A125724 Results 1-1 of 1 results found. %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 Search completed in 0.001 seconds