OFFSET
1,4
LINKS
Paul D. Hanna, Table of n, a(n), n=1..100.
EXAMPLE
G.f.: A(x) = x - x^2 + x^3 - 2*x^4 + 4*x^5 - 10*x^6 + 26*x^7 +...
The table of coefficients in the iterations of g.f. A(x) begins:
[1, 0, 0, 0, 0, 0, 0, 0, 0,...];
[1, -1, 1, -2, 4, -10, 26, -72, 205, -590,...];
[1, -2, 4, -10, 28, -86, 282, -973, 3489, -12881,...];
[1, -3, 9, -30, 108, -414, 1668, -7003, 30426, -136026,...];
[1, -4, 16, -68, 304, -1420, 6884, -34462, 177414, -935984,...];
[1, -5, 25, -130, 700, -3890, 22230, -130250, 780475, -4772240,...];
[1, -6, 36, -222, 1404, -9090, 60126, -405587, 2785695, -19453005,...];
[1, -7, 49, -350, 2548, -18886, 142352, -1089893, 8467424,...];
[1, -8, 64, -520, 4288, -35864, 304008, -2609868, 22675676,...];
[1, -9, 81, -738, 6804, -63450, 598194, -5698812, 54834489,...];
[1, -10, 100, -1010, 10300, -106030, 1101410, -11541225,...]; ...
in which the antidiagonal sums yield [1,1,0,0,0,0,0,...].
PROG
(PARI) {a(n)=local(A=[1], F, G, s); for(j=2, n, A=concat(A, 0); F=x*Ser(A); for(k=1, #A, s=0; G=x; for(i=1, k, G=subst(G, x, F); s+=Vec(G)[k-i+1])); A[ #A]=-s); A[n]}
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Aug 14 2010
STATUS
approved