OFFSET
0,5
COMMENTS
Related to A102573: letting T(q,r) be the coefficient of n^(r+1) in the polynomial 2^(q-n)/n times Sum_{k=0..n} binomial(n,k)*k^q, then A009453(x) equals (-1)^(x+1) times the real part of Sum_{k=0..x-1} T(x,k)*i^k, where i is the imaginary unit. See Mathematica code below. [John M. Campbell, Nov 17 2011]
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
FORMULA
E.g.f. sin(log(1 + tanh(x))).
a(n) = Sum_{m=0..(n-1)/2} (Sum_{r=2*m+1..n} (stirling1(r,2*m+1) * Sum_{k=r..n} binomial(k-1,r-1)*k!*2^(n-k)*stirling2(n,k)*(-1)^(r+k+m)/r!)). - Vladimir Kruchinin, Jun 21 2011
MATHEMATICA
poly[q_] := 2^(q - n)/n*FunctionExpand[Sum[Binomial[n, k]*k^q, {k, 0, n}]]; T[q_, r_] := First[Take[CoefficientList[poly[q], n], {r + 1, r + 1}]]; Table[Re[Sum[T[x, k]*I^k, {k, 0, x - 1}]]*(-1)^(x + 1), {x, 0, 23}] (* John M. Campbell, Nov 17 2011 *)
With[{nn = 30}, CoefficientList[Series[Sin[Log[1 + Tanh[x]]], {x, 0, nn}], x] Range[0, nn]!] (* Vincenzo Librandi, Feb 09 2014 *)
PROG
(Maxima) a(n):=sum(sum((stirling1(r, 2*m+1)*sum(binomial(k-1, r-1)*k!*2^(n-k)*stirling2(n, k)*(-1)^(r+k+m), k, r, n))/r!, r, 2*m+1, n), m, 0, (n-1)/2); /* Vladimir Kruchinin, Jun 21 2011 */
CROSSREFS
KEYWORD
sign,easy
AUTHOR
EXTENSIONS
Extended with signs by Olivier Gérard, Mar 15 1997
STATUS
approved