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.

Previous Showing 11-20 of 29 results. Next

A006414 Number of nonseparable toroidal tree-rooted maps with n + 2 edges and n + 1 vertices.

Original entry on oeis.org

1, 9, 40, 125, 315, 686, 1344, 2430, 4125, 6655, 10296, 15379, 22295, 31500, 43520, 58956, 78489, 102885, 133000, 169785, 214291, 267674, 331200, 406250, 494325, 597051, 716184, 853615, 1011375, 1191640, 1396736, 1629144, 1891505, 2186625, 2517480, 2887221
Offset: 0

Views

Author

Keywords

Comments

The number of faces is 1.
a(n) = K(Oa(2,3,n)), Kekulé numbers of certain benzenoid structures (see the Cyvin - Gutman reference).
Sequence of partial sums of A006322. - L. Edson Jeffery, Dec 13 2011
The sequence b(n) = a(n-2) with a(-1) = 0, for n >= 1, is b(n) = n^3*(n^2 - 1)/4!. It is obtained by comparing the result for the powers n^5 from Worpitzky's identity (see a formula in A000584) with the result obtained from the counting of degrees of freedom for the decomposition of a rank 5 tensor in n dimensions via the standard Young tableaux version with 5 boxes corresponding to the seven partitions of 5. The difference of the two versions gives: 10*(binomial(n+3, 5) + 3*binomial(n+2, 5) + binomial(n+1, 5)) = 5*n*(binomial(n+2, 4) + binomial(n+1, 4)) = 10*b(n). See the formula for a(n) below. - Wolfdieter Lang, Jul 18 2019

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988, p. 105, eq. (ii), and p. 186.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Differences of A006542 (C(n, 3)*C(n-1, 3)/4).

Programs

Formula

a(n) = (n+1)*(n+2)^3*(n+3)/24. - N. J. A. Sloane, Apr 02 2004
a(n) = (n+2)^3*((n+2)^2 - 1)/24. - Paul Richards, Mar 04 2007
G.f.: (1 + 3*x + x^2)/(1-x)^6. - Colin Barker, Feb 21 2012
a(n) = (Sum_{k=0..n+1} k*(n+1)*((n+1)^2 - k^2))/6 for n > 0, which is the sum of all areas of Pythagorean triangles with arms 2*k*(n+1) and (n+1)^2 - k^2 with hypotenuse k^2 + (n+1)^2. - J. M. Bergot, May 12 2014
a(n) = A143945(n+2)/8. - J. M. Bergot, Jun 14 2014
Sum_{n>=0} 1/a(n) = 30 - 24*zeta(3). - Jaume Oliver Lafont, Jul 09 2017
a(n) = binomial(n+5, 5) + 3*binomial(n+4, 5) + binomial(n+3, 5) = ((n+2)/2)*(binomial(n+4, 4) + binomial(n+3, 4)), for n >= 0. See a comment above on the sequence b(n) = a(n-2) = n^3*(n^2 - 1)/4!. - Wolfdieter Lang, Jul 19 2019
E.g.f.: (24 + 192*x + 276*x^2 + 124*x^3 + 20*x^4 + x^5)*exp(x)/4!. - G. C. Greubel, Sep 02 2019
Sum_{n>=0} (-1)^n/a(n) = 18*zeta(3) + 48*log(2) - 54. - Amiram Eldar, Jan 09 2022

Extensions

More terms from Robert Newstedt (Patternfinder(AT)webtv.net)
Name clarified by Andrew Howroyd, Apr 05 2021

A006857 a(n) = binomial(n+5,5) * binomial(n+5,4)/(n+5).

Original entry on oeis.org

1, 15, 105, 490, 1764, 5292, 13860, 32670, 70785, 143143, 273273, 496860, 866320, 1456560, 2372112, 3755844, 5799465, 8756055, 12954865, 18818646, 26883780, 37823500, 52474500, 71867250, 97260345, 130179231, 172459665, 226296280, 294296640, 379541184
Offset: 0

Views

Author

Keywords

Comments

Number of permutations of n+5 that avoid the pattern 132 and have exactly 4 descents.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
Partial sums of A114242. - Peter Bala, Sep 21 2007
Dimensions of certain Lie algebra (see reference for precise definition).

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (pp. 167-169, Table 10.5/II/1).
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 239.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

