cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A128386 Expansion of c(3*x^2)/(1-x*c(3*x^2)), c(x) the g.f. of A000108.

This page as a plain text file.
%I A128386 #28 Aug 06 2024 04:48:12
%S A128386 1,1,4,7,28,58,232,523,2092,4966,19864,48838,195352,492724,1970896,
%T A128386 5068915,20275660,52955950,211823800,560198962,2240795848,5987822380,
%U A128386 23951289520,64563867454,258255469816,701383563388,2805534253552
%N A128386 Expansion of c(3*x^2)/(1-x*c(3*x^2)), c(x) the g.f. of A000108.
%C A128386 Hankel transform is 3^C(n+1,2) = A047656(n+1).
%C A128386 Series reversion of x*(1+x)/(1+2*x+4*x^2).
%H A128386 G. C. Greubel, <a href="/A128386/b128386.txt">Table of n, a(n) for n = 0..1000</a>
%H A128386 Alin Bostan, <a href="https://citeseerx.ist.psu.edu/pdf/749aef4c6f3668e652b5074e5268346ccecc88c9">Computer Algebra for Lattice Path Combinatorics</a>, Slides, Séminaire de Combinatoire Ph. Flajolet, March 28 2013.
%H A128386 Alin Bostan, Andrew Elvey Price, Anthony John Guttmann, and Jean-Marie Maillard, <a href="https://arxiv.org/abs/2001.00393">Stieltjes moment sequences for pattern-avoiding permutations</a>, arXiv:2001.00393 [math.CO], 2020.
%F A128386 G.f.: (sqrt(1-12*x^2)+2*x-1)/(2*x*(1-4*x)).
%F A128386 a(n) = (1/(n+1))*Sum_{k=0..n+1} Sum_{j=0..k} C(n,k)*C(k,j)*C(2*n-2*k+j, n-2*k+j)*(-1)^(n+j)*2^(2*k-j).
%F A128386 a(n) = Sum_{k=0..floor(n/2)} C(n,n-k)*(n-2*k+1)*3^k/(n-k+1);
%F A128386 a(n) = Sum_{k=0..floor(n/2)} A009766(n-k,k)*3^k.
%F A128386 a(n) = Sum_{k=0..n} 3^k*A120730(n,n-k). - _Philippe Deléham_, Mar 03 2007
%F A128386 D-finite with recurrence (n+1)*a(n) - 4*(n+1)*a(n-1) + 12*(2-n)*a(n-2) + 48*(n-2)*a(n-3) = 0. - _R. J. Mathar_, Nov 14 2011
%t A128386 A120730[n_, k_]:= If[n>2*k, 0, Binomial[n,k]*(2*k-n+1)/(k+1)];
%t A128386 A126386[n_]:= Sum[3^k*A120730[n, n-k], {k,0,n}];
%t A128386 Table[A126386[n], {n,0,50}] (* _G. C. Greubel_, Nov 07 2022 *)
%o A128386 (Magma) R<x>:=PowerSeriesRing(Rationals(), 40); Coefficients(R!( (Sqrt(1-12*x^2)+2*x-1)/(2*x*(1-4*x)) )); // _G. C. Greubel_, Nov 07 2022
%o A128386 (SageMath)
%o A128386 def A120730(n, k): return 0 if (n>2*k) else binomial(n, k)*(2*k-n+1)/(k+1)
%o A128386 def A126386(n): return sum(3^k*A120730(n,n-k) for k in range(n+1))
%o A128386 [A126386(n) for n in range(51)] # _G. C. Greubel_, Nov 07 2022
%Y A128386 Cf. A000108, A009766, A047656, A120730.
%K A128386 easy,nonn
%O A128386 0,3
%A A128386 _Paul Barry_, Feb 28 2007