%I #32 Nov 24 2023 15:59:10
%S 1,1,2,3,12,3,20,35,280,252,2520,2310,27720,25740,24024,45045,720720,
%T 226893,4084080,775975,246341,235144,5173168,14872858,356948592,
%U 343219800,2974571600,2868336900,80313433200,77636318760,2329089562800
%N a(n) = floor((denominator of H(n))/n), where H(n) = Sum_{k=1..n} 1/k, the n-th harmonic number.
%C This is very similar to A027611, which is a different sequence. - _N. J. A. Sloane_, Nov 21 2008
%C Indices where a(n) differs from A027611 are terms of A074791. - _Gary Detlefs_, Sep 03 2011
%H Amiram Eldar, <a href="/A128438/b128438.txt">Table of n, a(n) for n = 1..2310</a>
%e The sequence denominator(H(n))/n begins 1, 1, 2, 3, 12, 10/3, 20, 35, 280, 252, 2520, 2310, ..., so the present sequence begins 1, 1, 2, 3, 12, 3, 20, 35, 280, 252, 2520, 2310, ...
%p H:=n->sum(1/k,k=1..n): a:=n->floor(denom(H(n))/n): seq(a(n),n=1..34); # _Emeric Deutsch_, Mar 25 2007
%t seq = {}; s = 0; Do[s += 1/n; AppendTo[seq, Floor[Denominator[s]/n]], {n, 1, 30}]; seq (* _Amiram Eldar_, Sep 18 2021 *)
%t Table[Floor[Denominator[HarmonicNumber[n]]/n],{n,40}] (* _Harvey P. Dale_, Nov 24 2023 *)
%o (Python)
%o from sympy import harmonic
%o def A128438(n): return harmonic(n).q//n # _Chai Wah Wu_, Sep 27 2021
%Y Cf. A128437, A002805, A027611, A074791.
%K nonn
%O 1,3
%A _Leroy Quet_, Mar 03 2007
%E More terms from _Emeric Deutsch_, Mar 25 2007
|