cp's OEIS Frontend

Number of edges in the join of two complete graphs of order n^2 and n, K_n^2 * K_n. - Roberto E. Martinez II, Jan 07 2002

The general formula for alternating sums of powers is in terms of the Swiss-Knife polynomials P(n,x) A153641 2^(-n-1)(P(n,1)-(-1)^k P(n,2k+1)). Thus a(k) = |2^(-5)(P(4,1)-(-1)^k P(4,2k+1))|. - Peter Luschny, Jul 12 2009

Define an infinite symmetric array by T(n,m) = n*(n-1) + m for 0 <= m <= n and T(n,m) = T(m,n), n >= 0. Then a(n) is the sum of terms in the top left (n+1) X (n+1) subarray: a(n) = Sum_{r=0..n} Sum_{c=0..n} T(r,c). - J. M. Bergot, Jul 05 2013

a(n) is the sum of all positive numbers less than A002378(n). - J. M. Bergot, Aug 30 2013

Except the first term, these are triangular numbers that remain triangular when divided by their index, e.g., 66 divided by 11 gives 6. - Waldemar Puszkarz, Sep 14 2017

a(n) is the semiperimeter of the unique primitive Pythagorean triple such that (a-b+c)/2 = T(n) = A000217(n). Its long leg and hypotenuse are consecutive natural numbers and the triple is (2*T(n) - 1, 2*T(n)*(T(n) - 1), 2*T(n)*(T(n) - 1) + 1). - Miguel-Ángel Pérez García-Ortega, May 27 2025

a(n) is divisible by A000537(n) if and only n is congruent to 1 mod 3 (see A016777) - Artur Jasinski, Oct 10 2007

This sequence is related to A000540 by a(n) = n*A000540(n) - Sum_{i=0..n-1} A000540(i). - Bruno Berselli, Apr 26 2010

Kekulé numbers for certain benzenoids.

6th column of the table of Narayana numbers A001263. - Zerinvary Lajos, Jun 18 2007

Sequence provided by binomial(n-1,m)*binomial(n,m)/(m+1) for m=5 and n>5 (these numbers are also called Runyon numbers, see T. Koshy in References). - Vincenzo Librandi, Sep 04 2014

These generalized unsigned Lah numbers are the instance L[2,1] of the Sheffer triangles called L[d,a], with integers d >= 1 and integers 0 <= a < d with gcd(d,a) = 1. The standard unsigned Lah numbers are L[1,0] = A271703.

The Sheffer structure of L[d,a] is ((1 - d*t)^(-2*a/d), t/(1 - d*t)). This follows from the defining property

risefac[d,a](x, n) = Sum_{m=0..n} L[d,a](n, m)*fallfac[d,a](x, m), where risefac[d,a](x, n):= Product_{0..n-1} (x + (a+d*j)) for n >= 1 and risefac[d,a](x, 0) := 1, and fallfac[d,a](x, n):= Product_{0..n-1} (x - (a+d*j)) = for n >= 1 and fallfac[d,a](x, 0) := 1. Such rising and falling factorials arise in the generalization of Stirling numbers of both kinds S2[d,a] and S1[d,a]. See the Peter Bala link under A143395 for these falling factorials called there [t;a,b,c]_n with t=x, a=d, b=0, c=a.

