A001720 -id:A001720 - cp's OEIS frontend

A162995 A scaled version of triangle A162990.

1, 3, 1, 12, 4, 1, 60, 20, 5, 1, 360, 120, 30, 6, 1, 2520, 840, 210, 42, 7, 1, 20160, 6720, 1680, 336, 56, 8, 1, 181440, 60480, 15120, 3024, 504, 72, 9, 1, 1814400, 604800, 151200, 30240, 5040, 720, 90, 10, 1

Offset: 1

Johannes W. Meijer, Jul 27 2009

We get this scaled version of triangle A162990 by dividing the coefficients in the left hand columns by their 'top-values' and then taking the square root.
T(n,k) = A173333(n+1,k+1), 1 <= k <= n. - Reinhard Zumkeller, Feb 19 2010
T(n,k) = A094587(n+1,k+1), 1 <= k <= n. - Reinhard Zumkeller, Jul 05 2012

			The first few rows of the triangle are:
[1]
[3, 1]
[12, 4, 1]
[60, 20, 5, 1]

Reinhard Zumkeller, Rows n = 1..150 of triangle, flattened

Cf. A094587.
A056542(n) equals the row sums for n>=1.
A001710, A001715, A001720, A001725, A001730, A049388, A049389, A049398, A051431 are related to the left hand columns.
A000012, A009056, A002378, A007531, A052762, A052787, A053625 and A159083 are related to the right hand columns.

a162995 n k = a162995_tabl !! (n-1) !! (k-1)
a162995_row n = a162995_tabl !! (n-1)
a162995_tabl = map fst $ iterate f ([1], 3)
   where f (row, i) = (map (* i) row ++ [1], i + 1)
-- Reinhard Zumkeller, Jul 04 2012

a := proc(n, m): (n+1)!/(m+1)! end: seq(seq(a(n, m), m=1..n), n=1..9); # Johannes W. Meijer, revised Nov 23 2012

Table[(n+1)!/(m+1)!, {n, 10}, {m, n}] (* Paolo Xausa, Mar 31 2024 *)

a(n,m) = (n+1)!/(m+1)! for n = 1, 2, 3, ..., and m = 1, 2, ..., n.

A167557 The lower left triangle of the ED1 array A167546.

Original entry on oeis.org

1, 1, 4, 2, 12, 32, 6, 48, 160, 384, 24, 240, 960, 2688, 6144, 120, 1440, 6720, 21504, 55296, 122880, 720, 10080, 53760, 193536, 552960, 1351680, 2949120, 5040, 80640, 483840, 1935360, 6082560, 16220160, 38338560, 82575360

Offset: 1

Johannes W. Meijer, Nov 10 2009

We discovered that the numbers that appear in the lower left triangle of the ED1 array A167546 (m <= n) behave in a regular way, see the formula below. This rather simple regularity doesn't show up in the upper right triangle of the ED1 array (m > n).

			The first few triangle rows are:
[1]
[1, 4]
[2, 12, 32]
[6, 48, 160, 384]
[24, 240, 960, 2688, 6144]
[120, 1440, 6720, 21504, 55296, 122880]

G. C. Greubel, Table of n, a(n) for the first 50 rows

A167546 is the ED1 array.
A047053 and A167558 are the first two right hand triangle columns.
A000142, 4*A001710 (n>=2), 32*A001720, 384*A001730, 6144*A049389, 122880*A051431 are the first six left hand triangle columns.
A167559 equals the row sums.

a := proc(n, m): 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)! end: seq(seq(a(n, m), m=1..n), n=1..8);  # Johannes W. Meijer, revised Nov 23 2012

Flatten[Table[(4^(m-1) (m-1)!(n+m-2)!)/(2m-2)!,{n,10},{m,n}]] (* Harvey P. Dale, Sep 29 2013 *)

a(n,m) = 4^(m-1)*(m-1)!*(n+m-2)!/(2*m-2)!.

