A062160 Square array T(n,k) = (n^k - (-1)^k)/(n+1), n >= 0, k >= 0, read by falling antidiagonals.
0, 1, 0, -1, 1, 0, 1, 0, 1, 0, -1, 1, 1, 1, 0, 1, 0, 3, 2, 1, 0, -1, 1, 5, 7, 3, 1, 0, 1, 0, 11, 20, 13, 4, 1, 0, -1, 1, 21, 61, 51, 21, 5, 1, 0, 1, 0, 43, 182, 205, 104, 31, 6, 1, 0, -1, 1, 85, 547, 819, 521, 185, 43, 7, 1, 0, 1, 0, 171, 1640, 3277, 2604, 1111, 300, 57, 8, 1, 0, -1, 1, 341, 4921, 13107, 13021, 6665, 2101, 455, 73, 9, 1, 0
Offset: 0
Examples
From _Seiichi Manyama_, Apr 12 2019: (Start) Square array begins: 0, 1, -1, 1, -1, 1, -1, 1, ... 0, 1, 0, 1, 0, 1, 0, 1, ... 0, 1, 1, 3, 5, 11, 21, 43, ... 0, 1, 2, 7, 20, 61, 182, 547, ... 0, 1, 3, 13, 51, 205, 819, 3277, ... 0, 1, 4, 21, 104, 521, 2604, 13021, ... 0, 1, 5, 31, 185, 1111, 6665, 39991, ... 0, 1, 6, 43, 300, 2101, 14706, 102943, ... (End)
Links
- Seiichi Manyama, Antidiagonals n = 0..139, flattened
- M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Crossrefs
Programs
-
Maple
seq(print(seq((n^k - (-1)^k)/(n+1), k = 0..10)), n = 0..10); # Peter Bala, May 31 2024
-
Mathematica
T[n_,k_]:=(n^k - (-1)^k)/(n+1); Join[{0},Table[Reverse[Table[T[n-k,k],{k,0,n}]],{n,12}]]//Flatten (* Stefano Spezia, Feb 20 2024 *)
Formula
T(n, k) = n^(k-1) - n^(k-2) + n^(k-3) - ... + (-1)^(k-1) = n^(k-1) - T(n, k-1) = n*T(n, k-1) - (-1)^k = (n - 1)*T(n, k-1) + n*T(n, k-2) = round[n^k/(n+1)] for n > 1.
T(n, k) = (-1)^(k+1) * resultant( n*x + 1, (x^k-1)/(x-1) ). - Max Alekseyev, Sep 28 2021
G.f. of row n: x/((1+x) * (1-n*x)). - Seiichi Manyama, Apr 12 2019
E.g.f. of row n: (exp(n*x) - exp(-x))/(n+1). - Stefano Spezia, Feb 20 2024
From Peter Bala, May 31 2024: (Start)
Binomial transform of the m-th row: Sum_{k = 0..n} binomial(n, k)*T(m, k) = (m + 1)^(n-1) for n >= 1.
Let R(m, x) denote the g.f. of the m-th row of the square array. Then R(m_1, x) o R(m_2, x) = R(m_1 + m_2 + m_1*m_2, x), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A109502.
T(m_1 + m_2 + m_1*m_2, k) = Sum_{i = 0..k} Sum_{j = i..k} binomial(k, i)* binomial(k-i, j-i)*T(m_1, j)*T(m_2, k-i). (End)
Comments