%I #7 Oct 22 2012 12:41:13
%S 1,1,1,1,9,1,1,36,36,1,1,100,419,100,1,1,225,2699,2699,225,1,1,441,
%T 12138,35052,12138,441,1,1,784,42865,286206,286206,42865,784,1,1,1296,
%U 127191,1696820,3932898,1696820,127191,1296,1,1,2025,330903,7958563
%N G.f.: A(x,y) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^4*y^k] * x^n/n ) = Sum_{n>=0,k=0..n} T(n,k)*x^n*y^k, as a triangle of coefficients T(n,k) read by rows.
%C Compare g.f. to that of the following triangle variants:
%C * Pascal's: exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)*y^k] * x^n/n );
%C * Narayana: exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2*y^k] * x^n/n );
%C * A181143: exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^3*y^k] * x^n/n );
%C * A218115: exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^5*y^k] * x^n/n );
%C * A218116: exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^6*y^k] * x^n/n ).
%e G.f.: A(x,y) = 1 + (1+y)*x + (1+9*y+y^2)*x^2 + (1+36*y+36*y^2+y^3)*x^3 + (1+100*y+419*y^2+100*y^3+y^4)*x^4 +...
%e The logarithm of the g.f. equals the series:
%e log(A(x,y)) = (1 + y)*x
%e + (1 + 2^4*y + y^2)*x^2/2
%e + (1 + 3^4*y + 3^4*y^2 + y^3)*x^3/3
%e + (1 + 4^4*y + 6^4*y^2 + 4^4*y^3 + y^4)*x^4/4
%e + (1 + 5^4*y + 10^4*y^2 + 10^4*y^3 + 5^4*y^4 + y^5)*x^5/5 +...
%e Triangle begins:
%e 1;
%e 1, 1;
%e 1, 9, 1;
%e 1, 36, 36, 1;
%e 1, 100, 419, 100, 1;
%e 1, 225, 2699, 2699, 225, 1;
%e 1, 441, 12138, 35052, 12138, 441, 1;
%e 1, 784, 42865, 286206, 286206, 42865, 784, 1;
%e 1, 1296, 127191, 1696820, 3932898, 1696820, 127191, 1296, 1;
%e 1, 2025, 330903, 7958563, 36955542, 36955542, 7958563, 330903, 2025, 1;
%e 1, 3025, 776688, 31205941, 261852055, 525079969, 261852055, 31205941, 776688, 3025, 1; ...
%e Note that column 1 forms the sum of cubes (A000537), and forms the squares of the triangular numbers.
%e Inverse binomial transform of columns begins:
%e [1];
%e [1, 8, 19, 18, 6];
%e [1, 35, 348, 1549, 3713, 5154, 4161, 1818, 333];
%e [1, 99, 2500, 27254, 161793, 589819, 1409579, 2282850, 2529900, 1893972, 917349, 259854, 32726]; ...
%o (PARI) {T(n,k)=polcoeff(polcoeff(exp(sum(m=1,n,sum(j=0,m,binomial(m,j)^4*y^j)*x^m/m)+O(x^(n+1))),n,x),k,y)}
%o for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print(""))
%Y Cf. A000537 (column 1), A166992 (row sums), A166898 (antidiagonal sums), A218140.
%Y Cf. variants: A001263 (Narayana), A181143, A218115, A218116.
%K nonn,tabl
%O 0,5
%A _Paul D. Hanna_, Oct 13 2010