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-7 of 7 results.

A007297 Number of connected graphs on n labeled nodes on a circle with straight-line edges that don't cross.

Original entry on oeis.org

1, 1, 4, 23, 156, 1162, 9192, 75819, 644908, 5616182, 49826712, 448771622, 4092553752, 37714212564, 350658882768, 3285490743987, 30989950019532, 294031964658430, 2804331954047160, 26870823304476690, 258548658860327880
Offset: 1

Views

Author

Keywords

Comments

Apart from the initial 1, reversion of g.f. for A162395 (squares with signs): see A263843.

Examples

			G.f. = x*(1 + x + 4*x^2 + 23*x^3 + 156*x^4 + 1162*x^5 + 9192*x^6 + 75819*x^7 + ...).
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A162395, A000290. 4th row of A107111. Row sums of A089434.
See A263843 for a variant.
Cf. A000108 (non-crossing set partitions), A001006, A001187, A054726 (non-crossing graphs), A054921, A099947, A194560, A293510, A323818, A324167, A324169, A324173.

Programs

  • Maple
    A007297:=proc(n) if n = 1 then 1 else add(binomial(3*n - 3, n + j)*binomial(j - 1, j - n + 1), j = n - 1 .. 2*n - 3)/(n - 1); fi; end;
  • Mathematica
    CoefficientList[ InverseSeries[ Series[(x-x^2)/(1+x)^3, {x, 0, 20}], x], x] // Rest (* From Jean-François Alcover, May 19 2011, after PARI prog. *)
    Table[Binomial[3n, 2n+1] Hypergeometric2F1[1-n, n, 2n+2, -1]/n, {n, 1, 20}] (* Vladimir Reshetnikov, Oct 25 2015 *)
  • PARI
    a(n)=if(n<0,0,polcoeff(serreverse((x-x^2)/(1+x)^3+O(x^(n+2))),n+1)) /* Ralf Stephan */

Formula

Apart from initial term, g.f. is the series reversion of (x-x^2)/(1+x)^3 (A162395). See A263843. - Vladimir Kruchinin, Feb 08 2013
G.f.: (g-z)/z, where g=-1/3+(2/3)*sqrt(1+9z)*sin((1/3)*arcsin((2+27z+54z^2)/2/(1+9*z)^(3/2))). - Emeric Deutsch, Dec 02 2002
a(n) = (1/n)*Sum_{k=0..n} binomial(3n, n-k-1)*binomial(n+k-1, k). - Paul Barry, May 11 2005
a(n) = 4^(n-1)*(Gamma(3*n/2-1)/Gamma(n/2+1)/Gamma(n) -Gamma((3*n-1)/2)/ Gamma( (n+1)/2)/Gamma(n+1)). - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = 4^n * binomial(3*n/2, n/2) / (9*n-6) - 4^(n-1) * binomial(3*(n-1)/2, (n-1)/2 ) / n. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
D-finite with recurrence: n*(n-1)*(3*n-4)*a(n) +36*(n-1)*a(n-1) -12*(3*n-8)*(3*n-1)*(3*n-7)*a(n-2)=0. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = (1/n)*Sum_{k=0..n} C(3n, k)*C(2n-k-2, n-1). - Paul Barry, Sep 27 2005
a(n) ~ (2-sqrt(3)) * 6^n * 3^(n/2) / (sqrt(2*Pi) * n^(3/2)). - Vaclav Kotesovec, Mar 17 2014
a(n) = binomial(3*n,2*n+1)*hypergeom([1-n,n], [2*n+2], -1)/n. - Vladimir Reshetnikov, Oct 25 2015
a(n) = 2*A078531(n) - A085614(n+1). - Vladimir Reshetnikov, Apr 24 2016

Extensions

Better description from Philippe Flajolet, Apr 20 2000
More terms from James Sellers, Aug 21 2000
Definition revised and initial a(1)=1 added by N. J. A. Sloane, Nov 05 2015 at the suggestion of Axel Boldt. Some of the formulas may now need to be adjusted slightly.

A088617 Triangle read by rows: T(n,k) = C(n+k,n)*C(n,k)/(k+1), for n >= 0, k = 0..n.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 6, 10, 5, 1, 10, 30, 35, 14, 1, 15, 70, 140, 126, 42, 1, 21, 140, 420, 630, 462, 132, 1, 28, 252, 1050, 2310, 2772, 1716, 429, 1, 36, 420, 2310, 6930, 12012, 12012, 6435, 1430, 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862
Offset: 0

Views

Author

N. J. A. Sloane, Nov 23 2003

