OFFSET
1,4
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
C. G. Bower, Transforms (2)
FORMULA
"CGK" (necklace, elements, unlabeled) transform of A032065 (shifted right one place).
PROG
(PARI)
BGK(v)={apply(p->subst(serlaplace(y^0*p + polcoeff(p, 1)), y, 1)/2, Vec(prod(k=1, #v, (1 + x^k*y + O(x*x^#v))^v[k])-1, -#v))}
CGK(v)={apply(p->subst(serlaplace(p/y), y, 1), Vec(prod(k=1, #v, (1 + x^k*y + O(x*x^#v))^v[k])-1, -#v))}
seq(n)={my(v=[1]); for(i=3, n, v=concat([1], BGK(v))); concat([1], CGK(v))} \\ Andrew Howroyd, Sep 20 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved