A110083 a(n+1) = Sum_{k=0..n} (n!/k!)*binomial(n,k)*a(k).

1, 1, 2, 8, 50, 442, 5212, 78664, 1472756, 33378740, 898227944, 28253387104, 1025373023848, 42467845178632, 1988513519453360, 104413376937507488, 6104596110052561808, 394921638012548722576, 28112685278602155590944, 2191142414957886078590080

Offset

0,3

Links

Alois P. Heinz, Table of n, a(n) for n = 0..307

Maple

a:= proc(n) option remember; `if`(n=0, 1, add(
      a(n-i)*binomial(n-1, i-1)^2*(i-1)!, i=1..n))
    end:
seq(a(n), n=0..20);  # Alois P. Heinz, Aug 13 2019

Mathematica

nmax=20; b = ConstantArray[0, nmax+2]; b[[1]]=1; Do[b[[n+2]] = Sum[n!/k!*Binomial[n, k]*b[[k+1]], {k, 0, n}], {n, 0, nmax}]; b (* Vaclav Kotesovec, Mar 02 2014 *)

Crossrefs

Cf. A001063.

Sequence in context: A050398 A135081 A296366 * A086922 A007128 A013085

Adjacent sequences: A110080 A110081 A110082 * A110084 A110085 A110086

Keyword

easy,nonn

Author

Vladeta Jovovic, Sep 04 2005

Status

approved