|
%I #7 Sep 08 2022 08:45:14
%S 1,15,232,3627,57016,899298,14216560,225110307,3568890328,56635884470,
%T 899474459280,14294357356110,227286593929136,3615608476770340,
%U 57538659207907552,915981394162628387,14586262906867731096
%N Row sums of triangle A097181, in which the n-th row polynomial R_n(y) is formed from the initial (n+1) terms of g.f. A097182(y)^(n+1), where R_n(1/2) = 8^n for all n>=0.
%H G. C. Greubel, <a href="/A097185/b097185.txt">Table of n, a(n) for n = 0..825</a>
%F G.f.: A(x) = 2/((1-16*x) + (1-16*x)^(7/8)).
%p seq(coeff(series(2/((1-16*x) + (1-16*x)^(7/8)), x, n+1), x, n), n = 0 ..30); # _G. C. Greubel_, Sep 17 2019
%t CoefficientList[Series[2/((1-16*x) +(1-16*x)^(7/8)), {x,0,30}], x] (* _G. C. Greubel_, Sep 17 2019 *)
%o (PARI) a(n)=polcoeff(2/((1-16*x)+(1-16*x+x*O(x^n))^(7/8)),n,x)
%o (Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( 2/((1-16*x) + (1-16*x)^(7/8)) )); // _G. C. Greubel_, Sep 17 2019
%o (Sage)
%o def A097185_list(prec):
%o P.<x> = PowerSeriesRing(QQ, prec)
%o return P(2/((1-16*x) + (1-16*x)^(7/8))).list()
%o A097185_list(30) # _G. C. Greubel_, Sep 17 2019
%Y Cf. A097181, A097182.
%K nonn
%O 0,2
%A _Paul D. Hanna_, Aug 03 2004
|
