A303486 -id:A303486 - cp's OEIS frontend

A303487 a(n) = n! * [x^n] 1/(1 - 4*x)^(n/4).

1, 1, 12, 231, 6144, 208845, 8648640, 422463195, 23781703680, 1515973484025, 107941254220800, 8491022274509775, 731304510986649600, 68444451854354701125, 6916953288171902976000, 750681472158682148959875, 87076954662428278259712000, 10751175443940144673035200625

Offset: 0

Ilya Gutkovskiy, Apr 24 2018

nonn

			a(1) = 1;
a(2) = 2*6 = 12;
a(3) = 3*7*11 = 231;
a(4) = 4*8*12*16 = 6144;
a(5) = 5*9*13*17*21 = 208845, etc.

Index entries for sequences related to factorial numbers

Column k=4 of A303489.

Cf. A000407, A001813, A007696, A008545, A034176, A034177, A047053, A051617, A051618, A051619, A051620, A051621, A051622, A113551, A303486, A303488.

Table[n! SeriesCoefficient[1/(1 - 4 x)^(n/4), {x, 0, n}], {n, 0, 17}]
Table[Product[4 k + n, {k, 0, n - 1}], {n, 0, 17}]
Table[4^n Pochhammer[n/4, n], {n, 0, 17}]

a(n) = Product_{k=0..n-1} (4*k + n).
a(n) = 4^n*Gamma(5*n/4)/Gamma(n/4).
a(n) ~ 5^(5*n/4-1/2)*n^n/exp(n).

A303489 Square array A(n,k), n >= 0, k >= 1, read by antidiagonals: A(n,k) = n! * [x^n] 1/(1 - k*x)^(n/k).

Original entry on oeis.org

1, 1, 1, 1, 1, 6, 1, 1, 8, 60, 1, 1, 10, 105, 840, 1, 1, 12, 162, 1920, 15120, 1, 1, 14, 231, 3640, 45045, 332640, 1, 1, 16, 312, 6144, 104720, 1290240, 8648640, 1, 1, 18, 405, 9576, 208845, 3674160, 43648605, 259459200, 1, 1, 20, 510, 14080, 375000, 8648640, 152152000, 1703116800, 8821612800

Offset: 0

Ilya Gutkovskiy, Apr 24 2018

			Square array begins:
      1,      1,       1,       1,       1,       1,  ...
      1,      1,       1,       1,       1,       1,  ...
      6,      8,      10,      12,      14,      16,  ...
     60,    105,     162,     231,     312,     405,  ...
    840,   1920,    3640,    6144,    9576,   14080,  ...
  15120,  45045,  104720,  208845,  375000,  623645,  ...
=========================================================
A(1,1) = 1;
A(2,1) = 2*3 = 6;
A(3,1) = 3*4*5 = 60;
A(4,1) = 4*5*6*7 = 840;
A(5,1) = 5*6*7*8*9 = 15120, etc.
...
A(1,2) = 1;
A(2,2) = 2*4 = 8;
A(3,2) = 3*5*7 = 105;
A(4,2) = 4*6*8*10 = 1920;
A(5,2) = 5*7*9*11*13 = 45045, etc.
...
A(1,3) = 1;
A(2,3) = 2*5 = 10;
A(3,3) = 3*6*9 = 162;
A(4,3) = 4*7*10*13 = 3640;
A(5,3) = 5*8*11*14*17 = 104720, etc.
...

Index entries for sequences related to factorial numbers

Columns k=1..5 give A000407, A113551, A303486, A303487, A303488.
Main diagonal gives A061711.
Cf. A008279, A131182, A256268, A265609.

Table[Function[k, n! SeriesCoefficient[1/(1 - k x)^(n/k), {x, 0, n}]][j - n + 1], {j, 0, 9}, {n, 0, j}] // Flatten
Table[Function[k, Product[k i + n, {i, 0, n - 1}]][j - n + 1], {j, 0, 9}, {n, 0, j}] // Flatten
Table[Function[k, k^n Pochhammer[n/k, n]][j - n + 1], {j, 0, 9}, {n, 0, j}] // Flatten

A(n,k) = Product_{j=0..n-1} (k*j + n).

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).

A303488 a(n) = n! * [x^n] 1/(1 - 5*x)^(n/5).

Original entry on oeis.org

1, 1, 14, 312, 9576, 375000, 17873856, 1004306688, 65006637696, 4763494479744, 389812500000000, 35237024762075136, 3487065897634615296, 374960171943074285568, 43532820293400237735936, 5427359437500000000000000, 723181462895975365595529216, 102563963819340862347122245632

