OFFSET
1,3
COMMENTS
Conjecture: nondecreasing sequence of odd numbers.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..2500 from Ivan Neretin)
MATHEMATICA
B[n_] := Plus @@ ((ks = Select[Range[n], GCD[n, #] == 1 &])^2)/Plus @@ ks; Select[B /@ Range[215], IntegerQ] (* Ivan Neretin, May 22 2015 *)
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)); }
list(lim) = {my(m); for(k = 1, lim, m = ah(k); if(denominator(m) == 1, print1(m, ", "))); } \\ Amiram Eldar, May 25 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 30 2010, Jul 31 2010
STATUS
approved
