OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
s = {}; sum1 = sum2 = 0; Do[sum1 += 1/j; sum2 += 1/j^2; If[Denominator[sum1^2] == Denominator[sum2], AppendTo[s, j]], {j, 1, 1000}]; s (* Amiram Eldar, Feb 18 2021 *)
PROG
(PARI) isok(k) = denominator(sum(j=1, k, 1/j)^2) == denominator(sum(j=1, k, 1/j^2)); \\ Michel Marcus, Feb 15 2021
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Mar 27 2002
STATUS
approved