The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A331538 Number of functions f:[n]->[n] such that there exists a k such that |f^(-1)(k)| = k. 1

%I #19 Feb 12 2021 12:06:00

%S 0,1,3,16,147,1756,25910,453594,9184091,211075288,5427652794,

%T 154380255250,4812088559014,163110595450466,5973198636395003,

%U 235010723141883563,9886231689434154971,442799642855527526848,21038043034795035118742,1056802542597653892224802,55962024535834950971809754

%N Number of functions f:[n]->[n] such that there exists a k such that |f^(-1)(k)| = k.

%D P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009.

%H Darij Grinberg, Marko Riedel, Markus Scheuer, et al., Math.StackExchange, <a href="https://math.stackexchange.com/questions/3509979/">Number of functions f:[n]->[n] such that there exists an i such that |f^(-1)(i)| = i.</a>

%F a(n) = n^n - n! * [z^n] Product_{k=1..n} (exp(z) - z^k/k!).

%F a(n) = n^n - n! * [z^n] Product_{k=1..n} (Sum_{q=0..k-1} z^q/q! + Sum_{q=k+1..n} z^q/q!).

%F a(n) = n^n - A331537(n).

%e For n = 0: a(0) = 0^0 - 0! [z^0] 1 = 0.

%e Functions from [2]->[2] are

%e * [1,1] - pre-images are [1,2] and [], no contribution

%e * [1,2] - pre-images are [1] and [2], pre-image of one has one element, one contribution

%e * [2,1] - pre-images are [2] and [1], pre-image of one has one element, one contribution

%e + [2,2] - pre-images are [] and [1,2], pre-image of two has two elements, one contribution

%e = total contributions is three.

%o (PARI) a(n)={n^n - n!*polcoef(prod(k=1, n, exp(x + O(x*x^n)) - x^k/k!), n)} \\ _Andrew Howroyd_, Jan 19 2020

%Y Cf. A000312, A331537.

%K nonn

%O 0,3

%A _Marko Riedel_, Jan 19 2020

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 26 05:37 EDT 2024. Contains 372807 sequences. (Running on oeis4.)