login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A080108 a(n) = Sum_{k=1..n} k^(n-k)*binomial(n-1,k-1). 21
1, 2, 6, 23, 104, 537, 3100, 19693, 136064, 1013345, 8076644, 68486013, 614797936, 5818490641, 57846681092, 602259154853, 6548439927680, 74180742421185, 873588590481988, 10674437936521069, 135097459659312176 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Row sums of triangle A154372. Example: a(3)=1+12+9+1=23. From A152818. - Paul Curtz, Jan 08 2009

Number of pointed set partitions of a pointed set k[1...k...n] with a prescribed point k. - Gus Wiseman, Sep 27 2015

With offset 0, a(n) is the number of partial functions (A000169) from [n]->[n] such that every element  in the domain of definition is mapped to a fixed point.  This implies a(n) is the number of idempotent partial functions Cf. A121337. - Geoffrey Critzer, Aug 07 2016

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..200

FORMULA

G.f.: Sum_{k>0} x^k/(1-k*x)^k.

E.g.f. (for offset 0): exp(x*(1+exp(x))). - Vladeta Jovovic, Aug 25 2003

a(n) = A185298(n)/n.

EXAMPLE

G.f. = x + 2*x^2 + 6*x^3 + 23*x^4 + 104*x^5 + 537*x^6 + 3100*x^7 + 19693*x^8 + ...

The a(4) = 23 pointed set partitions of 1[1 2 3 4] are 1[1[1 2 3 4]], 1[1[1] 2[2 3 4]], 1[1[1] 3[2 3 4]], 1[1[1] 4[2 3 4]], 1[1[1 2] 3[3 4]], 1[1[1 2] 4[3 4]], 1[1[1 3] 2[2 4]], 1[1[1 3] 4[2 4]], 1[1[1 4] 2[2 3]], 1[1[1 4] 3[2 3]], 1[1[1 2 3] 4[4]], 1[1[1 2 4] 3[3]], 1[1[1 3 4] 2[2]], 1[1[1] 2[2] 3[3 4]], 1[1[1] 2[2] 4[3 4]], 1[1[1] 2[2 3] 4[4]], 1[1[1] 2[2 4] 3[3]], 1[1[1] 3[3] 4[2 4]], 1[1[1] 3[2 3] 4[4]], 1[1[1 2] 3[3] 4[4]], 1[1[1 3] 2[2] 4[4]], 1[1[1 4] 2[2] 3[3]], 1[1[1] 2[2] 3[3] 4[4]].

MATHEMATICA

Table[Sum[k^(n-k) Binomial[n-1, k-1], {k, n}], {n, 30}] (* Harvey P. Dale, Aug 19 2012 *)

Table[SeriesCoefficient[Sum[x^k/(1-k*x)^k, {k, 0, n}], {x, 0, n}], {n, 1, 20}] (* Vaclav Kotesovec, Aug 06 2014 *)

CoefficientList[Series[E^(x*(1+E^x)), {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, Aug 06 2014 *)

PROG

(PARI) a(n)=sum(k=1, n, k^(n-k)*binomial(n-1, k-1)) \\ Anders Hellström, Sep 27 2015

CROSSREFS

First column of array A098697.

Cf. A185298, A262671.

Sequence in context: A137546 A137547 A137548 * A113227 A200406 A165489

Adjacent sequences:  A080105 A080106 A080107 * A080109 A080110 A080111

KEYWORD

nonn

AUTHOR

Vladeta Jovovic, Mar 15 2003

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 25 06:09 EST 2021. Contains 340416 sequences. (Running on oeis4.)