OFFSET
2,7
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 2..10000
MATHEMATICA
f[p_, e_] := p^(2^(-1 + Position[Reverse @ IntegerDigits[e, 2], _?(# == 1 &)])); isigma[1] = 1; isigma[n_] := Times @@ (Flatten@(f @@@ FactorInteger[n]) + 1);
infCoprimeQ[n1_, n2_] := Module[{g = GCD[n1, n2]}, If[g == 1, True, AllTrue[ FactorInteger[g][[;; , 1]], BitAnd @@ IntegerExponent[{n1, n2}, #] == 0 &]]];
a[n_] := Abs[Numerator[n^2 / isigma[n] - Sum[Boole[infCoprimeQ[j, n]], {j, 1, n-1}]]]; Array[a, 100, 2] (* Amiram Eldar, Mar 20 2025 *)
CROSSREFS
KEYWORD
nonn,frac,less
AUTHOR
Vladimir Shevelev, Feb 17 2011
STATUS
approved
