OFFSET
1,5
COMMENTS
FORMULA
T(n,k) = [x^k] Sum_{i>=0} F_i(x)^n / 2^(i+1) where F_0(x)=x, F_{n+1}(x) = F_n(x+x^2); a sum involving n-th powers of self-compositions of x+x^2 (cf. A122888).
EXAMPLE
Table begins:
1, 1, 2, 7, 34, 214, 1652, 15121, 160110, 1925442, 25924260, ...;
1, 2, 7, 34, 214, 1652, 15121, 160110, 1925442, 25924260, ...;
1, 3, 15, 94, 726, 6645, 70361, 846144, 11392530, 169785124, ...;
1, 4, 26, 200, 1831, 19388, 233154, 3139200, 46784118, ...;
1, 5, 40, 365, 3865, 46481, 625820, 9326720, 152426170, ...;
1, 6, 57, 602, 7239, 97470, 1452610, 23739936, 422171622, ...;
1, 7, 77, 924, 12439, 185388, 3029782, 53879148, 1035760670, ...;
1, 8, 100, 1344, 20026, 327296, 5820360, 111889248, 2312153223, ...;
1, 9, 126, 1875, 30636, 544824, 10473576, 216432783, 4784414985, ...;
1, 10, 155, 2530, 44980, 864712, 17868995, 395007850, 9301284465, ...;
Given that A122940 begins:
[1, 1, 4, 17, 106, 796, 7176, 75057, 894100, 11946906, ...],
demonstrate A122940(n)/n = Sum_{m=1..n} (-1)^(m-1)*T(m,n-m+1)/m
at n=4: A122940(4)/4 = 17/4 = 7/1 - 7/2 + 3/3 - 1/4;
at n=5: A122940(5)/5 = 106/5 = 34/1 - 34/2 + 15/3 - 4/4 + 1/5;
at n=6: A122940(6)/6 = 796/6 = 214/1 - 214/2 + 94/3 - 26/4 + 5/5 - 1/6.
PROG
(PARI) /* Get T(n, k) from H(n, ), the n-th self-composition of x+x^2: */
{H(n, p)=local(F=x+x^2, G=x+x*O(x^p)); if(n==0, G=x, for(i=1, n, G=subst(F, x, G)); G)}
{T(n, k)=round(polcoeff( sum(i=0, 6*n+100, H(i, k+n-1)^n/2^(i+1)), k+n-1))}
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Paul D. Hanna, Sep 25 2006
STATUS
approved