OFFSET
0,2
COMMENTS
The graphs considered here can have loops but not parallel edges.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..500
FORMULA
PROG
(PARI) \\ TreeGf gives gf of A000081.
TreeGf(N)={my(A=vector(N, j, 1)); for (n=1, N-1, A[n+1] = 1/n * sum(k=1, n, sumdiv(k, d, d*A[d]) * A[n-k+1] ) ); x*Ser(A)}
seq(n)={my(t=TreeGf(n)); my(g(e)=subst(t + O(x*x^(n\e)), x, x^e) + O(x*x^n)); Vec(1 + (sum(d=1, n, eulerphi(d)/d*log(1/(1-g(d)))) + ((1+g(1))^2/(1-g(2))-1)/2 + 2*g(1) - 2*g(1)^2 )/2) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Feb 02 2024
STATUS
approved