OFFSET
1,2
COMMENTS
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
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
E. R. Bobo, A sequence related to the harmonic series, College Math. J. 26 (1995), 308-310.
D. T. Clancy and S. J. Kifowit, A closer look at Bobo's sequence, College Math. J. 45 (2014), 199-206.
A. Sintamarian, A generalization of Euler's constant, Numer. Algor. 46 (2007), pp. 141-151.
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 *)
PROG
(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
CROSSREFS
KEYWORD
nonn
AUTHOR
David W. Wilson, Apr 14 2008
STATUS
approved