login
Expansion of e.g.f.: exp(x * (1 - x)).
15

%I #22 Jul 12 2024 21:59:55

%S 1,1,-1,-5,1,41,31,-461,-895,6481,22591,-107029,-604031,1964665,

%T 17669471,-37341149,-567425279,627491489,19919950975,-2669742629,

%U -759627879679,-652838174519,31251532771999,59976412450835,-1377594095061119,-4256461892701199

%N Expansion of e.g.f.: exp(x * (1 - x)).

%H Seiichi Manyama, <a href="/A293604/b293604.txt">Table of n, a(n) for n = 0..732</a>

%F a(n) = (-1)^n * A000321(n).

%F a(n) = a(n-1) - 2 * (n-1) * a(n-2) for n > 1.

%F E.g.f.: Product_{k>=1} (1 + x^k)^(mu(k)/k). - _Ilya Gutkovskiy_, May 23 2019

%F a(n) = Hermite(n, 1/2). - _G. C. Greubel_, Jul 12 2024

%t CoefficientList[Series[E^(x*(1-x)), {x, 0, 20}], x] * Range[0, 20]! (* _Vaclav Kotesovec_, Oct 13 2017 *)

%o (PARI) my(N=66, x='x+O('x^N)); Vec(serlaplace(exp(x*(1-x))))

%o (PARI) a(n) = polhermite(n, 1/2); \\ _Michel Marcus_, Oct 13 2017

%o (Magma)

%o R<x>:=PowerSeriesRing(Rationals(), 30);

%o Coefficients(R!(Laplace( Exp(x-x^2) ))); // _G. C. Greubel_, Jul 12 2024

%o (SageMath)

%o [hermite(n, 1/2) for n in range(31)] # _G. C. Greubel_, Jul 12 2024

%Y Cf. A000321, A111884, A293571, A293572, A293573.

%Y Sequences with e.g.f = exp(x + q*x^2): A158968 (q=-9), A158954 (q=-4), A362177 (q=-3), A362176 (q=-2), this sequence (q=-1), A000012 (q=0), A047974 (q=1), A115329 (q=2), A293720 (q=4).

%K sign

%O 0,4

%A _Seiichi Manyama_, Oct 12 2017