The expression binomial(m+n-1,n)^2-binomial(m+n,n+1)*binomial(m+n-2,n-1) for the values m = 2 through 14 produces the sequences A000012, A000217, A002415, A006542, A006857, A108679, A134288, A134289, A134290, A134291, A140925, A140935, A169937.
5th column of the table of Narayana numbers A001263.

Programs

  • Magma
    A006857:= func< n | Binomial(n+4,3)*Binomial(n+5,5)/4 >;
    [A006857(n): n in [0..40]]; // G. C. Greubel, Mar 12 2025
    
  • Maple
    a:=n->(n+1)*(n+2)^2*(n+3)^2*(n+4)^2*(n+5)/2880: seq(a(n),n=0..38); # Emeric Deutsch, Nov 18 2005
  • Mathematica
    Table[Binomial[n+5,5] * Binomial[n+5,4]/(n+5), {n, 0, 50}] (* T. D. Noe, May 29 2012 *)
    LinearRecurrence[{9,-36,84,-126,126,-84,36,-9,1},{1,15,105,490,1764,5292,13860,32670,70785},40] (* Harvey P. Dale, Oct 19 2024 *)
  • PARI
    a(n) = binomial(n+5,5) * binomial(n+5,4)/(n+5) \\ Charles R Greathouse IV, Jun 11 2015
    
  • PARI
    Vec((1+6*x+6*x^2+x^3)/(1-x)^9 + O(x^99)) \\ Altug Alkan, Sep 01 2016
    
  • SageMath
    def A006857(n): return binomial(n+4,3)*binomial(n+5,5)//4
    print([A006857(n) for n in range(41)]) # G. C. Greubel, Mar 12 2025

Formula

From - Vladeta Jovovic, Jan 29 2003: (Start)
a(n) = (n+4)!*(n+5)!/(2880*n!*(n+1)!).
E.g.f.: 1/2880*(2880 + 40320*x + 109440*x^2 + 105120*x^3 + 45000*x^4 + 9504*x^5 + 1016*x^6 + 52*x^7 + x^8)*exp(x). (End)
From Mike Zabrocki, Aug 26 2004: (Start)
a(n) = C(n+5,8) + 6*C(n+6,8) + 6*C(n+7,8) + C(n+8,8).
a(n) = C(n+4,4)*C(n+5,4)/5.
O.g.f.: (1 + 6*x + 6*x^2 + x^3)/(1-x)^9. (End)
From Wolfdieter Lang, Nov 13 2007: (Start)
a(n) = A001263(n+5,5).
Numerator polynomial of the g.f is the fourth row polynomial of the Narayana triangle. (End)
a(n)= C(n+4,4)^2 - C(n+4,3)*C(n+4,5). - Gary Detlefs, Dec 05 2011
a(n) = Product_{i=1..4} A002378(n+i)/A002378(i). - Bruno Berselli, Sep 01 2016
From Amiram Eldar, Oct 19 2020: (Start)
Sum_{n>=0} 1/a(n) = 25 * (79 - 8*Pi^2).
Sum_{n>=0} (-1)^n/a(n) = 595/3 - 20*Pi^2. (End)

Extensions

More terms from Vladeta Jovovic, Jan 29 2003
Better description from Mike Zabrocki, Aug 26 2004
New definition from N. J. A. Sloane, Aug 28 2010
Zabrocki formulas offset corrected by Gary Detlefs, Dec 05 2011

A132813 Triangle read by rows: A001263 * A127648 as infinite lower triangular matrices.

Original entry on oeis.org

1, 1, 2, 1, 6, 3, 1, 12, 18, 4, 1, 20, 60, 40, 5, 1, 30, 150, 200, 75, 6, 1, 42, 315, 700, 525, 126, 7, 1, 56, 588, 1960, 2450, 1176, 196, 8, 1, 72, 1008, 4704, 8820, 7056, 2352, 288, 9, 1, 90, 1620, 10080, 26460, 31752, 17640, 4320, 405, 10
Offset: 0

Views

Author

Gary W. Adamson, Sep 01 2007

Keywords

Comments

Also T(n,k) = binomial(n-1, k-1)*binomial(n, k-1), related to Narayana polynomials (see Sulanke reference). - Roger L. Bagula, Apr 09 2008
h-vector for cluster complex associated to the root system B_n. See p. 8, Athanasiadis and C. Savvidou. - Tom Copeland, Oct 19 2014

