OFFSET
1,2
PROG
(PARI) precision 1000 digits : for(n=1, 300, if(sqrt(denominator(sum(i=1, n, 1/i^3))/denominator(sum(i=1, n, 1/i))) == floor(sqrt(denominator(sum(i=1, n, 1/i^3))/denominator(sum(i=1, n, 1/i)))), print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 07 2002
STATUS
approved