A284861
Triangle read by rows: T(n, k) = S2[3,1](n, k)*k! with the Sheffer triangle S2[3,1] = (exp(x), exp(3*x) -1) given in A282629.
Original entry on oeis.org
1, 1, 3, 1, 15, 18, 1, 63, 216, 162, 1, 255, 1890, 3564, 1944, 1, 1023, 14760, 52650, 68040, 29160, 1, 4095, 109458, 659340, 1516320, 1487160, 524880, 1, 16383, 790776, 7578522, 27624240, 46539360, 36741600, 11022480, 1, 65535, 5633730, 82902204, 450057384, 1158993360, 1535798880, 1014068160, 264539520
Offset: 0
The triangle T(n, k) begins
n\k 0 1 2 3 4 5 6 7 ...
0: 1
1: 1 3
2: 1 15 18
3: 1 63 216 162
4: 1 255 1890 3564 1944
5: 1 1023 14760 52650 68040 29160
6: 1 4095 109458 659340 1516320 1487160 524880
7: 1 16383 790776 7578522 27624240 46539360 36741600 11022480
...
row n=8: 1 65535 5633730 82902204 450057384 1158993360 1535798880 1014068160 264539520,
row n=9: 1 262143 39829320 879725610 6845572440 25294754520 50042059200 54561276000 30951123840 7142567040,
row n=10: 1 1048575 280378098 9155719980 99549149040 507399658920 1406104706160 2251231315200 2083248720000 1035672220800 214277011200.
------------------------------------------------------------------
T(2, 1) = -1 + 4^2 = 15 = 2*A225117(2,2) + 1*A225117(2,1) = 2*1 + 1*13.
- P. Bala, Deformations of the Hadamard product of power series
- M. Dukes, C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
- Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.
Cf.
A000012,
A002450,
A016223,
A021874,
A032031,
A111577,
A122803,
A131689, |
A151919|,
A225117,
A282629,
A019538,
A145901.
-
Table[Sum[Binomial[k, m] (-1)^(k - m) (1 + 3m)^n, {m, 0, k}], {n, 0, 10}, {k, 0, n}]// Flatten (* Indranil Ghosh, Apr 09 2017 *)
-
for(n=0, 10, for(k=0, n, print1(sum(m=0, k, binomial(k, m) * (-1)^(k - m)*(1 + 3*m)^n),", "); ); print();) \\ Indranil Ghosh, Apr 09 2017
-
# Indranil Ghosh, Apr 09 2017
from sympy import binomial
for n in range(11):
print([sum([binomial(k, m)*(-1)**(k - m)*(1 + 3*m)**n for m in range(k + 1)]) for k in range(n + 1)])
A037961
a(n) = n^2*(n+1)*(n+3)!/48.
Original entry on oeis.org
0, 1, 30, 540, 8400, 126000, 1905120, 29635200, 479001600, 8083152000, 142702560000, 2637143308800, 50999300352000, 1031319184896000, 21785854970880000, 480178027929600000, 11029155770400768000
Offset: 0
- Identity (1.19) in H. W. Gould, Combinatorial Identities, Morgantown, 1972; page 3.
A135456
Number of surjections from an n-element set onto a seven-element set.
Original entry on oeis.org
5040, 141120, 2328480, 29635200, 322494480, 3162075840, 28805736960, 248619571200, 2060056318320, 16540688324160, 129568848121440, 995210916336000, 7524340159588560, 56163512390086080, 414847224363337920
Offset: 7
-
LinearRecurrence[{28, -322, 1960, -6769, 13132, -13068, 5040}, {5040, 141120, 2328480, 29635200, 322494480, 3162075840, 28805736960}, 25] (* G. C. Greubel, Oct 14 2016 *)
A002869
Largest number in n-th row of triangle A019538.
Original entry on oeis.org
1, 1, 2, 6, 36, 240, 1800, 16800, 191520, 2328480, 30240000, 479001600, 8083152000, 142702560000, 2731586457600, 59056027430400, 1320663933388800, 30575780537702400, 783699448602470400, 21234672840116736000, 591499300737945600000
Offset: 0
- 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).
- Danny Rorabaugh, Table of n, a(n) for n = 0..400 (first 251 terms from Reinhard Zumkeller)
- Jens Gulin and Kalle Åström, Alternative implementations of the Auxiliary Duplicating Permutation Invariant Training, Proc Work-in-Progress Papers at 14th Int'l Conf. Indoor Positioning Indoor Nav. (IPIN-WiP 2024). See p. 6.
- Victor Meally, Comparison of several sequences given in Motzkin's paper "Sorting numbers for cylinders...", letter to N. J. A. Sloane, N. D.
- T. S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176. [Annotated, scanned copy]
- OEIS Wiki, Sorting numbers
A000670 gives sum of terms in n-th row.
-
a002869 0 = 1
a002869 n = maximum $ a019538_row n
-- Reinhard Zumkeller, Dec 15 2013
-
f := proc(n) local t1, k; t1 := 0; for k to n do if t1 < A019538(n, k) then t1 := A019538(n, k) fi; od; t1; end;
-
A019538[n_, k_] := k!*StirlingS2[n, k]; f[0] = 1; f[n_] := Module[{t1, k}, t1 = 0; For[k = 1, k <= n, k++, If[t1 < A019538[n, k], t1 = A019538[n, k]]]; t1]; Table[f[n], {n, 0, 20}] (* Jean-François Alcover, Dec 26 2013, after Maple *)
-
def A002869(n):
return max(factorial(k)*stirling_number2(n,k) for k in range(1,n+1))
[A002869(i) for i in range(1, 20)] # Danny Rorabaugh, Oct 10 2015
A196347
Triangle T(n, k) read by rows, T(n, k) = n!*binomial(n, k).
Original entry on oeis.org
1, 1, 1, 2, 4, 2, 6, 18, 18, 6, 24, 96, 144, 96, 24, 120, 600, 1200, 1200, 600, 120, 720, 4320, 10800, 14400, 10800, 4320, 720, 5040, 35280, 105840, 176400, 176400, 105840, 35280, 5040, 40320, 322560, 1128960, 2257920, 2822400, 2257920, 1128960, 322560, 40320
Offset: 0
Triangle begins:
1;
1, 1;
2, 4, 2;
6, 18, 18, 6;
24, 96, 144, 96, 24;
120, 600, 1200, 1200, 600, 120;
...
- G. C. Greubel, Table of n, a(n) for n = 0..495
- P. Bala, Deformations of the Hadamard product of power series
- Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
- M. Dukes, C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
-
/* As triangle */ [[Factorial(n)*Binomial(n, k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 28 2015
-
Table[n!*Binomial[n, j], {n, 0, 30}, {j, 0, n}] (* G. C. Greubel, Sep 27 2015 *)
-
factorial(n)*binomial(n,k) # Danny Rorabaugh, Sep 27 2015
A225472
Triangle read by rows, k!*S_3(n, k) where S_m(n, k) are the Stirling-Frobenius subset numbers of order m; n >= 0, k >= 0.
Original entry on oeis.org
1, 2, 3, 4, 21, 18, 8, 117, 270, 162, 16, 609, 2862, 4212, 1944, 32, 3093, 26550, 72090, 77760, 29160, 64, 15561, 230958, 1031940, 1953720, 1662120, 524880, 128, 77997, 1941030, 13429962, 39735360, 57561840, 40415760, 11022480, 256, 390369, 15996222, 165198852
Offset: 0
[n\k][0, 1, 2, 3, 4, 5, 6 ]
[0] 1,
[1] 2, 3,
[2] 4, 21, 18,
[3] 8, 117, 270, 162,
[4] 16, 609, 2862, 4212, 1944,
[5] 32, 3093, 26550, 72090, 77760, 29160,
[6] 64, 15561, 230958, 1031940, 1953720, 1662120, 524880.
-
SF_SO := proc(n, k, m) option remember;
if n = 0 and k = 0 then return(1) fi;
if k > n or k < 0 then return(0) fi;
m*k*SF_SO(n-1, k-1, m) + (m*(k+1)-1)*SF_SO(n-1, k, m) end:
seq(print(seq(SF_SO(n, k, 3), k=0..n)), n = 0..5);
-
EulerianNumber[n_, k_, m_] := EulerianNumber[n, k, m] = (If[ n == 0, Return[If[k == 0, 1, 0]]]; Return[(m*(n-k)+m-1)*EulerianNumber[n-1, k-1, m] + (m*k+1)*EulerianNumber[n-1, k, m]]); SFSO[n_, k_, m_] := Sum[ EulerianNumber[n, j, m]*Binomial[j, n-k], {j, 0, n}]; Table[ SFSO[n, k, 3], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 29 2013, translated from Sage *)
-
@CachedFunction
def EulerianNumber(n, k, m) :
if n == 0: return 1 if k == 0 else 0
return (m*(n-k)+m-1)*EulerianNumber(n-1, k-1, m)+ (m*k+1)*EulerianNumber(n-1, k, m)
def SF_SO(n, k, m):
return add(EulerianNumber(n, j, m)*binomial(j, n - k) for j in (0..n))
for n in (0..6): [SF_SO(n, k, 3) for k in (0..n)]
A225473
Triangle read by rows, k!*S_4(n, k) where S_m(n, k) are the Stirling-Frobenius subset numbers of order m; n >= 0, k >= 0.
Original entry on oeis.org
1, 3, 4, 9, 40, 32, 27, 316, 672, 384, 81, 2320, 9920, 13824, 6144, 243, 16564, 127680, 326400, 337920, 122880, 729, 116920, 1536992, 6428160, 11642880, 9584640, 2949120, 2187, 821356, 17842272, 114866304, 324065280, 453304320, 309657600, 82575360, 6561
Offset: 0
[n\k][0, 1, 2, 3, 4, 5, 6 ]
[0] 1,
[1] 3, 4,
[2] 9, 40, 32,
[3] 27, 316, 672, 384,
[4] 81, 2320, 9920, 13824, 6144,
[5] 243, 16564, 127680, 326400, 337920, 122880,
[6] 729, 116920, 1536992, 6428160, 11642880, 9584640, 2949120.
-
SF_SO := proc(n, k, m) option remember;
if n = 0 and k = 0 then return(1) fi;
if k > n or k < 0 then return(0) fi;
m*k*SF_SO(n-1, k-1, m) + (m*(k+1)-1)*SF_SO(n-1, k, m) end:
seq(print(seq(SF_SO(n, k, 4), k=0..n)), n = 0..5);
-
EulerianNumber[n_, k_, m_] := EulerianNumber[n, k, m] = (If[ n == 0, Return[If[k == 0, 1, 0]]]; Return[(m*(n-k)+m-1)*EulerianNumber[n-1, k-1, m] + (m*k+1)*EulerianNumber[n-1, k, m]]); SFSO[n_, k_, m_] := Sum[ EulerianNumber[n, j, m]*Binomial[j, n-k], {j, 0, n}]; Table[ SFSO[n, k, 4], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 29 2013, translated from Sage *)
-
@CachedFunction
def EulerianNumber(n, k, m) :
if n == 0: return 1 if k == 0 else 0
return (m*(n-k)+m-1)*EulerianNumber(n-1, k-1, m)+ (m*k+1)*EulerianNumber(n-1, k, m)
def SF_SO(n, k, m):
return add(EulerianNumber(n, j, m)*binomial(j, n - k) for j in (0..n))
for n in (0..6): [SF_SO(n, k, 4) for k in (0..n)]
A281478
Central coefficients of Joffe's central differences of zero (assuming offset 0 and T(n,k) extended to 0 <= k <= n in A241171).
Original entry on oeis.org
1, 1, 126, 126720, 494053560, 5283068427000, 126301275727704000, 5896518025761483120000, 488276203972584492344880000, 66735969985432035804226510800000, 14236685931434801591697761172512160000, 4533351707244550464920840944132383960960000, 2077486542875366717627638783543223150778585600000
Offset: 0
A281479
Central coefficients of the polynomials defined in A278073.
Original entry on oeis.org
1, 1, 1364, 42771456, 10298900437056, 11287986820196486400, 41397337338743872194508800, 414528538783792919989135797964800, 9808376038359632185170127842947907993600, 492228239722024416239987973400425228541016064000
Offset: 0
A281480
Central coefficients of the polynomials defined in A278074.
Original entry on oeis.org
1, 1, 16510, 17651304000, 286988816206755000, 35284812773848049161035000, 21735699944364325706210750640600000, 51125456932397825107093888817556205542000000, 378603085421985456745667562645258531056443927230000000, 7641597761030055776217194099395682779700673105680593973250000000
Offset: 0
Comments