OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1749 from G. C. Greubel)
MATHEMATICA
f[n_] := 2 Plus @@ (Select[ Range@n, GCD[ #, n] == 1 &]^2)/(n*EulerPhi@n); Select[ Range@ 111, ! PrimeQ@# && ! IntegerQ@f@# &] (* Robert G. Wilson v, Aug 02 2010 *)
PROG
(PARI) isok(k) = if(isprime(k), 0, my(f = factor(k)); if(k == 1, 0, denominator(2*k/3 + (1/3) * prod(i = 1, #f~, 1 - f[i, 1])/eulerphi(f)) > 1)); \\ Amiram Eldar, May 25 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 30 2010
EXTENSIONS
More terms from Robert G. Wilson v, Aug 02 2010
STATUS
approved
