This site is supported by donations to The OEIS Foundation.

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A001600 Harmonic means of divisors of harmonic numbers. (Formerly M0609 N0220) 9
 1, 2, 3, 5, 6, 5, 8, 9, 11, 10, 7, 15, 15, 14, 17, 24, 24, 21, 13, 19, 27, 25, 29, 26, 44, 44, 29, 46, 39, 46, 27, 42, 47, 47, 54, 35, 41, 60, 51, 37, 48, 45, 49, 50, 49, 53, 77, 86, 86, 51, 96, 75, 70, 80, 99, 110, 81, 84, 13, 102, 82, 96, 114, 53, 108, 115, 105, 116, 91, 85, 105 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 COMMENTS Values of n*tau(n)/sigma(n) corresponding to terms of A001599, where tau(n) (A000005) is the number of divisors of n and sigma(n) is the sum of the divisors of n (A000203). REFERENCES N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence). N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence). LINKS R. J. Mathar, Table of n, a(n) for n = 1..937, extending the former b-file of T. D. Noe. Marco Abrate, Stefano Barbero, Umberto Cerruti, Nadir Murru, The Biharmonic mean, arXiv:1601.03081 [math.NT], 2016. G. L. Cohen, Email to N. J. A. Sloane, Apr. 1994 M. Garcia, On numbers with integral harmonic mean, Amer. Math. Monthly 61, (1954). 89-96. Takeshi Goto, All harmonic numbers less than 10^14 Takeshi Goto, Table of a(n) for n=1..937 O. Ore, On the averages of the divisors of a number, Amer. Math. Monthly, 55 (1948), 615-619. O. Ore, On the averages of the divisors of a number (annotated scanned copy) MATHEMATICA A001600 = Reap[Do[tau = DivisorSigma[0, n]; sigma = DivisorSigma[1, n]; h = n*tau/sigma; If[IntegerQ[h], Print[h]; Sow[h]], {n, 1, 90000000}]][[2, 1]](* Jean-François Alcover, May 11 2012 *) PROG (Haskell) a001600 n = a001600_list !! (n-1) a001600_list =    [numerator m | x <- [1..], let m = hm x, denominator m == 1] where    hm x = genericLength divs * recip (sum \$ map recip divs)           where divs = map fromIntegral \$ a027750_row x -- Reinhard Zumkeller, Apr 01 2014 (PARI) lista(nn) = for (n=1, nn, if (denominator(q=n*numdiv(n)/sigma(n)) == 1, print1(q, ", "))); \\ Michel Marcus, Jan 13 2016 CROSSREFS Cf. A001599, A090240 (sorted values). Sequence in context: A191665 A254105 A247548 * A175578 A316609 A000036 Adjacent sequences:  A001597 A001598 A001599 * A001601 A001602 A001603 KEYWORD nonn,nice,easy AUTHOR EXTENSIONS More terms from Matthew Conroy, Jan 15 2006 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 18 11:28 EST 2019. Contains 319271 sequences. (Running on oeis4.)