A034000 -id:A034000 - cp's OEIS frontend

A048966 A convolution triangle of numbers obtained from A025748.

1, 3, 1, 15, 6, 1, 90, 39, 9, 1, 594, 270, 72, 12, 1, 4158, 1953, 567, 114, 15, 1, 30294, 14580, 4482, 1008, 165, 18, 1, 227205, 111456, 35721, 8667, 1620, 225, 21, 1, 1741905, 867834, 287199, 73656, 15075, 2430, 294, 24, 1, 13586859, 6857136, 2328183, 623106, 136323, 24354, 3465, 372, 27, 1

Offset: 1

Wolfdieter Lang

A generalization of the Catalan triangle A033184.

			Triangle begins:
     1;
     3,    1;
    15,    6,    1;
    90,   39,    9,    1;
   594,  270,   72,   12,    1;
  4158, 1953,  567,  114,   15,    1;

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened
W. Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.

Cf. A034000, A049213, A049223, A049224. a(n, 1)= A025748(n), a(n, 1)= 3^(n-1)*2*A034000(n-1)/n!, n >= 2. Row sums = A025756.

a048966 n k = a048966_tabl !! (n-1) !! (k-1)
a048966_row n = a048966_tabl !! (n-1)
a048966_tabl = [1] : f 2 [1] where
   f x xs = ys : f (x + 1) ys where
     ys = map (flip div x) $ zipWith (+)
          (map (* 3) $ zipWith (*) (map (3 * (x - 1) -) [1..]) (xs ++ [0]))
          (zipWith (*) [1..] ([0] ++ xs))
-- Reinhard Zumkeller, Feb 19 2014

a[n_, m_] /; n >= m >= 1 := a[n, m] = 3*(3*(n-1) - m)*a[n-1, m]/n + m*a[n-1, m-1]/n; a[n_, m_] /; n < m := 0; a[n_, 0] = 0; a[1, 1] = 1; Table[a[n, m], {n, 1, 10}, {m, 1, n}] // Flatten (* Jean-François Alcover, Apr 26 2011, after given formula *)

a(n, m) = 3*(3*(n-1)-m)*a(n-1, m)/n + m*a(n-1, m-1)/n, n >= m >= 1; a(n, m) := 0, n

G.f. for m-th column: ((1-(1-9*x)^(1/3))/3)^m.

a(n,m) = m/n * sum(k=0..n-m, binomial(k,n-m-k) * 3^k*(-1)^(n-m-k) * binomial(n+k-1,n-1)). - Vladimir Kruchinin, Feb 08 2011

A051608 a(n) = (3*n+8)!!!/8!!!.

Original entry on oeis.org

1, 11, 154, 2618, 52360, 1204280, 31311280, 908027120, 29056867840, 1016990374400, 38645634227200, 1584471003315200, 69716724145868800, 3276686034855833600, 163834301742791680000, 8683217992367959040000, 486260207572605706240000, 28689352246783736668160000

Offset: 0

Wolfdieter Lang

Related to A008544(n+1) ((3*n+2)!!! triple factorials).

Row m=8 of the array A(4; m,n) := ((3*n+m)(!^3))/m(!^3), m >= 0, n >= 0.

G. C. Greubel, Table of n, a(n) for n = 0..377

Cf. A032031, A007559(n+1), A034000(n+1), A034001(n+1), A051604, A051605, A051606, A051607, A051609 (rows m=0..9).

Cf. A008544.

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-3*x)^(11/3))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018

