OFFSET
1,1
COMMENTS
Numbers k such that A362802(k) > 0.
Includes all the harmonic numbers (A001599) except for 1, since the set of their divisors has an integer harmonic mean (in this case the partition is into a single part).
This sequence is infinite. For example, if k is a term and p is a prime that does not divide k, then k*p is also a term.
EXAMPLE
12 is a term since its set of divisors, {1, 2, 3, 4, 6, 12} can be partitioned into 2 disjoint parts, {1, 2, 3, 6} and {4, 12}, whose harmonic means, 2 and 6, are both integers.
MATHEMATICA
harmQ[s_] := AllTrue[s, Length[#] > 1 && IntegerQ[HarmonicMean[#]] &]; q[n_] := Module[{d = Divisors[n], r}, r = ResourceFunction["SetPartitions"][d]; AnyTrue[r, harmQ]]; Do[If[q[n], Print[n]], {n, 1, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 04 2023
STATUS
approved