OFFSET
1,2
FORMULA
a(1)=1 and for n>1, a(n)=floor((e-1/2)*n!)
MATHEMATICA
Join[{a = 1}, Table[a = Ceiling[n*(a + 1/a)], {n, 2, 30}]]
PROG
(PARI) a(n)=if(n<2, 1, floor((exp(1)-1/2)*n!))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, May 06 2003
STATUS
approved