OFFSET
1,2
COMMENTS
a(14) > 10^12. - Giovanni Resta, Jun 02 2016
EXAMPLE
For n=2, 2*Sd(2) = 2*2 = 4 and Sd(2^2) = Sd(4) = 4.
For n=3, 2*Sd(27) = 2*9 = 18 and Sd(27^2) = Sd(729) = 18;
and 3*Sd(27) = 3*9 = 27 and Sd(27^3) = Sd(19683) = 27.
MAPLE
T:=proc(w) local x, y, z; x:=w; y:=0; for z from 1 to ilog10(x)+1 do
y:=y+(x mod 10); x:=trunc(x/10); od; y; end:
P:=proc(q) local a, j, k, n, ok, t; t:=1; for k from 1 to q do
for n from t to q do ok:=1; a:=T(n);
for j from 2 to k do if j*a<>T(n^j) then ok:=0; break; fi; od;
if ok=1 then t:=n; print(n); break; fi; od; od; end: P(10^20);
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Paolo P. Lava, Jun 01 2016
EXTENSIONS
a(12)-a(13) from Giovanni Resta, Jun 02 2016
STATUS
approved