|
%I M5168 #40 Jan 31 2022 01:32:07
%S 1,-24,852,-35744,1645794,-80415216,4094489992,-214888573248,
%T 11542515402255,-631467591949480,35063515239394764,
%U -1971043639046131296,111949770626330347638,-6414671157989386260432,370360217892318010055832,-21525284426246779936288192
%N Sequence of coefficients arising in connection with a rapidly converging series for Pi.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Georg Fischer, <a href="/A005149/b005149.txt">Table of n, a(n) for n = 1..100</a>
%H M. Newman and D. Shanks, <a href="http://dx.doi.org/10.1090/S0025-5718-1984-0725996-9">On a sequence arising in series for pi</a>, Math. Comp., 42 (1984), 199-217.
%H <a href="/index/Ph#Pi314">Index entries for sequences related to the number Pi</a>
%H <a href="/index/Res#revert">Index entries for reversions of series</a>
%F REVERT(A014103). a(n) = -(-1)^n * A195130(n).
%e x - 24*x^2 + 852*x^3 - 35744*x^4 + 1645794*x^5 - 80415216*x^6 + 4094489992*x^7 + ...
%t CoefficientList[InverseSeries[Series[x*Product[(1+x^k)^24, {k, 1, 16}], {x, 0, 16}], x], x] // Rest (* _Jean-François Alcover_, Mar 29 2011 *)
%t a[ n_] := If[ n < 1, 0, SeriesCoefficient[ InverseSeries[ Series[ q Product[ 1 + q^k, {k, n}]^24, {q, 0, n}], x], {x, 0, n}]] (* _Michael Somos_, Sep 11 2011 *)
%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. A005148, A014103, A195130.
%K sign,nice
%O 1,2
%A _Simon Plouffe_ and _N. J. A. Sloane_
%E Formula and more terms from _Michael Somos_, Nov 24 2001
%E Signs corrected, Dec 24 2001
|
