cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A353537 Numbers whose abundancy index is larger than Pi^2/6.

Original entry on oeis.org

4, 6, 8, 10, 12, 14, 16, 18, 20, 24, 28, 30, 32, 36, 40, 42, 44, 45, 48, 50, 52, 54, 56, 60, 63, 64, 66, 68, 70, 72, 75, 76, 78, 80, 84, 88, 90, 92, 96, 98, 100, 102, 104, 105, 108, 110, 112, 114, 116, 120, 124, 126, 128, 130, 132, 135, 136, 138, 140, 144, 148, 150
Offset: 1

Views

Author

Amiram Eldar, Apr 25 2022

Keywords

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.

Examples

			4 is a term since sigma(4)/4 = 7/4 = 1.75 > Pi^2/6 = 1.644...
		

References

  • Harold Davenport, Über numeri abundantes, Sitzungsberichte der Preußischen Akademie der Wissenschaften, phys.-math. Klasse, No. 6 (1933), pp. 830-837.

Crossrefs

Programs

  • Mathematica
    Select[Range[150], DivisorSigma[-1, #] > Pi^2/6 &]
  • PARI
    isok(k) = sigma(k)/k > Pi^2/6; \\ Michel Marcus, Apr 25 2022

A074920 Integers k such that phi(k) = 6k/Pi^2 rounded to nearest integer.

Original entry on oeis.org

1, 2, 3, 4, 33, 39, 99, 3237, 9711, 2069301, 6207903, 45502509, 24091595067
Offset: 1

Views

Author

Joseph L. Pe, Oct 01 2002

Keywords

Comments

The average value of phi(n) can be approximated for large n by 6n/Pi^2 (Tattersall, p. 162).

Examples

			phi(99) = 60 and 6*99/Pi^2 = 60.1848...., which rounds off to 60. Hence 99 is a term of the sequence.
		

References

  • James J. Tattersall, "Elementary Number Theory in Nine Chapters", Cambridge University Press, 2001.

Crossrefs

Cf. A059956 (6/Pi^2).

Programs

  • Mathematica
    s = Pi^2; l = {}; Do[ If[Abs[EulerPhi[n] - 6 n/ s] <= 0.5, l = Append[l, n]], {n, 1, 10^5}]; l

Extensions

a(10)-a(13) from Amiram Eldar, May 10 2019
Name edited by Michel Marcus, May 10 2019

A308045 Numbers k such that usigma(k) = round(zeta(2)/zeta(3)*k), where usigma(k) is the sum of unitary divisors of k (A034448).

Original entry on oeis.org

1, 2, 3, 4, 35, 44, 111, 123, 1105, 1900, 2920, 12452, 17889, 34200, 65067, 716148, 14134055, 179040201, 221709100, 221743300, 221766100, 221788900, 1120968741, 1272582040, 1441454511, 7339101375
Offset: 1

Views

Author

Amiram Eldar, May 10 2019

Keywords

Comments

The unitary version of A072355.
zeta(2)/zeta(3) is the asymptotic mean of the unitary abundancy index usigma(k)/k (A306633).
a(27) > 10^10.

Examples

			35 is in the sequence since usigma(35) = 48, and (zeta(2)/zeta(3)) * 35 = 47.895... has a round value of 48.
		

Crossrefs

Programs

  • Mathematica
    usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); meanAb = Zeta[2]/Zeta[3]; Select[Range[10^6], usigma[#] == Round[meanAb*#] &]
Showing 1-3 of 3 results.