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.

Showing 1-4 of 4 results.

A050351 Number of 3-level labeled linear rooted trees with n leaves.

Original entry on oeis.org

1, 1, 5, 37, 365, 4501, 66605, 1149877, 22687565, 503589781, 12420052205, 336947795317, 9972186170765, 319727684645461, 11039636939221805, 408406422098722357, 16116066766061589965, 675700891505466507541
Offset: 0

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Comments

Lists of lists of sets.

Examples

			G.f. = 1 + x + 5*x^2 + 37*x^3 + 365*x^4 + 4501*x^5 + 66605*x^6 + ...
		

References

  • T. S. Motzkin, Sorting numbers ...: for a link to an annotated scanned version of this paper see A000262.
  • T. S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176.

Crossrefs

Equals 1/2 * A004123(n) for n>0.

Programs

  • Maple
    with(combstruct); SeqSeqSetL := [T, {T=Sequence(S), S=Sequence(U,card >= 1), U=Set(Z,card >=1)},labeled];
  • Mathematica
    With[{nn=20},CoefficientList[Series[(2-E^x)/(3-2*E^x),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Feb 29 2012 *)
    a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ 1/(2 - 1/(2 - Exp[x])), {x, 0, n}]]; (* Michael Somos, Nov 28 2014 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( 1/(2 - 1/(2 - exp(x + x * O(x^n)))), n))};
    
  • PARI
    {a(n)=if(n==0, 1, (1/6)*round(suminf(k=1, k^n * (2/3)^k *1.)))} \\ Paul D. Hanna, Nov 28 2014
    
  • Sage
    A050351 = lambda n: sum(stirling_number2(n,k)*(2^(k-1))*factorial(k) for k in (0..n)) if n>0 else 1
    [A050351(n) for n in (0..17)] # Peter Luschny, Jan 18 2016

Formula

E.g.f.: (2-exp(x))/(3-2*exp(x)).
a(n) is asymptotic to (1/6)*n!/log(3/2)^(n+1). - Benoit Cloitre, Jan 30 2003
For m-level trees (m>1), e.g.f. is (m-1-(m-2)*e^x)/(m-(m-1)*e^x) and number of trees is 1/(m*(m-1))*sum(k>=0, (1-1/m)^k*k^n). Here m=3, so a(n)=(1/6)*sum(k>=0, (2/3)^k*k^n) (for n>0). - Benoit Cloitre, Jan 30 2003
a(n) = Sum_{k=1..n} Stirling2(n, k)*k!*2^(k-1). - Vladeta Jovovic, Sep 28 2003
Recurrence: a(n+1) = 1 + 2*Sum_{j=1..n} binomial(n+1, j)*a(j). - Jon Perry, Apr 25 2005
With p(n) = the number of integer partitions of n, p(i) = the number of parts of the i-th partition of n, d(i) = the number of different parts of the i-th partition of n, p(j, i) = the j-th part of the i-th partition of n, m(i, j) = multiplicity of the j-th part of the i-th partition of n, sum_{i=1}^{p(n)} = sum over i and prod_{j=1}^{d(i)} = product over j one has: a(n)=sum_{i=1}^{p(n)}(n!/(prod_{j=1}^{p(i)}p(i, j)!))*(p(i)!/(prod_{j=1}^{d(i)} m(i, j)!))*2^(p(i)-1). - Thomas Wieder, May 18 2005
Let f(x) = (1+x)*(1+2*x). Let D be the operator g(x) -> d/dx(f(x)*g(x)). Then for n>=1, a(n) = D^(n-1)(1) evaluated at x = 1/2. Compare with the result A000670(n) = D^(n-1)(1) at x = 0. See also A194649. - Peter Bala, Sep 05 2011
E.g.f.: 1 + x/(G(0)-3*x) where G(k)= x + k + 1 - x*(k+1)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Jul 11 2012
a(n) = (1/6) * Sum_{k>=1} k^n * (2/3)^k for n>0. - Paul D. Hanna, Nov 28 2014
E.g.f. A(x) satisfies 0 = 2 - A'(x) - 7*A(x) + 6*A(x)^2. - Michael Somos, Nov 28 2014

