|
%I #14 Sep 18 2019 11:10:13
%S 1,1,2,6,26,153,1159,10787,118990,1514889,21829045,350831856,
%T 6217711541,120414106210,2529280077269,57262796466087,
%U 1389919284470906,36003846648705324,991298773140316816,28907953199669176721,890060047425265014457,28852689090331738589427
%N Shifts left under "DIJ" (bracelet, indistinct, labeled) transform, with a(1) = 1.
%H Andrew Howroyd, <a href="/A032271/b032271.txt">Table of n, a(n) for n = 1..200</a>
%t m = 22; A[_] = 0;
%t Do[A[x_] = Integrate[1 + (1/2)(A[x] + A[x]^2/2 - Log[1 - A[x]]) + O[x]^m, x] // Normal, {m}];
%t CoefficientList[A[x], x] Range[0, m]! // Rest (* _Jean-François Alcover_, Sep 18 2019 *)
%o (PARI) DIJ(p)={(p + p^2/2 - log(1-p))/2}
%o seq(n)={my(p=O(x)); for(i=1, n, p=intformal(1+DIJ(p))); Vec(serlaplace(p))} \\ _Andrew Howroyd_, Sep 19 2018
%Y Cf. A032272, A032273.
%K nonn,eigen
%O 1,3
%A _Christian G. Bower_
%E Terms a(20) and beyond from _Andrew Howroyd_, Sep 19 2018
|
