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-2 of 2 results.

A072355 Numbers k such that sigma(k) = (Pi^2)*(k/6) rounded off (where 0.5 is rounded to 0).

Original entry on oeis.org

2, 4, 22, 63, 4202, 4246, 444886, 1161238, 9362914, 26996486, 545614671, 1640386293, 2242930954, 2243031802, 2243065418, 2243115842, 18000691527
Offset: 1

Views

Author

Joseph L. Pe, Jul 18 2002

Keywords

Comments

In 1838 Dirichlet showed that the average value of sigma(n) is (Pi^2)*(n/6) for large n (see Tattersall).

Examples

			sigma(444886) = 731808 = (Pi^2 * 444886)/6 rounded off; so 444886 is a term of the sequence.
		

References

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

Crossrefs

Cf. A013661 (Pi^2/6), A074920.

Programs

  • Mathematica
    Select[Range[10^6], Round[(Pi^2 * #)/6] == DivisorSigma[1, # ] &]

Extensions

More terms from Robert G. Wilson v, Jul 27 2002
a(10)-a(17) from Giovanni Resta, Apr 03 2017

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-2 of 2 results.