login
A376526
G.f. A(x) satisfies (4*A(x^9)/A(x) - 1)^3 = 28*A(x^3)^4/A(x)^4 - 1 + 4*x.
3
1, 1, 7, 101, 1182, 16303, 243618, 3794541, 61004863, 1004801055, 16865731633, 287434810434, 4960389320078, 86507430659736, 1522194102412522, 26991460147758051, 481824329994858705, 8651729550012408710, 156162244332280459203, 2831802125153075595781, 51565152204237125945253
OFFSET
0,3
COMMENTS
Compare to the following cubic modular identity of Jacobi theta functions,
(3*H(q^9)/H(q) - 1)^3 = 9*H(q^3)^4/H(q)^4 - 1,
which holds for H(q) = theta_2(q), H(q) = theta_3(q), and H(q) = theta_4(q) (see Mathworld link).
LINKS
Weisstein, Eric W., Modular Equation. From MathWorld -- A Wolfram Web Resource.
EXAMPLE
G.f.: A(x) = 1 + x + 7*x^2 + 101*x^3 + 1182*x^4 + 16303*x^5 + 243618*x^6 + 3794541*x^7 + 61004863*x^8 + 1004801055*x^9 + 16865731633*x^10 + ...
RELATED SERIES.
(4*A(x^9)/A(x) - 1)^3 = 27 - 108*x - 504*x^2 - 7840*x^3 - 73332*x^4 - 996240*x^5 - 14618520*x^6 + ...
28*A(x^3)^4 / A(x)^4 = 28 - 112*x - 504*x^2 - 7840*x^3 - 73332*x^4 - 996240*x^5 - 14618520*x^6 + ...
so that 28*A(x^3)^4/A(x)^4 - (4*A(x^9)/A(x) - 1)^3 = 1 - 4*x.
PROG
(PARI) {a(n) = my(V=[1], A); for(m=0, n, V = concat(V, 0); A = Ser(V);
V[#V] = polcoef( 28*subst(A, x, x^3)^4/A^4 - 1 + 4*x - (4*subst(A, x, x^9)/A - 1)^3, #V-1)/4 ); V[n+1]}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
Cf. A376525.
Sequence in context: A175345 A142358 A210684 * A367137 A367152 A357334
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 14 2024
STATUS
approved