OFFSET
0,8
COMMENTS
Table 3.1 in Hopkins thesis is the same below the main diagonal. - F. Chapoton, Sep 04 2025
LINKS
Alois P. Heinz, Antidiagonals n = 0..140, flattened
Noureddine Chair, Explicit Computations for the Intersection Numbers on Grassmannians, and on the Space of Holomorphic Maps from CP¹ into G<sub>r</sub>(C<sup>n</sup>), arXiv:hep-th/9808170, 1998. (Cf. Table 4).
E. Hopkins, Free Complexes over the Exterior Algebra with Small Homology, Ph.D. thesis, University of Nebraska, 2021. (Cf. Table 3.1.)
G. Kreweras, Sur les éventails de segments, Cahiers du Bureau Universitaire de Recherche Opérationelle, Cahier no. 15, Paris, 1970, pp. 3-41.
G. Kreweras, Sur les éventails de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #15 (1970), 3-41. [Annotated scanned copy]
Genki Shibukawa, New identities for some symmetric polynomials and their applications, arXiv:1907.00334 [math.CA], 2019.
FORMULA
Let F(n, k) = Sum_{j=0..(n-2)/2} (-1)^j*binomial(n-1-j, j+1)*F(n, k-2-2*j) for k > 0; F(n, 0) = 1 and F(n, k) = 0 if k < 0. Then A(n, k) = F(n+1, 2*k). See [Shibukawa] and A309896. - Peter Luschny, Aug 21 2019
For n >= k >= 2, A(n, k) = binomial(n+2*k,k-2)*(n^2+3*n-2*k+2)/(k*(k-1)). - F. Chapoton, Sep 05 2025
EXAMPLE
The first few antidiagonals are:
1;
1, 0;
1, 1, 0;
1, 2, 1, 0;
1, 3, 4, 1, 0;
1, 4, 8, 8, 1, 0;
1, 5, 13, 21, 16, 1, 0;
1, 6, 19, 40, 55, 32, 1, 0;
1, 7, 26, 66, 121, 144, 64, 1, 0;
...
Square array starts:
[0] 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
[1] 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
[2] 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, ...
[3] 1, 3, 8, 21, 55, 144, 377, 987, 2584, 6765, 17711, ...
[4] 1, 4, 13, 40, 121, 364, 1093, 3280, 9841, 29524, 88573, ...
[5] 1, 5, 19, 66, 221, 728, 2380, 7753, 25213, 81927, 266110, ...
[6] 1, 6, 26, 100, 364, 1288, 4488, 15504, 53296, 182688, 625184, ...
[7] 1, 7, 34, 143, 560, 2108, 7752, 28101, 100947, 360526, 1282735, ...
[8] 1, 8, 43, 196, 820, 3264, 12597, 47652, 177859, 657800, 2417416, ...
[9] 1, 9, 53, 260, 1156, 4845, 19551, 76912, 297275, 1134705, 4292145, ...
MAPLE
F:= proc(n) option remember;
`if`(n<2, 1, expand(F(n-1)-t*F(n-2)))
end:
A:= (n, k)-> coeff(series(1/F(n+1), t, k+1), t, k):
seq(seq(A(d-k, k), k=0..d), d=0..12); # Alois P. Heinz, Jul 04 2015
MATHEMATICA
F[n_] := F[n] = If[n<2, 1, Expand[F[n-1] - t*F[n-2]]]; A[n_, k_] := SeriesCoefficient[1/F[n+1], { t, 0, k}]; Table[A[d-k, k], {d, 0, 12}, {k, 0, d}] // Flatten (* Jean-François Alcover, Feb 17 2016, after Alois P. Heinz *)
PROG
(SageMath)
@cached_function
def F(n, k):
if k < 0: return 0
if k == 0: return 1
return sum((-1)^j*binomial(n-1-j, j+1)*F(n, k-2-2*j) for j in (0..(n-2)/2))
def A(n, k): return F(n+1, 2*k)
print([A(n-k, k) for n in (0..11) for k in (0..n)]) # Peter Luschny, Aug 21 2019
(Python) # The lower triangular array computed by F. Chapoton's formula:
from math import comb as binomial
def T(n: int, k: int) -> int:
if k < 0: return 0
if k == 0: return 1
if k == 1: return n
return (binomial(n + 2*k, k - 2) * (n**2 + 3*n - 2*k + 2)) // (k * (k - 1))
for n in range(9): print([T(n, k) for k in range(n+1)]) # Peter Luschny, Sep 06 2025
CROSSREFS
KEYWORD
AUTHOR
N. J. A. Sloane, Jul 03 2015
EXTENSIONS
More terms from Alois P. Heinz, Jul 04 2015
STATUS
approved
