A142979 -id:A142979 - cp's OEIS frontend

A002162 Decimal expansion of the natural logarithm of 2.

6, 9, 3, 1, 4, 7, 1, 8, 0, 5, 5, 9, 9, 4, 5, 3, 0, 9, 4, 1, 7, 2, 3, 2, 1, 2, 1, 4, 5, 8, 1, 7, 6, 5, 6, 8, 0, 7, 5, 5, 0, 0, 1, 3, 4, 3, 6, 0, 2, 5, 5, 2, 5, 4, 1, 2, 0, 6, 8, 0, 0, 0, 9, 4, 9, 3, 3, 9, 3, 6, 2, 1, 9, 6, 9, 6, 9, 4, 7, 1, 5, 6, 0, 5, 8, 6, 3, 3, 2, 6, 9, 9, 6, 4, 1, 8, 6, 8, 7

Offset: 0

N. J. A. Sloane

Newton calculated the first 16 terms of this sequence.
Area bounded by y = tan x, y = cot x, y = 0. - Clark Kimberling, Jun 26 2020
Choose four values independently and uniformly at random from the unit interval [0,1]. Sort them, and label them a,b,c,d from least to greatest (so that a b^2+c^2. - Akiva Weinberger, Dec 02 2024
Define the trihyperboloid to be the intersection of the three solid hyperboloids x^2+y^2-z^2<1, x^2-y^2+z^2<1, and -x^2+y^2+z^2<1. This fits perfectly within the cube [-1,1]^3. Then this is the ratio of the volume of the trihyperboloid to its bounding cube. - Akiva Weinberger, Dec 02 2024

			0.693147180559945309417232121458176568075500134360255254120680009493393...

G. Boros and V. H. Moll, Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals, Cambridge University Press, 2004.
Calvin C. Clawson, Mathematical Mysteries: The Beauty and Magic of Numbers, Springer, 2013. See p. 227.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 24, 250.
Steven R. Finch, Mathematical Constants, Cambridge, 2003, Sections 1.3.3, 2.21, 6.2, and 7.2.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 25 and appendix A, equations 25:14:3 and A:7:3 at pages 232, 670.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 29.

Harry J. Smith, Table of n, a(n) for n = 0..20000
D. H. Bailey and J. M. Borwein, Experimental Mathematics: Examples, Methods and Implications, Notices of the AMS, May 2005, Volume 52, Issue 5.
Peter Bala, New series for old functions
J. M. Borwein, P. B. Borwein, and K. Dilcher, Pi, Euler numbers and asymptotic expansions, Amer. Math. Monthly, 96 (1989), 681-687.
Paul Cooijmans, Odds.
X. Gourdon and P. Sebah, The logarithm constant:log(2), 2001.
X. Gourdon and P. Sebah, The logarithm constant:log(2), 2004. (Revised and extended version of above article.)
M. Kontsevich and D. Zagier, Periods, pp. 4-5.
Mathematical Reflections, Solution to Problem U376, Issue 4, 2016, p 17.
Isaac Newton, The method of fluxions and infinite series; with its application to the geometry of curve-lines, 1736; see p. 96.
Michael Penn, an alternating floor sum., YouTube video, 2020.
Michael Penn, A wonderful limit from the 2011 Virginia Tech Regional Math Competition, YouTube video, 2022.
Simon Plouffe, log(2), natural logarithm of 2 to 2000 places.
S. Ramanujan, Question 260, J. Ind. Math. Soc., III, p. 43.
Albert Stadler, Problem 3567, Crux Mathematicorum, Vol. 36 (Oct. 2010), p. 396; Oliver Geupel, Solution, Crux Mathematicorum, Vol. 37 (Oct. 2011), pp. 400-401.
D. W. Sweeney, On the computation of Euler's constant, Math. Comp., 17 (1963), 170-178.
Horace S. Uhler, Recalculation and extension of the modulus and of the logarithms of 2, 3, 5, 7 and 17, Proc. Nat. Acad. Sci. U. S. A. 26, (1940). 205-212.
Eric Weisstein's World of Mathematics, Natural Logarithm of 2, Masser-Gramain Constant, Logarithmic Integral, Dirichlet Eta Function.
Wikipedia, Natural logarithm of 2.
Wikipedia, Period (algebraic geometry).
Index entries for transcendental numbers.

Cf. A016730 (continued fraction), A002939, A008288, A142979, A142992.

RealDigits[N[Log[2],200]][[1]] (* Vladimir Joseph Stephan Orlovsky, Feb 21 2011 *)
RealDigits[Log[2],10,120][[1]] (* Harvey P. Dale, Jan 25 2024 *)

{ default(realprecision, 20080); x=10*log(2); for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b002162.txt", n, " ", d)); } \\ Harry J. Smith, Apr 21 2009

log(2) = Sum_{k>=1} 1/(k*2^k) = Sum_{j>=1} (-1)^(j+1)/j.
log(2) = Integral_{t=0..1} dt/(1+t).
log(2) = (2/3) * (1 + Sum_{k>=1} 2/((4*k)^3-4*k)) (Ramanujan).
log(2) = 4*Sum_{k>=0} (3-2*sqrt(2))^(2*k+1)/(2*k+1) (Y. Luke). - R. J. Mathar, Jul 13 2006
log(2) = 1 - (1/2)*Sum_{k>=1} 1/(k*(2*k+1)). - Jaume Oliver Lafont, Jan 06 2009, Jan 08 2009
log(2) = 4*Sum_{k>=0} 1/((4*k+1)*(4*k+2)*(4*k+3)). - Jaume Oliver Lafont, Jan 08 2009
log(2) = 7/12 + 24*Sum_{k>=1} 1/(A052787(k+4)*A000079(k)). - R. J. Mathar, Jan 23 2009
From Alexander R. Povolotsky, Jul 04 2009: (Start)
log(2) = (1/4)*(3 - Sum_{n>=1} 1/(n*(n+1)*(2*n+1))).
log(2) = (230166911/9240 - Sum_{k>=1} (1/2)^k*(11/k + 10/(k+1) + 9/(k+2) + 8/(k+3) + 7/(k+4) + 6/(k+5) - 6/(k+7) - 7/(k+8) - 8/(k+9) - 9/(k+10) - 10/(k+11)))/35917. (End)
log(2) = A052882/A000670. - Mats Granvik, Aug 10 2009
From log(1-x-x^2) at x=1/2, log(2) = (1/2)*Sum_{k>=1} L(k)/(k*2^k), where L(n) is the n-th Lucas number (A000032). - Jaume Oliver Lafont, Oct 24 2009
log(2) = Sum_{k>=1} 1/(cos(k*Pi/3)*k*2^k) (cf. A176900). - Jaume Oliver Lafont, Apr 29 2010
log(2) = (Sum_{n>=1} 1/(n^2*(n+1)^2*(2*n+1)) + 11)/16. - Alexander R. Povolotsky, Jan 13 2011
log(2) = ((Sum_{n>=1} (2*n+1)/(Sum_{k=1..n} k^2)^2)+396)/576. - Alexander R. Povolotsky, Jan 14 2011
From Alexander R. Povolotsky, Dec 16 2008: (Start)
log(2) = 105*(319/44100 - Sum_{n>=1} 1/(2*n*(2*n+1)*(2*n+3)*(2*n+5)*(2*n+7))).
log(2) = 319/420 - (3/2)*Sum_{n>=1} 1/(6*n^2+39*n+63). (End)
log(2) = Sum_{k>=1} A191907(2,k)/k. - Mats Granvik, Jun 19 2011
log(2) = Integral_{x=0..oo} 1/(1 + e^x) dx. - Jean-François Alcover, Mar 21 2013
log(2) = lim_{s->1} zeta(s)*(1-1/2^(s-1)). - Mats Granvik, Jun 18 2013
From Peter Bala, Dec 10 2013: (Start)
log(2) = 2*Sum_{n>=1} 1/( n*A008288(n-1,n-1)*A008288(n,n) ), a result due to Burnside.
log(2) = (1/3)*Sum_{n >= 0} (5*n+4)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(1/2)^n = (1/12)*Sum_{n >= 0} (28*n+17)/( (3*n+1)*(3*n+2)*C(3*n,n) )*(-1/4)^n.
log(2) = (3/16)*Sum_{n >= 0} (14*n+11)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(1/4)^n = (1/12)*Sum_{n >= 0} (34*n+25)/( (4*n+1)*(4*n+3)*C(4*n,2*n) )*(-1/18)^n. For more series of this type see the Bala link.
See A142979 for series acceleration formulas for log(2) obtained from the Mercator series log(2) = Sum_{n >= 1} (-1)^(n+1)/n. See A142992 for series for log(2) related to the root lattice C_n. (End)
log(2) = lim_{n->oo} Sum_{k=2^n..2^(n+1)-1} 1/k. - Richard R. Forberg, Aug 16 2014
From Peter Bala, Feb 03 2015: (Start)
log(2) = (2/3)*Sum_{k >= 0} 1/((2*k + 1)*9^k).
Define a pair of integer sequences A(n) = 9^n*(2*n + 1)!/n! and B(n) = A(n)*Sum_{k = 0..n} 1/((2*k + 1)*9^k). Both satisfy the same second-order recurrence equation u(n) = (40*n + 16)*u(n-1) - 36*(2*n - 1)^2*u(n-2). From this observation we obtain the continued fraction expansion log(2) = (2/3)*(1 + 2/(54 - 36*3^2/(96 - 36*5^2/(136 - ... - 36*(2*n - 1)^2/((40*n + 16) - ... ))))). Cf. A002391, A073000 and A105531 for similar expansions. (End)
log(2) = Sum_{n>=1} (Zeta(2*n)-1)/n. - Vaclav Kotesovec, Dec 11 2015
From Peter Bala, Oct 30 2016: (Start)
Asymptotic expansions:
for N even, log(2) - Sum_{k = 1..N/2} (-1)^(k-1)/k ~ (-1)^(N/2)*(1/N - 1/N^2 + 2/N^4 - 16/N^6 + 272/N^8 - ...), where the sequence of unsigned coefficients [1, 1, 2, 16, 272, ...] is A000182 with an extra initial term of 1. See Borwein et al., Theorem 1 (b);
for N odd, log(2) - Sum_{k = 1..(N-1)/2} (-1)^(k-1)/k ~ (-1)^((N-1)/2)*(1/N - 1/N^3 + 5/N^5 - 61/N^7 + 1385/N^9 - ...), by Borwein et al., Lemma 2 with f(x) := 1/(x + 1/2), h := 1/2 and then set x = (N - 1)/2, where the sequence of unsigned coefficients [1, 1, 5, 61, 1385, ...] is A000364. (End)
log(2) = lim_{n->oo} Sum_{k=1..n} sin(1/(n+k)). See Mathematical Reflections link. - Michel Marcus, Jan 07 2017
log(2) = Sum_{n>=1} A006519(n) / ((1 + 2^A006519(n)) * A000265(n) * (1 + A000265(n))). - Nicolas Nagel, Mar 19 2018
From Amiram Eldar, Jul 02 2020: (Start)
Equals Sum_{k>=2} zeta(k)/2^k.
Equals -Sum_{k>=2} log(1 - 1/k^2).
Equals Sum_{k>=1} 1/A002939(k).
Equals Integral_{x=0..Pi/3} tan(x) dx. (End)
log(2) = Integral_{x=0..Pi/2} (sec(x) - tan(x)) dx. - Clark Kimberling, Jul 08 2020
From Peter Bala, Nov 14 2020: (Start)
log(2) = Integral_{x = 0..1} (x - 1)/log(x) dx (Boros and Moll, p. 97).
log(2) = (1/2)*Integral_{x = 0..1} (x + 2)*(x - 1)^2/log(x)^2 dx.
log(2) = (1/4)*Integral_{x = 0..1} (x^2 + 3*x + 4)*(x - 1)^3/log(x)^3 dx. (End)
log(2) = 2*arcsinh(sqrt(2)/4) = 2*sqrt(2)*Sum_{n >= 0} (-1)^n*C(2*n,n)/ ((8*n+4)*32^n) = 3*Sum_{n >= 0} (-1)^n/((8*n+4)*(2^n)*C(2*n,n)). - Peter Bala, Jan 14 2022
log(2) = Integral_{x=0..oo} ( e^(-x) * (1-e^(-2x)) * (1-e^(-4x)) * (1-e^(-6x)) ) / ( x * (1-e^(-14x)) ) dx (see Crux Mathematicorum link). - Bernard Schott, Jul 11 2022
From Peter Bala, Oct 22 2023: (Start)
log(2) = 23/32 + 2!^3/16 * Sum_{n >= 1} (-1)^n * (n + 1)/(n*(n + 1)*(n + 2))^2 = 707/1024 - 4!^3/(16^2 * 2!^2) * Sum_{n >= 1} (-1)^n * (n + 2)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4))^2 = 42611/61440 + 6!^3/(16^3 * 3!^2) * Sum_{n >= 1} (-1)^n * (n + 3)/(n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*(n + 6))^2.
More generally, it appears that for k >= 0, log(2) = c(k) + (2*k)!^3/(16^k * k!^2) * Sum_{n >= 1} (-1)^(n+k+1) * (n + k)/(n*(n + 1)*...*(n + 2*k))^2 , where c(k) is a rational approximation to log(2). The first few values of c(k) are [0, 23/32, 707/1024, 42611/61440, 38154331/55050240, 76317139/110100480, 26863086823/38755368960, ...].
Let P(n,k) = n*(n + 1)*...*(n + k).
Conjecture: for k >= 0 and r even with r - 1 <= k, the series Sum_{n >= 1} (-1)^n * (d/dn)^r (P(n,k)) / (P(n,k)^2  = A(r,k)*log(2) + B(r,k), where A(r,k) and B(r,k) are both rational numbers. (End)
From Peter Bala, Nov 13 2023: (Start)
log(2) = 5/8 + (1/8)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)^2 / ( k*(k + 1) )^4
  = 257/384 + (3!^5/2^9)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)^2*(2*k + 5) / ( k*(k + 1)*(k + 2)*(k + 3) )^4
  = 267515/393216 + (5!^5/2^19)*Sum_{k >= 1} (-1)^(k+1) * (2*k + 1)*(2*k + 3)*(2*k + 5)^2*(2*k + 7)*(2*k + 9) / ( k*(k + 1)*(k + 2)*(k + 3)*(k + 4)*(k + 5) )^4
