login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A004081
a(n) = n-th positive integer such that only one integer lies between exp(s(m)) and exp(s(m+1)), where s(m) = 1 + 1/2 + 1/3 + . . . + 1/m.
3
4, 8, 13, 17, 22, 26, 31, 36, 40, 45, 49, 54, 58, 63, 68, 72, 77, 81, 86, 90, 95, 99, 104, 109, 113, 118, 122, 127, 131, 136, 141, 145, 150, 154, 159, 163, 168, 173, 177, 182, 186, 191, 195, 200
OFFSET
1,1
COMMENTS
Equivalently, this is a list of integers m for which there is exactly one integer k satisfying 1+1/2+...+1/m < log(k) < 1+1/2+...+1/(m+1). - Clark Kimberling May 30 2013.
LINKS
FORMULA
Conjecturally, a(n) = floor(n/(2 - exp(g)) - 1/2 + exp(g)/(24n)), where g is the Euler-Mascheroni constant. - David W. Cantrell (DWCantrell(AT)sigmaxi.net), Aug 11 2007
MATHEMATICA
Flatten[Position[Differences[Table[Floor[E^Sum[1/k, {k, 1, m}]], {m, 1, 500}]], 1]] (* Clark Kimberling, May 30 2013 *)
PROG
(PARI) lista(n) = {old = 1; expo = exp(old); for (i=2, n, new = old + 1/i; expn = exp(new); if (floor(expn)==ceil(expo), print1(i-1, ", ")); old = new; expo = expn; ); } \\ Michel Marcus, Mar 21 2013
CROSSREFS
Sequence in context: A311934 A311935 A311936 * A127264 A311937 A311938
KEYWORD
nonn
STATUS
approved