OFFSET
1,3
COMMENTS
MATHEMATICA
f[n_] := 2Plus @@ (Select[ Range@n, GCD[ #, n] == 1 &]^2)/(n*EulerPhi@n); f[1] = 1; Denominator@ Array[f, 80] (* Robert G. Wilson v, Jul 01 2010 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, May 31 2010
EXTENSIONS
More terms from Robert G. Wilson v, Jul 01 2010
STATUS
approved