log(2) = 3/4 - 1/128 * Sum_{k >= 0} (-1/16)^k * (10*k + 12)*binomial(2*k+2,k+1)/ ((k + 1)*(2*k + 3)). The terms of the series are O(1/(k^(3/2)*4^n)). (End)
log(2) = eta(1) is a period, where eta(x) is the Dirichlet eta function. - Andrea Pinos, Mar 19 2024
log(2) = K_{n>=0} (n^2 + [n=0])/1, where K is the Gauss notation for an infinite continued fraction. In the expanded form, log(2) = 1/(1 + 1/(1 + 4/(1 + 9/1 + 16/(1 + 25/(1 + ... (see Clawson at p. 227). - Stefano Spezia, Jul 01 2024
log(2) = lim_{n->oo} Sum_{k=1..n} 1/(n + k) = lim_{x->0} (2^x - 1)/x = lim_{x->0} (2^x - 2^(-x))/(2*x) (see Finch). - Stefano Spezia, Oct 19 2024
From Colin Linzer, Nov 08 2024: (Start)
log(2) = Integral_{t=0...oo} (1 - tanh(t)) dt.
log(2) = Integral_{t=0...1} arctanh(t) dt.
log(2) = (1/2) * Integral_{t=-1...1} |arctanh(t)| dt. (End)
log(2) = 1 + Sum_{n >= 1} (-1)^n/(n*(4*n^2 - 1)) = 1/2 + (1/2)*Sum_{n >= 1} 1/(n*(4*n^2 - 1)). - Peter Bala, Jan 07 2025
log(2) = Integral_{x=0..1} Integral_{y=0..1} 1/((1 - x*y)*(1 + x)*(1 + y)) dy dx. - Kritsada Moomuang, May 22 2025

A008288 Square array of Delannoy numbers D(i,j) (i >= 0, j >= 0) read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 5, 1, 1, 7, 13, 7, 1, 1, 9, 25, 25, 9, 1, 1, 11, 41, 63, 41, 11, 1, 1, 13, 61, 129, 129, 61, 13, 1, 1, 15, 85, 231, 321, 231, 85, 15, 1, 1, 17, 113, 377, 681, 681, 377, 113, 17, 1, 1, 19, 145, 575, 1289, 1683, 1289, 575, 145, 19, 1, 1, 21, 181, 833, 2241, 3653, 3653

Offset: 0

N. J. A. Sloane

In the Formula section, some contributors use T(n,k) = D(n-k, k) (for 0 <= k <= n), which is the triangular version of the square array (D(n,k): n,k >= 0). Conversely, D(n,k) = T(n+k,k) for n,k >= 0. - Petros Hadjicostas, Aug 05 2020
Also called the tribonacci triangle [Alladi and Hoggatt (1977)]. - N. J. A. Sloane, Mar 23 2014
D(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (1,0), (0,1), (1,1). - Joerg Arndt, Jul 01 2011 [Corrected by N. J. A. Sloane, May 30 2020]
Or, triangle read by rows of coefficients of polynomials P[n](x) defined by P[0] = 1, P[1] = x+1; for n >= 2, P[n] = (x+1)*P[n-1] + x*P[n-2].
D(n, k) is the number of k-matchings of a comb-like graph with n+k teeth. Example: D(1, 3) = 7 because the graph consisting of a horizontal path ABCD and the teeth Aa, Bb, Cc, Dd has seven 3-matchings: four triples of three teeth and the three triples {Aa, Bb, CD}, {Aa, Dd, BC}, {Cc, Dd, AB}. Also D(3, 1)=7, the 1-matchings of the same graph being the seven edges: {AB}, {BC}, {CD}, {Aa}, {Bb}, {Cc}, {Dd}. - Emeric Deutsch, Jul 01 2002
Sum of n-th antidiagonal of the array D is A000129(n+1). - Reinhard Zumkeller, Dec 03 2004 [Edited by Petros Hadjicostas, Aug 05 2020 so that the counting of antidiagonals of D starts at n = 0. That is, the sum of the terms in the n-th row of the triangles T is  A000129(n+1).]
The A-sequence for this Riordan type triangle (see one of Paul Barry's comments under Formula) is A112478 and the Z-sequence the trivial: {1, 0, 0, 0, ...}. See the W. Lang link under A006232 for Sheffer a- and z-sequences where also Riordan A- and Z-sequences are explained. This leads to the recurrence for the triangle given below. - Wolfdieter Lang, Jan 21 2008
The triangle or chess sums, see A180662 for their definitions, link the Delannoy numbers with twelve different sequences, see the crossrefs. All sums come in pairs due to the symmetrical nature of this triangle. The knight sums Kn14 and Kn15 have been added. It is remarkable that all knight sums are related to the tribonacci numbers, that is, A000073 and A001590, but none of the others. - Johannes W. Meijer, Sep 22 2010
This sequence, A008288, is jointly generated with A035607 as an array of coefficients of polynomials u(n,x): initially, u(1,x) = v(1,x) = 1; for n > 1, u(n,x) = x*u(n-1,x) + v(n-1) and v(n,x) = 2*x*u(n-1,x) + v(n-1,x). See the Mathematica section. - Clark Kimberling, Mar 09 2012
Row n, for n > 0, of Roger L. Bagula's triangle in the Example section shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(n)*x^n which is the numerator of the n-th convergent of the continued fraction [k, k, k, ...], where k = sqrt(x) + 1/sqrt(x); see A230000. - Clark Kimberling, Nov 13 2013
In an n-dimensional hypercube lattice, D(n,k) gives the number of nodes situated at a Minkowski (Manhattan) distance of k from a given node. In cellular automata theory, the cells at Manhattan distance k are called the von Neumann neighborhood of radius k. For k=1, see A005843. - Dmitry Zaitsev, Dec 10 2015
These numbers appear as the coefficients of series relating spherical and bispherical harmonics, in the solutions of Laplace's equation in 3D. [Majic 2019, Eq. 22] - Matt Majic, Nov 24 2019
From Peter Bala, Feb 19 2020: (Start)
The following remarks assume an offset of 1 in the row and column indices of the triangle.
The sequence of row polynomials T(n,x), beginning with T(1,x) = x, T(2,x) = x + x^2, T(3,x) = x + 3*x^2 + x^3, ..., is a strong divisibility sequence of polynomials in the ring Z[x]; that is, for all positive integers n and m, poly_gcd(T(n,x), T(m,x)) = T(gcd(n, m), x) - apply Norfleet (2005), Theorem 3. Consequently, the sequence (T(n,x): n >= 1) is a divisibility sequence in the polynomial ring Z[x]; that is, if n divides m then T(n,x) divides T(m,x) in Z[x].
Let S(x) = 1 + 2*x + 6*x^2 + 22*x^3 + ... denote the o.g.f. for the large Schröder numbers A006318. The power series (x*S(x))^n, n = 2, 3, 4, ..., can be expressed as a linear combination with polynomial coefficients of S(x) and 1: (x*S(x))^n = T(n-1,-x) - T(n,-x)*S(x). The result can be extended to negative integer n if we define T(0,x) = 0 and T(-n,x) = (-1)^(n+1) * T(n,x)/x^n. Cf. A115139.
[In the previous two paragraphs, D(n,x) was replaced with T(n,x) because the contributor is referring to the rows of the triangle T(n,k), not the rows of the array D(n,k). - Petros Hadjicostas, Aug 05 2020] (End)
Named after the French amateur mathematician Henri-Auguste Delannoy (1833-1915). - Amiram Eldar, Apr 15 2021
D(i,j) = D(j,i). With this and Dmitry Zaitsev's Dec 10 2015 comment, D(i,j) can be considered the number of points at L1 distance <= i in Z^j or the number of points at L1 distance <= j in Z^i from any given point. The rows and columns of D(i,j) are the crystal ball sequences on cubic lattices. See the first example below. The n-th term in the k-th crystal ball sequence can be considered the number of points at distance <= n from any point in a k-dimensional cubic lattice, or the number of points at distance <= k from any point in an n-dimensional cubic lattice. - Shel Kaphan, Jan 01 2023 and Jan 07 2023
Dimensions of hom spaces Hom(R^{(i)}, R^{(j)}) in the Delannoy category attached to the oligomorphic group of order preserving self-bijections of the real line. - Noah Snyder, Mar 22 2023

			The square array D(i,j) (i >= 0, j >= 0) begins:
  1, 1,  1,   1,   1,   1,    1,    1,    1,    1, ... = A000012
  1, 3,  5,   7,   9,  11,   13,   15,   17,   19, ... = A005408
  1, 5, 13,  25,  41,  61,   85,  113,  145,  181, ... = A001844
  1, 7, 25,  63, 129, 231,  377,  575,  833, 1159, ... = A001845
  1, 9, 41, 129, 321, 681, 1289, 2241, 3649, 5641, ... = A001846
  ...
For D(2,5) = 61, which is seen above in the row labeled A001844, we calculate the sum (9 + 11 + 41) of the 3 nearest terms above and/or to the left. - _Peter Munn_, Jan 01 2023
D(2,5) = 61 can also be obtained from the row labeled A005408 using a recurrence mentioned in the formula section:  D(2,5) = D(1,5) + 2*Sum_{k=0..4} D(1,k), so D(2,5) = 11 + 2*(1+3+5+7+9) = 11 + 2*25. - _Shel Kaphan_, Jan 01 2023
As a triangular array (on its side) this begins:
   0,   0,   0,   0,   1,   0,  11,   0, ...
   0,   0,   0,   1,   0,   9,   0,  61, ...
   0,   0,   1,   0,   7,   0,  41,   0, ...
   0,   1,   0,   5,   0,  25,   0, 129, ...
   1,   0,   3,   0,  13,   0,  63,   0, ...
   0,   1,   0,   5,   0,  25,   0, 129, ...
   0,   0,   1,   0,   7,   0,  41,   0, ...
   0,   0,   0,   1,   0,   9,   0,  61, ...
   0,   0,   0,   0,   1,   0,  11,   0, ...
   [Edited by _Shel Kaphan_, Jan 01 2023]
From _Roger L. Bagula_, Dec 09 2008: (Start)
As a triangle T(n,k) (with rows n >= 0 and columns k = 0..n), this begins:
   1;
   1,  1;
   1,  3,   1;
   1,  5,   5,   1;
   1,  7,  13,   7,    1;
   1,  9,  25,  25,    9,    1;
   1, 11,  41,  63,   41,   11,    1;
   1, 13,  61, 129,  129,   61,   13,   1;
   1, 15,  85, 231,  321,  231,   85,  15,   1;
   1, 17, 113, 377,  681,  681,  377, 113,  17,  1;
   1, 19, 145, 575, 1289, 1683, 1289, 575, 145, 19, 1;
   ... (End)
Triangle T(n,k) recurrence: 63 = T(6,3) = 25 + 13 + 25 = T(5,2) + T(4,2) + T(5,3).
Triangle T(n,k) recurrence with A-sequence A112478: 63 = T(6,3) = 1*25 + 2*25 - 2*9 + 6*1 (T entries from row n = 5 only). [Here the formula T(n,k) = Sum_{j=0..n-k} A112478(j) * T(n-1, k-1+j) is used with n = 6 and k = 3; i.e., T(6,3) = Sum_{j=0..3} A111478(j) * T(5, 2+j). - _Petros Hadjicostas_, Aug 05 2020]
From _Philippe Deléham_, Mar 29 2012: (Start)
Subtriangle of the triangle given by (1, 0, 1, -1, 0, 0, 0, ...) DELTA (0, 1, 0, 0, 0, ...) where DELTA is the operator defined in A084938:
   1;
   1,  0;
   1,  1,  0;
   1,  3,  1,  0;
   1,  5,  5,  1,  0;
   1,  7, 13,  7,  1,  0;
   1,  9, 25, 25,  9,  1, 0;
   1, 11, 41, 63, 41, 11, 1, 0;
   ...
Subtriangle of the triangle given by (0, 1, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, ...) where DELTA is the operator defined in A084938:
   1;
   0, 1;
   0, 1,  1;
   0, 1,  3,  1;
   0, 1,  5,  5,  1;
   0, 1,  7, 13,  7,  1;
   0, 1,  9, 25, 25,  9,  1;
   0, 1, 11, 41, 63, 41, 11, 1;
   ... (End)

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 593.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
L. Moser and W. Zayachkowski, Lattice paths with diagonal steps, Scripta Mathematica, 26 (1963), 223-229.
G. Picou, Note #2235, L'Intermédiaire des Mathématiciens, 8 (1901), page 281. - N. J. A. Sloane, Mar 02 2022
D. B. West, Combinatorial Mathematics, Cambridge, 2021, p. 28.

T. D. Noe, Table of n, a(n) for n = 0..5150
K. Alladi and V. E. Hoggatt Jr., On tribonacci numbers and related functions, Fibonacci Quart. 15 (1977), 42-45.
Said Amrouche and Hacène Belbachir, Asymmetric extension of Pascal-Dellanoy triangles, arXiv:2001.11665 [math.CO], 2020.
J.-M. Autebert et al., H.-A. Delannoy et les oeuvres posthumes d'Édouard Lucas, Gazette des Mathématiciens - no 95, Jan 2003 (in French).
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], 2017. See Sect. 2.3.
C. Banderier and S. Schwer, Why Delannoy numbers?, arXiv:math/0411128 [math.CO], 2004.
Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, 9 (2006), #06.2.4.
Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, 15 (2012), #12.8.2.
Paul Barry, On the Inverses of a Family of Pascal-Like Matrices Defined by Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.5.6.
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
Paul Barry, The Central Coefficients of a Family of Pascal-like Triangles and Colored Lattice Paths, J. Int. Seq., 22 (2019), #19.1.3.
Paul Barry, Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles, J. Int. Seq., 22 (2019), #19.5.8.
Paul Barry, Notes on Riordan arrays and lattice paths, arXiv:2504.09719 [math.CO], 2025. See pp. 3, 29.
Frédéric Bihan, Francisco Santos, and Pierre-Jean Spaenlehauer, A Polyhedral Method for Sparse Systems with many Positive Solutions, arXiv:1804.05683 [math.CO], 2018.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids, English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 37.
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I, Math. Zeit. 233 (2000), 1-19.
C. Carré, N. Debroux, M. Deneufchatel, J.-P. Dubernard et al., Dirichlet convolution and enumeration of pyramid polycubes, hal-00905889, 2013.
C. Carre, N. Debroux, M. Deneufchatel, J.-Ph. Dubernard, C. Hillariet, J.-G. Luque, and O. Mallet, Enumeration of Polycubes and Dirichlet Convolutions, J. Int. Seq. 18 (2015), #15.11.4.
J. S. Caughman et al., A note on lattice chains and Delannoy numbers, Discrete Math., 308 (2008), 2623-2628.
Swee Hong Chan, Igor Pak, and Greta Panova, Log-concavity in planar random walks, arXiv:2106.10640 [math.CO], 2021.
H. Delannoy, Emploi de l'échiquier pour la résolution de certains problèmes de probabilités, Association Française pour l'Avancement des Sciences, 24th session, 1895, pp. 70-90 (see the table given on p. 76).
Jerry Ray Dias, Properties and relationships of conjugated polyenes having a reciprocal eigenvalue spectrum - dendralene and radialene hydrocarbons, Croatica Chem. Acta, 77 (2004), 325-330.
M. Dziemianczuk, Generalizing Delannoy numbers via counting weighted lattice paths, INTEGERS, 13 (2013), #A54.
James East and Nicholas Ham, Lattice paths and submonoids of Z^2, arXiv:1811.05735 [math.CO], 2018.
Steven Edwards and William Griffiths, Generalizations of Delannoy and cross polytope numbers, Fib. Q., Vol. 55, No. 4 (2017), pp. 356-366.
Steven Edwards and William Griffiths, On Generalized Delannoy Numbers, J. Int. Seq., 23 (2020), #20.3.6.
R. Feria-Puron, H. Perez-Roses, and J. Ryan, Searching for Large Circulant Graphs, arXiv:1503.07357 [math.CO], 2015.
R. Feria-Purón, J. Ryan, and H. Pérez-Rosés, Searching for Large Multi-Loop Networks, Electronic Notes in Discrete Mathematics, 46 (2014), 233-240.
Nate Harman, Andrew Snowden, and Noah Snyder, The Delannoy Category, arxiv:2211.15392 [math.RT], 2023.
Rebecca Hartman-Baker, The Diffusion Equation Method for Global Optimization and Its Application to Magnetotelluric Geoprospecting, University of Illinois, Urbana-Champaign, 2005.
G. Hetyei, Shifted Jacobi polynomials and Delannoy numbers, arXiv:0909.5512 [math.CO], 2009.
G. Hetyei, Links we almost missed between Delannoy numbers and Legendre polynomials.
V. E. Hoggatt, Jr., Letters to N. J. A. Sloane, 1974-1975.
Kentaro Ihara, Yayoi Nakamura, and Shuji Yamamoto, Interpolant of truncated multiple zeta functions, arXiv:2407.20509 [math.NT], 2024. See p. 21.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
Milan Janjić and B. Petkovic, A Counting Function, arXiv:1301.4550 [math.CO], 2013. - _N. J. A. Sloane_, Feb 13 2013
Milan Janjić and B. Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014), #14.3.5.
Svante Janson, Patterns in random permutations avoiding some sets of multiple patterns, arXiv:1804.06071 [math.PR], 2018.
Shel Kaphan, Tables of Sequences Related to Delannoy Numbers and Cubic Lattice Coordination Numbers
Shel Kaphan, Illustration of a recurrence relation on the Delannoy numbers and their connection with geometry.
G. Kreweras, Sur les hiérarchies de segments, Cahiers Bureau Universitaire Recherche Opérationnelle, # 20, Inst. Statistiques, Univ. Paris, 1973, pp. 4-10.
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
G. Kreweras, Aires des chemins surdiagonaux et application à un problème économique, Cahiers du Bureau universitaire de recherche opérationnelle Série Recherche 24 (1976), 1-8. [Annotated scanned copy]
Eon Lee, Andrés R. Vindas-Meléndez, and Zhi Wang, Generalized snake posets, order polytopes, and lattice-point enumeration, arXiv:2411.18695 [math.CO], 2024. See p. 15.
Yi-Lin Lee, Off-diagonally symmetric domino tilings of the Aztec diamond of odd order, arXiv:2404.09057 [math.CO], 2024. See p. 20.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 1.
M. LLadser, Uniform formulas for coefficients of meromorphic functions, arXiv:math/0604152 [math.CO], 2006.
E. Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 174.
Matt Majic, Relationships between spherical and bispherical harmonics, and an electrostatic T-matrix for dimers, preprint, 2019, DOI:10.13140/RG.2.2.21203.12320.
J. W. Meijer, Famous numbers on a chessboard, Acta Nova, 4(4) (2010), 589-598.
Mirka Miller, Hebert Perez-Roses, and Joe Ryan, The Maximum Degree-and-Diameter-Bounded Subgraph in the Mesh, arXiv:1203.4069 [math.CO], 2012.
Alejandro H. Morales, Igor Pak, and Greta Panova, Hook formulas for skew shapes IV. Increasing tableaux and factorial Grothendieck polynomials, arXiv:2108.10140 [math.CO], 2021.
Lili Mu and Sai-nan Zheng, On the Total Positivity of Delannoy-Like Triangles, Journal of Integer Sequences, 20 (2017), #17.1.6.
M. Norfleet, Characterization of second-order strong divisibility sequences of polynomials, The Fibonacci Quarterly, 43(2) (2005), 166-169.
Richard L. Ollerton and Anthony G. Shannon, Some properties of generalized Pascal squares and triangles, Fib. Q., 36 (1998), 98-109. See Table 9.
L. Pachter and B. Sturmfels, The mathematics of phylogenomics, arXiv:math/0409132 [math.ST], 2004-2005.
R. Pemantle and M. C. Wilson, Asymptotics of multivariate sequences, I: smooth points of the singular variety, arXiv:math/0003192 [math.CO], 2000.
J. L. Ramirez and V. F. Sirvent, Incomplete Tribonacci Numbers and Polynomials, Journal of Integer Sequences, 17 (2014), #14.4.2. See Table 1. - _N. J. A. Sloane_, Mar 23 2014
Marko Razpet, A self-similarity structure generated by king's walk, Algebraic and topological methods in graph theory (Lake Bled, 1999). Discrete Math. 244(1-3) (2002), 423--433. MR1844050 (2002k:05022).
Shiva Samieinia, Digital straight line segments and curves, Licentiate Thesis. Stockholm University, Department of Mathematics, Report 2007:6.
Shiva Samieinia, The number of continuous curves in digital geometry, Port. Math. 67(1) (2010), 75-89.
Seunghyun Seo, The Catalan Threshold Arrangement, Journal of Integer Sequences, 20 (2017), #17.1.1.
Yuriy Shablya, Combinatorial Generation Algorithms for Some Lattice Paths Using the Method Based on AND/OR Trees, Algorithms (2023) Vol. 16, No. 6, 266.
M. Shattuck, Combinatorial identities for incomplete tribonacci polynomials, arXiv:1406.2755 [math.CO], 2014.
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Yi Wang, Zheng Sai-Nan, and Chen Xi, Analytic aspects of Delannoy numbers, Discrete Mathematics 342(8) (2019), 2270-2277.
Eric Weisstein's World of Mathematics, Delannoy Number.
Meng-Han Wu, Henryk A. Witek, Rafał Podeszwa, Clar Covers and Zhang-Zhang Polynomials of Zigzag and Armchair Carbon Nanotubes, MATCH Commun. Math. Comput. Chem. (2025) Vol. 93, 415-462. See p. 445.
Dmitry Zaitsev, k-neighborhood for Cellular Automata, arXiv:1605.08870 [cs.DM], 2016.
Liang Zhao and Fengyao Yan, Note on Total Positivity for a Class of Recursive Matrices, Journal of Integer Sequences, 19 (2016), #16.6.5.

Sums of antidiagonals: A000129 (Pell numbers).
Main diagonal: A001850 (central Delannoy numbers), which has further information and references.
A002002, A026002, and A190666 are +-k-diagonals for k=1, 2, 3 resp. - Shel Kaphan, Jan 01 2023
Rows 0..10: A000012, A005408, A001844, A001845, A001846, A001847, A001848, A001849, A008417, A008419, A008421.
See also A027618.
Cf. A059446.
Has same main diagonal as A064861. Different from A100936.
Cf. A101164, A101167, A128966, A131887, A131935.
Cf. A035607, A108625, A142979, A142992, A143007.
Read mod small primes: A211312, A211313, A211314, A211315.
Triangle sums (see the comments): A000129 (Row1); A056594 (Row2); A000073 (Kn11 & Kn21); A089068 (Kn12 & Kn22); A180668 (Kn13 & Kn23); A180669 (Kn14 & Kn24); A180670 (Kn15 & Kn25); A099463 (Kn3 & Kn4); A116404 (Fi1 & Fi2); A006498 (Ca1 & Ca2); A006498(3*n) (Ca3 & Ca4); A079972 (Gi1 & Gi2); A079972(4*n) (Gi3 & Gi4); A079973(3*n) (Ze1 & Ze2); A079973(2*n) (Ze3 & Ze4).
Cf. A102413, A128966. (D(n,1)) = A005843. Cf. A115139.

a008288 n k = a008288_tabl !! n !! k
a008288_row n = a008288_tabl !! n
a008288_tabl = map fst $ iterate
    (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                       zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [1, 1])
-- Reinhard Zumkeller, Jul 21 2013

A008288 := proc(n, k) option remember; if k = 0 then 1 elif n=k then 1 else procname(n-1, k-1) + procname(n-2, k-1) + procname(n-1, k) end if; end proc: seq(seq(A008288(n,k),k=0..n), n=0..10); # triangular indices n and k
P[0]:=1; P[1]:=x+1; for n from 2 to 12 do P[n]:=expand((x+1)*P[n-1]+x*P[n-2]); lprint(P[n]); lprint(seriestolist(series(P[n],x,200))); end do:

(* Next, A008288 jointly generated with A035607 *)
u[1, x_] := 1; v[1, x_] := 1; z = 16;
u[n_, x_] := x*u[n - 1, x] + v[n - 1, x];
v[n_, x_] := 2 x*u[n - 1, x] + v[n - 1, x];
Table[Expand[u[n, x]], {n, 1, z/2}]
Table[Expand[v[n, x]], {n, 1, z/2}]
cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
TableForm[cu]
Flatten[%]    (* A008288 *)
Table[Expand[v[n, x]], {n, 1, z}]
cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
TableForm[cv]
Flatten[%]    (* A035607 *)
(* Clark Kimberling, Mar 09 2012 *)
d[n_, k_] := Binomial[n+k, k]*Hypergeometric2F1[-k, -n, -n-k, -1]; A008288 = Flatten[Table[d[n-k, k], {n, 0, 12}, {k, 0, n}]] (* Jean-François Alcover, Apr 05 2012, after 3rd formula *)

from functools import cache
@cache
def delannoy_row(n: int) -> list[int]:
    if n == 0: return [1]
    if n == 1: return [1, 1]
    rov = delannoy_row(n - 2)
    row = delannoy_row(n - 1) + [1]
    for k in range(n - 1, 0, -1):
        row[k] += row[k - 1] + rov[k - 1]
    return row
for n in range(10): print(delannoy_row(n))  # Peter Luschny, Jul 30 2023

for k in range(8):  # seen as an array, read row by row
    a = lambda n: hypergeometric([-n, -k], [1], 2)
    print([simplify(a(n)) for n in range(11)]) # Peter Luschny, Nov 19 2014

D(n, 0) = 1 = D(0, n) for n >= 0; D(n, k) = D(n, k-1) + D(n-1, k-1) + D(n-1, k).
Bivariate o.g.f.: Sum_{n >= 0, k >= 0} D(n, k)*x^n*y^k = 1/(1 - x - y - x*y).
D(n, k) = Sum_{d = 0..min(n,k)} binomial(k, d)*binomial(n+k-d, k) = Sum_{d=0..min(n,k)} 2^d*binomial(n, d)*binomial(k, d). [Edited by Petros Hadjicostas, Aug 05 2020]
Seen as a triangle read by rows: T(n, 0) = T(n, n) = 1 for n >= 0 and T(n, k) = T(n-1, k-1) + T(n-2, k-1) + T(n-1, k), 0 < k < n and n > 1. - Reinhard Zumkeller, Dec 03 2004
Read as a number triangle, this is the Riordan array (1/(1-x), x(1+x)/(1-x)) with T(n, k) = Sum_{j=0..n-k} C(n-k, j) * C(k, j) * 2^j. - Paul Barry, Jul 18 2005
T(n,k) = Sum_{j=0..n-k} C(k,j)*C(n-j,k). - Paul Barry, May 21 2006
Let y^k(n) be the number of Khalimsky-continuous functions f from [0,n-1] to Z such that f(0) = 0 and f(n-1) = k. Then y^k(n) = D(i,j) for i = (1/2)*(n-1-k) and j = (1/2)*(n-1+k) where n-1+k belongs to 2Z. - Shiva Samieinia (shiva(AT)math.su.se), Oct 08 2007
Recurrence for triangle from A-sequence (see the Wolfdieter Lang comment above): T(n,k) = Sum_{j=0..n-k} A112478(j) * T(n-1, k-1+j), n >= 1, k >= 1. [For k > n, the sum is empty, in which case T(n,k) = 0.]
From Peter Bala, Jul 17 2008: (Start)
The n-th row of the square array is the crystal ball sequence for the product lattice A_1 x ... x A_1 (n copies). A035607 is the table of the associated coordination sequences for these lattices.
The polynomial p_n(x) := Sum {k = 0..n} 2^k * C(n,k) * C(x,k) = Sum_{k = 0..n} C(n,k) * C(x+k,n), whose values [p_n(0), p_n(1), p_n(2), ... ] give the n-th row of the square array, is the Ehrhart polynomial of the n-dimensional cross polytope (the hyperoctahedron) [Bump et al. (2000), Theorem 6].
The first few values are p_0(x) = 1, p_1(x) = 2*x + 1, p_2(x) = 2*x^2 + 2*x + 1 and p_3(x) = (4*x^3 + 6*x^2 + 8*x + 3)/3.
The reciprocity law p_n(m) = p_m(n) reflects the symmetry of the table.
The polynomial p_n(x) is the unique polynomial solution of the difference equation (x+1)*f(x+1) - x*f(x-1) = (2*n+1)*f(x), normalized so that f(0) = 1.
These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane; that is, the polynomials p_n(x-1), n = 1,2,3,..., satisfy a Riemann hypothesis [Bump et al. (2000), Theorem 4]. The o.g.f. for the p_n(x) is (1 + t)^x/(1 - t)^(x + 1) = 1 + (2*x + 1)*t + (2*x^2 + 2*x + 1)*t^2 + ... .
The square array of Delannoy numbers has a close connection with the constant log(2). The entries in the n-th row of the array occur in the series acceleration formula log(2) = (1 - 1/2 + 1/3 - ... + (-1)^(n+1)/n) + (-1)^n * Sum_{k>=1} (-1)^(k+1)/(k*D(n,k-1)*D(n,k)). [T(n,k) was replaced with D(n,k) in the formula to agree with the beginning of the paragraph. - Petros Hadjicostas, Aug 05 2020]
For example, the fourth row of the table (n = 3) gives the series log(2) = 1 - 1/2 + 1/3 - 1/(1*1*7) + 1/(2*7*25) - 1/(3*25*63) + 1/(4*63*129) - ... . See A142979 for further details.
Also the main diagonal entries (the central Delannoy numbers) give the series acceleration formula Sum_{n>=1} 1/(n*D(n-1,n-1)*D(n,n)) = (1/2)*log(2), a result due to Burnside. [T(n,n) was replaced here with D(n,n) to agree with the previous paragraphs. - Petros Hadjicostas, Aug 05 2020]
Similar relations hold between log(2) and the crystal ball sequences of the C_n lattices A142992. For corresponding results for the constants zeta(2) and zeta(3), involving the crystal ball sequences for root lattices of type A_n and A_n x A_n, see A108625 and A143007 respectively. (End)
From Peter Bala, Oct 28 2008: (Start)
Hilbert transform of Pascal's triangle A007318 (see A145905 for the definition of this term).
D(n+a,n) = P_n(a,0;3) for all integer a such that a >= -n, where P_n(a,0;x) is the Jacobi polynomial with parameters (a,0) [Hetyei]. The related formula A(n,k) = P_k(0,n-k;3) defines the table of asymmetric Delannoy numbers, essentially A049600. (End)
Seen as a triangle read by rows: T(n, k) = Hyper2F1([k-n, -k], [1], 2). - Peter Luschny, Aug 02 2014, Oct 13 2024.
From Peter Bala, Jun 25 2015: (Start)
O.g.f. for triangle T(n,k): A(z,t) = 1/(1 - (1 + t)*z - t*z^2) = 1 + (1 + t)*z + (1 + 3*t + t^2)*z^2 + (1 + 5*t + 5*t^2 + t^3)*z^3 + ....
1 + z*d/dz(A(z,t))/A(z,t) is the o.g.f. for A102413. (End)
E.g.f. for the n-th subdiagonal of T(n,k), n >= 0, equals exp(x)*P(n,x), where P(n,x) is the polynomial Sum_{k = 0..n} binomial(n,k)*(2*x)^k/k!. For example, the e.g.f. for the second subdiagonal is exp(x)*(1 + 4*x + 4*x^2/2) = 1 + 5*x + 13*x^2/2! + 25*x^3/3! + 41*x^4/4! + 61*x^5/5! + .... - Peter Bala, Mar 05 2017 [The n-th subdiagonal of triangle T(n,k) is the n-th row of array D(n,k).]
Let a_i(n) be multiplicative with a_i(p^e) = D(i, e), p prime and e >= 0, then Sum_{n > 0} a_i(n)/n^s = (zeta(s))^(2*i+1)/(zeta(2*s))^i for i >= 0. - Werner Schulte, Feb 14 2018
Seen as a triangle read by rows: T(n,k) = Sum_{i=0..k} binomial(n-i, i) * binomial(n-2*i, k-i) for 0 <= k <= n. - Werner Schulte, Jan 09 2019
Univariate generating function: Sum_{k >= 0} D(n,k)*z^k = (1 + z)^n/(1 - z)^(n+1). [Dziemianczuk (2013), Eq. 5.3] - Matt Majic, Nov 24 2019
(n+1)*D(n+1,k) = (2*k+1)*D(n,k) + n*D(n-1,k). [Majic (2019), Eq. 22] - Matt Majic, Nov 24 2019
For i, j >= 1, D(i,j) = D(i,j-1) + 2*Sum_{k=0..i-1} D(k,j-1), or, because D(i,j) = D(j,i), D(i,j) = D(i-1,j) + 2*Sum_{k=0..j-1} D(i-1,k). - Shel Kaphan, Jan 01 2023
Sum_{k=0..n} T(n,k)^2 = A026933(n). - R. J. Mathar, Nov 07 2023
Let S(x) = (1 - x - (1 - 6*x + x^2)^(1/2))/(2*x) denote the g.f. of the sequence of large Schröder numbers A006318. Read as a lower triangular array, the signed n-th row polynomial R(n, -x) = 1/sqrt(1 - 6*x + x^2) *( 1/S(x)^(n+1) + (x*S(x))^(n+1) ). For example, R(4, -x) = 1 - 7*x + 13*x^2 - 7*x^3 + x^4 = 1/sqrt(1 - 6*x + x^2) * ( 1/S(x)^5 + (x*S(x))^5 ). Cf. A102413. - Peter Bala, Aug 01 2024

Expanded description from Clark Kimberling, Jun 15 1997
Additional references from Sylviane R. Schwer (schwer(AT)lipn.univ-paris13.fr), Nov 28 2001
Changed the notation to make the formulas more precise. - N. J. A. Sloane, Jul 01 2002

A143003 a(0) = 0, a(1) = 1, a(n+1) = (2n+1)(n^2+n+5)a(n) - n^6a(n-1).

Original entry on oeis.org

0, 1, 21, 1091, 114520, 21298264, 6410456640, 2923097201856, 1920450126458880, 1747596822651334656, 2133806329230225408000, 3405545462439659704320000, 6950705677729940374290432000, 17807686090745585163974737920000

Offset: 0

Peter Bala, Jul 19 2008

This is the case m = 1 of the general recurrence a(0) = 0, a(1) = 1, a(n+1) = (2*n+1)*(n^2+n+2*m^2+2*m+1)*a(n) - n^6*a(n-1) (we suppress the dependence of a(n) on m), which arises when accelerating the convergence of the series Sum_{k>=1} 1/k^3 for Apery's constant zeta(3). For other cases see A066989 (m=0), A143004 (m=2), A143005 (m=3) and A143006 (m=4).
The solution to the general recurrence may be expressed as a sum: a(n) = n!^3*p_m(n)*Sum_{k = 1..n} 1/(k^3*p_m(k-1)*p_m(k)), where p_m(x) = Sum_{k = 0..n} C(2*k,k)^2*C(n+k,2*k)*C(x+k,2*k) is a polynomial in x of degree 2*m.
The first few are p_0(x) = 1, p_1(x) = 2*x^2 + 2*x + 1, p_2(x) = (3*x^4 + 6*x^3 + 9*x^2 + 6*x + 2)/2 and p_3(x) = (10*x^6 + 30*x^5 + 85*x^4 + 120*x^3 + 121*x^2 + 66*x + 18)/18. For fixed n, the sequence [p_n(k)]k>=0 is the crystal ball sequence for the product lattice A_n x A_n. See A143007 for the table of values [p_n(k)] n,k >= 0. Observe that [p_n(n)] n >= 0 is the sequence of Apery numbers A005259.
The reciprocity law p_m(n) = p_n(m) holds for nonnegative integers m and n. In particular we have p_m(1) = 2*m^2 + 2*m + 1 and p_m(2) = (3*m^4 + 6*m^3 + 9*m^2 + 6*m + 2)/2.
The polynomial p_m(x) is the unique polynomial solution of the difference equation (x+1)^3*f(x+1) + x^3*f(x-1) = (2*x+1)*(x^2+x+2*m^2+2*m+1)*f(x), normalized so that f(0) = 1. The reciprocity law now yields the Apery-like recursion m^3*p_m(x) + (m-1)^3*p_(m-2)(x) = (2*m-1)*(m^2-m+1+2*x^2+2*x)*p_(m-1)(x).
The polynomial functions p_m(x) have their zeros on the vertical line Re x = -1/2 in the complex plane; that is, the polynomials p_m(x-1), m = 1,2,3,..., satisfy a Riemann hypothesis (adapt the proof of the lemma on p. 4 of [BUMP et al.]).
The general recurrence in the first paragraph above has a second solution b(n) = n!^3*p_m(n) with initial conditions b(0) = 1, b(1) = 2*m^2+2*m+1. Hence the behavior of a(n) for large n is given by lim_{n -> infinity} a(n)/b(n) = Sum_{k>=1} 1/(k^3*p_m(k-1)*p_m(k)) = 1/((2*m^2+2*m+1) - 1^6/(3*(2*m^2+2*m+3) - 2^6/(5*(2*m^2+2*m+7) - 3^6/(7*(2*m^2+2*m+13) - ...)))) = Sum_{k>=1} 1/(m+k)^3. The final equality follows from a result of Ramanujan; see [Berndt, Chapter 12, Entry 32(iii)].
For the corresponding results for the constant zeta(2) see A142995. For corresponding results for the constant log(2) see A142979 and A142992.

Bruce C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag.

Seiichi Manyama, Table of n, a(n) for n = 0..181
D. Bump, K. Choi, P. Kurlberg and J. Vaaler, A local Riemann hypothesis, I, Math. Zeit. 233, (2000), 1-19.

Cf. A066989, A008288, A108625, A142979, A142992, A142995, A143004, A143005, A143006, A143007.

The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

p := n -> 2*n^2+2*n+1: a := n -> n!^3*p(n)*sum (1/(k^3*p(k-1)*p(k)), k = 1..n): seq(a(n), n = 0..14)

RecurrenceTable[{a[0]==0,a[1]==1,a[n+1]==(2n+1)(n^2+n+5)a[n]- n^6 a[n-1]}, a[n],{n,15}] (* Harvey P. Dale, Jun 20 2011 *)

a(n) = n!^3*p(n)*Sum_{k = 1..n} 1/(k^3*p(k-1)*p(k)), where p(n) = 2*n^2 + 2*n + 1 = A001844(n).
Recurrence: a(0) = 0, a(1) = 1, a(n+1) = (2*n+1)*(n^2+n+5)*a(n) - n^6*a(n-1).
The sequence b(n):= n!^3*p(n) satisfies the same recurrence with the initial conditions b(0) = 1, b(1) = 5. Hence we obtain the finite continued fraction expansion a(n)/b(n) = 1/(5 - 1^6/(21 - 2^6/(55 - 3^6/(119 - ... - (n-1)^6/((2*n-1)*(n^2-n+5)))))), for n >= 2. The behavior of a(n) for large n is given by lim_{n -> infinity} a(n)/b(n) = Sum_{k>=1} 1/(k^3*(4*k^4 + 1)) = 1/(5 - 1^6/(21 - 2^6/(55 - 3^6/(119 - ... - n^6/((2*n+1)*(n^2+n+5) - ...))))) = zeta(3) - 1, where the final equality follows from a result of Ramanujan; see [Berndt, Chapter 12, Entry 32(iii) at x = 1].

A024167 a(n) = n!*(1 - 1/2 + 1/3 - ... + c/n), where c = (-1)^(n+1).

Original entry on oeis.org

1, 1, 5, 14, 94, 444, 3828, 25584, 270576, 2342880, 29400480, 312888960, 4546558080, 57424792320, 948550176000, 13869128448000, 256697973504000, 4264876094976000, 87435019510272000, 1627055289796608000, 36601063093905408000, 754132445894209536000

Offset: 1

Clark Kimberling

Stirling transform of (-1)^n*a(n-1) = [0, 1, -1, 5, -14, 94, ...] is A000629(n-2) = [0, 1, 2, 6, 26, ...]. - Michael Somos, Mar 04 2004
Stirling transform of a(n) = [1, 1, 5, 14, 94, ...] is A052882(n) = [1, 2, 9, 52, 375, ...]. - Michael Somos, Mar 04 2004
a(n) is the number of n-permutations that have a cycle with length greater than n/2. - Geoffrey Critzer, May 28 2009
From Jens Voß, May 07 2010: (Start)
a(4n) is divisible by 6*n + 1 for all n >= 1; the quotient of a(4*n) and 6*n+1 is A177188(n).
a(4*n+3) is divisible by 6*n + 5 for all n >= 0; the quotient of a(4*n+3) and 6*n + 5 is A177174(n). (End)

			G.f. = x + x^2 + 5*x^3 + 14*x^4 + 94*x^5 + 444*x^6 + 3828*x^7 + 25584*x^8 + ...

Seiichi Manyama, Table of n, a(n) for n = 1..449

Cf. A024168, A000254, A054651, A094587, A142979, A177174, A177188.

a := n -> n!*(log(2) - (-1)^n*LerchPhi(-1, 1, n+1));
seq(simplify(a(n)), n=1..20); # Peter Luschny, Dec 27 2018

f[k_] := k (-1)^(k + 1)
t[n_] := Table[f[k], {k, 1, n}]
a[n_] := SymmetricPolynomial[n - 1, t[n]]
Table[a[n], {n, 1, 18}]    (* A024167 signed *)
(* Clark Kimberling, Dec 30 2011 *)
a[ n_] := If[ n < 0, 0, n! Sum[ -(-1)^k / k, {k, n}]]; (* Michael Somos, Nov 28 2013 *)
a[ n_] := If[ n < 0, 0, n! (PolyGamma[n + 1] - PolyGamma[(n + Mod[n, 2, 1]) / 2])]; (* Michael Somos, Nov 28 2013 *)
a[ n_] := If[ n < 1, 0, (-1)^Quotient[n, 2] SymmetricPolynomial[ n - 1, Table[ -(-1)^k k, {k, n}]]]; (* Michael Somos, Nov 28 2013 *)

{a(n) = if( n<0, 0, n! * polcoeff( log(1 + x + x * O(x^n)) / (1 - x), n))}; /* Michael Somos, Mar 02 2004 */

x='x+O('x^33); Vec(serlaplace(log(1+x)/(1-x))) \\ Joerg Arndt, Dec 27 2018

def A():
    a, b, n = 1, 1, 2
    yield(a)
    while True:
        yield(a)
        b, a = a, a + b * n * n
        n += 1
a = A(); print([next(a) for  in range(20)]) # _Peter Luschny, May 19 2020

E.g.f.: log(1 + x)/(1 - x). - Vladeta Jovovic, Aug 25 2002
a(n) = a(n-1) + a(n-2) * (n-1)^2, n > 1. - Michael Somos, Oct 29 2002
b(n) = n! satisfies the above recurrence with b(1) = 1, b(2) = 2. This gives the finite continued fraction expansion a(n)/n! = 1/(1 + 1^2/(1 + 2^2/(1 + 3^2/(1 + ... + (n-1)^2/1)))). Cf. A142979. - Peter Bala, Jul 17 2008
a(n) = A081358(n) - A092691(n). - Gary Detlefs, Jul 09 2010
E.g.f.: (x/(x-1))/G(0) where G(k) = -1 + (x-1)*k + x*(k+1)^2/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Aug 18 2012
a(n) ~ log(2)*n!. - Daniel Suteu, Dec 03 2016
a(n) = (1/2)*n!*((-1)^n*(digamma((n+1)/2) - digamma((n+2)/2)) + log(4)). - Daniel Suteu, Dec 03 2016
a(n) = n!*(log(2) - (-1)^n*LerchPhi(-1, 1, n+1)). - Peter Luschny, Dec 27 2018
a(n) = A054651(n,n-1). - Pontus von Brömssen, Oct 25 2020
a(n) = Sum_{k=0..n} (-1)^k*k!*A094587(n, k+1). - Mélika Tebni, Jun 20 2022
a(n) = n * a(n-1) - (-1)^n * (n-1)! for n > 1. - Werner Schulte, Oct 20 2024

More terms from Benoit Cloitre, Jan 27 2002

a(21)-a(22) from Pontus von Brömssen, Oct 25 2020

A142992 Square array, read by ascending antidiagonals, of the crystal ball sequences for the root lattices of type C_n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 9, 5, 1, 1, 19, 25, 7, 1, 1, 33, 85, 49, 9, 1, 1, 51, 225, 231, 81, 11, 1, 1, 73, 501, 833, 489, 121, 13, 1, 1, 99, 985, 2471, 2241, 891, 169, 15, 1, 1, 129, 1765, 6321, 8361, 4961, 1469, 225, 17, 1

Offset: 0

Peter Bala, Jul 18 2008

The lattice C_n consists of all integer lattice points v = (x_1,...,x_n) in Z^n such that the sum x_1 + ... + x_n is even. Let ||v|| = 1/2 * Sum_{i = 1..n} |x_i|; this defines a norm on C_n. The k-th term of the crystal ball sequence of C_n gives the number of lattice points v in C_n with ||v|| <= k [Bacher et al.]. The case n = 2 is illustrated in the Example section below.
This array has a remarkable relationship with the constant log(2). The row, column and (conjecturally) the diagonal entries of the array occur in series acceleration formulas for log(2) (see the Formula section below for some examples).
See A103884 for the table of coordination sequences of the C_n lattices. For the crystal ball sequences for the A_n and D_n lattices see A108625 and A108553 respectively. For the crystal ball sequences for the product lattices A_1 x ... x A_1(n copies) and A_n x A_n see A008288 and A143007 respectively.

			The square array begins
n\k|0...1....2.....3.....4......5
=================================
.0.|1...1....1.....1.....1......1
.1.|1...3....5.....7.....9.....11
.2.|1...9...25....49....81....121 A016754
.3.|1..19...85...231...489....891 A063496
.4.|1..33..225...833..2241...4961 A142993
.5.|1..51..501..2471..8361..22363 A142994
...
Triangular array begins
n\k|0...1...2...3...4...5
=========================
.0.|1
.1.|1...1
.2.|1...3...1
.3.|1...9...5...1
.4.|1..19..25...7...1
.5.|1..33..85..49...9...1
Case n = 2: The C_2 lattice consists of all integer lattice points v = (x,y) in Z x Z such that x + y is even, equipped with the taxicab type norm ||v|| = 1/2 * (|x| + |y|). There are 8 lattice points (marked with a 1 on the figure below) satisfying ||v|| = 1 and 16 lattice points (marked with a 2 on the figure) satisfying ||v|| = 2. Hence the crystal ball sequence for the C_2 lattice (row 2 of the table) begins 1, 1+8 = 9, 1+8+16 = 25, ... .
. . . . . . . . . . .
. . . . . 2 . . . . .
. . . . 2 . 2 . . . .
. . . 2 . 1 . 2 . . .
. . 2 . 1 . 1 . 2 . .
. 2 . 1 . 0 . 1 . 2 .
. . 2 . 1 . 1 . 2 . .
. . . 2 . 1 . 2 . . .
. . . . 2 . 2 . . . .
. . . . . 2 . . . . .
. . . . . . . . . . .

Seiichi Manyama, Rows n = 0..139, flattened
R. Bacher, P. de la Harpe and B. Venkov, Séries de croissance et séries d'Ehrhart associées aux réseaux de racines, C. R. Acad. Sci. Paris, 325 (Series 1) (1997), 1137-1142.
D. Bump, K. Choi, P. Kurlberg and J. Vaaler, A local Riemann hypothesis, I, Math. Zeit. 233, (2000), 1-19.
J. H. Conway and N. J. A. Sloane, Low dimensional lattices VII Coordination sequences, Proc. R. Soc. Lond., Ser. A, 453 (1997), 2369-2389.

Row 2-10 give A016754, A063496, A142993, A142994, A305549, A305721, A305722, A305723, A305724.

Cf. A008288, A103884, A108625, A142979, A143007.

with combinat: T := (n,k) -> add(binomial(2n,2i)*binomial(k+i,n),i = 0..n): for n from 0 to 9 do seq(T(n,k), k = 0..9) end do;

t[n_, k_] := Sum[ Binomial[2*n, 2*i]*Binomial[k+i, n], {i, 0, n}]; Table[t[n-k, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 06 2013 *)

T(n,k) = Sum_{i = 0..n} C(2*n,2*i)*C(k+i,n).
O.g.f. for row n: 1/(1-x)^(n+1) * Sum_{k = 0..n} C(2*n,2*k)*x^k = 1/(1-x) * T(n,(1+x)/(1-x)), where T(n,x) denotes the Chebyshev polynomial of the first kind.
O.g.f. for the array: 1/(1-x) * {(1-t) - x*(1+t)}/{(1-t)^2 - x*(1+t)^2} = (1+x+x^2+x^3+...) + (1+3*x+5*x^2+7*x^3+...)*t + (1+9*x+25*x^2+49*x^3+...)*t^2 + ... .
Row n of the array has the form [p_n(0),p_n(1),p_n(2),...], where the polynomial function p_n(x) = Sum_{k = 0..n} C(2*n,2*k)*C(x+k,n). The first few are p_0(x) = 1, p_1(x) = 2*x+1, p_2(x) = (2*x+1)^2, p_3(x) = (2*x+1)*(8*x^2+8*x+3)/3 and p_4(x) = (2*x+1)^2*(4*x^2+4*x+3)/3.
Alternative expressions for p_n(x) include p_n(x) = Sum_{k = 0..n} 2^(2*k)*n/(n+k)*C(n+k,2*k)*C(x,k) and p_n(x) = Sum_{k = 1..n} 2^(k-1)*C(n-1,k-1)*C(2*x+1,k).
The polynomials p_n(x) satisfy the 3-term recurrence relation n*p_n(x) = 2*(2*x+1)*p_(n-1)(x)+(n-2)*p(n-2)(x) for n >= 2; their generating function is 1/2*((1+t)/(1-t))^(2*x+1) = 1/2 + (2*x+1)*t + (2*x+1)^2*t^2 + (2*x+1)*(8*x^2+8*x+3)/3*t^3 + ... . Thus p_n(x) is, apart from a constant factor, the Meixner polynomial of the first kind M_n(2*x+1;b,c) at b = 0, c = -1. Compare with A142979.
The polynomial p_n(x) is the unique polynomial solution to the difference equation (2*x+1)*{f(x+1/2) - f(x-1/2)} = 2*n*f(x), normalized so that f(0) = 1. The function p_n(x) is also the unique polynomial solution to the difference equation (2*x+1)*{(x+1)*f(x+1) + x*f(x-1)} = ((2*x+1)^2 + 2*n^2)*f(x), normalized so that f(0) = 1.
The zeros of p_n(x) lie on the vertical line Re x = -1/2 in the complex plane, that is, the polynomials p_n(x-1), n = 1,2,3,..., satisfy a Riemann hypothesis (adapt the proof of the lemma on p.4 of [BUMP et al.]).
For n > 0, the entries in row n of the array occur in series acceleration formulas for log(2): 2*log(2) = 1 + (1/2 - 1/6 +...+(-1)^n/(n*(n-1))) + (-1)^(n+1)*Sum_{k >= 1} 1/(k*T(n,k-1)*T(n,k)). For example, the fourth row of the table (n = 3) gives 2*log(2) = 4/3 + 1/(1*1*19) + 1/(2*19*85) + 1/(3*85*231) + ... .
The corresponding result for column k is 2*log(2) = 1 + (1/(1*3) + 1/(2*3*5) +...+ 1/(k*(2*k-1)*(2k+1)) + (2*k+1)*Sum_{n >= 1} (-1)^(n+1)/(n*(n+1)*T(n,k)* T(n+1,k)).
For example, the third column of the table (k = 2) gives 2*log(2) = 41/30 + 5*(1/(1*2*5*25) - 1/(2*3*25*85) + 1/(3*4*85*225) - ... ).
For the main diagonal calculation suggests the result: 2*log(2) = 4/3 + Sum_{n >= 1} (-1)^(n+1)*(5*n+3)/(n*(n+1)*T(n,n)*T(n+1,n+1)).
Similar series acceleration formulas for log(2) come from the row, column and diagonal entries of the square array of Delannoy numbers, A008288 (which may viewed as the array of crystal ball sequences for the product lattices A_1 x...x A_1). For corresponding results for the constants zeta(2) and zeta(3) see A108625 and A143007 respectively.

A142983 a(1) = 1, a(2) = 2, a(n+2) = 2a(n+1) + (n + 1)(n + 2)*a(n).

Original entry on oeis.org

1, 2, 10, 44, 288, 1896, 15888, 137952, 1419840, 15255360, 186693120, 2387093760, 33898314240, 502247692800, 8123141376000, 136785729024000, 2483065912320000, 46822564905984000, 942853671825408000, 19678282007924736000, 435355106182520832000

Offset: 1

Peter Bala, Jul 17 2008

This is the case m = 1 of the general recurrence a(1) = 1, a(2) = 2*m, a(n+2) = 2*m*a(n+1) + (n + 1)*(n + 2)*a(n) (we suppress the dependence of a(n) on m), which arises when accelerating the convergence of the series 1/2 + 1/2*Sum_{k > 1} (-1)^(k+1)/(k*(k + 1)) = log(2). For other cases see A142984 (m = 2), A142985 (m = 3), A142986 (m = 4) and A142987 (m = 5).
The solution to the general recurrence may be expressed as a sum: a(n) = n!*p_m(n+1)*Sum_{k = 1..n} (-1)^(k+1)/(p_m(k)*p_m(k+1)), where p_m(x) = Sum_{k = 1..m} 2^(k-1)*C(m-1,k-1)*C(x,k) is the polynomial that gives the regular polytope numbers for the m-dimensional cross polytope as defined by [Kim] (see A142978). The first few values are p_1(x) = x, p_2(x) = x^2, p_3(x) = (2*x^3 + x)/3 and p_4(x) = (x^4 + 2*x^2)/3.
The polynomial p_m(x) is the unique polynomial solution of the difference equation x*(f(x+1) - f(x-1)) = 2*m*f(x), normalized so that f(1) = 1.
The o.g.f. for the p_m(x) is 1/2*((1 + t)/(1 - t))^x = 1/2 + x*t + x^2*t^2 + (2*x^3 + x)/3*t^3 + .... Thus p_m(x) is, apart from a constant factor, the Meixner polynomial of the first kind M_m(x;b,c) at b = 0, c = -1, also known as a Mittag-Leffler polynomial.
The general recurrence in the first paragraph above has a second solution b(n) = n!*p_m(n+1) with b(1) = 2*m, b(2) = m^2 + 2. Hence the behavior of a(n) for large n is given by Limit_{n-> oo} a(n)/b(n) = Sum_{k >= 1} (-1)^(k+1)/(p_m(k)*p_m(k+1)) = 1/((2*m) + 1*2/((2*m) + 2*3/((2*m) + 3*4/((2*m) + ... + n*(n + 1)/((2*m) + ...))))) = 1 + (-1)^(m+1) * (2*m)*(log(2) - (1 - 1/2 + 1/3 - ... + (-1)^(m+1)/m)), where the final equality follows by a result of Ramanujan (see [Berndt, Chapter 12, Entry 32(i)]).
See A142979, A142988 and A142992 for similar results. For corresponding results for Napier's constant e, the constant zeta(2) and Apery's constant zeta(3) refer to A000522, A142995 and A143003, respectively.

Bruce C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75.
Eric Weisstein's World of Mathematics, Meixner polynomial of the first kind.
Eric Weisstein's World of Mathematics, Mittag-Leffler polynomial.

Cf. A000522, A142979, A142984, A142985, A142986, A142987, A142988, A142992.

Cf. A002378.

a142983 n = a142983_list !! (n-1)
a142983_list = 1 : 2 : zipWith (+)
                       (map (* 2) $ tail a142983_list)
                       (zipWith (*) (drop 2 a002378_list) a142983_list)
-- Reinhard Zumkeller, Jul 17 2015

a := n -> (n+1)!*sum ((-1)^(k+1)/(k*(k+1)), k = 1..n): seq(a(n), n = 1..20);

Rest[CoefficientList[Series[(-x+2*Log[x+1])/(x-1)^2,{x,0,20}],x]*Range[0,20]!] (* Vaclav Kotesovec, Oct 21 2012 *)

a(n) = n!*p(n+1)*Sum_{k = 1..n} (-1)^(k+1)/(p(k)*p(k+1)), where p(n) = n.
Recurrence: a(1) = 1, a(2) = 2, a(n+2) = 2*a(n+1) + (n + 1)*(n + 2)*a(n).
The sequence b(n) := n!*p(n+1) satisfies the same recurrence with b(1) = 2 and b(2) = 6.
Hence we obtain the finite continued fraction expansion a(n)/b(n) = 1/(2 + 1*2/(2 + 2*3/(2 + 3*4/(2 + ... + (n - 1)*n/2)))), for n >= 2.
The behavior of a(n) for large n is given by Limit_{n -> oo} a(n)/b(n) = 1/(2 + 1*2/(2 + 2*3/(2 + 3*4/(2 + ... + n*(n+1)/(2 + ...))))) = Sum_{k >= 1} (-1)^(k+1)/(k*(k + 1)) = 2*log(2) - 1.
E.g.f.: (2*log(x+1)-x)/(x-1)^2. - Vaclav Kotesovec, Oct 21 2012

A060524 Triangle read by rows: T(n,k) = number of degree-n permutations with k odd cycles, k=0..n, n >= 0.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 5, 0, 1, 9, 0, 14, 0, 1, 0, 89, 0, 30, 0, 1, 225, 0, 439, 0, 55, 0, 1, 0, 3429, 0, 1519, 0, 91, 0, 1, 11025, 0, 24940, 0, 4214, 0, 140, 0, 1, 0, 230481, 0, 122156, 0, 10038, 0, 204, 0, 1, 893025, 0, 2250621, 0, 463490, 0, 21378, 0, 285, 0, 1, 0

Offset: 0

Vladeta Jovovic, Apr 01 2001

The row polynomials t(n,x):=Sum_{k=0..n} T(n,k)*x^k satisfy the recurrence relation t(n,x) = x*t(n-1,x) + ((n-1)^2)*t(n-2,x); t(-1,x)=0, t(0,x)=1. - Wolfdieter Lang, see above.
This is an example of a Sheffer triangle (coefficient triangle for Sheffer polynomials). In the umbral calculus (see the Roman reference given under A048854) s(n,x) := Sum_{k=0..n} T(n,k)*x^k would be called Sheffer polynomials for (1/cosh(t),tanh(t)), which translates to the e.g.f. for column number k>=0 given by (1/sqrt(1-x^2))*((arctanh(x))^k)/k!. The e.g.f. given below is rewritten in this Sheffer context as (1/sqrt(1-x^2))*exp(y*log(sqrt((1+x)/(1-x))))= (1/sqrt(1-x^2))*exp(y*arctanh(x)). The rows of the Jabotinsky type triangle |A049218| provide the coefficients of the associated polynomials. - Wolfdieter Lang, Feb 24 2005
The solution of the differential-difference relation f(n+1,x)= (d/dx)f(n,x) + (n^2)*f(n-1,x), n >= 1, with inputs f(0,x) and f(1,x) = (d/dx)f(0,x) is f(n,x) = t(n,d_x)*f(0,x), with the differential operator d_x:=d/dx and the row polynomials t(n,x) defined above. This problem appears in a computation of thermo field dynamics where f(0,x)=1/cosh(x). See the triangle A060081. - Wolfdieter Lang, Feb 24 2005
The inverse of the Sheffer matrix T with elements T(n,k) is the Sheffer matrix A060081. - Wolfdieter Lang, Jul 22 2005
T(n,k)=0 if n-k= 1(mod 2), else T(n,k) = sum of M2(n,p), p from {1,...,A000041(n)} restricted to partitions with exactly k odd parts and any nonnegative number of even parts. For the M2-multinomial numbers in A-St order see A036039(n,p). - Wolfdieter Lang, Aug 07 2007

			Triangle begins:
  [1],
  [0, 1],
  [1, 0, 1],
  [0, 5, 0, 1],
  [9, 0, 14, 0, 1],
  [0, 89, 0, 30, 0, 1],
  [225, 0, 439, 0, 55, 0, 1],
  [0, 3429, 0, 1519, 0, 91, 0, 1],
  [11025, 0, 24940, 0, 4214, 0, 140, 0, 1],
  [0, 230481, 0, 122156, 0, 10038, 0, 204, 0, 1],
  [893025, 0, 2250621, 0, 463490, 0, 21378, 0, 285, 0, 1],
  [0, 23941125, 0, 14466221, 0, 1467290, 0, 41778, 0, 385, 0, 1],
  ...
Signed version begins:
  [1],
  [0, 1],
  [-1, 0, 1],
  [0, -5, 0, 1],
  [9, 0, -14, 0, 1],
  [0, 89, 0, -30, 0, 1],
  [-225, 0, 439, 0, -55, 0, 1],
  [0, -3429, 0, 1519, 0, -91, 0, 1],
  ...
From _Peter Bala_, Feb 23 2024: (Start)
Maple can verify the following series for Pi:
Row 1 polynomial R(1, x) = x:
Pi = 3 + 4*Sum_{n >= 1} (-1)^(n+1)/((2*n + 1)*R(1, 2*n)*R(1, 2*n+2)).
Row 3 polynomial R(3, x) = 5*x + x^3:
(3/2)^2 * Pi = 7 + 4*(3^4)*Sum_{n >= 1} (-1)^(n+1)/((2*n + 1)*R(3, 2*n)*R(3, 2*n+2)).
Row 5 polynomial R(5, x) = 89*x + 30*x^3 + x^5:
((3*5)/(2*4))^2 * Pi = 11 + 4*(3*5)^4*Sum_{n >= 1} (-1)^(n+1)/((2*n + 1)*R(5, 2*n)*R(5, 2*n+2)). (End)

Alois P. Heinz, Rows n = 0..140, flattened
Marcelo Aguiar, Sarah Brauner, and Vic Reiner, Configuration spaces and peak representations, Sém. Lotharingien Comb., Proc. 36th Conf. Formal Power Series Alg. Comb. (2024) Vol. 91B, Art. No. 13. See p. 12.
Peter Bala, Meixner polynomials and Brouncker's continued fraction for Pi
Paul L. Butzer and Tom H. Koornwinder, Josef Meixner: His life and his orthogonal polynomials, Indagationes Mathematicae, Volume 30, Issue 1, January 2019, Pages 250-264.
Adel Hamdi and Jiang Zeng, Orthogonal polynomials and operator orderings, J. Math. Phys., 51:043506, (2010); arXiv:1006.0808 [math.CO], 2010.
Eric Weisstein's World of Mathematics, Meixner polynomial of the first kind.

Cf. A060338, A060523, A094368, A028353 (col 1), A103916 (col 2), A103917 (col 3), A103918 (col 4).

Cf. A111594 (associated Sheffer polynomials), A142979, A142983.

with(combinat):
b:= proc(n, i) option remember; expand(`if`(n=0, 1, `if`(i<1, 0,
      add(multinomial(n, n-i*j, i$j)*(i-1)!^j/j!*b(n-i*j, i-1)*
      `if`(irem(i, 2)=1, x^j, 1), j=0..n/i))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n$2)):
seq(T(n), n=0..12);  # Alois P. Heinz, Mar 09 2015
# alternative
A060524 := proc(n,k)
    option remember;
    if nR. J. Mathar, Jul 06 2023

nn = 6; Range[0, nn]! CoefficientList[
   Series[(1 - x^2)^(-1/2) ((1 + x)/(1 - x))^(y/2), {x, 0, nn}], {x, y}] // Grid  (* Geoffrey Critzer, Aug 28 2012 *)

E.g.f.: (1+x)^((y-1)/2)/(1-x)^((y+1)/2).
T(n, k) = T(n-1, k-1) + ((n-1)^2)*T(n-2, k); T(-1, k):=0, T(n, -1):=0, T(0, 0)=1, T(n, k)=0 if nWolfdieter Lang, see above.
The Meixner polynomials defined by S_0(x)=1, S_1(x) = x; S_{n+1}(x) = x*S_n(x) - n^2*S_{n-1}(x) give a signed version of this triangle (cf. A060338). - N. J. A. Sloane, May 30 2013
From Peter Bala, Apr 10 2024: (Start)
The n-th row polynomial R(n, x) satisfies
(4*n + 2)*R(n, x) = (x + 1)*R(n, x+2) - (x - 1)*R(n, x-2).
Series for Pi involving the row polynomials R(n, x): for n >= 0 there holds
((2*n + 1)!!/(2^n*n!))^2 * Pi = (4*n + 3) + 4*((2*n + 1)!!^4) * Sum_{k >= 1} (-1)^(k+1)/((2*k + 1)*R(2*n+1, 2*k)*R(2*n+1, 2*k+2)). Cf. A142979 and A142983.
R(2*n, 0) = A001147(n)^2 = A001818(n); R(2*n+1, 0) = 0.
R(n, 1) = n! = A000142(n).
R(2*n, 2) = (4*n + 1)*A001147(n)^2 = (4*n + 1)*((2*n)!/(2^n*n!))^2;
R(2*n+1, 2) = 2*A001447(n+1)^2 = 2*(2*n + 1)!^2/(n!^2*4^n).
R(n, 3) = (2*n + 1)*n! = A007680(n). (End)

cp's OEIS Frontend

A002162 Decimal expansion of the natural logarithm of 2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A008288 Square array of Delannoy numbers D(i,j) (i >= 0, j >= 0) read by antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Python

Sage

Formula

Extensions

A143003 a(0) = 0, a(1) = 1, a(n+1) = (2*n+1)*(n^2+n+5)*a(n) - n^6*a(n-1).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A024167 a(n) = n!*(1 - 1/2 + 1/3 - ... + c/n), where c = (-1)^(n+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A142992 Square array, read by ascending antidiagonals, of the crystal ball sequences for the root lattices of type C_n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A142983 a(1) = 1, a(2) = 2, a(n+2) = 2*a(n+1) + (n + 1)*(n + 2)*a(n).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

A143003 a(0) = 0, a(1) = 1, a(n+1) = (2n+1)(n^2+n+5)a(n) - n^6a(n-1).

A142983 a(1) = 1, a(2) = 2, a(n+2) = 2a(n+1) + (n + 1)(n + 2)*a(n).

A142995 a(0) = 0, a(1) = 1, a(n+1) = (2n^2 + 2n + 3)a(n) - n^4a(n-1), n >= 1.

A142999 a(0) = 0, a(1) = 1; for n > 1, a(n+1) = (2n + 1)a(n) + n^4*a(n-1).

A142988 a(1) = 1, a(2) = 3, a(n+2) = 3a(n+1) + (n + 1)(n + 3)*a(n).