cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 17 results. Next

A144422 a(n) = n!*A144416(n).

Original entry on oeis.org

1, 3, 62, 5052, 1087104, 487424520, 393702654960, 519740602925040, 1046019551260199040, 3046052768591313895680, 12322848899623787148556800, 67036205966744406283600262400, 477463744343745629221074425625600, 4354451197836949908463868034925593600
Offset: 0

Views

Author

David Applegate and N. J. A. Sloane, Dec 08 2008, corrected Jan 14 2009

Keywords

Comments

Number of placing of at least n objects but no more than 3*n distinct objects in n distinct boxes with the condition that in each box there is at least 1 object but no more than 3. - Adi Dani, May 04 2011

Crossrefs

Cf. A144416.

Programs

  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 3*n := t[n, k] = t[n - 1, k - 1] + (k - 1)*t[n - 1, k - 2] + (1/2)*(k - 1)*(k - 2)*t[n - 1, k - 3]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 3*n}]; Table[n!*a[n], {n, 0, 20}] (* Michael De Vlieger, Aug 04 2025, after Jean-François Alcover at A144416 *)

A001515 Bessel polynomial y_n(x) evaluated at x=1.

Original entry on oeis.org

1, 2, 7, 37, 266, 2431, 27007, 353522, 5329837, 90960751, 1733584106, 36496226977, 841146804577, 21065166341402, 569600638022431, 16539483668991901, 513293594376771362, 16955228098102446847, 593946277027962411007, 21992967478132711654106, 858319677924203716921141
Offset: 0

Views

Author

Keywords

Comments

For some applications it is better to start this sequence with an extra 1 at the beginning: 1, 1, 2, 37, 266, 2431, 27007, 353522, 5329837, ... (again with offset 0). This sequence now has its own entry - see A144301.
Number of partitions of {1,...,k}, n <= k <= 2n, into n blocks with no more than 2 elements per block. Restated, number of ways to use the elements of {1,...,k}, n <= k <= 2n, once each to form a collection of n sets, each having 1 or 2 elements. - Bob Proctor, Apr 18 2005, Jun 26 2006. E.g., for n=2 we get: (k=2): {1,2}; (k=3): {1,23}, {2,13}, {3,12}; (k=4): {12,34}, {13,24}, {14,23}, for a total of a(2) = 7 partitions.
Equivalently, number of sequences of n unlabeled items such that each item occurs just once or twice (cf. A105749). - David Applegate, Dec 08 2008
Numerator of (n+1)-th convergent to 1+tanh(1). - Benoit Cloitre, Dec 20 2002
The following Maple lines show how this sequence and A144505, A144498, A001514, A144513, A144506, A144514, A144507, A144301 are related.
f0:=proc(n) local k; add((n+k)!/((n-k)!*k!*2^k),k=0..n); end; [seq(f0(n),n=0..10)];
# that is this sequence
f1:=proc(n) local k; add((n+k+1)!/((n-k)!*k!*2^k),k=0..n); end; [seq(f1(n),n=0..10)];
# that is A144498
f2:=proc(n) local k; add((n+k+2)!/((n-k)!*k!*2^k),k=0..n); end; [seq(f2(n),n=0..10)];
# that is A144513; divided by 2 gives A001514
f3:=proc(n) local k; add((n+k+3)!/((n-k)!*k!*2^k),k=0..n); end; [seq(f3(n),n=0..10)];
# that is A144514; divided by 6 gives A144506
f4:=proc(n) local k; add((n+k+4)!/((n-k)!*k!*2^k),k=0..n); end; [seq(f4(n),n=0..10)];
# that divided by 24 gives A144507
a(n) is also the numerator of the continued fraction sequence beginning with 2 followed by 3 and the remaining odd numbers: [2,3,5,7,9,11,13,...]. - Gil Broussard, Oct 07 2009
Also, number of scenarios in the Gift Exchange Game when a gift can be stolen at most once. - N. J. A. Sloane, Jan 25 2017

Examples

			The first few Bessel polynomials are (cf. A001497, A001498):
  y_0 = 1
  y_1 = 1 +   x
  y_2 = 1 + 3*x +  3*x^2
  y_3 = 1 + 6*x + 15*x^2 + 15*x^3, etc.
