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-8 of 8 results.

A177453 Partial sums of A001863.

Original entry on oeis.org

0, 1, 5, 31, 267, 3027, 42599, 715191, 13942995, 309522515, 7707841015, 212783127799, 6449579387715, 212939326904131, 7606688596589431, 292321288041079671, 12025358303201356019, 527265684696785414387
Offset: 1

Views

Author

Jonathan Vos Post, May 09 2010

Keywords

Comments

Partial sums of normalized total height of rooted trees with n nodes. The subsequence of primes in the partial sums begins: 5, 31, no more through a(15).

Examples

			a(5) = 0 + 1 + 4 + 26 + 236 = 267 = 3 * 89.
		

Crossrefs

Programs

  • Maple
    A001863 := proc(n) if n = 1 then 0; else add( (n-2)!*n^k/k!,k=0..n-2) ; end if; end proc:
    A177453 := proc(n) add(A001863(i),i=0..n) ; end proc: seq(A177453(n),n=1..20) ; # R. J. Mathar, May 28 2010
  • Mathematica
    Accumulate[Table[Sum[(n-2)! n^k/k!,{k,0,n-2}],{n,20}]] (* Harvey P. Dale, Jun 19 2016 *)
  • Python
    from math import comb
    def A177453(n): return sum(((sum(comb(i,k)*(i-k)**(i-k)*k**k for k in range(1,(i+1>>1)))<<1) + (0 if i&1 else comb(i,m:=i>>1)*m**i))//i//(i-1) for i in range(2,n+1)) # Chai Wah Wu, Apr 25-26 2023

Formula

a(n) = Sum_{i=1..n} A001863(i).

Extensions

Extended by R. J. Mathar, May 28 2010

A001865 Number of connected functions on n labeled nodes.

Original entry on oeis.org

1, 3, 17, 142, 1569, 21576, 355081, 6805296, 148869153, 3660215680, 99920609601, 2998836525312, 98139640241473, 3478081490967552, 132705415800984825, 5423640496274200576, 236389784118231290049, 10944997108429625524224, 536484538620663729658993
Offset: 1

Views

Author

Keywords

Comments

If one randomly selects a ball from an urn containing n different balls, with replacement, until exactly one ball has been selected twice, the probability that that ball was also the first ball selected once is a(n)/n^n. See also A000435. - Matthew Vandermast, Jun 15 2004
a(n) equals the permanent of the (n-1) X (n-1) matrix with n+1's along the main diagonal and 1's everywhere else. - John M. Campbell, Apr 20 2012

References

  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 1, p. 112.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nuernberg, Jul 27 1994
  • 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

a(n) = A000435(n) + n^(n-1). See also A063169.
Column k=1 of A060281.

Programs

  • Maple
    spec := [B, {A=Prod(Z,Set(A)), B=Cycle(A)}, labeled]; [seq(combstruct[count](spec,size=n), n=0..20)];
    seq(simplify(GAMMA(n,n)*exp(n)),n=1..20); # Vladeta Jovovic, Jul 21 2005
  • Mathematica
    t=Sum[n^(n-1)x^n/n!,{n,1,20}];
    Range[0,20]! CoefficientList[Series[Log[1/(1-t)]+1,{x,0,20}],x] (* Geoffrey Critzer, Mar 12 2011 *)
    f[n_] := Sum[n! n^(n - k - 1)/(n - k)!, {k, n}]; Array[f, 18] (* Robert G. Wilson v *)
    a[n_] := Exp[n]*Gamma[n, n]; Table[a[n] // FunctionExpand, {n, 1, 18}] (* Jean-François Alcover, May 13 2013, after Vladeta Jovovic *)
  • PARI
    a(n)=if(n<0,0,n!*sum(k=1,n,n^(n-k-1)/(n-k)!))
    
  • PARI
    a(n)=(1/n)*sum(k=1,n,binomial(n,k)*(n-k)^(n-k)*k^k) \\ Paul D. Hanna, Jul 04 2013
    
  • PARI
    N=20; x='x+O('x^N); Vec(serlaplace(log(sum(k=0, N, (k*x)^k/k!)))) \\ Seiichi Manyama, May 27 2019
    
  • Python
    from math import comb
    def A001865(n): return ((sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n))//n + n**(n-1) # Chai Wah Wu, Apr 25-26 2023

Formula

a(n) = Sum_{k=1..n} n!*n^(n-k-1) / (n-k)!.
E.g.f.: -log(1+LambertW(-x)). - Vladeta Jovovic, Apr 11 2001
E.g.f. satisfies 0=2y'^4+2y''^2-y'''y'-y''y'^2. - Michael Somos, Aug 23 2003
Integral representation in terms of the incomplete Gamma function: a(n) = exp(n+1)*Gamma(n+1,n+1) = exp(n+1)*Integral_{x=n+1..oo} x^n exp(-x) dx.
Asymptotics: sqrt(Pi*n/2)*n^(n-1). - N-E. Fahssi, Jan 25 2008, corrected by Vaclav Kotesovec, Nov 27 2012
a(n) = exp(1)*Integral_{x=1..oo} (n+x)^n*exp(-x) dx. - Gerald McGarvey, Apr 16 2008
a(n) = (1/n) * Sum_{k=1..n} C(n,k) * (n-k)^(n-k) * k^k. - Paul D. Hanna, Jul 04 2013
From Peter Bala, Jun 29 2016: (Start)
It appears that a(n) = (n-1)!*( e^n - Sum_{k >= 0} n^(n + k)/(n + k)! ) = (n-1)!*( e^n - Sum_{k >= 0} k^2*n^(n + k - 1)/(n + k)! ).
Note that (n-1)!*( e^n - Sum_{k >= 0} k^3*n^(n + k - 1)/(n + k)! ) also appears to be an integer sequence beginning [1, 5, 37, 370, 4681, 71736, 1292005, ...]. (End)
a(n) = Sum_{k=1..n} (n!/(n-k)!) * k^2 * n^(n-k-2). - Brian P Hawkins, Feb 07 2024

Extensions

More terms from James Sellers, May 23 2000

A000435 Normalized total height of all nodes in all rooted trees with n labeled nodes.

Original entry on oeis.org

0, 1, 8, 78, 944, 13800, 237432, 4708144, 105822432, 2660215680, 73983185000, 2255828154624, 74841555118992, 2684366717713408, 103512489775594200, 4270718991667353600, 187728592242564421568, 8759085548690928992256, 432357188322752488126152, 22510748754252398927872000
Offset: 1

Views

Author

Keywords

Comments

This is the sequence that started it all: the first sequence in the database!
The height h(V) of a node V in a rooted tree is its distance from the root. a(n) = Sum_{all nodes V in all n^(n-1) rooted trees on n nodes} h(V)/n.
In the trees which have [0, n-1] = (0, 1, ..., n-1) as their ordered set of nodes, the number of nodes at distance i from node 0 is f(n,i) = (n-1)...(n-i)(i+1)n^(j-1), 0 <= i < n-1, i+j = n-1 (and f(n,n-1) = (n-1)!): (n-1)...(n-i) counts the words coding the paths of length i from any node to 0, n^(j-1) counts the Pruefer codes of the rest, words build by iterated deletion of the greater node of degree 1 ... except the last one, (i+1), necessary pointing at the path. If g(n,i) = (n-1)...(n-i)n^j, i+j = n-1, f(n,i) = g(n,i) - g(n,i+1), g(n,i) = Sum_{k>=i} f(n,k), the sequence is Sum_{i=1..n-1} g(n,i). - Claude Lenormand (claude.lenormand(AT)free.fr), Jan 26 2001
If one randomly selects one ball from an urn containing n different balls, with replacement, until exactly one ball has been selected twice, the probability that this ball was also the second ball to be selected once is a(n)/n^n. See also A001865. - Matthew Vandermast, Jun 15 2004
a(n) is the number of connected endofunctions with no fixed points. - Geoffrey Critzer, Dec 13 2011
a(n) is the number of weakly connected simple digraphs on n labeled nodes where every node has out-degree 1. A digraph where all out-degrees are 1 can be called a functional digraph due to the correspondence with endofunctions. - Andrew Howroyd, Feb 06 2024

Examples

			For n = 3 there are 3^2 = 9 rooted labeled trees on 3 nodes, namely (with o denoting a node, O the root node):
   o
   |
   o     o   o
   |      \ /
   O       O
The first can be labeled in 6 ways and contains nodes at heights 1 and 2 above the root, so contributes 6*(1+2) = 18 to the total; the second can be labeled in 3 ways and contains 2 nodes at height 1 above the root, so contributes 3*2=6 to the total, giving 24 in all. Dividing by 3 we get a(3) = 24/3 = 8.
For n = 4 there are 4^3 = 64 rooted labeled trees on 4 nodes, namely (with o denoting a node, O the root node):
   o
   |
   o     o        o   o
   |     |         \ /
   o     o   o      o     o o o
   |      \ /       |      \|/
   O       O        O       O
  (1)     (2)      (3)     (4)
Tree (1) can be labeled in 24 ways and contains nodes at heights 1, 2, 3 above the root, so contributes 24*(1+2+3) = 144 to the total;
tree (2) can be labeled in 24 ways and contains nodes at heights 1, 1, 2 above the root, so contributes 24*(1+1+2) = 96 to the total;
tree (3) can be labeled in 12 ways and contains nodes at heights 1, 2, 2 above the root, so contributes 12*(1+2+2) = 60 to the total;
tree (4) can be labeled in 4 ways and contains nodes at heights 1, 1, 1 above the root, so contributes 4*(1+1+1) = 12 to the total;
giving 312 in all. Dividing by 4 we get a(4) = 312/4 = 78.
		

References

  • 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

Cf. A001863, A001864, A001854, A002862 (unlabeled version), A234953, A259334.
Column k=1 of A350452.

Programs

  • Maple
    A000435 := n-> (n-1)!*add (n^k/k!, k=0..n-2);
    seq(simplify((n-1)*GAMMA(n-1,n)*exp(n)),n=1..20); # Vladeta Jovovic, Jul 21 2005
  • Mathematica
    f[n_] := (n - 1)! Sum [n^k/k!, {k, 0, n - 2}]; Array[f, 18] (* Robert G. Wilson v, Aug 10 2010 *)
    nx = 18; Rest[ Range[0, nx]! CoefficientList[ Series[ LambertW[-x] - Log[1 + LambertW[-x]], {x, 0, nx}], x]] (* Robert G. Wilson v, Apr 13 2013 *)
  • PARI
    x='x+O('x^30); concat(0, Vec(serlaplace(lambertw(-x)-log(1+lambertw(-x))))) \\ Altug Alkan, Sep 05 2018
    
  • PARI
    A000435(n)=(n-1)*A001863(n) \\ M. F. Hasler, Dec 10 2018
    
  • Python
    from math import comb
    def A000435(n): return ((sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n))//n # Chai Wah Wu, Apr 25-26 2023

