OFFSET
1,2
COMMENTS
The sequence u is define recursively by u(n) = u(n-1) + u(n-2)/(n-2), with u(1) = 0 and u(2) = 1. Let d(n) = a(n+1) - a(n). It appears that d(n) is in {1,2} for n >= 1, that d(n+1) - d(n) is in {-1,0,1}, and that similar bounds hold for higher differences.
REFERENCES
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, p. 19.
LINKS
Clark Kimberling, Table of n, a(n) for n = 1..1000
EXAMPLE
Approximations for the first few terms |(n+1)/u(n+1) - e| and 1/n^n are shown here:
n ... |(n+1)/u(n+1)-e| .. 1/n^n
1 ... 0.7182818285 ...... 1
2 ... 0.28171817 ........ 0.25
3 ... 0.051615161 ....... 0.037037
4 ... 0.0089908988 ...... 0.00390625
5 ... 0.0013006963 ...... 0.00032000
a(2) = 3 because |4/u(4) - e| < 1/2^2 < |3/u(3) - e|.
MATHEMATICA
$RecursionLimit = Infinity; $MaxExtraPrecision = Infinity;
z = 500; u[1] = 0; u[2] = 1; u[n_] := u[n] = u[n - 1] + u[n - 2]/(n - 2);
f[n_] := f[n] = Select[Range[z], Abs[(# + 1)/u[# + 1] - E] < n^-n &, 1];
u = Flatten[Table[f[n], {n, 1, z}]] (* A247914 *)
w = Differences[u]
f1 = Flatten[Position[w, 1]] (* A247915 *)
f2 = Flatten[Position[w, 2]] (* A247916 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Sep 27 2014
STATUS
approved