OFFSET
0,5
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..1000
Andrey Zabolotskiy, Python script for the sequence
FORMULA
G.f.: Sum_{k>=0} x^(2*k)*(Q(k,x)^2 + Q(k,x^2))*(1 + x*P(k,x))/2, where P(x,k) = 1/Product_{i=1..k} (1-x^i) and Q(x,k) = 1/Product_{i=1..k-1} (1-x^i)^(k-i). - Andrew Howroyd, Feb 06 2025
EXAMPLE
From Andrey Zabolotskiy, Nov 21 2016: (Start)
Three trees that are counted in A000055(10) but not in a(10):
(1)
o o-o-o
| |
o----o
| |
o o-o-o
(2)
o-o-o
|
o-o-o-o
|
o-o-o
(3)
o-o-o-o-o-o-o
|
o-o-o
(End)
PROG
(PARI) seq(n)={my(s=1+x, p=1+O(x^n), p2=p, q=p, q2=p); for(k=1, n\2, q*=p^2; q2*=p2; p /= 1-x^k; p2 /= 1-x^(2*k); s+=x^(2*k)*(q+q2)*(1+x*p)/2); Vec(s+O(x*x^n))} \\ Andrew Howroyd, Feb 06 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Gabriel Burns, Oct 31 2016
EXTENSIONS
Corrections and more terms from Andrey Zabolotskiy, Nov 21 2016
a(24) onwards from Andrew Howroyd, Feb 06 2025
STATUS
approved
