cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A094816 Triangle read by rows: T(n,k) are the coefficients of Charlier polynomials: A046716 transposed, for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 8, 6, 1, 1, 24, 29, 10, 1, 1, 89, 145, 75, 15, 1, 1, 415, 814, 545, 160, 21, 1, 1, 2372, 5243, 4179, 1575, 301, 28, 1, 1, 16072, 38618, 34860, 15659, 3836, 518, 36, 1, 1, 125673, 321690, 318926, 163191, 47775, 8274, 834, 45, 1, 1, 1112083, 2995011
Offset: 0

Views

Author

Philippe Deléham, Jun 12 2004

Keywords

Comments

The a-sequence for this Sheffer matrix is A027641(n)/A027642(n) (Bernoulli numbers) and the z-sequence is A130189(n)/ A130190(n). See the W. Lang link.
Take the lower triangular matrix in A049020 and invert it, then read by rows! - N. J. A. Sloane, Feb 07 2009
Exponential Riordan array [exp(x), log(1/(1-x))]. Equal to A007318*A132393. - Paul Barry, Apr 23 2009
A signed version of the triangle appears in [Gessel]. - Peter Bala, Aug 31 2012
T(n,k) is the number of permutations over all subsets of {1,2,...,n} (Cf. A000522) that have exactly k cycles. T(3,2) = 6: We permute the elements of the subsets {1,2}, {1,3}, {2,3}. Each has one permutation with 2 cycles. We permute the elements of {1,2,3} and there are three permutations that have 2 cycles. 3*1 + 1*3 = 6. - Geoffrey Critzer, Feb 24 2013
From Wolfdieter Lang, Jul 28 2017: (Start)
In Chihara's book the row polynomials (with rising powers) are the Charlier polynomials (-1)^n*C^(a)_n(-x), with a = -1, n >= 0. See p. 170, eq. (1.4).
In Ismail's book the present Charlier polynomials are denoted by C_n(-x;a=1) on p. 177, eq. (6.1.25). (End)
The triangle T(n,k) is a representative of the parametric family of triangles T(m,n,k), whose columns are the coefficients of the standard expansion of the function f(x) = (-log(1-x))^(k)*exp(-m*x)/k! for the case m=-1. See A381082. - Igor Victorovich Statsenko, Feb 14 2025

Examples

			From _Paul Barry_, Apr 23 2009: (Start)
Triangle begins
  1;
  1,     1;
  1,     3,     1;
  1,     8,     6,     1;
  1,    24,    29,    10,     1;
  1,    89,   145,    75,    15,    1;
  1,   415,   814,   545,   160,   21,   1;
  1,  2372,  5243,  4179,  1575,  301,  28,  1;
  1, 16072, 38618, 34860, 15659, 3836, 518, 36, 1;
Production matrix is
  1, 1;
  0, 2, 1;
  0, 1, 3,  1;
  0, 1, 3,  4,  1;
  0, 1, 4,  6,  5,  1;
  0, 1, 5, 10, 10,  6,  1;
  0, 1, 6, 15, 20, 15,  7,  1;
  0, 1, 7, 21, 35, 35, 21,  8, 1;
  0, 1, 8, 28, 56, 70, 56, 28, 9, 1; (End)
		

References

  • T. S. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, New York, London, Paris, 1978, Ch. VI, 1., pp. 170-172.
  • Classical and Quantum Orthogonal Polynomials in One Variable, Cambridge University Press, 2005, EMA, Vol. 98, p. 177.

Crossrefs

Columns k=0..4 give A000012, A002104, A381021, A381022, A381023.
Diagonals: A000012, A000217.
Row sums A000522, alternating row sums A024000.
KummerU(-n,1-n-x,z): this sequence (z=1), |A137346| (z=2), A327997 (z=3).

