OFFSET
0,2
COMMENTS
a(n) = sum over all multinomials M2(2*(n+1),k), k from {1..p(2*(n+1))} restricted to partitions with exactly two odd and any nonnegative number even parts. p(2*(n+1)) = A000041(2*(n+1)) (partition numbers) and for the M2-multinomial numbers in A-St order see A036039(2*(n+1),k). - Wolfdieter Lang, Aug 07 2007
FORMULA
E.g.f. (with alternating zeros): A(x) = (d^2/dx^2)a(x) with a(x):=(1/(sqrt(1-x^2))*(log(sqrt((1+x)/(1-x))))^2)/2!.
a(n) ~ log(2*n)^2 * 2^(2*n) * n^(2*n + 2) / (exp(2*n)) * (1 + (2*gamma + 6*log(2))/log(2*n) + (gamma^2 + 6*gamma*log(2) + 9*log(2)^2 - Pi^2/2) / log(2*n)^2), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Jul 21 2019
EXAMPLE
Multinomial representation for a(2): partitions of 2*3=6 with two odd parts: (1,5) with A-St position k=2; (3^2) with k=4; (1^2,4) with k=5; (1,2,3) with k=6 and (1^2,2^2) with k=9. The M2 numbers for these partitions are 144, 40, 90, 120, 45, adding up to 439 = a(2).
MATHEMATICA
nmax = 20; Table[(CoefficientList[Series[(4 + 8*x*Log[(1 + x)/(1 - x)] + (1/2 + x^2)*Log[(1 + x)/(1 - x)]^2)/(4*(1 - x^2)^(5/2)), {x, 0, 2*nmax}], x]*Range[0, 2*nmax]!)[[2*n + 1]], {n, 0, nmax}] (* Vaclav Kotesovec, Jul 21 2019 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Feb 24 2005
STATUS
approved