login
A249936
G.f. satisfies: A(x) = Sum_{n>=0} x^n * Sum_{k=0..n} x^k * {[x^k] A(x)^(3*n)}.
3
1, 1, 4, 7, 49, 85, 565, 1016, 7350, 13195, 95564, 173088, 1275166, 2317299, 17166633, 31318180, 233446920, 427044441, 3195627268, 5859732238, 43994595982, 80828819606, 608411734754, 1119676074515, 8445900399153, 15565569301793, 117619945650292, 217042442850592, 1642503844147711, 3034229201413924
OFFSET
0,3
COMMENTS
Here [x^k] A(x)^(3*n) denotes the coefficient of x^k in A(x)^(3*n).
FORMULA
G.f. satisfies: A(x) = (1 + 3*x^2*G'(x^2)/G(x^2)) / (1 - x*G(x^2)^3), where A(x) = G(x/A(x)^3) and G(x) = A(x*G(x)^3) = ( (1/x)*Series_Reversion(x/A(x)^3) )^(1/3).
EXAMPLE
G.f.: A(x) = 1 + x + 4*x^2 + 7*x^3 + 49*x^4 + 85*x^5 + 565*x^6 + 1016*x^7 +...
Related expansions:
A(x)^3 = 1 + 3*x + 15*x^2 + 46*x^3 + 249*x^4 + 786*x^5 + 3907*x^6 +...
A(x)^6 = 1 + 6*x + 39*x^2 + 182*x^3 + 999*x^4 + 4446*x^5 + 22116*x^6 +...
A(x)^9 = 1 + 9*x + 72*x^2 + 435*x^3 + 2655*x^4 + 14247*x^5 + 77145*x^6 +...
A(x)^12 = 1 + 12*x + 114*x^2 + 832*x^3 + 5703*x^4 + 35076*x^5 + 208630*x^6 +...
GENERATING METHOD.
The initial terms, k=0..n, of the (3*n)-th power of g.f. A(x) begin:
n=0: [1];
n=3: [1, 3];
n=6: [1, 6, 39];
n=9: [1, 9, 72, 435];
n=12: [1, 12, 114, 832, 5703];
n=15: [1, 15, 165, 1400, 10710, 73683];
n=18: [1, 18, 225, 2166, 18324, 138924, 982281];
n=21: [1, 21, 294, 3157, 29274, 242004, 1847629, 13193967];
n=24: [1, 24, 372, 4400, 44370, 396720, 3251592, 24825192, 179183031]; ...
from which the antidiagonal sums form this sequence:
a(0) = 1;
a(1) = 1;
a(2) = 1 + 3 = 4;
a(3) = 1 + 6 = 7;
a(4) = 1 + 9 + 39 = 49;
a(5) = 1 + 12 + 72 = 85;
a(6) = 1 + 15 + 114 + 435 = 565;
a(7) = 1 + 18 + 165 + 832 = 1016; ...
ALTERNATE GENERATING METHOD.
Define G(x) such that G(x) = A(x*G(x)^3) = ( (1/x)*Series_Reversion(x/A(x)^3) )^(1/3):
G(x) = 1 + x + 7*x^2 + 55*x^3 + 548*x^4 + 5748*x^5 + 64357*x^6 + 745777*x^7 + 8903267*x^8 + 108597640*x^9 + 1348164588*x^10 +...
then A(x) = (1 + 3*x^2*G'(x^2)/G(x^2)) / (1 - x*G(x^2)^3).
Note that 1 + 3*x^2*G'(x^2)/G(x^2) begins:
1 + 3*x^2 + 39*x^4 + 435*x^6 + 5703*x^8 + 73683*x^10 + 982281*x^12 +...
where the coefficients form the main diagonal of the above triangle.
PROG
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0, n, x^m*sum(k=0, m, x^k*polcoeff((A+x*O(x^m))^(3*m), k))+x*O(x^n))); polcoeff(A, n)}
for(n=0, 40, print1(a(n), ", "))
(PARI) /* ALTERNATE GENERATING METHOD (faster) */
{a(n)=local(A=1+x, G=1); for(i=0, #binary(n)+1, G=(1/x*serreverse(x/A^3 +x^2*O(x^n)))^(1/3); A=(1+3*x^2*subst(G'/G, x, x^2))/(1-x*subst(G^3, x, x^2))); polcoeff(A, n)}
for(n=0, 40, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 24 2014
STATUS
approved