|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
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
|
|
|
FORMULA
|
|
|
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.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|