G.f. = 1 + 2*x + 7*x^2 + 37*x^3 + 266*x^4 + 2431*x^5 + 27007*x^6 + 353522*x^7 + ...
		

References

  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
  • 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).

Crossrefs

See A144301 for other formulas and comments.
Row sums of Bessel triangle A001497 as well as of A001498.
Partial sums: A105748.
First differences: A144498.
Replace "sets" with "lists" in comment: A001517.
The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are this sequence, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Haskell
    a001515 = sum . a001497_row -- Reinhard Zumkeller, Nov 24 2014
    
  • Magma
    [(&+[Binomial(n+j, 2*j)*Catalan(j)*Factorial(j+1)/2^j: j in [0..n]]): n in [0..30]]; // G. C. Greubel, Sep 26 2023
    
  • Maple
    A001515 := proc(n) option remember; if n=0 then 1 elif n=1 then 2 else (2*n-1)*A001515(n-1)+A001515(n-2); fi; end;
    A001515:=proc(n) local k; add( (n+k)!/((n-k)!*k!*2^k),k=0..n); end;
    A001515:= n-> hypergeom( [n+1,-n],[],-1/2);
    bessel := proc(n,x) add(binomial(n+k,2*k)*(2*k)!*x^k/(k!*2^k),k=0..n); end;
  • Mathematica
    RecurrenceTable[{a[0]==1,a[1]==2,a[n]==(2n-1)a[n-1]+a[n-2]},a[n], {n,25}] (* Harvey P. Dale, Jun 18 2011 *)
    Table[Sum[BellY[n+1, k, (2 Range[n+1] - 3)!!], {k, n+1}], {n, 0, 20}] (* Vladimir Reshetnikov, Nov 09 2016 *)
  • PARI
    {a(n) = if( n<0, n = -1 - n); sum( k=0, n, (2*n - k)! / (k! * (n-k)!) * 2^(k-n))} /* Michael Somos, Apr 08 2012 */
    
  • SageMath
    [sum(binomial(n+j,2*j)*binomial(2*j,j)*factorial(j)//2^j for j in range(n+1)) for n in range(31)] # G. C. Greubel, Sep 26 2023

Formula

The following formulas can all be found in (or are easily derived from formulas in) Grosswald's book.
D-finite with recurrence: a(0) = 1, a(1) = 2; thereafter a(n) = (2*n-1)*a(n-1) + a(n-2).
E.g.f.: exp(1-sqrt(1-2*x))/sqrt(1-2*x).
a(n) = Sum_{ k = 0..n } binomial(n+k,2*k)*(2*k)!/(k!*2^k).
Equivalently, a(n) = Sum_{ k = 0..n } (n+k)!/((n-k)!*k!*2^k) = Sum_{ k = n..2n } k!/((2n-k)!*(k-n)!*2^(k-n)).
a(n) = Hypergeometric2F0( [n+1, -n] ; - ; -1/2).
a(n) = A105749(n)/n!.
a(n) ~ exp(1)*(2n)!/(n!*2^n) as n -> oo. [See Grosswald, p. 124]
a(n) = A144301(n+1).
G.f.: 1/(1-x-x/(1-x-2*x/(1-x-3*x/(1-x-4*x/(1-x-5*x/(1-.... (continued fraction). - Paul Barry, Feb 08 2009
From Michael Somos, Apr 08 2012: (Start)
a(-1 - n) = a(n).
(a(n+1) + a(n+2))^2 = a(n)*a(n+2) + a(n+1)*a(n+3) for all integer n. (End)
G.f.: 1/G(0) where G(k) = 1 - x - x*(2*k+1)/(1 - x - 2*x*(k+1)/G(k+1)); (continued fraction). - Sergei N. Gladkovskii, Oct 05 2012
E.g.f.: E(0)/(2*sqrt(1-2*x)), where E(k) = 1 + 1/(1 - 2*x/(2*x + (k+1)*(1+sqrt(1-2*x))/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 23 2013
G.f.: T(0)/(1-x), where T(k) = 1 - (k+1)*x/((k+1)*x - (1-x)^2/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 15 2013
a(n) = (2*BesselI(1/2, 1)+BesselI(3/2, 1))*BesselK(n+1/2, 1). - Jean-François Alcover, Feb 03 2014
a(n) = exp(1)*sqrt(2/Pi)*BesselK(1/2+n,1). - Gerry Martens, Jul 22 2015
From Peter Bala, Apr 14 2017: (Start)
a(n) = (1/n!)*Integral_{x = 0..inf} exp(-x)*x^n*(1 + x/2)^n dx.
E.g.f.: d/dx( exp(x*c(x/2)) ) = 1 + 2*x + 7*x^2/2! + 37*x^3/3! + ..., where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. (End)
From G. C. Greubel, Aug 16 2017: (Start)
a(n) = (1/2)_{n} * 2^n * hypergeometric1f1(-n; -2*n; 2).
G.f.: (1/(1-t))*hypergeometric2f0(1, 1/2; -; 2*t/(1-t)^2). (End)

Extensions

Extensively edited by N. J. A. Sloane, Dec 07 2008

A144508 a(n) = total number of partitions of [1, 2, ..., k] into exactly n blocks, each of size 1, 2, 3 or 4, for 0 <= k <= 4n.

Original entry on oeis.org

1, 4, 121, 18252, 7958726, 7528988476, 13130817809439, 38001495237579931, 169490425291053577442, 1102725620990181693266071, 10030550674270068548738783696, 123317200510025161580777179001154, 1993320784474917266370637900936051186, 41401645296339316791633672053851083955147
Offset: 0

Views

Author

David Applegate and N. J. A. Sloane, Dec 15 2008

Keywords

Comments

Also, number of scenarios in the Gift Exchange Game when a gift can be stolen at most 3 times. - N. J. A. Sloane, Jan 25 2017

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 4*n := t[n, k] = t[n - 1, k - 1] + (k - 1)*t[n - 1, k - 2] + (1/2)*(k - 1)*(k - 2)*t[n - 1, k - 3] + (1/6)*(k - 1)*(k - 2)*(k - 3)*t[n - 1, k - 4]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 4*n}]; Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 4*n, i!*polcoef(sum(j=1, 4, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

A144385 Triangle read by rows: T(n,k) is the number of partitions of [1, 2, ..., k] into exactly n blocks, each of size 1, 2 or 3 (n >= 0, 0 <= k <= 3n).

Original entry on oeis.org

1, 0, 1, 1, 1, 0, 0, 1, 3, 7, 10, 10, 0, 0, 0, 1, 6, 25, 75, 175, 280, 280, 0, 0, 0, 0, 1, 10, 65, 315, 1225, 3780, 9100, 15400, 15400, 0, 0, 0, 0, 0, 1, 15, 140, 980, 5565, 26145, 102025, 323400, 800800, 1401400, 1401400, 0, 0, 0, 0, 0, 0, 1, 21, 266, 2520, 19425, 125895, 695695, 3273270, 12962950, 42042000, 106506400, 190590400, 190590400
Offset: 0

Views

Author

David Applegate and N. J. A. Sloane, Dec 07 2008, Dec 17 2008

Keywords

Comments

Row n has 3n+1 entries.

Examples

			Triangle begins:
[1]
[0, 1, 1, 1]
[0, 0, 1, 3, 7, 10, 10]
[0, 0, 0, 1, 6, 25, 75, 175, 280, 280]
[0, 0, 0, 0, 1, 10, 65, 315, 1225, 3780, 9100, 15400, 15400]
[0, 0, 0, 0, 0, 1, 15, 140, 980, 5565, 26145, 102025, 323400, 800800, 1401400, 1401400]
		

Crossrefs

See A144399, A144402, A144417, A111246 for other versions of this triangle.
Column sums give A001680, row sums give A144416. Taking last nonzero entry in each row gives A025035.
Diagonals include A000217, A001296, A027778, A144516; also A025035.
A generalization of the triangle in A144331 (and in several other entries).
Cf. A144643.

Programs

  • Maple
    T := proc(n, k)
    option remember;
    if n = k then 1;
    elif k < n then 0;
    elif n < 1 then 0;
    else T(n - 1, k - 1) + (k - 1)*T(n - 1, k - 2) + 1/2*(k - 1)*(k - 2)*T(n - 1, k - 3);
    end if;
    end proc;
    for n from 0 to 12 do lprint([seq(T(n,k),k=0..3*n)]); od:
  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 3*n := t[n, k] = t[n-1, k-1] + (k-1)*t[n-1, k-2] + (1/2)*(k-1)*(k-2)*t[n-1, k-3]; t[, ] = 0; Table[t[n, k], {n, 0, 12}, {k, 0, 3*n}] // Flatten (* Jean-François Alcover, Jan 14 2014 *)

Formula

T(n, k) = T(n - 1, k - 1) + (k - 1)*T(n - 1, k - 2) + (1/2)*(k - 1)*(k - 2)*T(n - 1, k - 3).
E.g.f.: Sum_{ n >= 0, k >= 0 } T(n, k) y^n x^k / k! = exp( y*(x+x^2/2+x^3/6) ). That is, the coefficient of y^n is the e.g.f. for row n. E.g. the e.g.f. for row 2 is (1/2)*(x+x^2/2+x^3/6)^2 = 1*x^2/2! + 3*x^3/3! + 7*x^4/4! + 10*x^5/5! + 10*x^6/6!.

A144509 a(n) = total number of partitions of [1, 2, ..., k] into exactly n blocks, each of size 1, 2, ..., 5, for 0 <= k <= 5n.

Original entry on oeis.org

1, 5, 456, 405408, 1495388159, 15467641899285, 361207016885536095, 16557834064546698285496, 1350410785161120363519024709, 182141025850694258874753732988078, 38395944834298393218465758049745918098, 12093097322244029427838390643054170162054258, 5485321312184901565806045962453632525844948020084
Offset: 0

Views

Author

David Applegate and N. J. A. Sloane, Dec 15 2008

Keywords

Comments

Also, number of scenarios in the Gift Exchange Game when a gift can be stolen at most 4 times. - N. J. A. Sloane, Jan 25 2017

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 5*n := t[n, k] = Sum[(1/j!)*Product[k - m, {m, 1, j}]*t[n - 1, k - j - 1], {j, 0, 4}]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 5*n}]; Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 5*n, i!*polcoef(sum(j=1, 5, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

A149187 a(n) = total number of partitions of [1, 2, ..., k] into exactly n blocks, each of size 1, 2, ..., 6, for 0 <= k <= 6n.

Original entry on oeis.org

1, 6, 1709, 9268549, 295887993624, 34155922905682979, 10893033763705794846727, 8064519699524417149584982475, 12261371699318896159811165091392898, 34949877647533654983311522321749656046802, 174047342897498341701547082125166096889157924610, 1431472607165249058159939223685478666695036430843693596
Offset: 0

Views

Author

N. J. A. Sloane, May 13 2009

Keywords

Comments

Also, number of scenarios in the Gift Exchange Game when a gift can be stolen at most 5 times. - N. J. A. Sloane, Jan 25 2017

Crossrefs

Cf. A144512.
The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i, t) option remember; `if`(t*i add(b(k, 6, n), k=0..6*n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Sep 17 2015
  • Mathematica
    multinomial[n_, k_List] := n!/Times @@ (k!); b[n_, i_, t_] := b[n, i, t] = If[t*i < n, 0, If[n == 0, If[t == 0, 1, 0], Sum[b[n-i*j, i-1, t-j]* multinomial[n, Prepend[Array[i&, j], n-i*j]]/j!, {j, 0, Min[t, n/i]}]]]; a[n_] := Sum[b[k, 6, n], {k, 0, 6*n}];  Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Dec 06 2016 after Alois P. Heinz *)
  • PARI
    {a(n) = sum(i=n, 6*n, i!*polcoef(sum(j=1, 6, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

A281358 Number of scenarios in the Gift Exchange Game when a gift can be stolen at most 6 times.

Original entry on oeis.org

1, 7, 6427, 216864652, 60790021361170, 79397199549271412737, 350521520018942991464535019, 4247805448772073978048752721163278, 122022975450467092259059357046375920848764, 7449370563518425038119522091529589590475534631830
Offset: 0

Views

Author

N. J. A. Sloane, Jan 25 2017

Keywords

Comments

The result from the recurrence has been confirmed up to a(63) by using an optimized version of equation (23) in the Applegate-Sloane paper. - Lars Blomberg, Feb 01 2017

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i, t) option remember; `if`(t*i add(b(k, 7, n), k=0..7*n):
    seq(a(n), n=0..12);  # Alois P. Heinz, Feb 01 2017
  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 7*n := t[n, k] = Sum[(1/j!)*Product[k - m, {m, 1, j}]*t[n - 1, k - j - 1], {j, 0, 6}]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 7*n}]; Table[a[n], {n, 0, 10}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 7*n, i!*polcoef(sum(j=1, 7, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

Extensions

More terms from Lars Blomberg, Feb 01 2017

A281359 Number of scenarios in the Gift Exchange Game when a gift can be stolen at most 7 times.

Original entry on oeis.org

1, 8, 24301, 5165454442, 12845435390707724, 191739533381111401455478, 11834912423104188943497126664597, 2371013832433361706367594400829713564440, 1299618941291522676629215597535104557826094801396, 1716119248126070756229849154290399886241087778087554633612
Offset: 0

Views

Author

N. J. A. Sloane, Jan 25 2017

Keywords

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i, t) option remember; `if`(t*i add(b(k, 8, n), k=0..8*n):
    seq(a(n), n=0..12);  # Alois P. Heinz, Feb 01 2017
  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 8*n := t[n, k] = Sum[(1/j!)*Product[k - m, {m, 1, j}]*t[n - 1, k - j - 1], {j, 0, 7}]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 8*n}]; Table[a[n], {n, 0, 10}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 8*n, i!*polcoef(sum(j=1, 8, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

A281360 Number of scenarios in the Gift Exchange Game when a gift can be stolen at most 8 times.

Original entry on oeis.org

1, 9, 92368, 124762262630, 2774049143394729653, 476872353039366288373555323, 414678423576860263798348331987688320, 1383884737648788823775562903922773021277571568, 14584126149704606223764458141727351569547933381159988406, 419715170056359079715862408734598208208707081189266290220651371206
Offset: 0

Views

Author

N. J. A. Sloane, Jan 25 2017

Keywords

Comments

More than the usual number of terms are shown in the DATA field because there are the initial values needed for one of the recurrences.

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i, t) option remember; `if`(t*i add(b(k, 9, n), k=0..9*n):
    seq(a(n), n=0..12);  # Alois P. Heinz, Feb 01 2017
  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 9*n := t[n, k] = Sum[(1/j!)*Product[k - m, {m, 1, j}]*t[n - 1, k - j - 1], {j, 0, 8}]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 9*n}]; Table[a[n], {n, 0, 10}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 9*n, i!*polcoef(sum(j=1, 9, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019

A281361 Number of scenarios in the Gift Exchange Game when a gift can be stolen at most 9 times.

Original entry on oeis.org

1, 10, 352705, 3047235458767, 609542744597785306189, 1214103036523322674154687139158, 14963835327495031822418126706099787884130, 836883118002221273912672040462907783367741190535388, 170589804359366329173961838612841486616626580885839826818966688, 107640669875812795238625627484701500354901860426640161278022882392148747562, 185260259482556646382994900799988470488841686941141661692183483756531004879305895810561
Offset: 0

Views

Author

N. J. A. Sloane, Jan 25 2017

Keywords

Comments

More than the usual number of terms are shown in the DATA field because there are the initial values needed for one of the recurrences.

Crossrefs

The gift scenarios sequences when a gift can be stolen at most s times, for s = 1..9, are A001515, A144416, A144508, A144509, A149187, A281358, A281359, A281360, A281361.

Programs

  • Maple
    with(combinat):
    b:= proc(n, i, t) option remember; `if`(t*i add(b(k, 10, n), k=0..10*n):
    seq(a(n), n=0..12);  # Alois P. Heinz, Feb 01 2017
  • Mathematica
    t[n_, n_] = 1; t[n_ /; n >= 0, k_] /; 0 <= k <= 10*n := t[n, k] = Sum[(1/j!)*Product[k - m, {m, 1, j}]*t[n - 1, k - j - 1], {j, 0, 9}]; t[, ] = 0; a[n_] := Sum[t[n, k], {k, 0, 10*n}]; Table[a[n], {n, 0, 11}] (* Jean-François Alcover, Feb 18 2017 *)
  • PARI
    {a(n) = sum(i=n, 10*n, i!*polcoef(sum(j=1, 10, x^j/j!)^n, i))/n!} \\ Seiichi Manyama, May 22 2019
Showing 1-10 of 17 results. Next