OFFSET
0,3
COMMENTS
For n > 0, a(n) = number of endofunctions of [n] mapping some x<>1 to 1. - Len Smiley, Nov 15 2001 (Endofunction interpretation from a(n) = n*(n^(n-1) - (n-1)^(n-1)).)
LINKS
Harry J. Smith, Table of n, a(n) for n = 0..100
D. Callan, A Bijection between Marked Trees
Leonard Smiley, Problem 10781, Amer. Math. Monthly, 107, Feb. 2000, p. 176.
FORMULA
a(n) = n*A055869(n-1).
Limit_{n -> oo} ( a(n)/a(n-1) - a(n-1)/a(n-2) ) -> e.
E.g.f.: (1-x)/(1-T), where T=T(x) is Euler's tree function (see A000169). The e.g.f. for n > 0 terms only (applicable to endofunctions) is (T - x)/(1 - T). - Len Smiley, Dec 10 2001
MAPLE
f := n-> n*sum(binomial(n-1, j-1)*(n-1)^(n-j), j=2..n);
g := n-> n^n -n*(n-1)^(n-1);
h := n-> sum(binomial(n, j)*j^(j-1)*(n-j)^(n-j), j=2..n);
k := n-> sum(binomial(n, j-1)*(j-1)^(j-1)*(n-j)^(n-j), j=2..n); # then a(n)=f(n)=g(n)=h(n)=k(n)
MATHEMATICA
Join[{1, 0}, Table[n^n-n*(n-1)^(n-1), {n, 2, 20}]] (* Harvey P. Dale, Nov 16 2012 *)
PROG
(PARI) { for (n=0, 100, write("b060226.txt", n, " ", n^n - n*(n - 1)^(n - 1)); ) } \\ Harry J. Smith, Jul 03 2009
(Haskell)
a060226 0 = 1
a060226 n = a000312 n - n * a000312 (n - 1)
-- Reinhard Zumkeller, Aug 27 2012
(Magma)
A060226:= func< n | n^n - n*(n-1)^(n-1) >;
[A060226(n): n in [0..30]]; // G. C. Greubel, Nov 03 2024
(SageMath)
def A060226(n): return n^n - n*(n-1)^(n-1)
[A060226(n) for n in range(31)] # G. C. Greubel, Nov 03 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Henry Bottomley, Jul 12 2001
STATUS
approved