OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Digit.
FORMULA
A055014(a(n)) | a(n).
EXAMPLE
21321 is a term since 2^5 + 1^5 + 3^5 + 2^5 + 1^5 = 309 and 21321 = 69*309.
54748 is a term since 5^5 + 4^5 + 7^5 + 4^5 + 8^5 = 54748.
MAPLE
with(numtheory):for n from 1 to 200000 do:l:=evalf(floor(ilog10(n))+1): n0:=n:indic:=0:s5:=0:for m from 1 to l do:q:=n0:u:=irem(q, 10):v:=iquo(q, 10):n0:=v :s5:=s5+u^5: od:if irem(n, s5)=0 then print (n):else fi:od:
MATHEMATICA
Select[Range[10^4], Divisible[#, Plus @@ (IntegerDigits[#]^5)] &] (* Amiram Eldar, Jan 31 2021 *)
PROG
(PARI) is_A169666(n)=!(n%sum(i=1, #n=Vecsmall(Str(n)), (n[i]-48)^5))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Apr 05 2010
EXTENSIONS
Corrected and edited by D. S. McNeil, Nov 20 2010
More terms from Amiram Eldar, Jan 31 2021
STATUS
approved