login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
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 (list; graph; refs; listen; history; text; internal format)
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

Clark Kimberling, Table of n, a(n) for n = 1..400

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

Adjacent sequences:  A004078 A004079 A004080 * A004082 A004083 A004084

KEYWORD

nonn

AUTHOR

Clark Kimberling

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 19 13:31 EDT 2021. Contains 347563 sequences. (Running on oeis4.)