OFFSET
1,1
COMMENTS
The average of the divisors of n is equal to sigma(n)/tau(n).
Up to 5*10^12, there are only 3 terms for which the mean is not an integer, namely 254641594575, 280895287491 and 328966666100.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..6934 (terms < 5*10^12)
EXAMPLE
91 is a term since the average of the divisors of 91 and 92 is the same. Indeed, (1+7+13+91)/4 = (1+2+4+23+46+92)/6.
MATHEMATICA
av[n_] := DivisorSigma[1, n]/DivisorSigma[0, n]; Select[Range[10^5], av[#] == av[# + 1] &]
SequencePosition[Table[DivisorSigma[1, n]/DivisorSigma[0, n], {n, 450000}], {x_, x_}][[All, 1]] (* Harvey P. Dale, Jun 01 2022 *)
PROG
(Python)
from sympy import divisors
from fractions import Fraction
def aupto(limit):
alst, prev_divavg = [], 1
for n in range(2, limit+2):
divs = divisors(n)
divavg = Fraction(sum(divs), len(divs))
if divavg == prev_divavg: alst.append(n-1)
prev_divavg = divavg
return alst
print(aupto(440013)) # Michael S. Branicky, May 14 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Resta, Feb 25 2014
STATUS
approved