login
A100433
Bisection of A005349.
2
2, 4, 6, 8, 10, 18, 21, 27, 36, 42, 48, 54, 63, 72, 81, 90, 102, 110, 112, 117, 126, 133, 140, 150, 153, 162, 180, 192, 198, 201, 207, 210, 220, 224, 228, 234, 243, 252, 264, 270, 285, 300, 308, 315, 322, 330, 336, 351, 364, 372, 378, 396, 400, 405, 408, 414
OFFSET
1,1
LINKS
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:= func< n | A005349[2*n] >;
[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)]
def A100433(n): return A005349[2*n-1]
[A100433(n) for n in range(1, 151)] # G. C. Greubel, Apr 09 2023
CROSSREFS
Cf. A005349.
Sequence in context: A078106 A353764 A111082 * A254343 A373574 A157006
KEYWORD
nonn,easy,base
AUTHOR
N. J. A. Sloane, Nov 20 2004
EXTENSIONS
More terms from Emeric Deutsch, Dec 16 2004
STATUS
approved