A030267
Compose the natural numbers with themselves, A(x) = B(B(x)) where B(x) = x/(1-x)^2 is the generating function for natural numbers.
Original entry on oeis.org
1, 4, 14, 46, 145, 444, 1331, 3926, 11434, 32960, 94211, 267384, 754309, 2116936, 5914310, 16458034, 45638101, 126159156, 347769719, 956238170, 2623278946, 7181512964, 19622668679, 53522804976, 145753273225, 396323283724, 1076167858046, 2918447861686
Offset: 1
From _Petros Hadjicostas_, Jun 24 2019: (Start)
Recall that with m-color compositions, a part of size m may be colored with one of m colors.
We have a(1) = 1 because we only have one colored composition, namely 1_1, that has only 1 part.
We have a(2) = 4 because we have the following colored compositions of n = 2: 2_1, 2_2, 1_1 + 1_1; hence, a(2) = 1 + 1 + 2 = 4.
We have a(3) = 14 because we have the following colored compositions of n = 3: 3_1, 3_2, 3_3, 1_1 + 2_1, 1_1 + 2_2, 2_1 + 1_1, 2_2 + 1_1, 1_1 + 1_1 + 1_1; hence, a(3) = 1 + 1 + 1 + 2 + 2 + 2 + 2 + 3 = 14.
We have a(14) = 46 because we have the following colored compositions of n = 4:
(i) 4_1, 4_2, 4_3, 4_4; with a total of 4 parts.
(ii) 1_1 + 3_1, 1_1 + 3_2, 1_1 + 3_3, 3_1 + 1_1, 3_2 + 1_1, 3_3 + 1_1, 2_1 + 2_1, 2_1 + 2_2, 2_2 + 2_1, 2_2 + 2_2; with a total of 2 x 10 = 20 parts.
(iii) 1_1 + 1_1 + 2_1, 1_1 + 1_1 + 2_2, 1_1 + 2_1 + 1_1, 1_1 + 2_2 + 1_1, 2_1 + 1_1 + 1_1, 2_2 + 1_1 + 1_1; with a total of 3 x 6 = 18 parts.
(iv) 1_1 + 1_1 + 1_1 + 1_1; with a total of 4 parts.
Hence, a(4) = 4 + 20 + 18 + 4 = 46.
(End)
- R. P. Grimaldi, Compositions and the alternate Fibonacci numbers, Congressus Numerantium, 186, 2007, 81-96.
- T. D. Noe, Table of n, a(n) for n = 1..200
- A. K. Agarwal, n-colour compositions, Indian J. Pure Appl. Math. 31 (11) (2000), 1421-1427.
- E. Barcucci, A. Del Lungo, S. Fezzi, and R. Pinzani, Nondecreasing Dyck paths and q-Fibonacci numbers, Discrete Math., 170 (1997), 211-217.
- C. G. Bower, Transforms (2).
- R. X. F. Chen and L. W. Shapiro, On sequences G(n) satisfying G(n)=(d+2)G(n-1)-G(n-2), J. Integer Seq. 10 (2007), Article #07.8.1; see Proposition 17.
- Éva Czabarka, Rigoberto Flórez, and Leandro Junes, Some Enumerations on Non-Decreasing Dyck Paths, Electron. J. Combin., 22(1) (2015), #P1.3.
- Éva Czabarka, Rigoberto Flórez, and Leandro Junes, A Discrete Convolution on the Generalized Hosoya Triangle, J. Integer Seq., 18 (2015), Article #15.1.6.
- Éva Czabarka, Rigoberto Flórez, Leandro Junes, and José L. Ramírez, Enumerations of peaks and valleys on non-decreasing Dyck paths, Discrete Math. 341 (10) (2018), 2789-2807.
- A. Denise and R. Simion, Two combinatorial statistics on Dyck paths, Discrete Math., 137 (1995), 155-176.
- Rigoberto Flórez, Leandro Junes, Luisa M. Montoya, and José L. Ramírez, Counting Subwords in Non-Decreasing Dyck Paths, Journal of Integer Sequences, Vol. 28 (2025), Article 25.1.6. See pp. 6, 14-15, 17, 19.
- Meghann Moriah Gibson, Combinatorics of compositions, Master of Science, Georgia Southern University, 2017.
- Meghann Moriah Gibson, Daniel Gray, and Hua Wang, Combinatorics of n-color compositions, Discrete Mathematics 341 (2018), 3209-3226.
- Milan Janjic, Hessenberg Matrices and Integer Sequences, J. Integer Seq. 13 (2010), Article #10.7.8.
- N. J. A. Sloane, Transforms.
- Index entries for two-way infinite sequences
- Index entries for linear recurrences with constant coefficients, signature (6,-11,6,-1).
-
with(combinat): L[0]:=2: L[1]:=1: for n from 2 to 60 do L[n]:=L[n-1] +L[n-2] end do: seq(2*fibonacci(2*n)*1/5+(1/5)*n*L[2*n],n=1..30); # Emeric Deutsch, Jul 21 2008
-
Table[Sum[k Binomial[n+k-1,2k-1],{k,n}],{n,30}] (* or *) LinearRecurrence[ {6,-11,6,-1},{1,4,14,46},30] (* Harvey P. Dale, Aug 01 2011 *)
-
a(n)=(2*n*fibonacci(2*n+1)+(2-n)*fibonacci(2*n))/5
Name clarified using a comment of the author by
Peter Luschny, Aug 03 2019
A007852
Number of antichains in rooted plane trees on n nodes.
Original entry on oeis.org
1, 2, 7, 29, 131, 625, 3099, 15818, 82595, 439259, 2371632, 12967707, 71669167, 399751019, 2247488837, 12723799989, 72474333715, 415046380767, 2388355096446, 13803034008095, 80082677184820, 466263828731640, 2723428895205210, 15954063529603565, 93711351580424391
Offset: 1
- O. Bodini, A. Genitrini and F. Peschanski, Enumeration and Random Generation of Concurrent Computations In proc. 23rd International Meeting on Probabilistic, Combinatorial and Asymptotic Methods for the Analysis of Algorithms (AofA'12), Discrete Mathematics and Theoretical Computer Science, pp 83-96, 2012.
- O. Bodini, A. Genitrini and F. Peschanski, A Quantitative Study of Pure Parallel Processes, Preprint, 2013.
- O. Bodini, A. Genitrini, and F. Peschanski, A Quantitative Study of Pure Parallel Processes, arXiv preprint arXiv:1407.1873 [cs.PL], 2014.
- G.-S. Cheon, H. Kim, and L. W. Shapiro, Mutation effects in ordered trees, arXiv preprint arXiv:1410.1249 [math.CO], 2014.
- V. Crismale, M. E. Griseta, and J. Wysoczański, Weakly Monotone Fock Space and Monotone Convolution of the Wigner Law, J Theor Probab 33, 268-294,2020.
- Martin Klazar, Twelve countings with rooted plane trees, European Journal of Combinatorics 18 (1997), 195-210; Addendum, 18 (1997), 739-740.
- Frank Ruskey, Listing and Counting Subtrees of a Tree, SIAM J. Computing, 10 (1981) 141-150.
- Index entries for sequences related to rooted trees
- Index entries for reversions of series
-
Rest[CoefficientList[Series[(1-(1-Sqrt[1-4*x])/2-Sqrt[1-5*x-(1-Sqrt[1-4*x])/2])/2, {x, 0, 20}], x]] (* Vaclav Kotesovec, Mar 08 2014 *)
-
a(n):=sum(binomial(2*i+1,i)*binomial(2*n-1,n-i-1),i,0,n)/((2*n-1)); /* Vladimir Kruchinin, Jun 09 2014 */
-
N = 33; x = 'x + O('x^N);
B = (1-sqrt(1-4*x))/2;
gf = (1-B-sqrt(1-5*x-B))/2;
v = Vec(gf)
\\ Joerg Arndt, Mar 14 2013
-
def a(n):
l = [0,1,2,7]
if n < 4:
return l[n]
for i in range(n-3):
l[i%4] = ( (-500*i+2000*i**3)*l[i%4]+(120-220*i-1380*i**2-920*i**3)*l[(i+1)%4]+(-1488-1626*i-387*i**2+21*i**3)*l[(i+2)%4]+(1088*i+1104+351*i**2+37*i**3)*l[(i+3)%4] ) // (+42*i**2+146*i+168+4*i**3)
return l[i%4]
# Antoine Genitrini, Mar 14 2013
A363308
Expansion of g.f. C(x*C(x)^3), where C(x) = 1 + x*C(x)^2 is the g.f. of the Catalan numbers (A000108).
Original entry on oeis.org
1, 1, 5, 26, 141, 790, 4542, 26668, 159333, 966038, 5930678, 36801660, 230491410, 1455283172, 9253674120, 59209786992, 380961295445, 2463303690790, 15998687418030, 104325569140156, 682768883525830, 4483232450501492, 29527005540912660, 195006621974036808
Offset: 0
G.f.: A(x) = 1 + x + 5*x^2 + 26*x^3 + 141*x^4 + 790*x^5 + 4542*x^6 + 26668*x^7 + 159333*x^8 + 966038*x^9 + 5930678*x^10 + ...
such that A(x) = C(x*C(x)^3), where
C(x) = 1 + x + 2*x^2 + 5*x^3 + 14*x^4 + 42*x^5 + ... + A000108(n)*x^n + ...
x*C(x)^3 = x + 3*x^2 + 9*x^3 + 28*x^4 + 90*x^5 + ... + A000245(n)*x^n + ...
Note that x*C(x)^3 = (C(x) - 1)*(1-x)/x - 1.
Also, the g.f. of related sequence A033296 begins
B(x) = 1 + x + 6*x^2 + 42*x^3 + 326*x^4 + 2706*x^5 + 23526*x^6 + ...
where A(x) = B(x/A(x)), B(x) = A(x*B(x)) = C(x*B(x)*C(x*B(x))^3).
-
{a(n) = if(n==0,1, sum(k=1,n, 3*k* binomial(2*k+1,k) * binomial(2*n+k,n-k) / ((2*k+1)*(2*n+k)) ) )}
for(n=0, 30, print1(a(n), ", "))
-
/* G.f. A(x) = C(x*C(x)^3), where C(x) = 1 + x*C(x)^2 */
{a(n) = my(C = (1 - sqrt(1 - 4*x +x^2*O(x^n)))/(2*x)); polcoeff( subst(C, x, x*C^3), n)}
for(n=0, 30, print1(a(n), ", "))
A153293
G.f.: A(x) = F(x*F(x)^3) = F(F(x)-1) where F(x) = 1 + x*F(x)^3 is the g.f. of A001764.
Original entry on oeis.org
1, 1, 6, 42, 317, 2508, 20517, 172180, 1474689, 12843768, 113444721, 1014062898, 9158151426, 83449247979, 766340138037, 7085966319858, 65919413472834, 616559331247512, 5794778945023698, 54700034442193302, 518375457403431600
Offset: 0
G.f.: A(x) = F(x*F(x)^3) = 1 + x + 6*x^2 + 42*x^3 + 317*x^4 +... where
F(x) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 +...
F(x)^2 = 1 + 2*x + 7*x^2 + 30*x^3 + 143*x^4 + 728*x^5 + 3876*x^6 +...
F(x)^3 = 1 + 3*x + 12*x^2 + 55*x^3 + 273*x^4 + 1428*x^5 + 7752*x^6 +...
-
S:= (1/2)*GAMMA(n+1/3)*GAMMA(n+2/3)*hypergeom([4/3, 5/3, -n+1], [5/2, 2*n+2], -27/4)*27^n*sqrt(3)/(Pi*GAMMA(2*n+2)):
1, seq(simplify(S),n=1..40); # Robert Israel, Dec 26 2017
-
F[x_] = 1 + InverseSeries[x/(1 + x)^3 + O[x]^21];
CoefficientList[F[F[x] - 1], x] (* Jean-François Alcover, Nov 02 2019 *)
-
{a(n)=if(n==0,1,sum(k=0,n,binomial(3*k+1,k)/(3*k+1)*binomial(3*(n-k)+3*k,n-k)*3*k/(3*(n-k)+3*k)))}
A153295
G.f.: A(x) = F(x*G(x)^2) where F(x) = G(x/F(x)) = 1 + x*F(x)^2 is the g.f. of A000108 (Catalan) and G(x) = F(x*G(x)) = 1 + x*G(x)^3 is the g.f. of A001764.
Original entry on oeis.org
1, 1, 4, 20, 110, 638, 3828, 23515, 146972, 930869, 5958094, 38462190, 250054804, 1635421543, 10750864640, 70987129653, 470542935654, 3129729034478, 20880459397920, 139689406647522, 936832986074664, 6297064070279195
Offset: 0
G.f.: A(x) = F(x*G(x)^2) = 1 + x + 4*x^2 + 20*x^3 + 110*x^4 +... where
F(x) = 1 + x + 2*x^2 + 5*x^3 + 14*x^4 + 42*x^5 + 132*x^6 +...
F(x)^2 = 1 + 2*x + 5*x^2 + 14*x^3 + 42*x^4 + 132*x^5 + 429*x^6 +...
G(x) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 +...
G(x)^2 = 1 + 2*x + 7*x^2 + 30*x^3 + 143*x^4 + 728*x^5 + 3876*x^6 +...
G(x)^3 = 1 + 3*x + 12*x^2 + 55*x^3 + 273*x^4 + 1428*x^5 +...
A(x)^2 = 1 + 2*x + 9*x^2 + 48*x^3 + 276*x^4 + 1656*x^5 +...
G(x)^2*A(x)^2 = 1 + 4*x + 20*x^2 + 110*x^3 + 638*x^4 +...
-
{a(n)=if(n==0,1,sum(k=0,n,binomial(2*k+1,k)/(2*k+1)*binomial(3*(n-k)+2*k,n-k)*2*k/(3*(n-k)+2*k)))}
A243585
Expansion of x*log'(C(C(x)-1)-1), C(x) = (1-sqrt(1-4*x))/(2*x).
Original entry on oeis.org
1, 4, 20, 106, 580, 3244, 18446, 106250, 618340, 3628600, 21438820, 127377980, 760346350, 4556473276, 27396081950, 165189725326, 998492094244, 6048338850560, 36706629690824, 223139239595840, 1358475322091620
Offset: 0
-
CoefficientList[Series[1/(Sqrt[(1-4*x)*(2*Sqrt[1-4*x]+5*x-2)/x]), {x, 0, 20}], x] (* Vaclav Kotesovec, Jun 08 2014 *)
A243585[n_] := Binomial[2 n, n] Hypergeometric2F1[1/2, -n, n + 1, -4];
Table[A243585[n], {n, 0, 20}] (* Peter Luschny, Aug 04 2019 *)
-
a(n):=sum(binomial(2*k,k)*binomial(2*n,n-k),k,0,n);
Showing 1-6 of 6 results.
Comments