Programs

  • Maple
    A094816 := (n,k) -> (-1)^(n-k)*add(binomial(-j-1,-n-1)*Stirling1(j,k), j=0..n):
    seq(seq(A094816(n, k), k=0..n), n=0..9); # Peter Luschny, Apr 10 2016
  • Mathematica
    nn=10;f[list_]:=Select[list,#>0&];Map[f,Range[0,nn]!CoefficientList[Series[ Exp[x]/(1-x)^y,{x,0,nn}],{x,y}]]//Grid  (* Geoffrey Critzer, Feb 24 2013 *)
    Flatten[Table[(-1)^(n-k) Sum[Binomial[-j-1,-n-1] StirlingS1[j,k],{j,0,n}], {n,0,9},{k,0,n}]] (* Peter Luschny, Apr 10 2016 *)
    p[n_] := HypergeometricU[-n, 1 - n - x, 1];
    Table[CoefficientList[p[n], x], {n,0,9}] // Flatten (* Peter Luschny, Oct 27 2019 *)
  • PARI
    {T(n, k)= local(A); if( k<0 || k>n, 0, A = x * O(x^n); polcoeff( n! * polcoeff( exp(x + A) / (1 - x + A)^y, n), k))} /* Michael Somos, Nov 19 2006 */
    
  • Sage
    def a_row(n):
        s = sum(binomial(n,k)*rising_factorial(x,k) for k in (0..n))
        return expand(s).list()
    [a_row(n) for n in (0..9)] # Peter Luschny, Jun 28 2019

Formula

E.g.f.: exp(t)/(1-t)^x = Sum_{n>=0} C(x,n)*t^n/n!.
Sum_{k = 0..n} T(n, k)*x^k = C(x, n), Charlier polynomials; C(x, n)= A024000(n), A000012(n), A000522(n), A001339(n), A082030(n), A095000(n), A095177(n), A096307(n), A096341(n), A095722(n), A095740(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Feb 27 2013
T(n+1, k) = (n+1)*T(n, k) + T(n, k-1) - n*T(n-1, k) with T(0, 0) = 1, T(0, k) = 0 if k>0, T(n, k) = 0 if k<0.
PS*A008275*PS as infinite lower triangular matrices, where PS is a triangle with PS(n, k) = (-1)^k*A007318(n, k). PS = 1/PS. - Gerald McGarvey, Aug 20 2009
T(n,k) = (-1)^(n-k)*Sum_{j=0..n} C(-j-1, -n-1)*S1(j, k) where S1 are the signed Stirling numbers of the first kind. - Peter Luschny, Apr 10 2016
Absolute values T(n,k) of triangle (-1)^(n+k) T(n,k) where row n gives coefficients of x^k, 0 <= k <= n, in expansion of Sum_{k=0..n} binomial(n,k) (-1)^(n-k) x^{(k)}, where x^{(k)} := Product_{i=0..k-1} (x-i), k >= 1, and x^{(0)} := 1, the falling factorial powers. - Daniel Forgues, Oct 13 2019
From Peter Bala, Oct 23 2019: (Start)
The n-th row polynomial is
R(n, x) = Sum_{k = 0..n} (-1)^k*binomial(n, k)*k! * binomial(-x, k).
These polynomials occur in series acceleration formulas for the constant
1/e = n! * Sum_{k >= 0} (-1)^k/(k!*R(n,k)*R(n,k+1)), n >= 0. (cf. A068985, A094816 and A137346). (End)
R(n, x) = KummerU[-n, 1 - n - x, 1]. - Peter Luschny, Oct 27 2019
Sum_{j=0..m} (-1)^(m-j) * Bell(n+j) * T(m,j) = m! * Sum_{k=0..n} binomial(k,m) * Stirling2(n,k). - Vaclav Kotesovec, Aug 06 2021
From Natalia L. Skirrow, Jun 11 2025: (Start)
G.f.: 2F0(1,y;x/(1-x)) / (1-x).
Polynomial for the n-th row is R(n,y) = 2F0(-n,y;-1).
Falling g.f. for n-th row: Sum_{k=0..n} a(n,k)*(y)_k = [x^0] 2F0(1,-n;-1/x) * (1+log(1/(1-x)))^y = [x^n] e^x * Gamma(n+1,x) * (1+log(1/(1-x)))^y, where (y)_k = y!/(y-k)! denotes the falling factorial. (End)

A269953 Triangle read by rows: T(n, k) = Sum_{j=0..n} binomial(-j-1, -n-1)*S1(j, k) where S1 are the Stirling cycle numbers A132393.

Original entry on oeis.org

1, -1, 1, 1, -1, 1, -1, 2, 0, 1, 1, 0, 5, 2, 1, -1, 9, 15, 15, 5, 1, 1, 35, 94, 85, 40, 9, 1, -1, 230, 595, 609, 315, 91, 14, 1, 1, 1624, 4458, 4844, 2779, 924, 182, 20, 1, -1, 13209, 37590, 43238, 26817, 9975, 2310, 330, 27, 1
Offset: 0

Views

Author

Peter Luschny, Apr 12 2016

Keywords

Comments

Replacing the Stirling cycle numbers in the definition by the Stirling set numbers leads to A105794.
From Wolfdieter Lang, Jun 19 2017: (Start)
The triangle t(n, k) = (-1)^(n-k)*T(n, k) is the matrix product of P = A007318 (Pascal) and s1 = A048994 (signed Stirling1). This is Sheffer (exp(t), log(1+t)).
The present triangle T is therefore the Sheffer triangle (exp(-t), -log(1-t)). Note that P is Sheffer (exp(t), t) (of the Appell type). (End)
The triangle T(n,k) is a representative of the parametric family of triangles T(m,n,k), whose columns are the coefficients of the standard expansion of the function f(x) = (-log(1-x))^(k)*exp(-m*x)/k! for the case m=1. See A381082. - Igor Victorovich Statsenko, Feb 14 2025

Examples

			Triangle starts:
   1;
  -1,  1;
   1, -1,  1;
  -1,  2,  0,  1;
   1,  0,  5,  2,  1;
  -1,  9, 15, 15,  5,  1;
   1, 35, 94, 85, 40,  9,  1.
		

Crossrefs

Columns k=0..4 give A033999, A002741, A381064, A381065, A381066.
Cf. A000166 (row sums), A080956 (diag n,n-1).
KummerU(-n,1-n-x,z): this sequence (z=-1), A094816 (z=1), |A137346| (z=2), A327997 (z=3).

Programs

  • Maple
    A269953 := (n,k) -> add(binomial(-j-1,-n-1)*abs(Stirling1(j,k)), j=0..n):
    seq(print(seq(A269953(n, k), k=0..n)), n=0..9);
    # Alternative:
    egf := exp(-t)*(1-t)^(-x): ser := series(egf, t, 12): p := n -> coeff(ser, t, n):
    seq(n!*seq(coeff(p(n), x, k), k=0..n), n=0..9); # Peter Luschny, Oct 28 2019
  • Mathematica
    Flatten[Table[Sum[Binomial[-j-1,-n-1] Abs[StirlingS1[j,k]], {j,0,n}], {n,0,9},{k,0,n}]]
    (* Or: *)
    p [n_] := HypergeometricU[-n, 1 - n - x, -1];
    Table[CoefficientList[p[n], x], {n, 0, 9}] (* Peter Luschny, Oct 28 2019 *)

Formula

From Wolfdieter Lang, Jun 19 2017: (Start)
E.g.f. of row polynomials R(n, x) = Sum_{k=0..n} T(n,k)*x^k: exp(-t)/(1 - t)^x.
E.g.f. of column k sequence: exp(-x)*(-log(1-x))^k/k!, k >= 0. (End)
From Peter Bala, Oct 26 2019: (Start)
Let R(n, x) = (-1)^n*Sum_{k >= 0} binomial(n,k)*k!* binomial(-x,k) the n-th row polynomial of this triangle.
R(n, x) = c_n(-x;-1), where c_n(x;a) denotes the n-th Poisson Charlier polynomial.
The series representation e = Sum_{k >= 0} 1/k! is the case n = 0 of the more general result e = n!*Sum_{k >= 0} 1/(k!*R(n,k)*R(n,k+1)), n = 0,2,3,4,.... (End)
R(n, x) = KummerU(-n, 1-n-x, -1). - Peter Luschny, Oct 28 2019
Showing 1-2 of 2 results.