Keywords

Comments

Row sums: A006318 (Schroeder numbers). Essentially same as triangle A060693 transposed.
T(n,k) is number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0), having k U's. E.g., T(2,1)=3 because we have UHD, UDH and HUD. - Emeric Deutsch, Dec 06 2003
Little Schroeder numbers A001003 have a(n) = Sum_{k=0..n} A088617(n,k)*(-1)^(n-k)*2^k. - Paul Barry, May 24 2005
Conjecture: The expected number of U's in a Schroeder n-path is asymptotically Sqrt[1/2]*n for large n. - David Callan, Jul 25 2008
T(n, k) is also the number of order-preserving and order-decreasing partial transformations (of an n-chain) of width k (width(alpha) = |Dom(alpha)|). - Abdullahi Umar, Oct 02 2008
The antidiagonals of this lower triangular matrix are the rows of A055151. - Tom Copeland, Jun 17 2015

Examples

			Triangle begins:
  [0] 1;
  [1] 1,  1;
  [2] 1,  3,   2;
  [3] 1,  6,  10,    5;
  [4] 1, 10,  30,   35,    14;
  [5] 1, 15,  70,  140,   126,    42;
  [6] 1, 21, 140,  420,   630,   462,   132;
  [7] 1, 28, 252, 1050,  2310,  2772,  1716,   429;
  [8] 1, 36, 420, 2310,  6930, 12012, 12012,  6435,  1430;
  [9] 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862;
		

References

  • Charles Jordan, Calculus of Finite Differences, Chelsea 1965, p. 449.

Crossrefs