Formula

a(n) = (n-1)! * Sum_{k=0..n-2} n^k/k!.
a(n) = A001864(n)/n.
E.g.f.: LambertW(-x) - log(1+LambertW(-x)). - Vladeta Jovovic, Apr 10 2001
a(n) = A001865(n) - n^(n-1).
a(n) = A001865(n) - A000169(n). - Geoffrey Critzer, Dec 13 2011
a(n) ~ sqrt(Pi/2)*n^(n-1/2). - Vaclav Kotesovec, Aug 07 2013
a(n)/A001854(n) ~ 1/2 [See Renyi-Szekeres, (4.7)]. Also a(n) = Sum_{k=0..n-1} k*A259334(n,k). - David desJardins, Jan 20 2017
a(n) = (n-1)*A001863(n). - M. F. Hasler, Dec 10 2018

Extensions

Additional references from Valery A. Liskovets
Editorial changes by N. J. A. Sloane, Feb 03 2012
Edited by M. F. Hasler, Dec 10 2018

A001864 Total height of rooted trees with n labeled nodes.

Original entry on oeis.org

0, 2, 24, 312, 4720, 82800, 1662024, 37665152, 952401888, 26602156800, 813815035000, 27069937855488, 972940216546896, 37581134047987712, 1552687346633913000, 68331503866677657600, 3191386068123595166656, 157663539876436721860608
Offset: 1

