A100432 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

A100432

Bisection of A005349.

1, 3, 5, 7, 9, 12, 20, 24, 30, 40, 45, 50, 60, 70, 80, 84, 100, 108, 111, 114, 120, 132, 135, 144, 152, 156, 171, 190, 195, 200, 204, 209, 216, 222, 225, 230, 240, 247, 261, 266, 280, 288, 306, 312, 320, 324, 333, 342, 360, 370, 375, 392, 399, 402, 407, 410, 420

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

LINKS

G. C. Greubel, Table of n, a(n) for n = 1..10000

FORMULA

a(n) = A005349(2*n-1). - 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-1], j=1..58); # Emeric Deutsch, Dec 16 2004

MATHEMATICA

Select[Range[1000], Divisible[#, Total[IntegerDigits[#]]] &][[1;; ;; 2]] (* G. C. Greubel, Apr 09 2023 *)

PROG

(Magma)

A005349:=[n: n in [1..10000] | n mod &+Intseq(n) eq 0];

A100432:= func< n | A005349[2*n-1] >;

[A100432(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)]

def A100432(n): return A005349[2*n-2]

[A100432(n) for n in range(1, 151)] # G. C. Greubel, Apr 09 2023

CROSSREFS