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”).

A247778
Least k such that e - (1 + 1/k)^k < 1/n.
7
1, 2, 4, 5, 6, 8, 9, 10, 12, 13, 15, 16, 17, 19, 20, 21, 23, 24, 25, 27, 28, 29, 31, 32, 34, 35, 36, 38, 39, 40, 42, 43, 44, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 63, 65, 66, 68, 69, 70, 72, 73, 74, 76, 77, 78, 80, 81, 82, 84, 85, 87, 88, 89
OFFSET
1,2
COMMENTS
a(n+1) - a(n) is in {1,2} for n >= 1.
LINKS
EXAMPLE
The values of e - (1 + 1/k)^k for k = 1..8 are approximately 0.718282, 0.468282, 0.347911, 0.276876, 0.229962, 0.196655, 0.171782, 0.152497, so that the first 6 terms of A247778 are 1,2,4,5,6,8.
MATHEMATICA
z = 600; f[n_] := f[n] = Select[Range[z], E - (1 + 1/#)^# < 1/n &, 1];
u = Flatten[Table[f[n], {n, 1, z}]] (*A247778*)
d1 = Flatten[Position[Differences[u], 1]] (*A247779*)
d2 = Flatten[Position[Differences[u], 2]] (*A247780*)
PROG
(PARI) a(n) = {k = 1; while ((exp(1) - (1 + 1/k)^k) >= 1/n, k++); k; } \\ Michel Marcus, Sep 25 2014
CROSSREFS
Sequence in context: A292636 A079709 A285496 * A003511 A059567 A006594
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Sep 23 2014
STATUS
approved