A032271 Shifts left under "DIJ" (bracelet, indistinct, labeled) transform, with a(1) = 1.

1, 1, 2, 6, 26, 153, 1159, 10787, 118990, 1514889, 21829045, 350831856, 6217711541, 120414106210, 2529280077269, 57262796466087, 1389919284470906, 36003846648705324, 991298773140316816, 28907953199669176721, 890060047425265014457, 28852689090331738589427

Offset

1,3

Links

Andrew Howroyd, Table of n, a(n) for n = 1..200

Mathematica

m = 22; A[_] = 0;
Do[A[x_] = Integrate[1 + (1/2)(A[x] + A[x]^2/2 - Log[1 - A[x]]) + O[x]^m, x] // Normal, {m}];
CoefficientList[A[x], x] Range[0, m]! // Rest (* Jean-François Alcover, Sep 18 2019 *)

Prog

(PARI) DIJ(p)={(p + p^2/2 - log(1-p))/2}
seq(n)={my(p=O(x)); for(i=1, n, p=intformal(1+DIJ(p))); Vec(serlaplace(p))} \\ Andrew Howroyd, Sep 19 2018

Crossrefs

Cf. A032272, A032273.

Sequence in context: A372236 A003659 A159602 * A205502 A107104 A141761

Adjacent sequences: A032268 A032269 A032270 * A032272 A032273 A032274

Keyword

nonn,eigen

Author

Christian G. Bower

Extensions

Terms a(20) and beyond from Andrew Howroyd, Sep 19 2018

Status

approved