OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = 8 * (A076274(n-1) + 1) for n > 3, as proved by Lawrence Sze. - Ralf Stephan, Nov 16 2004
EXAMPLE
For n=2896: d=10 divisors, r=1440 coprimes, u=1447 unrelated or n - u = r + d - 1 = 1449 related numbers to n; thus abs(1449 - 1447) = 2.
MATHEMATICA
Do[r=EulerPhi[n]; d=DivisorSigma[0, n]; u=n-r-d+1; df=2*u-n; If[Equal[Abs[df], 2], Print[n(*, {d, r, u}*)]], {n, 1, 3000}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 07 2003
STATUS
approved