OFFSET
1,1
FORMULA
Positive integers n such that A195860(n)=17.
EXAMPLE
210^1=210 is multiple of Sum_digits(210)=3
210^2=44100 is multiple of Sum_digits(44100)=9
...
210^16=14305686902419853283210000000000000000 is a multiple of Sum_digits(210^16)=90
while
210^17=3004194249508169189474100000000000000000 is not multiple of Sum_digits(210^17)=99
MAPLE
readlib(log10); P:=proc(n, m) local a, i, k, w, x, ok; for i from 1 by 1 to n do a:=simplify(log10(i)); if not (trunc(a)=a) then ok:=1; x:=1; while ok=1 do w:=0; k:=i^x; while k>0 do w:=w+k-(trunc(k/10)*10); k:=trunc(k/10); od; if trunc(i^x/w)=i^x/w then x:=x+1; else if x-1=m then print(i); fi; ok:=0; fi; od; fi; od; end: P(45000, 16);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava and Giorgio Balzarotti, Nov 26 2007
EXTENSIONS
Terms a(7) onward from Max Alekseyev, Sep 24 2011
STATUS
approved