Views

Author

Keywords

Comments

a(n) is the total number of nonrecurrent elements mapped into a recurrent element in all functions f:{1,2,...,n}->{1,2,...,n}. a(n) = Sum_{k=1..n-1} A216971(n,k)*k. - Geoffrey Critzer, Jan 01 2013
a(n) is the sum of the lengths of all cycles over all functions f:{1,2,...,n}->{1,2,...,n}. Fixed points are taken to have length zero. a(n) = Sum_{k=2..n} A066324(n,k)*(k-1). - Geoffrey Critzer, Aug 19 2013

References

  • 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

Programs

  • Maple
    A001864 := proc(n) local k; add(n!*n^k/k!, k=0..n-2); end;
  • Mathematica
    Table[Sum[Binomial[n,k](n-k)^(n-k) k^k,{k,1,n-1}],{n,20}] (* Harvey P. Dale, Oct 10 2011 *)
    a[n_] := n*(n-1)*Exp[n]*Gamma[n-1, n] // Round; Table[a[n], {n, 1, 18}]  (* Jean-François Alcover, Jun 24 2013 *)
  • PARI
    a(n)=sum(k=1,n-1,binomial(n,k)*(n-k)^(n-k)*k^k)
    
  • Python
    from math import comb
    def A001864(n): return (sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n) # Chai Wah Wu, Apr 25-26 2023

