|
%I #8 Sep 21 2018 02:20:58
%S 1,2,2,3,4,7,17,39,102,258,651,1620,4028,9956,24588,60742,150108,
%T 371689,921822,2291257,5706077,14240518,35607189,89202758,223857382,
%U 562716179,1416686294,3571809027,9017514253,22794726763
%N "DHK" (bracelet, identity, unlabeled) transform of 2,1,1,1,...
%H Andrew Howroyd, <a href="/A032252/b032252.txt">Table of n, a(n) for n = 0..500</a>
%H C. G. Bower, <a href="/transforms2.html">Transforms (2)</a>
%o (PARI)
%o 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}
%o seq(n)={Vec(1 + DHK(x+x/(1-x) + O(x*x^n), n))} \\ _Andrew Howroyd_, Sep 21 2018
%K nonn
%O 0,2
%A _Christian G. Bower_
%E a(0)=1 prepended by _Andrew Howroyd_, Sep 21 2018
|
