OFFSET
1,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A005349(2*n). - G. C. Greubel, Apr 09 2023
MAPLE
s:=proc(n) local N:N:=convert(n, base, 10):sum(N[j], j=1..nops(N)) end:p:=proc(n) if floor(n/s(n))=n/s(n) then n else fi end: A:=[seq(p(n), n=1..440)]: seq(A[2*j], j=1..58); # Emeric Deutsch, Dec 16 2004
MATHEMATICA
Select[Range[1000], Divisible[#, Total[IntegerDigits[#]]] &][[2;; ;; 2]] (* G. C. Greubel, Apr 09 2023 *)
PROG
(Magma)
A005349:=[n: n in [1..10000] | n mod &+Intseq(n) eq 0];
[A100433(n): n in [1..150]]; // G. C. Greubel, Apr 09 2023
(SageMath)
A005349=[n for n in (1..10^4) if sum(n.digits(base=10)).divides(n)]
[A100433(n) for n in range(1, 151)] # G. C. Greubel, Apr 09 2023
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
N. J. A. Sloane, Nov 20 2004
EXTENSIONS
More terms from Emeric Deutsch, Dec 16 2004
STATUS
approved