OFFSET
0,3
FORMULA
a(n)=2*sum(k=1..2*n, sum(r=0..2*n-k, (stirling1(r,k)*sum(j=1..r+k, ((sum(i=0..(j-1)/2, (j-2*i)^(2*n)*binomial(j,i)))*(-1)^(r+k+n-j)*binomial(r+k,j))/2^j))/(r)!)), n>0, a(0)=1.
MATHEMATICA
With[{nn=40}, Take[CoefficientList[Series[Cos[x]^(Cos[x]-1), {x, 0, nn}], x] Range[0, nn]!, {1, -1, 2}]] (* Harvey P. Dale, Mar 19 2018 *)
PROG
(Maxima)
a(n):=2*sum(sum((stirling1(r, k)*sum(((sum((j-2*i)^(2*n)*binomial(j, i), i, 0, (j-1)/2))*(-1)^(r+k+n-j)*binomial(r+k, j))/2^j, j, 1, r+k))/(r)!, r, 0, 2*n-k), k, 1, 2*n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, Jun 03 2011
STATUS
approved