Programs

  • Magma
    [[Binomial(n+k,n)*Binomial(n,k)/(k+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jun 18 2015
    
  • Maple
    R := n -> simplify(hypergeom([-n, n + 1], [2], -x)):
    Trow := n -> seq(coeff(R(n, x), x, k), k = 0..n):
    seq(print(Trow(n)), n = 0..9); # Peter Luschny, Apr 26 2022
  • Mathematica
    Table[Binomial[n+k, n] Binomial[n, k]/(k+1), {n,0,10}, {k,0,n}]//Flatten (* Michael De Vlieger, Aug 10 2017 *)
  • PARI
    {T(n, k)= if(k+1, binomial(n+k, n)*binomial(n, k)/(k+1))}
    
  • SageMath
    flatten([[binomial(n+k, 2*k)*catalan_number(k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 22 2022

Formula

Triangle T(n, k) read by rows; given by [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...] DELTA [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
T(n, k) = A085478(n, k)*A000108(k); A000108 = Catalan numbers. - Philippe Deléham, Dec 05 2003
Sum_{k=0..n} T(n, k)*x^k*(1-x)^(n-k) = A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11. - Philippe Deléham, Aug 18 2005
Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A107841(n), A080243(n), A000007(n), A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Oct 18 2007
O.g.f. (with initial 1 excluded) is the series reversion with respect to x of (1-t*x)*x/(1+x). Cf. A062991 and A089434. - Peter Bala, Jul 31 2012
G.f.: 1 + (1 - x - T(0))/y, where T(k) = 1 - x*(1+y)/( 1 - x*y/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 03 2013
From Peter Bala, Jul 20 2015: (Start)
O.g.f. A(x,t) = ( 1 - x - sqrt((1 - x)^2 - 4*x*t) )/(2*x*t) = 1 + (1 + t)*x + (1 + 3*t + 2*t^2)*x^2 + ....
1 + x*(dA(x,t)/dx)/A(x,t) = 1 + (1 + t)*x + (1 + 4*t + 3*t^2)*x^2 + ... is the o.g.f. for A123160.
For n >= 1, the n-th row polynomial equals (1 + t)/(n+1)*Jacobi_P(n-1,1,1,2*t+1). Removing a factor of 1 + t from the row polynomials gives the row polynomials of A033282. (End)
From Tom Copeland, Jan 22 2016: (Start)
The o.g.f. G(x,t) = {1 - (2t+1) x - sqrt[1 - (2t+1) 2x + x^2]}/2x = (t + t^2) x + (t + 3t^2 + 2t^3) x^2 + (t + 6t^2 + 10t^3 + 5t^3) x^3 + ... generating shifted rows of this entry, excluding the first, was given in my 2008 formulas for A033282 with an o.g.f. f1(x,t) = G(x,t)/(1+t) for A033282. Simple transformations presented there of f1(x,t) are related to A060693 and A001263, the Narayana numbers. See also A086810.
The inverse of G(x,t) is essentially given in A033282 by x1, the inverse of f1(x,t): Ginv(x,t) = x [1/(t+x) - 1/(1+t+x)] = [((1+t) - t) / (t(1+t))] x - [((1+t)^2 - t^2) / (t(1+t))^2] x^2 + [((1+t)^3 - t^3) / (t(1+t))^3] x^3 - ... . The coefficients in t of Ginv(xt,t) are the o.g.f.s of the diagonals of the Pascal triangle A007318 with signed rows and an extra initial column of ones. The numerators give the row o.g.f.s of signed A074909.
Rows of A088617 are shifted columns of A107131, whose reversed rows are the Motzkin polynomials of A055151, related to A011973. The diagonals of A055151 give the rows of A088671, and the antidiagonals (top to bottom) of A088617 give the rows of A107131 and reversed rows of A055151. The diagonals of A107131 give the columns of A055151. The antidiagonals of A088617 (bottom to top) give the rows of A055151.
(End)
T(n, k) = [x^k] hypergeom([-n, 1 + n], [2], -x). - Peter Luschny, Apr 26 2022

A062991 Coefficient triangle for certain polynomials N(2; n,x) (rising powers of x).

Original entry on oeis.org

1, 2, -1, 5, -6, 2, 14, -28, 20, -5, 42, -120, 135, -70, 14, 132, -495, 770, -616, 252, -42, 429, -2002, 4004, -4368, 2730, -924, 132, 1430, -8008, 19656, -27300, 23100, -11880, 3432, -429, 4862, -31824, 92820, -157080, 168300, -116688, 51051, -12870, 1430
Offset: 0

Views

Author

Wolfdieter Lang, Jul 12 2001

Keywords

Comments

The g.f. for the sequence of column m of triangle A009766(n,m) (or Catalan A033184(n,n-m) diagonals) is N(2; m-1,x)*(x^m)/(1-x)^(m+1), m >= 1, with N(2; n,x) = Sum_{k=0..n} T(n,k)*x^k.
For k=0..1 the column sequences give A000108(n+1) (Catalan), -A002694. The row sums give A000012 (powers of 1) and (unsigned) A062992.
Another version of [1, 1, 1, 1, 1, 1, 1, 1, ...] DELTA [0, -1, -1, -1, -1, -1, -1, -1, ...] = 1; 1, 0; 2, -1, 0; 5, -6, 2, 0; 14, -28, 20, -5, 0; 42, -120, 135, -70, 14, 0; ... where DELTA is Deléham's operator defined in A084938.
The positive triangle is |T(n,k)| = binomial(2*n+2, n-k)*binomial(n+k, k)/(n+1). - Paul Barry, May 11 2005

Examples

			The triangle N2 = {a(n,k)} begins:
n\k      0       1      2       3       4       5      6       7     8     9
----------------------------------------------------------------------------
0:       1
1:       2      -1
2:       5      -6      2
3:      14     -28     20      -5
4:      42    -120    135     -70      14
5:     132    -495    770    -616     252     -42
6:     429   -2002   4004   -4368    2730    -924    132
7:    1430   -8008  19656  -27300   23100  -11880   3432    -429
8:    4862  -31824  92820 -157080  168300 -116688  51051  -12870  1430
9:   16796 -125970 426360 -852720 1108536 -969969 570570 -217360 48620 -4862
... formatted by _Wolfdieter Lang_, Jan 20 2020
N(2; 2, x)= 5 - 6*x + 2*x^2.
		

Crossrefs

For an unsigned version see Borel's triangle, A234950.
Sums include: A000012 (row), A000079 (diagonal), A064062 (signed row), A071356 (signed diagonal).

Programs

  • Magma
    A062991:= func< n,k | (-1)^k*Binomial(2*n+2,n-k)*Binomial(n+k,k)/(n+1) >;
    [A062991(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 27 2024
    
  • Mathematica
    T[n_, k_] := 2 (-1)^k Binomial[2n+1, n] (n-k+1) Binomial[n+1, k]/((k+n+1)(k+n+2));
    Table[T[n, k], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, Sep 19 2018 *)
  • SageMath
    def A062991(n,k): return (-1)^k*binomial(2*n+2,n-k)*binomial(n+k,k)/(n+1)
    flatten([[A062991(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Sep 27 2024

Formula

T(n, k) = [x^k] N(2; n, x) with N(2; n, x) = (N(2; n-1, x) - A000108(n)*(1-x)^(n+1))/x, N(2; 0, x) = 1.
T(n, k) = T(n-1, k+1) + (-1)^k*binomial(n+1, k+1)*binomial(2*n+1, n)/(2*n+1) if k=0, .., (n-2); T(n, k) = (-1)^k*binomial(n+1, k+1)*binomial(2*n+1, n)/(2*n+1) if k=(n-1) or n; else 0.
O.g.f. (with offset 1) is the series reversion w.r.t. x of x*(1+x*t)/(1+x)^2. If R(n,t) denotes the n-th row polynomial of this triangle then R(n,1-t) is the n-th row polynomial of A009766. Cf. A089434. - Peter Bala, Jul 15 2012
From G. C. Greubel, Sep 27 2024: (Start)
Sum_{k=0..n} T(n, k) = A000012(n).
Sum_{k=0..n} (-1)^k*T(n, k) = A064062(n+1).
Sum_{k=0..floor(n/2)} T(n-k, k) = A000079(n).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = A071356(n). (End)

A108410 Triangle T(n,k) read by rows: number of 12312-avoiding matchings on [2n] with exactly k crossings (n >= 1, 0 <= k <= n-1).

Original entry on oeis.org

1, 2, 1, 5, 5, 2, 14, 21, 15, 5, 42, 84, 84, 49, 14, 132, 330, 420, 336, 168, 42, 429, 1287, 1980, 1980, 1350, 594, 132, 1430, 5005, 9009, 10725, 9075, 5445, 2145, 429, 4862, 19448, 40040, 55055, 55055, 40898, 22022, 7865, 1430, 16796, 75582
Offset: 1

Views

Author

Ralf Stephan, Jun 03 2005

Keywords

Examples

			Triangle begins
     1;
     2,     1;
     5,     5,     2;
    14,    21,    15,     5;
    42,    84,    84,    49,    14;
   132,   330,   420,   336,   168,    42;
   429,  1287,  1980,  1980,  1350,   594,   132;
  1430,  5005,  9009, 10725,  9075,  5445,  2145,  429;
  4862, 19448, 40040, 55055, 55055, 40898, 22022, 7865, 1430;
		

Crossrefs

Left-hand columns include A000108 and A002054. Right-hand columns include A000108 and A007851+1. Row sums are A001764. A089434.

Programs

  • Maple
    T:=(n,k)->binomial(n-1+k,n-1)*binomial(2*n-k,n+1)/n: for n from 1 to 10 do seq(T(n,k),k=0..n-1) od; # yields sequence in triangular form - Emeric Deutsch, Dec 19 2006
  • Mathematica
    T[n_, k_] := Binomial[n + k - 1, n - 1]*Binomial[2*n - k, n + 1]/n;
    Table[T[n, k], {n, 1, 10}, {k, 0, n - 1}] // Flatten (* Jean-François Alcover, Nov 11 2017, after Emeric Deutsch *)
  • PARI
    T(n,k) = binomial(n-1+k, n-1)*binomial(2*n-k, n+1)/n; \\ Andrew Howroyd, Nov 06 2017

Formula

T(n, k) = Sum_{i=n..2*n-1} (-1)^(n+k+i)/i*C(i, n)*C(3*n, i+1+n)*C(i-n, k).
T(n,k) = C(n-1+k,n-1)*C(2*n-k,n+1)/n, (0 <= k <= n-1). [Chen et al.] - Emeric Deutsch, Dec 19 2006
O.g.f. equals the series reversion with respect to x of x*(1 + x*(1 - t))/(1 + x)^3. If R(n,t) is the n-th row polynomial of this triangle then R(n, 1+t) is the n-th row polynomial of A089434. - Peter Bala, Jul 15 2012

A045742 Number of interior faces in all noncrossing connected graphs on n nodes on a circle.

Original entry on oeis.org

0, 1, 13, 141, 1456, 14778, 149031, 1499773, 15089932, 151927854, 1531242362, 15451614738, 156114597744, 1579223536788, 15993825704427, 162159485143581, 1645827425223220, 16720488433727910, 170023231905932790
Offset: 2

Views

Author

Keywords

Crossrefs

Cf. A089434.

Programs

  • Maple
    A045742 := proc(n)
            binomial(3*n-3,n-3)*hypergeom([n, 3 - n], [2*n + 1], -1) ;
            simplify(%) ;
    end proc: # R. J. Mathar, Mar 27 2012
  • Mathematica
    Rest[Table[Sum[(k Binomial[n+k-2,k]Binomial[3n-3,n-2-k])/(n-1),{k,n-2}], {n,20}]] (* Harvey P. Dale, Nov 29 2011 *)
  • PARI
    a(n) = if(n>1, sum(k=1, n-2, k*binomial(n+k-2, k)*binomial(3*n-3, n-2-k))/(n-1)); \\ Andrew Howroyd, Nov 12 2017

Formula

a(n) = Sum_{k=1..n-2} k*binomial(n+k-2, k)*binomial(3*n-3, n-2-k)/(n-1).
a(n) = Sum_{k=1..n-2} k*A089434(n,k). - Andrew Howroyd, Nov 12 2017
a(n) ~ (9 - 5*sqrt(3)) * 2^(n - 5/2) * 3^((3*n-4)/2) / sqrt(Pi*n). - Vaclav Kotesovec, Dec 10 2021
Conjecture D-finite with recurrence n*(n-1)*(523*n-1993)*a(n) -6*(n-1)*(759*n^2-6019*n+12790)*a(n-1) -12*(n-2)*(4707*n^2-17937*n+6260)*a(n-2) +72*(3*n-10)*(759*n-2224)*(3*n-11)*a(n-3)=0. - R. J. Mathar, Jul 26 2022

A089433 Number of noncrossing connected graphs on n nodes having exactly two interior faces.

Original entry on oeis.org

2, 30, 315, 2856, 23940, 191268, 1480050, 11196900, 83304936, 611931320, 4450217772, 32104210320, 230080173960, 1639890119016, 11634355574100, 82216112723640, 579022013389050, 4065827626164150, 28475852003986695
Offset: 4

Views

Author

Emeric Deutsch, Dec 28 2003

Keywords

Examples

			a(4)=2 because the only connected graphs on the nodes A,B,C,D having exactly two interior faces are {AB,BC,CD,DA,AC} and {AB,BC,CD,DA,BD}.
		

Crossrefs

Column k=2 of A089434.
Cf. A007297.

Programs

  • PARI
    a(n) = n*binomial(3*n-3, n-4)/2; \\ Michel Marcus, Oct 26 2015

Formula

a(n) = n*binomial(3n-3, n-4)/2.
D-finite with recurrence -2*(2*n+1)*(n-4)*a(n) +3*(3*n-4)*(3*n-5)*a(n-1)=0. - R. J. Mathar, Jul 26 2022

A127537 Triangle read by rows: T(n,k) (n >= 2, 1 <= k <= 2n-3) is the number of non-crossing connected graphs on n nodes on a circle, having k edges. Rows are indexed 2,3,4,...; columns are indexed 0,1,2,....

Original entry on oeis.org

1, 0, 3, 1, 0, 0, 12, 9, 2, 0, 0, 0, 55, 66, 30, 5, 0, 0, 0, 0, 273, 455, 315, 105, 14, 0, 0, 0, 0, 0, 1428, 3060, 2856, 1428, 378, 42, 0, 0, 0, 0, 0, 0, 7752, 20349, 23940, 15960, 6300, 1386, 132, 0, 0, 0, 0, 0, 0, 0, 43263, 134596, 191268, 159390, 83490, 27324, 5148, 429
Offset: 2

Views

Author

Emeric Deutsch, Jan 24 2007

Keywords

Comments

Row n contains 2n-3 terms, the first n-2 of which are equal to 0.
T(n,n-1) = A001764(n-1). T(n,2n-3) = A000108(n-2) (the Catalan numbers).
T(n,k) = A089434(n,k+1-n).
Sum_{k=n-1..2n-3} k*T(n,k) = A045741(n).
Sum_{n=k..2k-2} T(n,k) = A065065(k).

Examples

			Triangle starts:
  1;
  0,  3,  1;
  0,  0, 12,  9,  2;
  0,  0,  0, 55, 66, 30,  5;
		

Crossrefs

Programs

  • Maple
    T:=(n,k)->binomial(3*n-3,n+k)*binomial(k-1,k-n+1)/(n-1): for n from 2 to 10 do seq(T(n,k),k=1..2*n-3) od; # yields sequence in triangular form
  • Mathematica
    T[n_, k_] := Binomial[3n - 3, n + k] Binomial[k - 1, k - n + 1]/(n - 1);
    Table[T[n, k], {n, 2, 10}, {k, 1, 2n - 3}] // Flatten (* Jean-François Alcover, Jul 29 2018 *)

Formula

T(n,k) = C(3n-3,n+k)C(k-1,k-n+1)/(n-1) (n >= 2, 0 <= k <= 2n-3).
G.f.: G=G(t,z) satisfies tG^3 + tG^2 - z(1+2t)G + z^2*(1+t) = 0.

Extensions

Keyword tabl changed to tabf by Michel Marcus, Apr 09 2013
Showing 1-7 of 7 results.