%I #21 Mar 12 2021 22:24:43
%S 1,8,28,64,134,288,568,1024,1809,3152,5316,8704,13990,22208,34696,
%T 53248,80724,121240,180068,264448,384940,556064,796760,1132544,
%U 1598789,2243056,3127360,4333568,5971922,8188096,11170160,15163392,20491033
%N McKay-Thompson series of class 12D for the Monster group.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%D D. Ford, J. McKay and S. P. Norton, More on replicable functions, Commun. Algebra 22, No. 13, 5175-5193 (1994).
%H Seiichi Manyama, <a href="/A101127/b101127.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Mat#McKay_Thompson">Index entries for McKay-Thompson series for Monster simple group</a>
%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/InfiniteProduct.html">Infinite Product</a>
%F Expansion of q^(1/3) * (eta(q^2)^2 / (eta(q) * eta(q^4)))^8 in powers of q.
%F Euler transform of period 4 sequence [8, -8, 8, 0, ...].
%F Given g.f. A(x), B(q) = A(q^3) / q satisfies 0 = f(B(q), B(q^2)) where f(u, v) = u*v*(u^3+v^3) -(u*v)^3 + 15*(u*v)^2 - 32*u*v + 16.
%F G.f.: (Product_{k>0} (1 + x^(2*k-1)))^8.
%F A007259(n) = (-1)^n * a(n). Convolution square of A112160.
%F a(n) ~ exp(2*Pi*sqrt(n/3)) / (2 * 3^(1/4) * n^(3/4)). - _Vaclav Kotesovec_, Aug 27 2015
%F Expansion of chi(x)^8 in powers of x where chi() is a Ramanujan theta function. - _Michael Somos_, Sep 12 2017
%F G.f.: exp(8*Sum_{k>=1} x^k/(k*(1 - (-x)^k))). - _Ilya Gutkovskiy_, Jun 07 2018
%e T12D = 1/q + 8*q^2 + 28*q^5 + 64*q^8 + 134*q^11 + 288*q^14 + 568*q^17 + ...
%t nmax = 50; CoefficientList[Series[Product[(1 + x^(2*k+1))^8, {k, 0, nmax}], {x, 0, nmax}], x] (* _Vaclav Kotesovec_, Aug 27 2015 *)
%t a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2]^8, {x, 0, n}]; (* _Michael Somos_, Sep 12 2017 *)
%o (PARI) {a(n) = my(A); if(n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 / (eta(x + A) * eta(x^4 + A)))^8, n))};
%o (PARI) {a(n) = my(A); if(n<0, 0, A = x * O(x^n); polcoeff( prod(k=1, (n+1)\2, 1 + x^(2*k-1), 1 + A)^8, n))};
%Y cf. A007259, A112160.
%K nonn
%O 0,2
%A _Michael Somos_, Dec 02 2004
|