OFFSET
0,4
FORMULA
a(n)=2*sum(m=1..n/2, (sum(j=0..m, binomial(2*j,j)*4^(m-2*j)*sum(i=0..j, (i-j)^(2*m)* binomial(2*j,i)*(-1)^(m+j-i))))*sum(r=2*m..n, (stirling1(r,2*m)*sum(k=r..n, binomial(k-1,r-1)*k!*2^(n-k)*stirling2(n,k)*(-1)^(r+k)))/r!)), n>0, a(0)=1.
MATHEMATICA
With[{nn=20}, CoefficientList[Series[Sec[Log[1+Tanh[x]]], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Oct 06 2011 *)
PROG
(Maxima)
a(n):=if n=0 then 1 else 2*sum((sum(binomial(2*j, j)*4^(m-2*j)*sum((i-j)^(2*m)*binomial(2*j, i)*(-1)^(m+j-i), i, 0, j), j, 0, m))*sum((stirling1(r, 2*m)*sum(binomial(k-1, r-1)*k!*2^(n-k)*stirling2(n, k)*(-1)^(r+k), k, r, n))/r!, r, 2*m, n), m, 1, n/2);
CROSSREFS
KEYWORD
sign
AUTHOR
Vladimir Kruchinin, Jun 21 2011
STATUS
approved