A094646 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -2, 1, 2, -3, 1, 0, 2, -3, 1, 0, 2, -1, -2, 1, 0, 4, 0, -5, 0, 1, 0, 12, 4, -15, -5, 3, 1, 0, 48, 28, -56, -35, 7, 7, 1, 0, 240, 188, -252, -231, 0, 42, 12, 1, 0, 1440, 1368, -1324, -1638, -231, 252, 114, 18, 1, 0, 10080, 11016, -7900, -12790, -3255, 1533, 1050, 240, 25, 1

Offset: 0

Vladeta Jovovic, May 17 2004

Triangle T(n,k), 0 <= k <= n, read by rows, given by [ -2, 1, -1, 2, 0, 3, 1, 4, 2, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 23 2006
From Wolfdieter Lang, Jun 23 2011: (Start)
The row polynomials s(n,x):=Sum_{k=0..n} T(n,k)*x^k satisfy risefac(x-2,n)=s(n,x), with the rising factorials risefac(x-2,n):=Product_{j=0..n-1} (x-2+j), n >= 1, risefac(x-2,0)=1. Compare this with the formula risefac(x,n)=|S1|(n,x), with the row polynomials |S1|(n,x) of A132393 (unsigned Stirling1).
This is the third triangle of an a-family of Sheffer arrays, call them |S1|(a), with e.g.f. of the row polynomials |S1|(a;x;z) = ((1-z)^a)*exp(-x*log(1-z)). In the notation showing the column e.g.f.s this is Sheffer ((1-z)^a,-log(1-z)). In the umbral notation (see the Roman reference, given under A094645) this is called Sheffer for (exp(a*t),1-exp(-t)). For a=0 this becomes the unsigned Stirling1 triangle |S1|(0) = A132393 with row polynomials |S1|(0;n,x) =: s1(n,x).
E.g.f. column number k (with leading zeros): ((1-x)^a)*((-log(1-x))^k)/k!, k >= 0.
E.g.f. for row sums is (1-x)^(a-1), and the e.g.f. for the alternating row sums is (1-x)^(a+1).
Row polynomial recurrence:
  |S1|(a;n,x)=(x+(n-1-a))*|S1|(a;n-1,m), |S1|(a;0,x)=1.
Meixner identity (see the reference under A060338):
  |S1|(a;n,x) - |S1|(a;n,x-1) = n*|S1|(a;n-1,x), n >= 1,
Also (from the corollary 3.7.2 on p. 50 of the Roman reference): |S1|(a;n,x) = (x-a)*|S1|(a;n-1,x+1), n >= 1.
Recurrence: |S1|(a;n,k) = |S1|(a;n-1,k-1) + (n-(a+1))*|S1|(a;n-1,k); |S1|(a;n,k)=0 if n < m, |S1|(a;n,-1)=0, |S1|(a;0,0)=1.
Connection to |Stirling1|=|S1|(0):
  |S1|(a;n,k) = Sum_{p=0..a} |S1|(a;a,p)*abs(Stirling1(n-a,k-p)), n >= a.
The exponential convolution identity is
  |S1|(a;n,x+y) = Sum_{k=0..n} binomial(n,k)*|S1|(a;k,y)*s1(n-k,x), n >= 0, with symmetry x <-> y.
The Sheffer a- and z-sequences are (see the W. Lang link under A006232): Sha(a;n)=A164555(n)/A027642(n) (independent of a) with e.g.f. x/(1-exp(-x)), and the z-sequence has e.g.f. (exp(a*x)-1)/(exp(-x)-1).
The inverse Sheffer matrix has e.g.f. exp(a*z)*exp(x*(1-exp(-z))), in short notation (exp(a*z),1-exp(-z)),
  (or in umbral notation ((1-t)^a,-log(1-t))).
(End)

			Triangle begins
   1;
  -2,  1;
   2, -3,  1;
   0,  2, -3,  1;
   0,  2, -1, -2,  1;
   0,  4,  0, -5,  0,  1;
   ...
risefac(x-2,3) = (x-2)*(x-1)*x = 2*x-3*x^2+x^3.
-1 = T(4,2) = T(3,1) + 1*T(3,2) =  2 + (-3).
T(4,3) = 2*abs(S1(2,3)) - 3*abs(S1(2,2)) + 1*abs(S1(2,1)) = 2*0 - 3*1 + 1*1 = -2.

Cf. A049444, A049458, A094645.

A094646_row := n -> seq((-1)^(n-k)*coeff(expand(pochhammer(x-n+3, n)), x, k), k=0..n): seq(print(A094646_row(n)), n = 0..6); # Peter Luschny, May 16 2013

Flatten[ Table[ CoefficientList[ Pochhammer[x-2, n], x], {n, 0, 10}]] (* Jean-François Alcover, Sep 26 2011 *)

E.g.f.: (1-y)^(2-x).
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000142(n+1), A001710(n+2), A001715(n+3), A001720(n+4), A001725(n+5), A001730(n+6), A049388(n), A049389(n), A049398(n), A051431(n) for x = 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 respectively. - Philippe Deléham, Nov 13 2007
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then |T(n,i)| = |f(n,i,-2)|, for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
From Wolfdieter Lang, Jun 23 2011: (Start)
risefac(x-2,n) = Sum_{k=0..n} T(n,k)*x^k, n >= 0, with the rising factorials (see a comment above).
Recurrence: T(n,k) = T(n-1,k-1) + (n-3)*T(n-1,k); T(n,k)=0 if n < m, T(n,-1)=0, T(0,0)=1.
T(n,k) = 2*abs(S1(n-2,k)) - 3*abs(S1(n-2,k-1)) + abs(S1(n-2,k-2)), n >= 2, with S1(n,k) = Stirling1(n,k) = A048994(n,k).
E.g.f. column number k (with leading zeros):
((1-x)^2)*((-log(1-x))^k)/k!, k >= 0.
E.g.f. for row sums is 1-x, i.e., [1,-1,0,0,...],
  and the e.g.f. for the alternating row sums is (1-x)^3. i.e., [1,-3,3,1,0,0,...]. (End)

A144355 Partition number array, called M31(5), related to A049353(n,m)= |S1(5;n,m)| (generalized Stirling triangle).

Original entry on oeis.org

1, 5, 1, 30, 15, 1, 210, 120, 75, 30, 1, 1680, 1050, 1500, 300, 375, 50, 1, 15120, 10080, 15750, 9000, 3150, 9000, 1875, 600, 1125, 75, 1, 151200, 105840, 176400, 220500, 35280, 110250, 63000, 78750, 7350, 31500, 13125, 1050, 2625, 105, 1, 1663200, 1209600, 2116800

Offset: 1

Wolfdieter Lang Oct 09 2008, Oct 28 2008

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k) =: M31(5;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42,...].
Fifth member (K=5) in the family M31(K) of partition number arrays.
If M31(5;n,k) is summed over those k with fixed number of parts m one obtains the unsigned triangle |S1(5)|:= A049353.

			[1];[5,1];[30,15,1];[210,120,75,30,1];[1680,1050,1500,300,375,50,1];...
