|
%I #21 Oct 24 2015 23:28:12
%S 1,24,852,35744,1645794,80415216,4094489992,214888573248,
%T 11542515402255,631467591949480,35063515239394764,1971043639046131296,
%U 111949770626330347638,6414671157989386260432,370360217892318010055832
%N Series inversion of A100130.
%H Michael Somos, <a href="/A195130/b195130.txt">Table of n, a(n) for n = 1..16</a>
%F A005149(n) = -(-1)^n * a(n).
%e G.f. = q + 24*q^2 + 852*q^3 + 35744*q^4 + 1645794*q^5 + 80415216*q^6 + ...
%t a[ n_] := If[ n < 1, 0, SeriesCoefficient[ InverseSeries[ Series[ q / QPochhammer[ -q, q^2]^24, {q, 0, n}], q], {q, 0, n}]];
%o (PARI) {a(n) = if( n<1, 0, polcoeff( serreverse( x * prod( k=1, n, 1 + (-x)^k, 1 + x * O(x^n) )^24), n))};
%Y Cf. A005149, A100130.
%K nonn
%O 1,2
%A _Michael Somos_, Sep 11 2011
|
