OFFSET
1,3
EXAMPLE
G.f.: A(x) = x - x^2 + 2*x^3 - 6*x^4 + 20*x^5 - 80*x^6 + 348*x^7 +...
Let F^n(x) denote the n-th iteration of F(x) = x+x^2 with F^0(x)=x,
then the table of coefficients in A(F^n(x)), n>=0, begins:
[1, -1, 2, -6, 20, -80, 348, -1778, 9892, -64392, ...];
[1, 0, 0, -1, 2, -14, 44, -348, 1476, -14148, 73920, ...];
[1, 1, 0, -1, -2, -10, -24, -231, -654, -9276, -32456, ...];
[1, 2, 2, 0, -6, -26, -108, -570, -3216, -22622, -162596, ...];
[1, 3, 6, 8, 0, -54, -324, -1776, -10594, -71702, -540448, ...];
[1, 4, 12, 29, 50, 0, -616, -4846, -32686, -228926, -1749972, ...];
[1, 5, 20, 69, 202, 436, 0, -8629, -84140, -680298, -5508864, ...];
[1, 6, 30, 134, 538, 1880, 4912, 0, -143442, -1672428, -15821492, ...];
[1, 7, 42, 230, 1164, 5404, 22108, 68098, 0, -2762748, -37526484, ...];
[1, 8, 56, 363, 2210, 12646, 67092, 315784, 1122952, 0, -60534272, ..];
[1, 9, 72, 539, 3830, 25930, 166520, 997581, 5322126, 21488640, 0, ..]; ...
in which the main diagonal equals all zeros after the initial '1';
the lower triangular portion of the above table forms triangle A187005.
PROG
(PARI) {ITERATE(F, n, p)=local(G=x); for(i=1, n, G=subst(F, x, G+x*O(x^p))); G}
{a(n)=local(A=[1]); for(i=1, n, A=concat(A, 0); A[#A]=-Vec(subst(x*Ser(A), x, ITERATE(x+x^2, i, #A)))[#A]); A[n]}
CROSSREFS
KEYWORD
sign,changed
AUTHOR
Paul D. Hanna, Mar 02 2011
STATUS
approved