OFFSET
0,4
FORMULA
A132616(n+1) = [x^n] (1+x)^(n^2) / A(x) for n>=0.
EXAMPLE
G.f.: A(x) = 1 + 4*x^3 + 90*x^4 + 2448*x^5 + 79716*x^6 + 3058740*x^7 + 135637242*x^8 + 6835557984*x^9 + 386119895256*x^10 + 24170805494868*x^11 + ...
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in (1+x)^((n-1)^2) / A(x) begins:
n=0: [1, 1, 0, -4, -94, -2538, -82148, -3137720, ...];
n=1: [1, 0, 0, -4, -90, -2448, -79700, -3058020, ...];
n=2: [1, 1, 0, -4, -94, -2538, -82148, -3137720, ...];
n=3: [1, 4, 6, 0, -105, -2832, -90048, -3391872, ...];
n=4: [1, 9, 36, 80, 0, -3276, -105224, -3871476, ...];
n=5: [1, 16, 120, 556, 1666, 0, -123900, -4673220, ...];
n=6: [1, 25, 300, 2296, 12460, 47232, 0, -5561820, ...];
n=7: [1, 36, 630, 7136, 58671, 368784, 1694704, 0, ...];
n=8: [1, 49, 1176, 18420, 211590, 1895322, 13604628, 73552752, 0, ...]; ...
in which the main diagonal is all zeros after the initial term, illustrating that [x^n] (1+x)^((n-1)^2) / A(x) = 0 for n>0.
RELATED SEQUENCES.
The secondary diagonal in the above table that begins
[1, 1, 6, 80, 1666, 47232, 1694704, 73552752, ...]
Related triangular matrix T = A132615 begins:
1;
1, 1;
1, 1, 1;
6, 3, 1, 1;
80, 25, 5, 1, 1;
1666, 378, 56, 7, 1, 1;
47232, 8460, 1020, 99, 9, 1, 1;
1694704, 252087, 26015, 2134, 154, 11, 1, 1;
73552752, 9392890, 855478, 61919, 3848, 221, 13, 1, 1; ...
in which row n equals row (n-1) of T^(2*n-1) followed by '1' for n > 0.
PROG
(PARI) {a(n) = my(A=[1]); for(i=1, n, A=concat(A, 0); m=#A; A[m] = Vec( (1+x +x*O(x^m))^((m-2)^2)/Ser(A) )[m] ); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, May 08 2018
STATUS
approved