a(4,3)= 75 = 3*|S1(5;2,1)|^2. The relevant partition of 4 is (2^2).

W. Lang, First 10 rows of the array and more.
W. Lang, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) 09.3.3.

A049378 (row sums).

A144354 (M31(4) array), A144356 (M31(6) array).

a(n,k)=(n!/product(e(n,k,j)!*j!^(e(n,k,j),j=1..n))*product(|S1(5;j,1)|^e(n,k,j),j=1..n)= M3(n,k)*product(|S1(5;j,1)|^e(n,k,j),j=1..n) with |S1(5;n,1)|= A001720(n+3) = (n+3)!/4!, n>=1 and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. M3(n,k)=A036040.

A303675 Triangle read by rows: coefficients in the sum of odd powers as expressed by Faulhaber's theorem, T(n, k) for n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 6, 1, 120, 30, 1, 5040, 1680, 126, 1, 362880, 151200, 17640, 510, 1, 39916800, 19958400, 3160080, 168960, 2046, 1, 6227020800, 3632428800, 726485760, 57657600, 1561560, 8190, 1, 1307674368000, 871782912000, 210680870400, 22313491200, 988107120, 14217840, 32766, 1

Offset: 1

Kolosov Petro, May 08 2018

T(n,k) are the coefficients in an identity due to Faulhaber: Sum_{j=0..n} j^(2*m-1) = Sum_{k=1..m} T(m,k) binomial(n+k, 2*k). See the Knuth reference, page 10.

More explicitly, Faulhaber's theorem asserts that, given integers n >= 0, m >= 1 and odd, Sum_{k=1..n} k^m = Sum_{k=1..(m+1)/2} C(n+k,n-k)*[(1/k)*Sum_{j=0..k-1} (-1)^j*C(2*k,j)*(k-j)^(m+1)]. The coefficients T(m, k) are indicated by square brackets. Sums similar to this inner part are A304330, A304334, A304336; however, these triangles are (0,0)-based and lead to equivalent but slightly more systematic representations. - Peter Luschny, May 12 2018

			The triangle begins (see the Knuth reference p. 10):
         1;
         6,          1;
       120,         30,         1;
      5040,       1680,       126,        1;
    362880,     151200,     17640,      510,       1;
  39916800,   19958400,   3160080,   168960,    2046,    1;
6227020800, 3632428800, 726485760, 57657600, 1561560, 8190, 1;
.
Let S(n, m) = Sum_{j=1..n} j^m. Faulhaber's formula gives for m = 7 (m odd!):
F(n, 7) = 5040*C(n+4, 8) + 1680*C(n+3, 6) + 126*C(n+2, 4) + C(n+1, 2).
Faulhaber's theorem asserts that for all n >= 1 S(n, 7) = F(n, 7).
If n = 43 the common value is 1600620805036.

John H. Conway and Richard Guy, The Book of Numbers, Springer (1996), p. 107.

Donald E. Knuth, Johann Faulhaber and Sums of Powers, arXiv:9207222 [math.CA], 1992.
Petro Kolosov, Polynomial identities involving central factorial numbers, GitHub, 2024. See pp. 3, 6.

First column is a bisection of A000142, second column is a bisection of A001720.
Row sums give A100868.
Cf. A008955, A008957, A036969 and A304330, A304334, A304336.

T := proc(n,k) local m; m := n-k;
2*(2*m+1)!*add((-1)^(j+m)*(j+1)^(2*n)/((j+m+2)!*(m-j)!), j=0..m) end:
seq(seq(T(n, k), k=1..n), n=1..8); # Peter Luschny, May 09 2018

(* After Peter Luschny's above formula. *)
T[n_, k_] := (1/(n-k+1))*Sum[(-1)^j*Binomial[2*(n-k+1), j]*((n-k+1) - j)^(2*n), {j, 0, n-k+1}]; Column[Table[T[n, k], {n, 1, 10}, {k, 1, n}], Center]

def A303675(n, k): return factorial(2*(n-k)+1)*A008957(n, k)
for n in (1..7): print([A303675(n, k) for k in (1..n)]) # Peter Luschny, May 10 2018

T(n, k) = (2*(n-k)+1)!*A008957(n, k), n >= 1, 1 <= k <= n.

T(n, k) = (1/m)*Sum_{j=0..m} (-1)^j*binomial(2*m,j)*(m-j)^(2*n) where m = n-k+1. - Peter Luschny, May 09 2018

New name by Peter Luschny, May 10 2018

A062261 Fourth (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 32, 720, 14400, 277200, 5322240, 103783680, 2075673600, 42810768000, 913296384000, 20183850086400, 462393656524800, 10981849342464000, 270322445352960000, 6893222356500480000, 181981070211612672000

Offset: 0

Wolfdieter Lang, Jun 19 2001

Harry J. Smith, Table of n, a(n) for n=0..100
Index entries for sequences related to Laguerre polynomials

Cf. A001720, A062199, A062260.

[Factorial(n+3)*Binomial(n+7,7)/6: n in [0..30]]; // G. C. Greubel, May 13 2018

Table[(n+3)!*Binomial[n+7, 7]/3!, {n, 0, 30}] (* G. C. Greubel, May 13 2018 *)

{ f=2; for (n=0, 100, f*=n + 3; write("b062261.txt", n, " ", f*binomial(n + 7, 7)/6) ) } \\ Harry J. Smith, Aug 03 2009

E.g.f.: (1+21*x+63*x^2+35*x^3)/(1-x)^11.
a(n) = A062140(n+3, 3).
a(n) = (n+3)!*binomial(n+7, 7)/3!.
If we define f(n,i,x) =  Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-3)=(-1)^(n-1)*f(n,3,-8), (n>=3). - Milan Janjic, Mar 01 2009

A144890 Partition number array, called M31hat(5).

Original entry on oeis.org

1, 5, 1, 30, 5, 1, 210, 30, 25, 5, 1, 1680, 210, 150, 30, 25, 5, 1, 15120, 1680, 1050, 900, 210, 150, 125, 30, 25, 5, 1, 151200, 15120, 8400, 6300, 1680, 1050, 900, 750, 210, 150, 125, 30, 25, 5, 1, 1663200, 151200, 75600, 50400, 44100, 15120, 8400, 6300, 5250, 4500

Offset: 1

Wolfdieter Lang Oct 09 2008, Oct 28 2008

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k) =: M31hat(5;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42,...].
Fourth member (K=5) in the family M31hat(K) of partition number arrays.
If M31hat(5;n,k) is summed over those k with fixed number of parts m one obtains the unsigned triangle S1hat(5):= A144891.

			[1];[5,1];[30,5,1];[210,30,25,5,1];[1680,210,150,30,25,5,1];...
a(4,3)= 25 = |S1(5;2,1)|^2. The relevant partition of 4 is (2^2).

W. Lang, First 10 rows of the array and more.
W. Lang, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) 09.3.3.

