|
%I #14 Sep 11 2019 06:26:05
%S 1,3,12,57,471,6159,105987,2207847,53715891,1493742279,46730999667,
%T 1624397276487,62111517207411,2590839950614599,117076879593900147,
%U 5697526788322329927,297074131510779228531,16522357697381918530119,976355872514689231565427
%N "DIJ" (bracelet, indistinct, labeled) transform of 3,3,3,3,...
%H Andrew Howroyd, <a href="/A032268/b032268.txt">Table of n, a(n) for n = 0..200</a>
%H C. G. Bower, <a href="/transforms2.html">Transforms (2)</a>
%F E.g.f.: 1 + (g(x) + g(x)^2/2 - log(1-g(x)))/2 where g(x) = 3*(exp(x) - 1). - _Andrew Howroyd_, Sep 12 2018
%t m = 19; (1 + (g[x] + g[x]^2/2 - Log[1 - g[x]])/2 /. g[x] -> 3*(Exp[x] - 1) + O[x]^m // CoefficientList[#, x]&)*Range[0, m-1]! (* _Jean-François Alcover_, Sep 11 2019 *)
%o (PARI) seq(n)={my(p=3*(exp(x + O(x*x^n))-1)); Vec(1 + serlaplace(p + p^2/2 - log(1-p))/2)} \\ _Andrew Howroyd_, Sep 12 2018
%Y Cf. A032266, A032267, A032269, A032270.
%K nonn
%O 0,2
%A _Christian G. Bower_
%E a(0)=1 prepended and terms a(17) and beyond from _Andrew Howroyd_, Sep 12 2018
|
