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!)
A209319 Number of functions f:{1,2,...,n}->{1,2,...,n} whose cycle lengths are <= 2. 3
1, 1, 4, 25, 218, 2451, 33832, 554527, 10535100, 227790505, 5525843696, 148673435769, 4394818486456, 141611317636075, 4940870266568160, 185595910032346111, 7468517348971708688, 320562141349559055633, 14619577651630443611200, 706025600924216704982425 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

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

FORMULA

E.g.f.: exp(T(x) + T(x)^2/2) = A(T(x)) where A(x) is the e.g.f. for A000085 and T(x) is the e.g.f. for A000169.

a(n) ~ 2*exp(3/2)*n^(n-1). - Vaclav Kotesovec, Sep 30 2013

EXAMPLE

a(3) = 25 because there are 27 functions from {1,2,3} into itself but 2 of these have cycle length of 3: 2,3,1, and 3,1,2.

MAPLE

T:= -LambertW(-x):

egf:= exp(T + T^2/2):

a:= n-> n!*coeff(series(egf, x, n+1), x, n):

seq(a(n), n=0..20);  # Alois P. Heinz, Jan 19 2013

MATHEMATICA

nn=20; t=Sum[n^(n-1)x^n/n!, {n, 1, nn}]; Range[0, nn]!CoefficientList[Series[Exp[Sum[t^i/i, {i, 1, 2}]], {x, 0, nn}], x]

CROSSREFS

Column k=2 of A246522.

Sequence in context: A218826 A060911 A060912 * A273996 A246530 A268163

Adjacent sequences:  A209316 A209317 A209318 * A209320 A209321 A209322

KEYWORD

nonn

AUTHOR

Geoffrey Critzer, Jan 19 2013

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 May 18 22:26 EDT 2021. Contains 344004 sequences. (Running on oeis4.)