%I #24 Sep 27 2021 11:24:57
%S 0,1,5,1,17,9,83,201,2089,2341,551,2861,64913,90653,114677,274399,
%T 5385503,2022061,42503239,9276623,3338549,3573693,87368107,276977179,
%U 7281378067,7624597867,71595952403,74464289303,2239777822987
%N Numerator of the fractional part of a harmonic number.
%H Clark Kimberling, <a href="/A104174/b104174.txt">Table of n, a(n) for n = 1..1000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BookStackingProblem.html">Book Stacking Problem.</a>
%F a(n) = numerator(frac(Sum_{k=1..n} 1/k)). [edited by _Michel Marcus_, Sep 27 2021]
%t h[n_] := Sum[1/k, {k, 1, n}]
%t Table[Numerator[FractionalPart[h[n]]], {n, 1, 30}]
%t (* _Clark Kimberling_, Aug 13 2012 *)
%t FractionalPart[HarmonicNumber[Range[30]]]//Numerator (* _Harvey P. Dale_, Jul 28 2019 *)
%o (Python)
%o from sympy import harmonic
%o def A104174(n): return (lambda x: x.p % x.q)(harmonic(n)) # _Chai Wah Wu_, Sep 26 2021
%o (PARI) a(n) = numerator(frac(sum(k=1, n, 1/k))); \\ _Michel Marcus_, Sep 27 2021
%Y Cf. A064169, A075135.
%K nonn
%O 1,3
%A Georg Haass (geha5001(AT)stud.uni-saarland.de), Mar 10 2005
|