OFFSET
1,2
COMMENTS
All terms for n>1 are even.
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = denominator((n+1)*H(n)^2-(2*n+1)*H(n)+2*n), where H(n) is the n-th harmonic number.
a(n) = denominator((S(n)*H(n)^2+(2*n-2*S(n)+1)*H(n) - 2*n)/(H(n) - 1)), where S(n) = the n-th partial sum of H(n).
EXAMPLE
The squares of the first three harmonic numbers are 1, 9/4, 121/36 which sum to 119/18 so a(3) = 18.
MAPLE
H2:= n-> add(harmonic(k)^2, k = 1..n): seq(denom(H2(n)), n=1..25);
PROG
(PARI) a(n) = denominator(sum(k=1, n, sum(i=1, k, 1/i)^2)); \\ Michel Marcus, Apr 07 2025
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Gary Detlefs, Apr 05 2025
STATUS
approved