Formula

a(n) = n*A000435(n).
E.g.f: (LambertW(-x)/(1+LambertW(-x)))^2. - Vladeta Jovovic, Apr 10 2001
a(n) = Sum_{k=1..n-1} binomial(n, k)*(n-k)^(n-k)*k^k. - Benoit Cloitre, Mar 22 2003
a(n) ~ sqrt(Pi/2)*n^(n+1/2). - Vaclav Kotesovec, Aug 07 2013
a(n) = n! * Sum_{k=0..n-2} n^k/k!. - Jianing Song, Aug 08 2022

A001866 Number of connected graphs with n nodes and n edges.

Original entry on oeis.org

0, 0, 1, 24, 936, 56640, 4968000, 598328640, 94916183040, 19200422062080, 4826695329792000, 1476585999504000000, 540272647694971699200, 233019960215154829516800, 117009251702203840384204800, 67680314823703303654732800000, 44677678066673631080900198400000
Offset: 0

Views

Author

Keywords

Comments

Or number of n X n (0,1) matrices with two 1's in each row the permanent of which equals to 2. Note that, if (0,1) matrix with two 1's in each row has positive permanent, then it is equal to a power of 2. - Vladimir Shevelev, Mar 25 2010

References

  • V. S. Shevelev, On the permanent of the stochastic (0,1)-matrices with equal row sums, Izvestia Vuzov of the North-Caucasus region, Nature sciences 1 (1997), 21-38 (in Russian). - Vladimir Shevelev, Mar 25 2010
  • 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

Cf. A174586.

Programs

  • Mathematica
    Join[{0}, Table[(n!^2*n^(n - 1)/2)*Sum[n^(-k)/(n - k)!, {k, 2, n}], {n, 20}]] (* T. D. Noe, Aug 10 2012 *)