Cf. A144892 (row sums).

Cf. A144885 (M31hat(4) array). A144891 (S1hat(5)).

a(n,k) = product(|S1(5;j,1)|^e(n,k,j),j=1..n) with |S1(5;n,1)| = A049353(n,1) = A001720(n+3) = [1,5,30,210,1680,...] = (n+3)!/4!, n>=1 and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n.

A153271 Triangle T(n, k) = Product_{j=0..k} (j*n + prime(m)), with T(n, 0) = prime(m) and m = 3, read by rows.

Original entry on oeis.org

5, 5, 30, 5, 35, 315, 5, 40, 440, 6160, 5, 45, 585, 9945, 208845, 5, 50, 750, 15000, 375000, 11250000, 5, 55, 935, 21505, 623645, 21827575, 894930575, 5, 60, 1140, 29640, 978120, 39124800, 1838865600, 99298742400, 5, 65, 1365, 39585, 1464645, 65909025, 3493178325, 213083877825, 14702787569925

Offset: 0

Roger L. Bagula, Dec 22 2008

Row sums are {5, 35, 355, 6645, 219425, 11640805, 917404295, 101177741765, 14919432040765, 2839006665525525, 677815000136926955, ...}.

			Triangle begins as:
  5;
  5, 30;
  5, 35, 315;
  5, 40, 440,  6160;
  5, 45, 585,  9945, 208845;
  5, 50, 750, 15000, 375000, 11250000;
  5, 55, 935, 21505, 623645, 21827575, 894930575;

