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

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A101313 Number of painted forests - exactly one of its trees is painted - on labeled vertex set [n]. 4
1, 3, 12, 68, 525, 5262, 65674, 987408, 17426565, 353759300, 8127640224, 208600774032, 5917247520457, 183872561612040, 6212370268252950, 226762373954676608, 8893485959056048521, 372980176625914811568, 16656844860594186642100, 789196576594282265505600 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..100

FORMULA

a(n) = f(n) + SUM{((n-i)^(n-i-2))*C((n-1), i)*a(i):i=1, 2, ..(n-1)}, where f(n)=number of forests on labeled vertex set [n], A001858.

Exponential convolution of A000272 and A001858: a(n) = Sum_{k=1..n} binomial(n, k)*k^(k-2)*A001858(n-k). E.g.f.: B(x)*exp(B(x)), where B(x) is e.g.f. for A000272. - Vladeta Jovovic, May 24 2005

a(n) = Sum_{m=1..n} A105599(n,m)*m. -  Geoffrey Critzer, Nov 04 2012

EXAMPLE

a(5) = 291 + (16*4*1)+(3*6*3)+(1*4*12)+(1*1*68) = 525.

MAPLE

B:= n-> exp(add(k^(k-2) *x^k/k!, k = 1..n )): b:= n-> coeff(series(B(n), x, n+1), x, n)*n!: a:= n-> add(binomial(n, k) *k^(k-2) *b(n-k), k=1..n): seq(a(n), n=1..25);  # Alois P. Heinz, Sep 10 2008

MATHEMATICA

nn=20; t=Sum[n^(n-1)x^n/n!, {n, 1, nn}]; Drop[Range[0, nn]!CoefficientList[Series[D[Exp[y(t-t^2/2)], y]/.y->1, {x, 0, nn}], x], 1]  (* Geoffrey Critzer, Nov 04 2012 *)

CROSSREFS

Sequence in context: A296979 A004127 A058115 * A257605 A265886 A295762

Adjacent sequences:  A101310 A101311 A101312 * A101314 A101315 A101316

KEYWORD

nonn

AUTHOR

Joseph G. Moser (jmoser(AT)wcupa.edu), Jan 26 2005

EXTENSIONS

More terms from Alois P. Heinz, Sep 10 2008

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 20 14:02 EST 2020. Contains 331094 sequences. (Running on oeis4.)