OFFSET
1,1
COMMENTS
The abundancy index of a number k is sigma(k)/k, where sigma is the sum of divisors function (A000203).
Pi^2/6 (A013661) is the asymptotic mean of the abundancy indices of the positive integers.
The least odd term is 45 and the least term that is coprime to 6 is 25025.
Davenport (1933) proved that sigma(k)/k possesses a continuous distribution function and that the asymptotic density of numbers with abundancy index that is larger than x exists for all x > 1 and is a continuous function of x. Therefore, this sequence has an asymptotic density.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 4, 41, 436, 4258, 42928, 428557, 4286145, 42864566, 428585795, 4286368677, 42861854540, ... Apparently, the asymptotic density is 0.4286... which means that the distribution of the abundancy indices is skewed with a positive nonparametric skew.
REFERENCES
Harold Davenport, Über numeri abundantes, Sitzungsberichte der Preußischen Akademie der Wissenschaften, phys.-math. Klasse, No. 6 (1933), pp. 830-837.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Wikipedia, Nonparametric skew.
EXAMPLE
4 is a term since sigma(4)/4 = 7/4 = 1.75 > Pi^2/6 = 1.644...
MATHEMATICA
Select[Range[150], DivisorSigma[-1, #] > Pi^2/6 &]
PROG
(PARI) isok(k) = sigma(k)/k > Pi^2/6; \\ Michel Marcus, Apr 25 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 25 2022
STATUS
approved