OFFSET
0,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..1000
EXAMPLE
a(2)=5 because the following compounds are possible:
| | | | | | | | |
-O-O-C-C- -O-C-O-C- -O-C-C-O- -C-O-O-C- -O-C-O-
| | | | | | | | |
-C-
|
PROG
(PARI) \\ here R(n) gives g.f. of A000598.
R(n)={my(g=O(x)); for(n=0, n, g = 1 + x*(g^3/6 + subst(g, x, x^2)*g/2 + subst(g, x, x^3)/3) + O(x*x^n)); g}
seq(n)={my(A=O(x*x^n), p=R(n), p2=subst(p, x, x^2) + A, q=(p^2+p2)/2, q2=subst(q, x, x^2) + A); Vec((p^2/(1 - x^2*q^2) + p2/(1 - x^2*q2))*(1 + x*q)/2)} \\ Andrew Howroyd, Feb 17 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Erich Friedman, Feb 17 2025
EXTENSIONS
a(10) onwards from Andrew Howroyd, Feb 17 2025
STATUS
approved
