OFFSET
1,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..2310
FORMULA
G.f. A(x) (for fractions) satisfies: A(x) = -log(1 - x)/(1 - x) - Sum_{k>=2} A(x^k)/k. - Ilya Gutkovskiy, Mar 31 2020
MATHEMATICA
Table[s=0; Do[If[GCD[i, n]==1, s=s+1/i], {i, n}]; Denominator[s], {n, 1, 35}]
Table[Denominator[Total[1/Select[Range[n], GCD[n, #]==1&]]], {n, 40}] (* Harvey P. Dale, Jun 07 2020 *)
PROG
(PARI) for(n=1, 40, print1(denominator(sum(k=1, n, if(gcd(k, n)==1, 1/k))), ", "))
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Sharon Sela (sharonsela(AT)hotmail.com), Apr 12 2002
EXTENSIONS
More terms from Jason Earls, Apr 14 2002
STATUS
approved