OFFSET
1,3
COMMENTS
The terms are of course the denominators of the fraction "in smallest terms", otherwise said: a(n) = n*d/gcd(n*d, n - d), which is unambiguous also for n = 1 and n = 2 where n - d = 0.
LINKS
M. F. Hasler, Table of n, a(n) for n = 1..10000 (terms a(3..10^4) from Michel Marcus), Apr 17 2022
EXAMPLE
The number n = 1 has d = 1 divisors, so (n-d)/(n*d) = 0/1 has denominator a(1) = 1.
The number n = 2 has d = 2 divisors, so (n-d)/(n*d) = 0/4 = 0/1 has denominator a(2) = 1 when written in smallest terms.
The number n = 3 has d = 2 divisors, so (n-d)/(n*d) = 1/6 has denominator a(3) = 6.
The number n = 4 has d = 3 divisors, so (n-d)/(n*d) = 1/12 has denominator a(4) = 12.
The number n = 6 has d = 4 divisors, so (n-d)/(n*d) = 2/24 = 1/12 has denominator a(6) = 12.
MATHEMATICA
a[n_] := Numerator[n*(d = DivisorSigma[0, n])/(n - d)]; Array[a, 100, 3] (* Amiram Eldar, Mar 18 2022 *)
PROG
(PARI) A352482(n, d=numdiv(n))=denominator((n-d)/(n*d))
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Michel Marcus, Mar 18 2022
EXTENSIONS
Edited and extended to offset 1 by M. F. Hasler, Apr 17 2022
STATUS
approved