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.

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