OFFSET
1,2
COMMENTS
Lim_{n-> infinity} (a(n)/A195506(n)) = 2*Zeta(3) [L. Euler].
Sum_{k = 1..n} H(k)/k^2 is an example of a multiple harmonic (star) sum. Euler's result Sum_{k = 1..inf} H(k)/k^2 = 2*zeta(3) was the first evaluation of a multiple zeta star value. - Peter Bala, Jan 31 2019
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..768
Leonhard Euler, Meditationes circa singulare serierum genus, Novi. Comm. Acad. Sci. Petropolitanae, 20 (1775), 140-186.
FORMULA
From Peter Bala, Jan 31 2019: (Start)
Let S(n) = Sum_{k = 1..n} H(k)/k^2. Then
S(n) = 1 + (1 + 1/2^3)*(n-1)/(n+1) + (1/2^3 + 1/3^3)*(n-1)*(n-2)/((n+1)*(n+2)) + (1/3^3 + 1/4^3)*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ...
S(n)/n = 1 + (1/2^4 - 1)*(n-1)/(n+1) + (1/3^4 - 1/2^4)*(n-1)*(n-2)/((n+1)*(n+2)) + (1/4^4 - 1/3^4)*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ...
For odd n >= 3, 1/2*S((n-1)/2) = (n-1)/(n+1) + 1/2^3*(n-1)*(n-3)/((n+1)*(n+3)) + 1/3^3*(n-1)*(n-3)*(n-5)/((n+1)*(n+3)*(n+5)) + ....
EXAMPLE
a(2) = 11 because 1 + (1 + 1/2)/2^2 = 11/8.
The first few fractions are 1, 11/8, 341/216, 2953/1728, 388853/216000, 403553/216000, 142339079/74088000, 1163882707/592704000, ... = A195505/A195506. - Petros Hadjicostas, May 06 2020
MATHEMATICA
s = 0; Table[s = s + HarmonicNumber[n]/n^2; Numerator[s], {n, 20}] (* T. D. Noe, Sep 20 2011 *)
PROG
(PARI) H(n) = sum(k=1, n, 1/k);
a(n) = numerator(sum(k=1, n, H(k)/k^2)); \\ Michel Marcus, May 07 2020
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
Franz Vrabec, Sep 19 2011
STATUS
approved