OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
PROG
(PARI) isok(k) = {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, Jul 31 2010
EXTENSIONS
More terms from Amiram Eldar, May 25 2025
STATUS
approved
