OFFSET
0,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..500
C. G. Bower, Transforms (2)
PROG
(PARI)
DHK(p, n)={my(q=((1+p)^2/(1-subst(p, x, x^2))-1)/2); p + (p^2-subst(p, x, x^2))/2 + sum(d=1, n, moebius(d)*(log(subst(1/(1+O(x*x^(n\d))-p), x, x^d))/d - subst(q + O(x*x^(n\d)), x, x^d)))/2}
seq(n)={Vec(1 + DHK(x+x/(1-x) + O(x*x^n), n))} \\ Andrew Howroyd, Sep 21 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(0)=1 prepended by Andrew Howroyd, Sep 21 2018
STATUS
approved