A050352 Number of 4-level labeled linear rooted trees with n leaves.

Original entry on oeis.org

1, 1, 7, 73, 1015, 17641, 367927, 8952553, 248956855, 7788499561, 270732878647, 10351919533033, 431806658432695, 19512813265643881, 949587798053709367, 49512355251796513513, 2753726282896986372535, 162725978752448205162601
Offset: 0

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Crossrefs

Equals 1/3 * A032033(n) for n>0.

Programs

  • Magma
    [1] cat [(&+[3^(j-1)*Factorial(j)*StirlingSecond(n,j): j in [1..n]]) : n in [0..20]]; // G. C. Greubel, Jun 08 2020
    
  • Maple
    seq(coeff(series( (3-2*exp(x))/(4-3*exp(x)), x, n+1)*n!, x, n), n = 0..20); # G. C. Greubel, Jun 08 2020
  • Mathematica
    With[{nn=20}, CoefficientList[Series[(3-2Exp[x])/(4-3Exp[x]),{x,0,nn}], x]*Range[0,nn]!] (* Harvey P. Dale, Aug 16 2012 *)
  • PARI
    a(n)=n!*if(n<0,0,polcoeff((3-2*exp(x))/(4-3*exp(x))+O(x^(n+1)),n))
    
  • PARI
    {a(n)=if(n==0, 1, (1/12)*round(suminf(k=1, k^n * (3/4)^k *1.)))} \\ Paul D. Hanna, Nov 28 2014
    
  • Sage
    [1]+[sum( 3^(j-1)*factorial(j)*stirling_number2(n,j) for j in (1..n)) for n in (1..20)] # G. C. Greubel, Jun 08 2020

Formula

E.g.f.: (3 - 2*exp(x))/(4 - 3*exp(x)).
a(n) is asymptotic to (1/12)*n!/log(4/3)^(n+1). - Benoit Cloitre, Jan 30 2003
For m-level trees (m>1), e.g.f. is (m-1-(m-2)*e^x)/(m-(m-1)*e^x) and number of trees is 1/(m*(m-1))*sum(k>=0, (1-1/m)^k*k^n). Here m=4, so a(n)=(1/12)*sum(k>=0, (3/4)^k*k^n) (for n>0). - Benoit Cloitre, Jan 30 2003
Let f(x) = (1+x)*(1+2*x). Let D be the operator g(x) -> d/dx(f(x)*g(x)). Then for n>=1, a(n) = D^(n-1)(1) evaluated at x = 1. Compare with the result A000670(n) = D^(n-1)(1) at x = 0. See also A194649. - Peter Bala, Sep 05 2011
E.g.f.: 1 + x/(G(0)-4*x) where G(k)= x + k + 1 - x*(k+1)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Jul 11 2012
a(n) = (1/12) * Sum_{k>=1} k^n * (3/4)^k for n>0. - Paul D. Hanna, Nov 28 2014
a(n) = Sum_{k=1..n} Stirling2(n, k) * k! * 3^(k-1). - Paul D. Hanna, Nov 28 2014, after Vladeta Jovovic in A050351
a(n) = 1 + 3 * Sum_{k=1..n-1} binomial(n,k) * a(k). - Ilya Gutkovskiy, Jun 08 2020

A050353 Number of 5-level labeled linear rooted trees with n leaves.

Original entry on oeis.org

1, 1, 9, 121, 2169, 48601, 1306809, 40994521, 1469709369, 59277466201, 2656472295609, 130952452264921, 7042235448544569, 410269802967187801, 25740278881968596409, 1730295054262416751321, 124066865052334175027769
Offset: 0

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Crossrefs

