login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A181805
Number of divisors of A181804(n) that are highly composite (A002182).
5
1, 2, 3, 3, 5, 6, 6, 7, 6, 7, 8, 8, 8, 10, 11, 14, 9, 9, 12, 14, 19, 15, 20, 21, 21, 20, 15, 22, 22, 22, 21, 23, 22, 17, 23, 23, 23, 24, 25, 24, 25, 23, 23, 25, 28, 25, 27, 27, 31, 22, 27, 26, 30, 18, 29, 25, 32, 33, 28, 29, 28, 35, 25, 33, 34, 31, 31, 38, 37
OFFSET
1,2
COMMENTS
a(n) = maximal number of members of A002182 that have a least common multiple of A181804(n). Also, a(n) = length of row A181804(n) in triangles A181802 and A181803.
4, 13 and 16 are the first three positive integers that appear nowhere in this sequence (and, therefore, nowhere in A181801). It would be interesting to know whether there are others.
LINKS
Eric Weisstein's World of Mathematics, Highly composite number.
FORMULA
a(n) = A181801(A181804(n)).
EXAMPLE
A181804(10) = 72 has exactly seven divisors that are members of A002182 (namely, 1, 2, 4, 6, 12, 24 and 36). Hence, a(10) = 7.
MATHEMATICA
seq[max_] := Module[{hcn = {}, hcnmax, d, dm = 0, s = {1}}, Do[d = DivisorSigma[0, n]; If[d > dm, dm = d; AppendTo[hcn, n]], {n, 1, max}]; hcnmax = hcn[[-1]]; Do[s = Union[Join[s, Select[LCM[hcn[[k]], s], # <= hcnmax &]]], {k, 2, Length[hcn]}]; Do[s[[k]] = Count[hcn, _?(Divisible[s[[k]], #] &)], {k, 1, Length[s]}]; s]; seq[300000] (* Amiram Eldar, Jun 23 2023 *)
CROSSREFS
A181806(m) is the m-th member of A181804 such that the value of a(n) increases to a record. See also A181807.
Sequence in context: A296834 A242642 A178041 * A369450 A212010 A366418
KEYWORD
nonn
AUTHOR
Matthew Vandermast, Nov 27 2010
EXTENSIONS
More terms from Amiram Eldar, Jun 23 2023
STATUS
approved