Offset: 0

Ilya Gutkovskiy, Apr 24 2018

nonn

			a(1) = 1;
a(2) = 2*7 = 14;
a(3) = 3*8*13 = 312;
a(4) = 4*9*14*19 = 9576;
a(5) = 5*10*15*20*25 = 375000, etc.

Index entries for sequences related to factorial numbers

Column k=5 of A303489.

Cf. A008546, A008548, A034300, A034301, A034323, A034325, A047055, A047056, A051687, A051688, A051689, A051690, A051691, A052562, A113551, A303486, A303487.

Table[n! SeriesCoefficient[1/(1 - 5 x)^(n/5), {x, 0, n}], {n, 0, 17}]
Table[Product[5 k + n, {k, 0, n - 1}], {n, 0, 17}]
Table[5^n Pochhammer[n/5, n], {n, 0, 17}]

a(n) = Product_{k=0..n-1} (5*k + n).
a(n) = 5^n*Gamma(6*n/5)/Gamma(n/5).
a(n) ~ 6^(6*n/5-1/2)*n^n/exp(n).

A384168 a(n) = 3^n * n! * binomial(4n/3,n) Sum_{k=1..n} 1/(n+3*k).

Original entry on oeis.org

1, 13, 234, 5566, 165944, 5966136, 251491120, 12169996912, 665146831680, 40530954643840, 2724842629685120, 200361647815660800, 15997170878205905920, 1378271357428552115200, 127459020533529062246400, 12593128815600367187507200, 1323895109721239722075136000

Offset: 1

Seiichi Manyama, May 21 2025

nonn

Cf. A098118, A384167, A384169.

Cf. A303486.

a(n) = sum(k=0, n, k*(n+3)^(k-1)*3^(n-k)*abs(stirling(n, k, 1)));

a(n) = Sum_{k=0..n} k * (n+3)^(k-1) * 3^(n-k) * |Stirling1(n,k)|.

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

A384256 a(n) = Product_{k=0..n-1} (n+3*k+1).

Original entry on oeis.org

1, 2, 18, 280, 6160, 174960, 6086080, 250490240, 11904278400, 641472832000, 38645634227200, 2573895458534400, 187787322731008000, 14894027431162880000, 1275931456704672768000, 117412145664335441920000, 11550258696757088788480000, 1209613643310990696210432000

Offset: 0

Seiichi Manyama, May 23 2025

Cf. A303486, A384257.

PARI
```
a(n) = prod(k=0, n-1, n+3*k+1);
```

def a(n): return 3^n*rising_factorial((n+1)/3, n)

a(n) = 3^n * RisingFactorial((n+1)/3,n).

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

A384257 a(n) = Product_{k=0..n-1} (n+3*k+2).

Original entry on oeis.org

1, 3, 28, 440, 9720, 276640, 9634240, 396809280, 18866848000, 1016990374400, 61283225203200, 4082333102848000, 297880548623257600, 23628360309345792000, 2024347339040266240000, 186294495108985303040000, 18327479444105919639552000, 1919453757320555804508160000

Offset: 0

Seiichi Manyama, May 23 2025

Cf. A303486, A384256.

PARI
```
a(n) = prod(k=0, n-1, n+3*k+2);
```

def a(n): return 3^n*rising_factorial((n+2)/3, n)

a(n) = 3^n * RisingFactorial((n+2)/3,n).

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

cp's OEIS Frontend

A303487 a(n) = n! * [x^n] 1/(1 - 4*x)^(n/4).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A303489 Square array A(n,k), n >= 0, k >= 1, read by antidiagonals: A(n,k) = n! * [x^n] 1/(1 - k*x)^(n/k).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

SageMath

Formula

A303488 a(n) = n! * [x^n] 1/(1 - 5*x)^(n/5).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A384168 a(n) = 3^n * n! * binomial(4*n/3,n) * Sum_{k=1..n} 1/(n+3*k).

Views

Author

Keywords

Crossrefs

Programs

PARI

Formula

A384256 a(n) = Product_{k=0..n-1} (n+3*k+1).

Views

Author

Keywords

Crossrefs

Programs

PARI

Sage

Formula

A384257 a(n) = Product_{k=0..n-1} (n+3*k+2).

Views

Author

Keywords

Crossrefs

Programs

PARI

Sage

Formula

A384168 a(n) = 3^n * n! * binomial(4n/3,n) Sum_{k=1..n} 1/(n+3*k).