OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..3465 from Carmine Suriano)
FORMULA
A003132[a(n)] | a(n). - R. J. Mathar, Feb 25 2007
EXAMPLE
a(100) = 4131 since 4^2+1^2+3^2+1^2=27 divides 4131. - Carmine Suriano, May 04 2013
MAPLE
isA034087 := proc(n) if n mod A003132(n) = 0 then true ; else false ; end if ; end proc:
for n from 1 to 1800 do if isA034087(n) then printf("%d ", n) ; end if ; end do ; # R. J. Mathar, Feb 25 2007
MATHEMATICA
Select[Range[1500], Divisible[#, Plus @@ (IntegerDigits[#]^2)] &] (* Amiram Eldar, Jan 31 2021 *)
PROG
(PARI) isok(m) = !(m % norml2(digits(m))); \\ Michel Marcus, Jan 31 2021
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
STATUS
approved