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.

A079526 a(n) = floor( exp(H_n)*log(H_n) ) - sigma(n).

Original entry on oeis.org

-1, -2, -1, -2, 2, -2, 4, 0, 4, 2, 10, -3, 13, 6, 9, 4, 20, 2, 23, 4, 17, 16, 31, -3, 29, 21, 26, 13, 42, 3, 46, 18, 36, 32, 41, 1, 57, 38, 45, 14, 65, 14, 69, 32, 41, 51, 78, 5, 75, 42, 66, 43, 90, 27, 78, 33, 76, 70, 103, -2, 107, 76, 71, 51, 98, 41, 120, 65, 98, 53
Offset: 1

Views

Author

N. J. A. Sloane, Jan 22 2003

Keywords

Comments

M. Kaneko has shown that the Riemann hypothesis is equivalent to the assertion that a(n) > 0 for n > 60.

Crossrefs

H_n = sum of harmonic series (see A002387), sigma(n) = A000203.

Programs

  • Magma
    [Floor(Exp(HarmonicNumber(n))*Log(HarmonicNumber(n))) - DivisorSigma(1,n): n in [1..80]]; // G. C. Greubel, Jan 15 2019
    
  • Mathematica
    f[n_] := Floor[Exp[HarmonicNumber[n]]Log[HarmonicNumber[n]]] - DivisorSigma[1, n]; Array[f, 70] (* Robert G. Wilson v, Dec 17 2016 *)
  • PARI
    {h(n) = sum(k=1, n, 1/k)};
    vector(80, n, floor( exp(h(n))*log(h(n))) - sigma(n,1) ) \\ G. C. Greubel, Jan 15 2019
    
  • Sage
    [floor(exp(harmonic_number(n))*log(harmonic_number(n))) - sigma(n,1) for n in (1..80)] # G. C. Greubel, Jan 15 2019

A306348 Numbers k such that exp(H_k)*log(H_k) <= sigma(k), where H_k is the harmonic number.

Original entry on oeis.org

1, 2, 3, 4, 6, 12, 24, 60
Offset: 1

Views

Author

Seiichi Manyama, Feb 09 2019

Keywords

Comments

If the Riemann hypothesis is true, there are no more terms.

Examples

			Let b(n) = exp(H_{a(n)})*log(H_{a(n)}).
n | a(n) |    b(n)    | sigma(a(n))
--+------+------------+-------------
1 |   1  |   0        |      1
2 |   2  |   1.817... |      3
3 |   3  |   3.791... |      4
4 |   4  |   5.894... |      7
5 |   6  |  10.384... |     12
6 |  12  |  25.218... |     28
7 |  24  |  57.981... |     60
8 |  60  | 166.296... |    168
		

Crossrefs

Programs

  • Mathematica
    For[k = 1, True, k++, If[Exp[HarmonicNumber[k]] Log[HarmonicNumber[k]] <= DivisorSigma[1, k], Print[k]]] (* Jean-François Alcover, Feb 14 2019 *)
Showing 1-2 of 2 results.