Examples

			First few rows of the triangle are:
  1;
  1,  2;
  1,  6,   3;
  1, 12,  18,   4;
  1, 20,  60,  40,   5;
  1, 30, 150, 200,  75,   6;
  1, 42, 315, 700, 525, 126, 7;
  ...
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), this sequence (c=2), A062196 (c=3), A062145 (c=4), A062264 (c=5), A062190 (c=6).
Columns: A000012 (k=0), A002378 (k=1), A006011 (k=2), 4*A006542 (k=3), 5*A006857 (k=4), 6*A108679 (k=5), 7*A134288 (k=6), 8*A134289 (k=7), 9*A134290 (k=8), 10*A134291 (k=9).
Diagonals: A000027 (k=n), A002411 (k=n-1), A004302 (k=n-2), A108647 (k=n-3), A134287 (k=n-4).
Main diagonal: A000894.
Sums: (-1)^floor((n+1)/2)*A001405 (signed row), A001700 (row), A203611 (diagonal).
Cf. A103371 (mirrored).

Programs

  • GAP
    Flat(List([0..10],n->List([0..n], k->(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1)))); # Muniru A Asiru, Feb 26 2019
    
  • Haskell
    a132813 n k = a132813_tabl !! n !! k
    a132813_row n = a132813_tabl !! n
    a132813_tabl = zipWith (zipWith (*)) a007318_tabl $ tail a007318_tabl
    -- Reinhard Zumkeller, Apr 04 2014
    
  • Magma
    /* triangle */ [[(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Oct 19 2014
    
  • Maple
    P := (n, x) -> hypergeom([1-n, -n], [1], x): for n from 1 to 9 do PolynomialTools:-CoefficientList(simplify(P(n,x)),x) od; # Peter Luschny, Nov 26 2014
  • Mathematica
    T[n_,k_]=Binomial[n-1,k-1]*Binomial[n,k-1]; Table[Table[T[n,k],{k,1,n}],{n,1,11}]; Flatten[%] (* Roger L. Bagula, Apr 09 2008 *)
    P[n_, x_] := HypergeometricPFQ[{1-n, -n}, {1}, x]; Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten (* Jean-François Alcover, Nov 27 2014, after Peter Luschny *)
  • PARI
    tabl(nn) = {for (n = 1, nn, for (k = 1, n, print1(binomial(n-1, k-1)*binomial(n, k-1) , ", ");););} \\ Michel Marcus, Feb 12 2014
    
  • SageMath
    def A132813(n,k): return binomial(n,k)*binomial(n+1,k)
    print(flatten([[A132813(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 12 2025

Formula

T(n,k) = (k+1)*binomial(n+1,k+1)*binomial(n+1,k)/(n+1), n >= k >= 0.
From Roger L. Bagula, May 14 2010: (Start)
T(n, m) = coefficients(p(x,n)), where
p(x,n) = (1-x)^(2*n)*Sum_{k >= 0} binomial(k+n-1, k)*binomial(n+k, k)*x^k,
or p(x,n) = (1-x)^(2*n)*Hypergeometric2F1([n, n+1], [1], x). (End)
T(n,k) = binomial(n,k) * binomial(n+1,k). - Reinhard Zumkeller, Apr 04 2014
These are the coefficients of the polynomials Hypergeometric2F1([1-n,-n], [1], x). - Peter Luschny, Nov 26 2014
G.f.: A(x,y) = A281260(x,y)/(1-A281260(x,y))/x. - Vladimir Kruchinin, Oct 10 2020

A124428 Triangle, read by rows: T(n,k) = binomial(floor(n/2),k)*binomial(floor((n+1)/2),k).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 4, 1, 1, 6, 3, 1, 9, 9, 1, 1, 12, 18, 4, 1, 16, 36, 16, 1, 1, 20, 60, 40, 5, 1, 25, 100, 100, 25, 1, 1, 30, 150, 200, 75, 6, 1, 36, 225, 400, 225, 36, 1, 1, 42, 315, 700, 525, 126, 7, 1, 49, 441, 1225, 1225, 441, 49, 1, 1, 56, 588, 1960, 2450, 1176, 196, 8
Offset: 0

Views

Author

Paul D. Hanna, Oct 31 2006

Keywords

Comments

Row sums form A001405, the central binomial coefficients: C(n,floor(n/2)). The eigenvector of this triangle is A124430.
T(n,k) is the number of dispersed Dyck paths of length n (i.e., Motzkin paths of length n with no (1,0) steps at positive heights) having k peaks. Example: T(5,2)=3 because, denoting U=(1,1), D=(1,-1), H=(1,0), we have HUDUD, UDHUD, and UDUDH. - Emeric Deutsch, Jun 01 2011
From Emeric Deutsch, Jan 18 2013: (Start)
T(n,k) is the number of Dyck prefixes of length n having k peaks. Example: T(5,2)=3 because we have (UD)(UD)U, (UD)U(UD), and U(UD)(UD); the peaks are shown between parentheses.
T(n,k) is the number of Dyck prefixes of length n having k ascents and descents of length >= 2. Example: T(5,2)=3 because we have (UU)(DD)U, (UU)D(UU), and (UUU)(DD); the ascents and descents of length >= 2 are shown between parentheses. (End)
T(n,k) is the number of noncrossing partitions of [n] having n-k blocks, such that the nontrivial blocks are of type {a,b}, with a < = n/2 and b > n/2. Such partitions have k nontrivial blocks, uniquely determined by the choice of k first elements among floor(n/2) elements, and the choice of k second elements among floor((n+1)/2) elements. Indeed, by planarity, any two blocs {a,b} and {c,d} satisfy a < c iff b > d. - Francesca Aicardi Nov 03 2022

Examples

			Triangle begins:
  1;
  1;
  1,   1;
  1,   2;
  1,   4,   1;
  1,   6,   3;
  1,   9,   9,   1;
  1,  12,  18,   4;
  1,  16,  36,  16,   1;
  1,  20,  60,  40,   5;
  1,  25, 100, 100,  25,   1;
  1,  30, 150, 200,  75,   6;
  1,  36, 225, 400, 225,  36,   1; ...
		

Crossrefs

Cf. A001405 (row sums), A056953, A026003, A124429 (antidiagonal sums), A124430 (eigenvector), A191521.
Columns = A002378, A006011, A006542, etc.

Programs

  • Magma
    [[Binomial(Floor(n/2), k)*Binomial(Floor((n+1)/2),k): k in [0..Floor(n/2)]]: n in [0..15]]; // G. C. Greubel, Feb 24 2019
    
  • Mathematica
    Table[Binomial[Floor[n/2], k]*Binomial[Floor[(n+1)/2], k], {n, 0, 15}, {k, 0, Floor[n/2]}]//Flatten (* G. C. Greubel, Feb 24 2019 *)
  • PARI
    T(n,k)=binomial(n\2,k)*binomial((n+1)\2,k)
    
  • Sage
    [[binomial(floor(n/2),k)*binomial(floor((n+1)/2),k) for k in (0..floor(n/2))] for n in (0..15)] # G. C. Greubel, Feb 24 2019

Formula

A056953(n) = Sum_{k=0..floor(n/2)} k!*T(n,k).
A026003(n) = Sum_{k=0..floor(n/2)} 2^k*T(n,k).

A108679 a(n) = (n+1)*(n+2)^2*(n+3)^2*(n+4)^2*(n+5)^2*(n+6)/86400.

Original entry on oeis.org

1, 21, 196, 1176, 5292, 19404, 60984, 169884, 429429, 1002001, 2186184, 4504864, 8836464, 16604784, 30046752, 52581816, 89311761, 147685461, 238369516, 376372920, 582481900, 885069900, 1322357400, 1945206900, 2820550005, 4035556161, 5702666256, 7965629056
Offset: 0

Views

Author

Emeric Deutsch, Jun 17 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.
6th column of the table of Narayana numbers A001263. - Zerinvary Lajos, Jun 18 2007
Sequence provided by binomial(n-1,m)*binomial(n,m)/(m+1) for m=5 and n>5 (these numbers are also called Runyon numbers, see T. Koshy in References). - Vincenzo Librandi, Sep 04 2014

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p. 232, # 1).
  • T. Koshy, Catalan Numbers with Applications, Oxford University Press, 2009, p. 7.
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; Prop. 8.4, case n=7. - N. J. A. Sloane, Aug 28 2010

Crossrefs

Cf. A001263, A002378, A005585, A006542, A006857 (sequences having a similar structure), A288876.

Programs

  • Magma
    [Binomial(n-1,5)*Binomial(n,5)/6: n in [6..45]]; // Vincenzo Librandi, Sep 04 2014
    
  • Maple
    a:=n->(n+1)*(n+2)^2*(n+3)^2*(n+4)^2*(n+5)^2*(n+6)/86400: seq(a(n),n=0..40);
  • Mathematica
    Table[(n + 1) (n + 2)^2 (n + 3)^2 (n + 4)^2 (n + 5)^2 (n + 6)/86400,{n, 0, 50}]  (* Harvey P. Dale, Mar 13 2011 *)
  • PARI
    Vec((1+10*x+20*x^2+10*x^3+x^4)/(1-x)^11 + O(x^99)) \\ Altug Alkan, Sep 02 2016
    
  • SageMath
    def A108679(n): return binomial(n+5,5)*binomial(n+6,5)//6
    print([A108679(n) for n in range(41)]) # G. C. Greubel, Mar 12 2025

Formula

a(n) = binomial(n+5, 5)*binomial(n+6, 5)/6 = binomial(n+6, 6)*binomial(n+6, 5)/(n+6).
a(n) = A001263(n+6,6).
G.f.: (1 + 10*x + 20*x^2 + 10*x^3 + x^4)/(1 - x)^11. Numerator polynomial is the fifth row polynomial of the Narayana triangle.
a(n) = binomial(n+5,5)^2 - binomial(n+5,4)*binomial(n+5,6). - Gary Detlefs, Dec 05 2011
a(n) = Product_{i=1..5} A002378(n+i)/A002378(i). - Bruno Berselli, Sep 01 2016
From Amiram Eldar, Oct 19 2020: (Start)
Sum_{n>=0} 1/a(n) = 27637/2 - 1400*Pi^2.
Sum_{n>=0} (-1)^n/a(n) = 2560*log(2) - 3547/2. (End)
a(n) = (A005585(n+2)^2 - A288876(n+1))/24. - Yasser Arath Chavez Reyes, Aug 19 2024

A131198 Triangle T(n,k), 0 <= k <= n, read by rows, given by [1,0,1,0,1,0,1,0,...] DELTA [0,1,0,1,0,1,0,1,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 6, 1, 0, 1, 10, 20, 10, 1, 0, 1, 15, 50, 50, 15, 1, 0, 1, 21, 105, 175, 105, 21, 1, 0, 1, 28, 196, 490, 490, 196, 28, 1, 0, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1, 0, 1, 45, 540, 2520, 5292, 5292, 2520, 540, 45, 1, 0
Offset: 0

Views

Author

Philippe Deléham, Oct 20 2007

Keywords

Comments

Mirror image of triangle A090181, another version of triangle of Narayana (A001263).
Equals A133336*A130595 as infinite lower triangular matrices. - Philippe Deléham, Oct 23 2007

Examples

			Triangle begins:
  1;
  1,  0;
  1,  1,   0;
  1,  3,   1,   0;
  1,  6,   6,   1,   0;
  1, 10,  20,  10,   1,   0;
  1, 15,  50,  50,  15,   1,  0;
  1, 21, 105, 175, 105,  21,  1, 0;
  1, 28, 196, 490, 490, 196, 28, 1, 0; ...
		

Crossrefs

Programs

  • Magma
    [[n le 0 select 1 else (n-k)*Binomial(n,k)^2/(n*(k+1)): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Feb 06 2018
  • Maple
    T := (n,k) -> `if`(n=0, 0^n, binomial(n,k)^2*(n-k)/(n*(k+1)));
    seq(print(seq(T(n,k), k=0..n)), n=0..5); # Peter Luschny, Jun 08 2014
    R := n -> simplify(hypergeom([1 - n, -n], [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[If[n == 0, 1, (n-k)*Binomial[n,k]^2/(n*(k+1))], {n,0,10}, {k,0,n}] //Flatten (* G. C. Greubel, Feb 06 2018 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(if(n==0,1, (n-k)*binomial(n,k)^2/(n* (k+1))), ", "))) \\ G. C. Greubel, Feb 06 2018
    

Formula

Sum_{k=0..n} T(n,k)*x^k = A000012(n), A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 respectively.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000007(n), A000108(n), A006318(n), A047891(n+1), A082298(n), A082301(n), A082302(n), A082305(n), A082366(n), A082367(n), for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 23 2007
Sum_{k=0..floor(n/2)} T(n-k,k) = A004148(n). - Philippe Deléham, Nov 06 2007
T(2*n,n) = A125558(n). - Philippe Deléham, Nov 16 2011
T(n, k) = [x^k] hypergeom([1 - n, -n], [2], x). - Peter Luschny, Apr 26 2022

A107891 a(n) = (n+1)*(n+2)^2*(n+3)^2*(n+4)*(3n^2 + 15n + 20)/2880.

Original entry on oeis.org

1, 19, 155, 805, 3136, 9996, 27468, 67320, 150645, 313027, 611611, 1134497, 2012920, 3436720, 5673648, 9093096, 14194881, 21643755, 32310355, 47319349, 68105576, 96479020, 134699500, 185562000, 252493605, 339663051, 452103939
Offset: 0

Views

Author

Emeric Deutsch, Jun 12 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.
Partial sums of A114239. First differences of A047819. - Peter Bala, Sep 21 2007

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see pp. 167, 187 and p. 105 eq. (iii) for k=2 and m=5).

Crossrefs

Programs

  • Maple
    a:=n->(1/2880)*(n+1)*(n+2)^2*(n+3)^2*(n+4)*(3*n^2+15*n+20): seq(a(n),n=0..32);
  • Mathematica
    Table[((1+n) (2+n)^2 (3+n)^2 (4+n) (20+3 n (5+n)))/2880,{n,0,40}] (* or *) LinearRecurrence[{9,-36,84,-126,126,-84,36,-9,1},{1,19,155,805,3136,9996,27468,67320,150645},40] (* Harvey P. Dale, Dec 10 2021 *)

Formula

a(n-2) = (1/8) * Sum_{1 <= x_1, x_2 <= n} (x_1*x_2)^2*(det V(x_1,x_2))^2 = 1/8*sum {1 <= i,j <= n} (i*j*(i-j))^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
G.f.: (1+10*x+20*x^2+10*x^3+x^4)/(1-x)^9. - Colin Barker, Feb 08 2012
a(n) = (A000330(n+2)*A000538(n+2) - (A000537(n+2))^2)/4. - J. M. Bergot, Sep 17 2013
Sum_{n>=0} 1/a(n) = 17095/4 - 240*Pi^2 - 162*sqrt(15)*Pi*tanh(sqrt(5/3)*Pi/2). - Amiram Eldar, May 29 2022

A114242 a(n) = (n+1)(n+2)^2*(n+3)^2*(n+4)(2n+5)/720.

Original entry on oeis.org

1, 14, 90, 385, 1274, 3528, 8568, 18810, 38115, 72358, 130130, 223587, 369460, 590240, 915552, 1383732, 2043621, 2956590, 4198810, 5863781, 8065134, 10939720, 14651000, 19392750, 25393095, 32918886, 42280434, 53836615, 68000360
Offset: 0

Views

Author

Emeric Deutsch, Nov 18 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.
Partial sums of A114244. First differences of A006857. - Peter Bala, Sep 21 2007

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (pp. 167-169, Table 10.5/II/2 and p. 105, eq. (ii) K(Ob(2,4,n))).

Crossrefs

Programs

  • Maple
    a:=n->(n+1)*(n+2)^2*(n+3)^2*(n+4)*(2*n+5)/720: seq(a(n),n=0..30);
  • Mathematica
    Table[((n+1)(n+2)^2 (n+3)^2 (n+4)(2n+5))/720,{n,0,30}] (* or *) LinearRecurrence[ {8,-28,56,-70,56,-28,8,-1},{1,14,90,385,1274,3528,8568,18810},30] (* Harvey P. Dale, Aug 21 2013 *)
  • PARI
    a(n)=(n+1)*(n+2)^2*(n+3)^2*(n+4)*(2*n+5)/720 \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: (1+x)(1 + 5x + x^2)/(1-x)^8.
a(n-2) = (1/6) * Sum_{1 <= x_1, x_2 <= n} (x_1)^2*x_2*(det V(x_1,x_2))^2 = (1/6)*Sum_{1 <= i,j <= n} i^2*j*(i-j)^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
a(n) = 8*a(n-1) - 28*a(n-2) + 56*a(n-3) - 70*a(n-4) + 56*a(n-5) - 28*a(n-6) + 8*a(n-7) - a(n-8). - Harvey P. Dale, Aug 21 2013
From Amiram Eldar, May 29 2022: (Start)
Sum_{n>=0} 1/a(n) = 3550 - 5120*log(2).
Sum_{n>=0} (-1)^n/a(n) = 3430 - 1280*Pi + 60*Pi^2. (End)

A359363 Triangle read by rows. The coefficients of the Baxter polynomials p(0, x) = 1 and p(n, x) = x*hypergeom([-1 - n, -n, 1 - n], [2, 3], -x) for n >= 1.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 10, 10, 1, 0, 1, 20, 50, 20, 1, 0, 1, 35, 175, 175, 35, 1, 0, 1, 56, 490, 980, 490, 56, 1, 0, 1, 84, 1176, 4116, 4116, 1176, 84, 1, 0, 1, 120, 2520, 14112, 24696, 14112, 2520, 120, 1, 0, 1, 165, 4950, 41580, 116424, 116424, 41580, 4950, 165, 1
Offset: 0

Views

Author

Peter Luschny, Dec 28 2022

Keywords

Comments

This triangle is a member of a family of Pascal-like triangles. Let T(n, k, m) = sf(m)*F(n - 1) / (F(k - 1)*F(n - k)) if k > 0 and otherwise k^n, where F(n) = Product_{j=0..m} (n + j)! and sf(m) are the superfactorials A000178. The case m = 2 gives this triangle, some other cases are given in the crossreferences. See also A342889 for a related representation of generalized binomial coefficients.

Examples

			Triangle T(n, k) starts:
[0] 1
[1] 0, 1
[2] 0, 1,   1
[3] 0, 1,   4,    1
[4] 0, 1,  10,   10,     1
[5] 0, 1,  20,   50,    20,     1
[6] 0, 1,  35,  175,   175,    35,     1
[7] 0, 1,  56,  490,   980,   490,    56,    1
[8] 0, 1,  84, 1176,  4116,  4116,  1176,   84,   1
[9] 0, 1, 120, 2520, 14112, 24696, 14112, 2520, 120, 1
.
Let p = (p1, p2,..., pn) denote a permutation of {1, 2,..., n}. The pair (p(i), p(i+1)) is a 'rise' if p(i) < p(i+1). Additionally a conventional rise is counted at the beginning of p.
T(n, k) is the number of Baxter permutations of {1,2,...,n} with k rises. For example for n = 4, [T(n, k) for k = 0..n] = [0, 1, 10, 10, 1]. The permutations, with preceding number of rises, are:
.
1 [4, 3, 2, 1],  3 [2, 3, 4, 1],  2 [3, 4, 2, 1],  3 [2, 3, 1, 4],
2 [3, 2, 4, 1],  3 [2, 1, 3, 4],  2 [3, 2, 1, 4],  3 [1, 3, 4, 2],
2 [2, 4, 3, 1],  3 [1, 3, 2, 4],  2 [4, 2, 3, 1],  3 [3, 4, 1, 2],
2 [2, 1, 4, 3],  3 [3, 1, 2, 4],  2 [4, 2, 1, 3],  3 [1, 2, 4, 3],
2 [1, 4, 3, 2],  3 [1, 4, 2, 3],  2 [4, 1, 3, 2],  3 [4, 1, 2, 3],
2 [4, 3, 1, 2],  4 [1, 2, 3, 4].
		

Crossrefs

Special cases of the general formula: A097805 (m = 0), (0,1)-Pascal triangle; A090181 (m = 1), triangle of Narayana; this triangle (m = 2); A056940 (m = 3), with 1,0,0...; A056941 (m = 4), with 1,0,0...; A142465 (m = 5), with 1,0,0....
Variant: A056939. Diagonals: A000292, A006542, A047819.

Programs

  • Maple
    p := (n, x) -> ifelse(n = 0, 1, x*hypergeom([-1-n, -n, 1-n], [2, 3], -x)):
    seq(seq(coeff(simplify(p(n, x)), x, k), k = 0..n), n = 0..10);
    # Alternative:
    T := proc(n, k) local F; F := n -> n!*(n+1)!*(n+2)!;
    ifelse(k = 0, k^n, 2*F(n-1)/(F(k-1)*F(n-k))) end:
    for n from 0 to 9 do seq(T(n, k), k = 0..n) od;
  • PARI
    C=binomial;
    T(n, k) = if(n==0 && k==0, 1, ( C(n+1,k-1) * C(n+1,k) * C(n+1,k+1) ) / ( C(n+1,1) * C(n+1,2) ) );
    for(n=0,10,for(k=0,n,print1(T(n,k),", "));print());
    \\ Joerg Arndt, Jan 04 2024
    
  • Python
    from functools import cache
    from math import factorial
    @cache
    def A359363Row(n: int) -> list[int]:
        @cache
        def F(n: int): return factorial(n) ** 3 * ((n+1) * (n+1) * (n+2))
        if n == 0: return [1]
        return [0] + [(2*F(n-1))//(F(k-1) * F(n-k)) for k in range(1, n+1)]
    for n in range(0, 10): print(A359363Row(n))
    # Peter Luschny, Jan 04 2024
  • SageMath
    def A359363(n):
        if n == 0: return SR(1)
        h = x*hypergeometric([-1 - n, -n, 1 - n], [2, 3], -x)
        return h.series(x, n + 1).polynomial(SR)
    for n in range(10): print(A359363(n).list())
    def PolyA359363(n, t): return Integer(A359363(n)(x=t).n())
    # Peter Luschny, Jan 04 2024
    

Formula

T(n, k) = [x^k] p(n, x).
T(n, k) = 2*F(n-1)/(F(k-1)*F(n-k)) for k > 0 where F(n) = n!*(n+1)!*(n+2)!.
p(n, 1) = A001181(n), i.e. the Baxter numbers are the values of the Baxter polynomials at x = 1.
(-1)^(n + 1)*p(2*n + 1, -1) = A217800(n) .

A123095 Sum of first n 11th powers.

Original entry on oeis.org

0, 1, 2049, 179196, 4373500, 53201625, 415998681, 2393325424, 10983260016, 42364319625, 142364319625, 427675990236, 1170684360924, 2962844754961, 7012409924625, 15662165784000, 33254351828416, 67526248136049, 131794658215281, 248284917113500, 453084917113500
Offset: 0

Views

Author

Zerinvary Lajos, Sep 27 2006

Keywords

Crossrefs

Sequences of the form Sum_{j=0..n} j^m : A000217 (m=1), A000330 (m=2), A000537 (m=3), A000538 (m=4), A000539 (m=5), A000540 (m=6), A000541 (m=7), A000542 (m=8), A007487 (m=9), A023002 (m=10), this sequence (m=11), A123094 (m=12), A181134 (m=13).

Programs

  • Magma
    [(&+[j^11: j in [0..n]]): n in [0..30]]; // G. C. Greubel, Jul 21 2021
    
  • Maple
    [seq(add(i^11, i=1..n), n=0..20)];
    a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=a[n-1]+n^11 od: seq(a[n], n=0..13); # Zerinvary Lajos, Feb 22 2008
  • Mathematica
    Table[Sum[k^11, {k, n}], {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2008 *)
    Accumulate[Range[0,20]^11] (* Harvey P. Dale, Sep 17 2021 *)
  • Python
    A123095_list, m = [0], [39916800, -199584000, 419126400, -479001600, 322494480, -129230640, 29607600, -3498000, 171006, -2046, 1, 0 , 0]
    for _ in range(10**2):
        for i in range(12):
            m[i+1]+= m[i]
        A123095_list.append(m[-1]) # Chai Wah Wu, Nov 05 2014
    
  • Sage
    [(bernoulli_polynomial(n+1, 12) - bernoulli(12))/12  for n in (0..30)] # G. C. Greubel, Jul 21 2021

Formula

a(n) = n*A023002(n) - Sum_{i=0..n-1} A023002(i). - Bruno Berselli, Apr 27 2010
a(n) = n^2*(n+1)^2*(2*n^8 +8*n^7 +4*n^6 -16*n^5 -5*n^4 +26*n^3 -3*n^2 -20*n +10)/24. - Bruno Berselli, Oct 03 2010
G.f.: x*(x^10 +2036*x^9 +152637*x^8 +2203488*x^7 +9738114*x^6 +15724248*x^5 +9738114*x^4 +2203488*x^3 +152637*x^2 +2036*x +1)/(1-x)^13. - Colin Barker, May 27 2012
a(n) = (-1)*Sum_{j=1..11} j*Stirling1(n+1,n+1-j)*Stirling2(n+11-j,n). - Mircea Merca, Jan 25 2014
a(n) = 1728*A006542(n+2)^2 + 216*A288876(n-2) + 96*A006542(n+2) + A000537(n). - Yasser Arath Chavez Reyes, May 25 2024
Previous Showing 11-20 of 29 results. Next