OFFSET
1,1
COMMENTS
For corresponding values of numbers h+1 see A179880. - Jaroslav Krizek, Jul 31 2010
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A179880(n) - 1. - Jaroslav Krizek, Jul 31 2010
PROG
(PARI) ah(n) = {my(f = factor(n)); if(n == 1, 1, 2*n/3 + (1/3) * prod(i = 1, #f~, 1 - f[i, 1])/eulerphi(f)); }
isok(k) = {my(ah1 = ah(k), ah2 = ah(k+1)); ah1 == ah2 && denominator(ah1) > 1; } \\ Amiram Eldar, May 24 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 30 2010
EXTENSIONS
a(36) corrected and more terms added by Amiram Eldar, May 24 2025
STATUS
approved