With[{nn = 30}, CoefficientList[Series[1/(1 - 3*x)^(11/3), {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Aug 15 2018 *)

x='x+O('x^30); Vec(serlaplace(1/(1-3*x)^(11/3))) \\ G. C. Greubel, Aug 15 2018

a(n) = ((3*n+8)(!^3))/8(!^3).
E.g.f.: 1/(1-3*x)^(11/3).
Sum_{n>=0} 1/a(n) = 1 + 9*(9*e)^(1/3)*(Gamma(11/3) - Gamma(11/3, 1/3)). - Amiram Eldar, Dec 23 2022

A112333 An invertible triangle of ratios of triple factorials.

Original entry on oeis.org

1, 2, 1, 10, 5, 1, 80, 40, 8, 1, 880, 440, 88, 11, 1, 12320, 6160, 1232, 154, 14, 1, 209440, 104720, 20944, 2618, 238, 17, 1, 4188800, 2094400, 418880, 52360, 4760, 340, 20, 1, 96342400, 48171200, 9634240, 1204280, 109480, 7820, 460, 23, 1, 2504902400

Offset: 0

Paul Barry, Sep 04 2005

First column is A008544. Second column is A034000. Third column is A051605. As a square array read by antidiagonals, columns have e.g.f. (1/(1-3x)^(2/3)) * (1/(1-3x))^k.

			Triangle begins
      1;
      2,    1;
     10,    5,    1;
     80,   40,    8,   1;
    880,  440,   88,  11,  1;
  12320, 6160, 1232, 154, 14, 1;
Inverse triangle A112334 begins
   1;
  -2,  1;
   0, -5,  1;
   0,  0, -8,   1;
   0,  0,  0, -11,   1;
   0,  0,  0,   0, -14,   1;
   0,  0,  0,   0,   0, -17, 1;

nmax:=8: for n from 0 to nmax do for k from 0 to n do if k<=n then T(n, k) := mul(3*k1-1, k1=1..n)/ mul(3*j-1, j=1..k) else T(n, k) := 0: fi: od: od: for n from 0 to nmax do seq(T(n, k), k=0..n) od: seq(seq(T(n, k), k=0..n), n=0..nmax); # Johannes W. Meijer, Jul 04 2011, revised Nov 23 2012

Number triangle T(n, k)=if(k<=n, Product{k=1..n, 3k-1}/Product{j=1..k, 3j-1}, 0); T(n, k)=if(k<=n, 3^(n-k)*(n-1/3)!/(k-1/3)!, 0).

A371077 Square array read by ascending antidiagonals: A(n, k) = 3^n*Pochhammer(k/3, n).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 4, 2, 1, 0, 28, 10, 3, 1, 0, 280, 80, 18, 4, 1, 0, 3640, 880, 162, 28, 5, 1, 0, 58240, 12320, 1944, 280, 40, 6, 1, 0, 1106560, 209440, 29160, 3640, 440, 54, 7, 1, 0, 24344320, 4188800, 524880, 58240, 6160, 648, 70, 8, 1

Offset: 0

Werner Schulte and Peter Luschny, Mar 10 2024

			The array starts:
  [0] 1,    1,     1,     1,     1,      1,      1,      1,      1, ...
  [1] 0,    1,     2,     3,     4,      5,      6,      7,      8, ...
  [2] 0,    4,    10,    18,    28,     40,     54,     70,     88, ...
  [3] 0,   28,    80,   162,   280,    440,    648,    910,   1232, ...
  [4] 0,  280,   880,  1944,  3640,   6160,   9720,  14560,  20944, ...
  [5] 0, 3640, 12320, 29160, 58240, 104720, 174960, 276640, 418880, ...
.
Seen as the triangle T(n, k) = A(n - k, k):
  [0] 1;
  [1] 0,       1;
  [2] 0,       1,      1;
  [3] 0,       4,      2,     1;
  [4] 0,      28,     10,     3,    1;
  [5] 0,     280,     80,    18,    4,   1;
  [6] 0,    3640,    880,   162,   28,   5,  1;
  [7] 0,   58240,  12320,  1944,  280,  40,  6, 1;
  [8] 0, 1106560, 209440, 29160, 3640, 440, 54, 7, 1;
.
Illustrating the LU decomposition of A:
    / 1                \   / 1 1 1 1 1 ... \   / 1   1   1    1    1 ... \
    | 0   1            |   |   1 2 3 4 ... |   | 0   1   2    3    4 ... |
    | 0   4   2        | * |     1 3 6 ... | = | 0   4  10   18   28 ... |
    | 0  28  24   6    |   |       1 4 ... |   | 0  28  80  162  280 ... |
    | 0 280 320 144 24 |   |         1 ... |   | 0 280 880 1944 3640 ... |
    | . . .            |   | . . .         |   | . . .                   |

Paolo Xausa, Table of n, a(n) for n = 0..11324 (first 150 antidiagonals, flattened).

Family m^n*Pochhammer(k/m, n): A094587 (m=1), A370419 (m=2), this sequence (m=3), A370915 (m=4).
Rows: A000012, A001477, A028552.
Columns: A000007, A007559, A008544, A032031, A007559 (shifted), A034000, A034001, A051604, A051605, A051606, A051607, A051608, A051609.
Cf. A303486 (main diagonal), A371079 (row sums of triangle), A371076, A371080.

A := (n, k) -> 3^n*pochhammer(k/3, n):
A := (n, k) -> local j; mul(3*j + k, j = 0..n-1):
# Read by antidiagonals:
T := (n, k) -> A(n - k, k): seq(seq(T(n, k), k = 0..n), n = 0..9);
seq(lprint([n], seq(T(n, k), k = 0..n)), n = 0..9);
# Using the generating polynomials of the rows:
P := (n, x) -> local k; add(Stirling1(n, k)*(-3)^(n - k)*x^k, k=0..n):
seq(lprint([n], seq(P(n, k), k = 0..9)), n = 0..5);
# Using the exponential generating functions of the columns:
EGFcol := proc(k, len) local egf, ser, n; egf := (1 - 3*x)^(-k/3);
ser := series(egf, x, len+2): seq(n!*coeff(ser, x, n), n = 0..len) end:
seq(lprint([k], EGFcol(k, 8)), k = 0..6);
# As a matrix product:
with(LinearAlgebra):
L := Matrix(7, 7, (n, k) -> A371076(n - 1,  k - 1)):
U := Matrix(7, 7, (n, k) -> binomial(n - 1, k - 1)):
MatrixMatrixMultiply(L, Transpose(U));

Table[3^(n-k)*Pochhammer[k/3, n-k], {n, 0, 10}, {k, 0, n}] (* Paolo Xausa, Mar 14 2024 *)

def A(n, k): return 3**n * rising_factorial(k/3, n)
def A(n, k): return (-3)**n * falling_factorial(-k/3, n)

A(n, k) = Product_{j=0..n-1} (3*j + k).
A(n, k) = A(n+1, k-3) / (k - 3) for k > 3.
A(n, k) = Sum_{j=0..n} Stirling1(n, j)*(-3)^(n - j)* k^j.
A(n, k) = k! * [x^k] (exp(x) * p(n, x)), where p(n, x) are the row polynomials of A371080.
E.g.f. of column k: (1 - 3*t)^(-k/3).
E.g.f. of row n: exp(x) * (Sum_{k=0..n} A371076(n, k) * x^k / (k!)).
Sum_{n>=0, k>=0} A(n, k) * x^k * t^n / (n!) = 1/(1 - x/(1 - 3*t)^(1/3)).
Sum_{n>=0, k>=0} A(n, k) * x^k * t^n /(n! * k!) = exp(x/(1 - 3*t)^(1/3)).
The LU decomposition of this array is given by the upper triangular matrix U which is the transpose of A007318 and the lower triangular matrix L = A371076, i.e., A(n, k) = Sum_{i=0..k} A371076(n, i) * binomial(k, i).

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

A081408 a(n) = (n+1)*a(n-5), with a(0)=a(1)=a(2)=a(3)=a(4)=1.

Original entry on oeis.org

1, 1, 1, 1, 1, 6, 7, 8, 9, 10, 66, 84, 104, 126, 150, 1056, 1428, 1872, 2394, 3000, 22176, 31416, 43056, 57456, 75000, 576576, 848232, 1205568, 1666224, 2250000, 17873856, 27143424, 39783744, 56651616, 78750000, 643458816, 1004306688, 1511782272

Offset: 0

Labos Elemer, Apr 01 2003

nonn

Quintic factorial sequences are generated by single 5-order recursion and appear in unified form.

			A008548, A034323, A034300, A034301, A034325 sequences are combed together as A081408(5n+r) with r=0,1,2,3,4.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000

Cf. A001147, A002866, A034001, A007599, A034000, A007696, A000407, A034176, A034177, A008548, A034323, A034300, A034301, A034325 [double, triple, quartic, quintic, factorial subsequences], generated together in A081405-A081408.

a:=[1,1,1,1,1];; for n in [6..40] do a[n]:=n*a[n-5]; od; a; # G. C. Greubel, Aug 15 2019

a081407 n = a081408_list !! n
a081407_list = 1 : 1 : 1 : 1 : zipWith (*) [5..] a081407_list
-- Reinhard Zumkeller, Jan 05 2012

[n le 5 select 1 else n*Self(n-5): n in [1..40]]; // G. C. Greubel, Aug 15 2019

a[0]=a[1]=a[2]=a[3]=a[4]=1; a[x_]:= (x+1)*a[x-5]; Table[a[n], {n, 40}]

m=30; v=concat([1,1,1,1,1], vector(m-5)); for(n=6, m, v[n]=n*v[n-5] ); v \\ G. C. Greubel, Aug 15 2019

def a(n):
    if (n<5): return 1
    else: return (n+1)*a(n-5)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 15 2019

Previous Showing 11-20 of 21 results. Next

cp's OEIS Frontend

A051606 a(n) = (3*n+6)!!!/6!!!, related to A032031 ((3*n)!!! triple factorials).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A051609 a(n) = (3*n+9)!!!/9!!!, related to A032031 ((3*n)!!! triple factorials).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A303486 a(n) = n! * [x^n] 1/(1 - 3*x)^(n/3).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A051607 a(n) = (3*n+7)!!!/7!!!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A048966 A convolution triangle of numbers obtained from A025748.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A051608 a(n) = (3*n+8)!!!/8!!!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A112333 An invertible triangle of ratios of triple factorials.

Views

Author

Keywords

Comments

Examples

Programs

A051606 a(n) = (3n+6)!!!/6!!!, related to A032031 ((3n)!!! triple factorials).

A051609 a(n) = (3n+9)!!!/9!!!, related to A032031 ((3n)!!! triple factorials).