login
G.f. A(x) satisfies (1 + 9*x*A(x)^3/A(x^9)^3)^3 = 1 + 27*x*A(x)^6/A(x^3)^6.
1

%I #13 Oct 23 2024 01:06:45

%S 1,3,27,321,4257,61155,924543,14492790,233421669,3839329521,

%T 64218149487,1088970739038,18678139900776,323477653749255,

%U 5648683685011425,99348712270884792,1758325597698951072,31292300465578800315,559640987226716260068,10052876829984340893552,181296262179794893505427

%N G.f. A(x) satisfies (1 + 9*x*A(x)^3/A(x^9)^3)^3 = 1 + 27*x*A(x)^6/A(x^3)^6.

%C Formula (2.6) given in the Borwein link implies the identity

%C (1 + 9*x*P(x)^3/P(x^9)^3)^3 = 1 + 27*x*P(x)^12/P(x^3)^12

%C for the partition function P(x) = 1/Product_{n>=1} (1 - x^n) (A000041).

%C For fixed integer m, the g.f. G(x) that satisfies (1 + 9*x*G(x)^3/G(x^9)^3)^3 = 1 + 27*x*G(x)^m/G(x^3)^m forms an integer power series in x only when m = 12 (the partition function P(x)), m = 6 (the g.f. A(x) of this sequence), or m = 0.

%H Paul D. Hanna, <a href="/A376529/b376529.txt">Table of n, a(n) for n = 0..730</a>

%H J. M. Borwein, P. B. Borwein and F. Garvan, <a href="http://dx.doi.org/10.1090/S0002-9947-1994-1243610-6">Some Cubic Modular Identities of Ramanujan</a>, Trans. Amer. Math. Soc. 343 (1994), p. 39, formula (2.6).

%e G.f.: A(x) = 1 + 3*x + 27*x^2 + 321*x^3 + 4257*x^4 + 61155*x^5 + 924543*x^6 + 14492790*x^7 + 233421669*x^8 + 3839329521*x^9 + ...

%e RELATED SERIES.

%e (1 + 9*x*A(x)^3/A(x^9)^3)^3 = 1 + 27*x + 486*x^2 + 8019*x^3 + 131706*x^4 + 2182626*x^5 + 36571014*x^6 + 619550127*x^7 + ...

%e which also equals 1 + 27*x*A(x)^6/A(x^3)^6.

%e A(x)^3/A(x^9)^3 = 1 + 9*x + 108*x^2 + 1476*x^3 + 21465*x^4 + 327321*x^5 + 5163804*x^6 + 83575719*x^7 + 1379999187*x^8 + ...

%e A(x)^6/A(x^3)^6 = 1 + 18*x + 297*x^2 + 4878*x^3 + 80838*x^4 + 1354482*x^5 + 22946301*x^6 + 392705496*x^7 + 6782271912*x^8 + ...

%o (PARI) {a(n) = my(V=[1], A); for(m=0, n, V = concat(V, 0); A = Ser(V);

%o V[#V] = polcoef( (1 + 9*x*Ser(A)^3/subst(Ser(A),x,x^9)^3 )^3 - 1 - 27*x*Ser(A)^6/subst(Ser(A),x,x^3)^6, #V)/81 ); V[n+1]}

%o for(n=0, 30, print1(a(n), ", "))

%K nonn,changed

%O 0,2

%A _Paul D. Hanna_, Oct 16 2024