Equals 1/4 * A094417(n) for n>0.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 20); Coefficients(R!(Laplace( (4-3*Exp(x))/(5-4*Exp(x)) ))); // G. C. Greubel, Jun 08 2020
    
  • Maple
    seq(coeff(series( (4-3*exp(x))/(5-4*exp(x)), x, n+1)*n!, x, n), n = 0..20); # G. C. Greubel, Jun 08 2020
  • Mathematica
    max = 16; f[x_] := (4-3*E^x) / (5-4*E^x); CoefficientList[ Series[ f[x], {x, 0, max}], x]*Range[0, max]! (* Jean-François Alcover, Nov 14 2011, after g.f. *)
  • PARI
    a(n)=n!*if(n<0,0,polcoeff((4-3*exp(x))/(5-4*exp(x))+O(x^(n+1)),n))
    
  • PARI
    {a(n)=if(n==0,1,(1/20)*round(suminf(k=1, k^n * (4/5)^k *1.)))} \\ Paul D. Hanna, Nov 28 2014
    
  • Sage
    [1]+[sum( 4^(j-1)*factorial(j)*stirling_number2(n,j) for j in (1..n)) for n in (1..20)] # G. C. Greubel, Jun 08 2020

Formula

E.g.f.: (4 - 3*exp(x))/(5 - 4*exp(x)).
a(n) is asymptotic to (1/20)*n!/log(5/4)^(n+1). More generally if m>1, the number of m-level labeled linear rooted trees with n leaves is asymptotic to n!/log(m/(m-1))^(n+1)/(m^2-m). - Benoit Cloitre, Jan 30 2003
For m-level trees (m>1), e.g.f. is (m-1-(m-2)*e^x)/(m-(m-1)*e^x) and number of trees is 1/(m*(m-1))*sum(k>=0, (1-1/m)^k*k^n). Here m=5, so a(n)=(1/20)*sum(k>=0, (4/5)^k*k^n) (for n>0). - Benoit Cloitre, Jan 30 2003
Let f(x) = (1+x)*(1+2*x). Let D be the operator g(x) -> d/dx(f(x)*g(x)). Then for n>=1, a(n) = D^(n-1)(1) evaluated at x = 3/2. Compare with the result A000670(n) = D^(n-1)(1) at x = 0. See also A194649. - Peter Bala, Sep 05 2011
E.g.f.: 1 + x/(G(0)-5*x) where G(k)= x + k + 1 - x*(k+1)/G(k+1); (continued fraction Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Jul 11 2012
a(n) = (1/20) * Sum_{k>=1} k^n * (4/5)^k for n>0. - Paul D. Hanna, Nov 28 2014
a(n) = Sum_{k=1..n} Stirling2(n, k) * k! * 4^(k-1). - Paul D. Hanna, Nov 28 2014, after Vladeta Jovovic in A050351
a(n) = 1 + 4 * Sum_{k=1..n-1} binomial(n,k) * a(k). - Ilya Gutkovskiy, Jun 08 2020

A098384 Triangle read by rows of coefficients used to generate diagonals of ordered factorizations as displayed in A098348.

Original entry on oeis.org

1, 3, 2, 13, 18, 8, 75, 158, 144, 48, 541, 1530, 2120, 1440, 384, 4683, 16622, 30960, 31920, 17280, 3840
Offset: 0

Views

Author

Alford Arnold, Sep 06 2004

Keywords

Comments

Note that the table includes the well-known sequence (A000165) discussed by Gordon on pages 636-645 of AMM 106 (1999).

Examples

			The table begins:
1
3 2
13 18 8
75 158 144 48
541 1530 2120 1440 384
The binomial transform of (13,18,8) yields 13,31,57,91,...
The binomial transform of 13,31,57,91,... yields 13,44,132,368,... A098385
		

Crossrefs

Formula

From Peter Bala, Apr 20 2012: (Start)
The following formulas are all conjectural:
T(n,k) = 2^k*sum {i = k+1..n+1} binomial(i,k+1)*(i-1)!*Stirling2(n+1,i) = 1/(k+1)*A194649(n+1,k).
Recurrence equation:
T(n,k) = 2*k*T(n-1,k-1) + 3*(k+1)*T(n-1,k) + (k+2)*T(n-1,k+1).
E.g.f.: exp(x)/((2-exp(x))*(2*t+2-(2*t+1)*exp(x))) = 1 + (3+2*t)*x + (13+18*t+8*t^2)*x^2/2! + ....
Column n generating function: 2^n*exp(x)*(1-exp(x))^n/(exp(x)-2)^(n+2) for n >= 0.
(End)
Showing 1-4 of 4 results.