OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..300
EXAMPLE
12 is a term since the harmonic means of the divisors of 12 and 13 are 18/7 and 13/7, respectively, and both have the denominator 7.
MATHEMATICA
dh[n_] := Denominator[DivisorSigma[0, n]/DivisorSigma[-1, n]]; Select[Range[10^6], dh[#] == dh[# + 1] &]
PROG
(PARI) f(n) = my(d=divisors(n)); denominator(#d/sum(k=1, #d, 1/d[k])); \\ A099378
isok(k) = f(k) == f(k+1); \\ Michel Marcus, Oct 20 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 17 2021
STATUS
approved