In matrix notation: L[d,a] = S1phat[d,a].S2hat[d,a] with the unsigned scaled Stirling1 and the scaled Stirling2 generalizations with Sheffer structures S1phat[d,a] = ((1 - d*t)^(-a/d), -(1/d)*(log(1 - d*t))) and S2hat[d,a] = (exp(a*t), (1/d)*(exp(d*t) - 1). See, e.g., S1phat[2,1] = A028338 and S2hat[2,1] = A039755.

The a- and z-sequences for these Sheffer matrices have e.g.f.s 1 + d*t and ((1 + d*t)/t)*(1 - (1 + d*t)^(-2*a/d)), respectively. See a W. Lang link under A006232 for these types of sequences.

E.g.f. of row polynomials R[d,a](n, x) := Sum_{m=0..n} L[d,a](n, m)*x^m

(1 - d*x)^(-2*a/d)*exp(t*x/(1 - d*x)) (this is the e.g.f. for the triangle).

E.g.f. of column m: (1 - d*t)^(-2*a/d)*(t/(1 - d*t))^m/m, m >= 0.

Meixner type identity for (monic) row polynomials: (D_x/(1 + d*D_x)) * R[d,a](n, x) = n*R[d,a](n-1, x), n >= 1, with R[d,a](0, x) = 1. The series in the differentiations D_x = d/dx terminates.

General Sheffer recurrence for row polynomials (see the Roman reference, p. 50, Corollary 3.7.2, rewritten for the present Sheffer notation):

R[d,a](n, x) = [(2*a+x)*1 + 2*d*(a + x)*D_x + d^2*x*(D_x)^2]*R[d,a](n-1, x), n >= 1, with R[d,a](0, x) = 1.

The inverse matrix L^(-1)[d,a] is Sheffer (g[d,a](-t), -f[d,a](-t)) with L[d,a] Sheffer (g[d,a](t), f[d,a](t)) from above. This means (see the column e.g.f. of Sheffer matrices) that L^(-1)[d,a](n, m) = (-1)^(n-m)*L[d,a](n, m). Therefore, the recurrence relations can easily be rewritten for L^(-1)[d,a] by replacing a -> -a and d -> -d.

fallfac[d,a](x, n) = Sum_{m=0..n} L^(-1)[d,a](n, m)*risefac[d,a](x, m), n >= 0.

From Wolfdieter Lang, Aug 12 2017: (Start)

The Sheffer row polynomials R[d,a](n, x) belong to the Boas-Buck class and satisfy therefore the Boas-Buck identity (see the reference, and we use the notation of Rainville, Theorem 50, p. 141, adapted to an exponential generating function) (E_x - n*1)*R[d,a](n, x) = - n*(2*a*1 + d*E_x) * Sum_{k=0..n-1} d^k*R(d,a;n-1-k,x)/(n-1-k)!, with E_x = x*d/dx (Euler operator).

This implies a recurrence for the sequence of column m: L[d,a](n, m) = (n!*(2*a + d*m)/(n-m))*Sum_{p=0..n-1-m} d^p*L[d,a](n-1-p, m)/(n-1-p)!, for n > m>=0, and input L[d,a](m, m) = 1. For the present [d,a] = [2,1] instance see the formula and example sections. (End)

From Wolfdieter Lang, Sep 14 2017: (Start)

The diagonal sequences are 2^D*D!*(binomial(m+D, m))^2, m >= 0, for D >= 0 (main diagonal D = 0). From the o.g.f.s obtained via Lagrange's theorem. See the second W. Lang link below for the general Sheffer case.

The o.g.f. of the diagonal D sequence is 2^D*D!*Sum_{m=0..D} A008459(D, m)*x^m /(1- x)^(2*D + 1), D >= 0. (End)

It appears that this is also the matrix square of unsigned triangle of coefficients of Laguerre polynomials n!*L_n(x), abs(A021009(n, k)). - Ali Pourzand, Mar 10 2025 [This observation is correct. - Peter Luschny, Mar 10 2025]

The row sums of this triangle are b(n) = A005259(n), for n >= 0. This sequence b was used in R. Apéry's 1979 proof of the irrationality of Zeta(3). See A005259 for references and links.

Row polynomials R(n, x) := Sum_{k=0..n} T(n, k)*x^k = hypergeometric([-n, -n, n+1, n+1], [1, 1, 1], x), hence b(n) = hypergeometric([-n, -n, n+1, n+1], [1, 1, 1], 1) (see the formula in A005259 given by K. A. Penson. This is the solution to Exercise 2.14 of the Koepf reference given there, p. 29).

The length of row n of this irregular triangle is 2*n+1.

A008459 gives the squares of the entries of Pascal's triangle A007318.

The e.g.f. of the (n+1)-th diagonal sequence of the squares of Pascal's triangle (A008459) is EDP2(n, t) = Sum_{m=0..n} binomial(n+m, m)^2*t^m/m!, for n >= 0. It turns out to be EDP2(n, t) = exp(t)*Sum_{k=0..2*n} T(n, k)*t^k/k!.

This has been computed from the corresponding o.g.f.: GDP2(n, x) = Sum_{m=0..n} binomial(n+m, m)^2*x^m, which is GDP2(n, x) = Sum_{m=0..n} binomial(n,m)^2*x^m / (1 - x)^(2*n+1) (see the triangle A008459, and comments in A288876 on how to compute these o.g.f.s). To obtain the e.g.f.s from the o.g.f.s the formulas (23) - (25) of the W. Lang link given in A060187 have been used.