G. C. Greubel, Rows n = 0..100 of triangle, flattened

Cf. A153271 (m=2), this sequence (m=3), A153272 (m=4).
Sequences related to m values:
m=2: A001147, A001710, A008545, A032031, A047056, A144739, A144758.
m=3: A001720, A007696, A008543, A034000, A051577, A052562, A147585, A147625, A147629.

m:=3;
function T(n,k)
  if k eq 0 then return NthPrime(m);
  else return (&*[j*n + NthPrime(m): j in [0..k]]);
  end if; return T; end function;
[T(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Dec 03 2019

m:=3; seq(seq(`if`(k=0, ithprime(m), mul(j*n + ithprime(m), j=0..k)), k=0..n), n=0..10); # G. C. Greubel, Dec 03 2019

T[n_, k_, m_]:= If[k==0, Prime[m], Product[j*n + Prime[m], {j,0,k}]];
Table[T[n,k,3], {n,0,10}, {k,0,n}]//Flatten

T(n,k) = my(m=3); if(k==0, prime(m), prod(j=0,k, j*n + prime(m)) ); \\ G. C. Greubel, Dec 03 2019

def T(n, k):
    m=3
    if (k==0): return nth_prime(m)
    else: return product(j*n + nth_prime(m) for j in (0..k))
[[T(n, k) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Dec 03 2019

T(n, k) = Product_{j=0..k} (j*n + prime(m)), with T(n, 0) = prime(m) and m = 3.

Edited by G. C. Greubel, Dec 03 2019

A172455 The case S(6,-4,-1) of the family of self-convolutive recurrences studied by Martin and Kearney.

Original entry on oeis.org

1, 7, 84, 1463, 33936, 990542, 34938624, 1445713003, 68639375616, 3676366634402, 219208706540544, 14397191399702118, 1032543050697424896, 80280469685284582812, 6725557192852592984064, 603931579625379293509683

Offset: 1

N. J. A. Sloane, Nov 20 2010

nonn

			G.f. = x + 7*x^2 + 84*x^3 + 1463*x^4 + 33936*x^5 + 990542*x^6 + 34938624*x^7 + ...
a(2) = 7 since (6*2 - 4) * a(2-1) - (a(1) * a(2-1)) = 7.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
R. J. Martin and M. J. Kearney, An exactly solvable self-convolutive recurrence, Aequat. Math., 80 (2010), 291-318. see p. 307.
R. J. Martin and M. J. Kearney, An exactly solvable self-convolutive recurrence, arXiv:1103.4936 [math.CO], 2011.
NIST Digital Library of Mathematical Functions, Airy Functions.
A. N. Stokes, Continued fraction solutions of the Riccati equation, Bull. Austral. Math. Soc. Vol. 25 (1982), 207-214.
Eric Weisstein's World of Mathematics, Airy Functions, contains the definitions of Ai(x), Bi(x).

Cf. A000079 S(1,1,-1), A000108 S(0,0,1), A000142 S(1,-1,0), A000244 S(2,1,-2), A000351 S(4,1,-4), A000400 S(5,1,-5), A000420 S(6,1,-6), A000698 S(2,-3,1), A001710 S(1,1,0), A001715 S(1,2,0), A001720 S(1,3,0), A001725 S(1,4,0), A001730 S(1,5,0), A003319 S(1,-2,1), A005411 S(2,-4,1), A005412 S(2,-2,1), A006012 S(-1,2,2), A006318 S(0,1,1), A047891 S(0,2,1), A049388 S(1,6,0), A051604 S(3,1,0), A051605 S(3,2,0), A051606 S(3,3,0), A051607 S(3,4,0), A051608 S(3,5,0), A051609 S(3,6,0), A051617 S(4,1,0), A051618 S(4,2,0), A051619 S(4,3,0), A051620 S(4,4,0), A051621 S(4,5,0), A051622 S(4,6,0), A051687 S(5,1,0), A051688 S(5,2,0), A051689 S(5,3,0), A051690 S(5,4,0), A051691 S(5,5,0), A053100 S(6,1,0), A053101 S(6,2,0), A053102 S(6,3,0), A053103 S(6,4,0), A053104 S(7,1,0), A053105 S(7,2,0), A053106 S(7,3,0), A062980 S(6,-8,1), A082298 S(0,3,1), A082301 S(0,4,1), A082302 S(0,5,1), A082305 S(0,6,1), A082366 S(0,7,1), A082367 S(0,8,1), A105523 S(0,-2,1), A107716 S(3,-4,1), A111529 S(1,-3,2), A111530 S(1,-4,3), A111531 S(1,-5,4), A111532 S(1,-6,5), A111533 S(1,-7,6), A111546 S(1,0,1), A111556 S(1,1,1), A143749 S(0,10,1), A146559 S(1,1,-2), A167872 S(2,-3,2), A172450 S(2,0,-1), A172485 S(-1,-2,3), A177354 S(1,2,1), A292186 S(4,-6,1), A292187 S(3, -5, 1).

a[1] = 1; a[n_]:= a[n] = (6*n-4)*a[n-1] - Sum[a[k]*a[n-k], {k, 1, n-1}]; Table[a[n], {n, 1, 20}] (* Vaclav Kotesovec, Jan 19 2015 *)

{a(n) = local(A); if( n<1, 0, A = vector(n); A[1] = 1; for( k=2, n, A[k] = (6 * k - 4) * A[k-1] - sum( j=1, k-1, A[j] * A[k-j])); A[n])} /* Michael Somos, Jul 24 2011 */

S(v1, v2, v3, N=16) = {
  my(a = vector(N)); a[1] = 1;
  for (n = 2, N, a[n] = (v1*n+v2)*a[n-1] + v3*sum(j=1,n-1,a[j]*a[n-j])); a;
};
S(6,-4,-1)
\\ test: y = x*Ser(S(6,-4,-1,201)); 6*x^2*y' == y^2 - (2*x-1)*y - x
\\ Gheorghe Coserea, May 12 2017

a(n) = (6*n - 4) * a(n-1) - Sum_{k=1..n-1} a(k) * a(n-k) if n>1. - Michael Somos, Jul 24 2011
G.f.: x / (1 - 7*x / (1 - 5*x / (1 - 13*x / (1 - 11*x / (1 - 19*x / (1 - 17*x / ... )))))). - Michael Somos, Jan 03 2013
a(n) = 3/(2*Pi^2)*int((4*x)^((3*n-1)/2)/(Ai'(x)^2+Bi'(x)^2), x=0..inf), where Ai'(x), Bi'(x) are the derivatives of the Airy functions. [Vladimir Reshetnikov, Sep 24 2013]
a(n) ~ 6^n * (n-1)! / (2*Pi) [Martin + Kearney, 2011, p.16]. - Vaclav Kotesovec, Jan 19 2015
6*x^2*y' = y^2 - (2*x-1)*y - x, where y(x) = Sum_{n>=1} a(n)*x^n. - Gheorghe Coserea, May 12 2017
G.f.: x/(1 - 2*x - 5*x/(1 - 7*x/(1 - 11*x/(1 - 13*x/(1 - ... - (6*n - 1)*x/(1 - (6*n + 1)*x/(1 - .... Cf. A062980. - Peter Bala, May 21 2017

A062262 Fifth (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 45, 1350, 34650, 831600, 19459440, 454053600, 10702692000, 256864608000, 6307453152000, 158947819430400, 4118193503424000, 109818493424640000, 3015784780968960000, 85303626661693440000

Offset: 0

Wolfdieter Lang, Jun 19 2001

Harry J. Smith, Table of n, a(n) for n = 0..100
Index entries for sequences related to Laguerre polynomials

Cf. A001720, A062199, A062260, A062261.

[Factorial(n+4)*Binomial(n+8,8)/24: n in [0..30]]; // G. C. Greubel, May 13 2018

Table[(n+4)!*Binomial[n+8, 8]/4!, {n, 0, 30}] (* G. C. Greubel, May 13 2018 *)
With[{nn=20},CoefficientList[Series[(1+32x+168x^2+224x^3+70x^4)/(1-x)^13,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Mar 18 2025 *)

{ f=6; for (n=0, 100, f*=n + 4; write("b062262.txt", n, " ", f*binomial(n + 8, 8)/24) ) } \\ Harry J. Smith, Aug 03 2009

E.g.f.: (1+32*x+168*x^2+224*x^3+70*x^4)/(1-x)^13.
a(n) = A062140(n+4, 4).
a(n) = (n+4)!*binomial(n+8, 8)/4!.
If we define f(n,i,x)= Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-4) = (-1)^n*f(n,4,-9), (n>=4). - Milan Janjic, Mar 01 2009

cp's OEIS Frontend

A162995 A scaled version of triangle A162990.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

A167557 The lower left triangle of the ED1 array A167546.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A094646 Generalized Stirling number triangle of first kind.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

A144355 Partition number array, called M31(5), related to A049353(n,m)= |S1(5;n,m)| (generalized Stirling triangle).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A303675 Triangle read by rows: coefficients in the sum of odd powers as expressed by Faulhaber's theorem, T(n, k) for n >= 1, 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

Extensions

A062261 Fourth (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A144890 Partition number array, called M31hat(5).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs