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.

A103762 a(n) = least k with Sum_{j = n..k} 1/j >= 1.

Original entry on oeis.org

1, 4, 7, 10, 12, 15, 18, 20, 23, 26, 29, 31, 34, 37, 39, 42, 45, 48, 50, 53, 56, 58, 61, 64, 67, 69, 72, 75, 77, 80, 83, 86, 88, 91, 94, 97, 99, 102, 105, 107, 110, 113, 116, 118, 121, 124, 126, 129, 132, 135, 137, 140, 143, 145, 148, 151, 154, 156, 159, 162
Offset: 1

Views

Author

David W. Wilson, Apr 14 2008

Keywords

Comments

a(n) = A136617(n) + n for n > 1. Also a(n) = A136616(n-1) + 1 for n > 1.
If you compare this to floor(e*n) = A022843, 2,5,8,10,13,16,..., it appears that floor(e*n)-a(n) = 1,1,1,0,1,1,1,1,1,1,0,..., initially consisting of 0's and 1's. The places where the 0's occur are 4, 11, 18, 25, 32, 36, 43, 50, 57, 64, 71, ... whose differences seem to be 4, 7 or 11.
There are some rather sharp estimates on this type of differences between harmonic numbers in Theorem 3.2 of the Sintamarian reference, which may help to uncover such a pattern. - R. J. Mathar, Apr 15 2008
a(n) = round(e*(n-1/2)) with the exception of the terms of A277603; at those values of n, a(n) = round(e*(n-1/2)) + 1. - Jon E. Schoenfield, Apr 03 2018

Crossrefs

Cf. A136616, A136617, A242679 (Bobo numbers).

Programs

  • Mathematica
    i = 0; s = 0; Table[While[s < 1, i++; s = s + 1/i]; s = s - 1/n; i, {n, 100}] (* T. D. Noe, Jun 26 2012 *)
  • PARI
    default(realprecision, 10^5); e=exp(1);
    a(n) = if(n<2, 1, floor(e*n+(1-e)/2+(e-1/e)/(24*n-12))); \\ Jinyuan Wang, Mar 06 2020

A242679 Bobo numbers: Numbers k with the property that floor(e*k) = least m with Sum_{j=k..m} 1/j > 1.

Original entry on oeis.org

4, 11, 18, 25, 32, 36, 43, 50, 57, 64, 71, 75, 82, 89, 96, 103, 114, 121, 128, 135, 142, 146, 153, 160, 167, 174, 185, 192, 199, 206, 213, 217, 224, 231, 238, 245, 256, 263, 270, 277, 284, 288, 295, 302, 309, 316, 327, 334, 341, 348, 355, 359, 366, 373, 380, 387, 398, 405, 412, 419, 426, 430, 437, 444, 451, 458, 469, 476, 483, 490, 497
Offset: 1

Views

Author

Steven J. Kifowit, May 20 2014

Keywords

Comments

These are the numbers n for which A103762(n) = floor(e*n).
If frac(e*n) > (e-1)/2, then n is a Bobo number, but not every Bobo number has this property. The exceptions are in A277603.
In Bobo's article (see Bobo link), the Bobo numbers through 2105 are listed. There is a typo: the number 143 is given in place of the correct number 142.
These numbers are mentioned in the comments associated with A103762. Differences between consecutive Bobo numbers are indeed 4, 7, or 11. An elementary proof is given in the Clancy/Kifowit link.

Crossrefs

Programs

  • PARI
    is(n)=my(e=exp(1),s); if(frac(e*n)>(e-1)/2, return(1)); s=sum(j=n,e*n\1-1, 1/j); s<=1 && s+e*n\1>1 \\ Charles R Greathouse IV, Sep 17 2016
Showing 1-2 of 2 results.