a(0)=0 and for n>=1, a(n)n!*sum(k=1..n-1,((sum(i=0,floor((k-1)/2),(k-2*i)^(n-k)*binomial(k,i)))*(-1)^((n-k)/2)*((-1)^(n-k)+1))/(2^k*(n-k)!)/k*(-1)^(k-1))+(-1)^(n-1)*(n-1)!. - Vladimir Kruchinin, Apr 23 2011
MATHEMATICA
With[{nn=30}, CoefficientList[Series[Log[1+Cos[x]x], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Nov 11 2011 *)
PROG
(Maxima)
a(n) := n! *sum(((sum((k-2*i)^(n-k)*binomial(k, i), i, 0, floor((k-1)/2)))*(-1)^((n-k)/2)*((-1)^(n-k)+1))/(2^k*(n-k)!)/k*(-1)^(k-1), k, 1, n-1)+(-1)^(n-1)*(n-1)!; /* Vladimir Kruchinin, Apr 23 2011 */
(PARI) my(x='x+O('x^30)); concat(0, Vec(serlaplace(log(1+x*cos(x))))) \\ Michel Marcus, Oct 29 2022