OFFSET
0,1
LINKS
Robert Israel, Table of n, a(n) for n = 0..898
FORMULA
From Robert Israel, Dec 22 2015: (Start)
a(2*k) = 2*k! + Sum_{j=0..k-1} k!/j! = 2*k! + k*e*Gamma(k,1).
a(2*k+1) = 2*(k+1)! + Sum_{j=0..k-1} (k+1)!/j! = 2*(k+1)! + k*(k+1)*e*Gamma(k,1).
a(n) ~ (e+2)*(ceiling(n/2))!. (End)
MAPLE
A[0]:= 2:
for n from 0 to 14 do
A[2*n+1]:= (n+1)*A[2*n];
A[2*n+2]:= (n+1)+A[2*n+1];
od:
seq(A[i], i=0..30); # Robert Israel, Dec 22 2015
MATHEMATICA
a = {2}; Do[If[EvenQ@ Length@ a, AppendTo[a, Floor[Length[a]/2] Last@ a],
AppendTo[a, Last@ a + Floor[Length[a] /2]]], {k, 27}]; Rest@ a (* Michael De Vlieger, Dec 22 2015 *)
PROG
(PARI) A019465(n, a=2)={for(i=2, n+1, if(bittest(i, 0), a+=i\2, a*=i\2)); a} \\ M. F. Hasler, Feb 25 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Edited by M. F. Hasler, Feb 25 2018
STATUS
approved