A097193 - OEIS

%I #15 Sep 08 2022 08:45:14

%S 1,12,204,3978,83538,1837836,41745132,970574319,22970258883,

%T 551286213192,13381219902024,327839887599588,8095123378420596,

%U 201221638263597672,5030540956589941800,126392341534322287725

%N G.f. A(x) satisfies A097191(x*A(x)) = A(x) and so equals the ratio of the g.f.s of any two adjacent diagonals of triangle A097190.

%H Vincenzo Librandi, <a href="/A097193/b097193.txt">Table of n, a(n) for n = 0..200</a>

%F G.f.: A(x) = (1-(1-27*x)^(1/9))/(3*x).

%F G.f.: A(x) = (1/x)*(series reversion of x/A097191(x)).

%F a(n) = A097192(n)/(n+1).

%F a(n) ~ 27^n / (Gamma(8/9) * n^(10/9)). - _Vaclav Kotesovec_, Feb 12 2014

%p seq(coeff(series((1-(1-27*x)^(1/9))/(3*x), x, n+2), x, n), n = 0 ..20); # _G. C. Greubel_, Sep 17 2019

%t CoefficientList[Series[(1-(1-27*x)^(1/9))/(3*x), {x, 0, 20}], x] (* _Vaclav Kotesovec_, Feb 12 2014 *)

%o (PARI) a(n)=polcoeff((1-(1-27*x+x^2*O(x^n))^(1/9))/(3*x),n,x)

%o (Magma) R<x>:=PowerSeriesRing(Rationals(), 20); Coefficients(R!( (1-(1-27*x)^(1/9))/(3*x) )); // _G. C. Greubel_, Sep 17 2019

%o (Sage)

%o def A097193_list(prec):

%o     P.<x> = PowerSeriesRing(QQ, prec)

%o     return P((1-(1-27*x)^(1/9))/(3*x)).list()

%o A097193_list(20) # _G. C. Greubel_, Sep 17 2019

%Y Cf. A097190, A097191, A097192.

%K nonn

%O 0,2

%A _Paul D. Hanna_, Aug 03 2004