OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..2000 from Colin Barker)
EXAMPLE
20 is a term because the divisors of 20 are [1,2,4,5,10,20] and 6 / (1/1 + 1/2 + 1/4 + 1/5 + 1/10 + 1/20) = 20/7.
MATHEMATICA
Select[Range[200], Numerator[HarmonicMean[Divisors[#]]]==#&] (* Harvey P. Dale, May 24 2017 *)
Select[Range[134], Numerator[DivisorSigma[0, #] * #/DivisorSigma[1, #]] == # &] (* Amiram Eldar, Mar 02 2020 *)
PROG
(PARI)
harmonicmean(v) = #v / sum(k=1, #v, 1/v[k])
s=[]; for(n=1, 500, if(numerator(harmonicmean(divisors(n)))==n, s=concat(s, n))); s
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Nov 12 2014
STATUS
approved