OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
PROG
(PARI) list(lim) = print1(1, ", "); for(k = 2, lim, if(!isprime(k), my(f = factor(k), b = 2*k/3 + (1/3) * prod(i = 1, #f~, 1 - f[i, 1])/eulerphi(f)); if(denominator(b) == 1, print1(b, ", ")))); \\ Amiram Eldar, May 26 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 30 2010, Jul 31 2010
EXTENSIONS
Incorrect formula removed by Amiram Eldar, May 26 2025
STATUS
approved