Formula

Explicit formula: a(n) = (n!^2*n^(n-1)/2)*Sum_{k=2..n} n^(-k)/(n-k)!; Recursion: a(2)=1, for n>=3, a(n) = n!*((n-1)!/2+Sum_{k=2..n-1} (-1)^(n+k+1)*k^(n-k)*binomial(n,k)*a(k)/k!). - Vladimir Shevelev, Mar 25 2010
a(n) ~ Pi * n^(2*n) / (2*exp(n)). - Vaclav Kotesovec, Nov 30 2017

A320064 The number of F_2 graphs on { 1, 2, ..., n } with a unique cycle of weight 1, which corresponds to the number of reflectable bases of the root system of type D_n.

Original entry on oeis.org

0, 1, 16, 312, 7552, 220800, 7597824, 301321216, 13545271296, 681015214080, 37879390720000, 2309968030334976, 153275504883695616, 10995166075754119168, 847974316241667686400, 69971459959477921382400, 6151490510604350965940224, 574035430519008722436489216, 56669921387839814123670994944
Offset: 1

Views

Author

Masaya Tomie, Oct 04 2018

Keywords

Crossrefs

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (&+[(&+[j^(j-1)*(2*x)^j/Factorial(j): j in [1..m+2]])^k/(4*k): k in [2..m+1]]) )); [0] cat [Factorial(n+1)*b[n]: n in [1..m-2]]; // G. C. Greubel, Dec 10 2018
    
  • Mathematica
    nmax = 20; Rest[CoefficientList[Series[Sum[1/(4*m)*(Sum[k^(k-1)*(2*x)^k/k!, {k, 1, nmax}])^m, {m, 2, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!] (* Vaclav Kotesovec, Oct 23 2018 *)
  • PARI
    seq(n)={Vec(serlaplace(sum(m=2, n, (sum(k=1, n, k^(k-1)*(2*x)^k/k!) + O(x^n))^m/(4*m))), -n)} \\ Andrew Howroyd, Nov 07 2018
    
  • PARI
    apply( A320064(n)=A001863(n)*(n-1)<<(n-2), [1..20]) \\ Defines the function A320064. The additional apply(...) provides a check and illustration. - M. F. Hasler, Dec 09 2018
    
  • Python
    from math import comb
    def A320064(n): return 0 if n<2 else ((sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n))//n<Chai Wah Wu, Apr 25-26 2023

Formula

E.g.f.: Sum_{m>=2} (1/(4*m)) (Sum_{k>=1} k^(k-1)*(2*x)^k/k!)^m.
a(n) = (n-1)*2^(n-2)*A001863(n). - M. F. Hasler, Dec 09 2018
a(n) = 2^(n-2)*A000435(n). - Chai Wah Wu, Apr 25 2023

A321233 a(n) is the number of reflectable bases of the root system of type D_n.

Original entry on oeis.org

0, 4, 128, 4992, 241664, 14131200, 972521472, 77138231296, 6935178903552, 697359579217920, 77576992194560000, 9461629052252061696, 1255632936007234486272, 180144800985155488448512, 27786422394606966747955200, 4585649599904345055716966400, 806288164205933489807717040128
Offset: 1

Views

Author

Masaya Tomie, Nov 01 2018

Keywords

Comments

The root systems of type D_n are only defined for n >= 4. See chapter 3 of the Humphreys reference. Sequence extended to n=1 using formula/recurrence.

References

  • J. E. Humphreys, Introduction to Lie algebras and representation theory, 2nd ed, Springer-Verlag, New York, 1972.

Crossrefs

Programs

  • Magma
    m:=25; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (&+[ (&+[ j^(j-1)*(4*x)^j/Factorial(j) :j in [1..m+3]])^k/(4*k) :k in [2..m+2]]) )); [0] cat [Factorial(n+1)*b[n]: n in [1..m-2]]; // G. C. Greubel, Dec 09 2018
    
  • Mathematica
    Rest[With[{m = 25}, CoefficientList[Series[Sum[Sum[j^(j - 1)*(4*x)^j/j!, {j, 1, m + 1}]^k/(4*k), {k, 2, m}], {x, 0, m}], x]*Range[0, m]!]] (* G. C. Greubel, Dec 09 2018 *)
  • PARI
    a(n)={n!*polcoef(sum(m=2, n, (sum(k=1, n, k^(k-1)*(4*x)^k/k!) + O(x^(n-m+2)))^m/(4*m)), n)} \\ Andrew Howroyd, Nov 01 2018
    
  • PARI
    A321233(n)=A001863(n)*(n-1)*4^(n-1) \\ M. F. Hasler, Dec 09 2018
    
  • Python
    from math import comb
    def A321233(n): return 0 if n<2 else ((sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n))//n<<(n-1<<1) # Chai Wah Wu, Apr 26 2023

Formula

E.g.f.: Sum_{m>=2} (1/(4*m)) (Sum_{k>=1} k^(k-1)*(4*x)^k/k!)^m.
a(n) = 2^n*A320064(n).
a(n) = (n-1)*4^(n-1)*A001863(n). - M. F. Hasler, Dec 09 2018

A134558 Array read by antidiagonals, a(n,k) = gamma(n+1,k)*e^k, where gamma(n,k) is the upper incomplete gamma function and e is the exponential constant 2.71828...

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 6, 5, 3, 1, 24, 16, 10, 4, 1, 120, 65, 38, 17, 5, 1, 720, 326, 168, 78, 26, 6, 1, 5040, 1957, 872, 393, 142, 37, 7, 1, 40320, 13700, 5296, 2208, 824, 236, 50, 8, 1, 362880, 109601, 37200, 13977, 5144, 1569, 366, 65, 9, 1, 3628800, 986410, 297856
Offset: 0

Views

Author

Ross La Haye, Jan 22 2008

Keywords

Examples

			Square array begins:
    1,    1,    1,     1,     1,     1,      1, ...
    1,    2,    3,     4,     5,     6,      7, ...
    2,    5,   10,    17,    26,    37,     50, ...
    6,   16,   38,    78,   142,   236,    366, ...
   24,   65,  168,   393,   824,  1569,   2760, ...
  120,  326,  872,  2208,  5144, 10970,  21576, ...
  720, 1957, 5296, 13977, 34960, 81445, 176112, ...
		

Crossrefs

Cf. a(n, 0) = A000142(n); a(n, 1) = A000522(n); a(n, 2) = A010842(n); a(n, 3) = A053486(n); a(n, 4) = A053487(n); a(n, 5) = A080954(n); a(n, 6) = A108869(n); a(1, k) = A000027(k+1); a(2, k) = A002522(k+1); a(n, n) = A063170(n); a(n, n+1) = A001865(n+1); a(n, n+2) = A001863(n+2).
Another version: A089258.
A transposed version: A080955.
Cf. A001113.

Programs

  • Mathematica
    T[n_,k_] := Gamma[n+1, k]*E^k; Table[T[n-k, k], {n, 0, 10}, {k, 0, n}] //Flatten (* Amiram Eldar, Jun 27 2020 *)

Formula

a(n,k) = gamma(n+1,k)*e^k = Sum_{m=0..n} m!*binomial(n,m)*k^(n-m).
a(n,k) = n*a(n-1,k) + k^n for n,k > 0.
E.g.f. (by columns) is e^(kx)/(1-x).
a(n,k) = the binomial transform by columns of a(n,k-1).
Conjecture: a(n,k) is the permanent of the n X n matrix with k+1 on the main diagonal and 1 elsewhere.

Extensions

More terms from Amiram Eldar, Jun 27 2020
Showing 1-8 of 8 results.