OFFSET
0,18
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..5150 (first 101 antidiagonals)
FORMULA
G.f.: HadamardSquare(Product_{k=1..n} 1/(1 - y*x^k))/(1 - y), where HadamardSquare(g) is the termwise product of the series g with itself.
EXAMPLE
Array begins:
================================================
n\k | 0 1 2 3 4 5 6 7 8 9 10 ...
-----+------------------------------------------
0 | 1 1 1 1 1 1 1 1 1 1 1 ...
1 | 0 0 1 1 1 1 1 1 1 1 1 ...
2 | 0 0 1 3 4 4 4 4 4 4 4 ...
3 | 0 0 1 3 6 8 9 9 9 9 9 ...
4 | 0 0 1 5 11 17 22 24 25 25 25 ...
5 | 0 0 1 5 13 23 33 41 46 48 49 ...
6 | 0 0 1 7 22 44 67 87 103 113 118 ...
7 | 0 0 1 7 24 54 92 130 163 189 207 ...
8 | 0 0 1 9 35 85 156 234 307 367 414 ...
9 | 0 0 1 9 39 107 214 344 478 598 697 ...
10 | 0 0 1 11 52 150 318 542 791 1031 1240 ...
...
The T(3,4) = 6 sequences are (-3,1,1,1), (-3,0,1,2), (-3,0,0,3), (-2,-1,1,2), (-2,-1,0,3), (-1,-1,-1,3).
PROG
(PARI)
Mtx(n, m=n)=my(y='y+O('y^(1+m)), g=1/prod(k=1, n, 1-y*x^k, 1 + O(x*x^n))); Mat([Vec(p+O('y^(1+m)), -m-1) | p<-Col(serconvol(g, g)/(1-y))]~)
{ my(A=Mtx(10, 10)); for(i=1, #A~, print(A[i, ])) }
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Andrew Howroyd, Oct 10 2025
STATUS
approved
