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-10 of 47 results. Next

A000096 a(n) = n*(n+3)/2.

Original entry on oeis.org

0, 2, 5, 9, 14, 20, 27, 35, 44, 54, 65, 77, 90, 104, 119, 135, 152, 170, 189, 209, 230, 252, 275, 299, 324, 350, 377, 405, 434, 464, 495, 527, 560, 594, 629, 665, 702, 740, 779, 819, 860, 902, 945, 989, 1034, 1080, 1127, 1175, 1224, 1274, 1325, 1377, 1430, 1484, 1539, 1595, 1652, 1710, 1769
Offset: 0

Views

Author

Keywords

Comments

For n >= 1, a(n) is the maximal number of pieces that can be obtained by cutting an annulus with n cuts. See illustration. - Robert G. Wilson v
n(n-3)/2 (n >= 3) is the number of diagonals of an n-gon. - Antreas P. Hatzipolakis (xpolakis(AT)otenet.gr)
n(n-3)/2 (n >= 4) is the degree of the third-smallest irreducible presentation of the symmetric group S_n (cf. James and Kerber, Appendix 1).
a(n) is also the multiplicity of the eigenvalue (-2) of the triangle graph Delta(n+1). (See p. 19 in Biggs.) - Felix Goldberg (felixg(AT)tx.technion.ac.il), Nov 25 2001
For n > 3, a(n-3) = dimension of the traveling salesman polytope T(n). - Benoit Cloitre, Aug 18 2002
Also counts quasi-dominoes (quasi-2-ominoes) on an n X n board. Cf. A094170-A094172. - Jon Wild, May 07 2004
Coefficient of x^2 in (1 + x + 2*x^2)^n. - Michael Somos, May 26 2004
a(n) is the number of "prime" n-dimensional polyominoes. A "prime" n-polyomino cannot be formed by connecting any other n-polyominoes except for the n-monomino and the n-monomino is not prime. E.g., for n=1, the 1-monomino is the line of length 1 and the only "prime" 1-polyominoes are the lines of length 2 and 3. This refers to "free" n-dimensional polyominoes, i.e., that can be rotated along any axis. - Bryan Jacobs (bryanjj(AT)gmail.com), Apr 01 2005
Solutions to the quadratic equation q(m, r) = (-3 +- sqrt(9 + 8(m - r))) / 2, where m - r is included in a(n). Let t(m) = the triangular number (A000217) less than some number k and r = k - t(m). If k is neither prime nor a power of two and m - r is included in A000096, then m - q(m, r) will produce a value that shares a divisor with k. - Andrew S. Plewe, Jun 18 2005
Sum_{k=2..n+1} 4/(k*(k+1)*(k-1)) = ((n+3)*n)/((n+2)*(n+1)). Numerator(Sum_{k=2..n+1} 4/(k*(k+1)*(k-1))) = (n+3)*n/2. - Alexander Adamchuk, Apr 11 2006
Number of rooted trees with n+3 nodes of valence 1, no nodes of valence 2 and exactly two other nodes. I.e., number of planted trees with n+2 leaves and exactly two branch points. - Theo Johnson-Freyd (theojf(AT)berkeley.edu), Jun 10 2007
If X is an n-set and Y a fixed 2-subset of X then a(n-2) is equal to the number of (n-2)-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
For n >= 1, a(n) is the number of distinct shuffles of the identity permutation on n+1 letters with the identity permutation on 2 letters (12). - Camillia Smith Barnes, Oct 04 2008
If s(n) is a sequence defined as s(1) = x, s(n) = kn + s(n-1) + p for n > 1, then s(n) = a(n-1)*k + (n-1)*p + x. - Gary Detlefs, Mar 04 2010
The only primes are a(1) = 2 and a(2) = 5. - Reinhard Zumkeller, Jul 18 2011
a(n) = m such that the (m+1)-th triangular number minus the m-th triangular number is the (n+1)-th triangular number: (m+1)(m+2)/2 - m(m+1)/2 = (n+1)(n+2)/2. - Zak Seidov, Jan 22 2012
For n >= 1, number of different values that Sum_{k=1..n} c(k)*k can take where the c(k) are 0 or 1. - Joerg Arndt, Jun 24 2012
On an n X n chessboard (n >= 2), the number of possible checkmate positions in the case of king and rook versus a lone king is 0, 16, 40, 72, 112, 160, 216, 280, 352, ..., which is 8*a(n-2). For a 4 X 4 board the number is 40. The number of positions possible was counted including all mirror images and rotations for all four sides of the board. - Jose Abutal, Nov 19 2013
If k = a(i-1) or k = a(i+1) and n = k + a(i), then C(n, k-1), C(n, k), C(n, k+1) are three consecutive binomial coefficients in arithmetic progression and these are all the solutions. There are no four consecutive binomial coefficients in arithmetic progression. - Michael Somos, Nov 11 2015
a(n-1) is also the number of independent components of a symmetric traceless tensor of rank 2 and dimension n >= 1. - Wolfdieter Lang, Dec 10 2015
Numbers k such that 8k + 9 is a square. - Juri-Stepan Gerasimov, Apr 05 2016
Let phi_(D,rho) be the average value of a generic degree D monic polynomial f when evaluated at the roots of the rho-th derivative of f, expressed as a polynomial in the averaged symmetric polynomials in the roots of f. [See the Wojnar et al. link] The "last" term of phi_(D,rho) is a multiple of the product of all roots of f; the coefficient is expressible as a polynomial h_D(N) in N:=D-rho. These polynomials are of the form h_D(N)= ((-1)^D/(D-1)!)*(D-N)*N^chi*g_D(N) where chi = (1 if D is odd, 0 if D is even) and g_D(N) is a monic polynomial of degree (D-2-chi). Then a(n) are the negated coefficients of the next to the highest order term in the polynomials N^chi*g_D(N), starting at D=3. - Gregory Gerard Wojnar, Jul 19 2017
For n >= 2, a(n) is the number of summations required to solve the linear regression of n variables (n-1 independent variables and 1 dependent variable). - Felipe Pedraza-Oropeza, Dec 07 2017
For n >= 2, a(n) is the number of sums required to solve the linear regression of n variables: 5 for two variables (sums of X, Y, X^2, Y^2, X*Y), 9 for 3 variables (sums of X1, X2, Y1, X1^2, X1*X2, X1*Y, X2^2, X2*Y, Y^2), and so on. - Felipe Pedraza-Oropeza, Jan 11 2018
a(n) is the area of a triangle with vertices at (n, n+1), ((n+1)*(n+2)/2, (n+2)*(n+3)/2), ((n+2)^2, (n+3)^2). - J. M. Bergot, Jan 25 2018
Number of terms less than 10^k: 1, 4, 13, 44, 140, 446, 1413, 4471, 14141, 44720, 141420, 447213, ... - Muniru A Asiru, Jan 25 2018
a(n) is also the number of irredundant sets in the (n+1)-path complement graph for n > 2. - Eric W. Weisstein, Apr 11 2018
a(n) is also the largest number k such that the largest Dyck path of the symmetric representation of sigma(k) has exactly n peaks, n >= 1. (Cf. A237593.) - Omar E. Pol, Sep 04 2018
For n > 0, a(n) is the number of facets of associahedra. Cf. A033282 and A126216 and their refinements A111785 and A133437 for related combinatorial and analytic constructs. See p. 40 of Hanson and Sha for a relation to projective spaces and string theory. - Tom Copeland, Jan 03 2021
For n > 0, a(n) is the number of bipartite graphs with 2n or 2n+1 edges, no isolated vertices, and a stable set of cardinality 2. - Christian Barrientos, Jun 13 2022
For n >= 2, a(n-2) is the number of permutations in S_n which are the product of two different transpositions of adjacent points. - Zbigniew Wojciechowski, Mar 31 2023
a(n) represents the optimal stop-number to achieve the highest running score for the Greedy Pig game with an (n-1)-sided die with a loss on a 1. The total at which one should stop is a(s-1), e.g. for a 6-sided die, one should pass the die at 20. See Sparks and Haran. - Nicholas Stefan Georgescu, Jun 09 2024

Examples

			G.f. = 2*x + 5*x^2 + 9*x^3 + 14*x^4 + 20*x^5 + 27*x^6 + 35*x^7 + 44*x^8 + 54*x^9 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), Table 22.7, p. 797.
  • Norman Biggs, Algebraic Graph Theory, 2nd ed. Cambridge University Press, 1993.
  • G. James and A. Kerber, The Representation Theory of the Symmetric Group, Encyclopedia of Maths. and its Appls., Vol. 16, Addison-Wesley, 1981, Reading, MA, U.S.A.
  • D. G. Kendall et al., Shape and Shape Theory, Wiley, 1999; see p. 4.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Complement of A007401. Column 2 of A145324. Column of triangle A014473, first skew subdiagonal of A033282, a diagonal of A079508.
Occurs as a diagonal in A074079/A074080, i.e., A074079(n+3, n) = A000096(n-1) for all n >= 2. Also A074092(n) = 2^n * A000096(n-1) after n >= 2.
Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488.
Similar sequences are listed in A316466.

Programs

Formula

G.f.: A(x) = x*(2-x)/(1-x)^3. a(n) = binomial(n+1, n-1) + binomial(n, n-1).
Connection with triangular numbers: a(n) = A000217(n+1) - 1.
a(n) = a(n-1) + n + 1. - Bryan Jacobs (bryanjj(AT)gmail.com), Apr 01 2005
a(n) = 2*t(n) - t(n-1) where t() are the triangular numbers, e.g., a(5) = 2*t(5) - t(4) = 2*15 - 10 = 20. - Jon Perry, Jul 23 2003
a(-3-n) = a(n). - Michael Somos, May 26 2004
2*a(n) = A008778(n) - A105163(n). - Creighton Dement, Apr 15 2005
a(n) = C(3+n, 2) - C(3+n, 1). - Zerinvary Lajos, Dec 09 2005
a(n) = A067550(n+1) / A067550(n). - Alexander Adamchuk, May 20 2006
a(n) = A126890(n,1) for n > 0. - Reinhard Zumkeller, Dec 30 2006
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Paul Curtz, Jan 02 2008
Starting (2, 5, 9, 14, ...) = binomial transform of (2, 3, 1, 0, 0, 0, ...). - Gary W. Adamson, Jul 03 2008
For n >= 0, a(n+2) = b(n+1) - b(n), where b(n) is the sequence A005586. - K.V.Iyer, Apr 27 2009
A002262(a(n)) = n. - Reinhard Zumkeller, May 20 2009
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=1, a(n-1)=coeff(charpoly(A,x),x^(n-2)). - Milan Janjic, Jul 08 2010
a(n) = Sum_{k=1..n} (k+1)!/k!. - Gary Detlefs, Aug 03 2010
a(n) = n(n+1)/2 + n = A000217(n) + n. - Zak Seidov, Jan 22 2012
E.g.f.: F(x) = 1/2*x*exp(x)*(x+4) satisfies the differential equation F''(x) - 2*F'(x) + F(x) = exp(x). - Peter Bala, Mar 14 2012
a(n) = binomial(n+3, 2) - (n+3). - Robert G. Wilson v, Mar 15 2012
a(n) = A181971(n+1, 2) for n > 0. - Reinhard Zumkeller, Jul 09 2012
a(n) = A214292(n+2, 1). - Reinhard Zumkeller, Jul 12 2012
G.f.: -U(0) where U(k) = 1 - 1/((1-x)^2 - x*(1-x)^4/(x*(1-x)^2 - 1/U(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Sep 27 2012
A023532(a(n)) = 0. - Reinhard Zumkeller, Dec 04 2012
a(n) = A014132(n,n) for n > 0. - Reinhard Zumkeller, Dec 12 2012
a(n-1) = (1/n!)*Sum_{j=0..n} binomial(n,j)*(-1)^(n-j)*j^n*(j-1). - Vladimir Kruchinin, Jun 06 2013
a(n) = 2n - floor(n/2) + floor(n^2/2). - Wesley Ivan Hurt, Jun 15 2013
a(n) = Sum_{i=2..n+1} i. - Wesley Ivan Hurt, Jun 28 2013
Sum_{n>0} 1/a(n) = 11/9. - Enrique Pérez Herrero, Nov 26 2013
a(n) = Sum_{i=1..n} (n - i + 2). - Wesley Ivan Hurt, Mar 31 2014
A023531(a(n)) = 1. - Reinhard Zumkeller, Feb 14 2015
For n > 0: a(n) = A101881(2*n-1). - Reinhard Zumkeller, Feb 20 2015
a(n) + a(n-1) = A008865(n+1) for all n in Z. - Michael Somos, Nov 11 2015
a(n+1) = A127672(4+n, n), n >= 0, where A127672 gives the coefficients of the Chebyshev C polynomials. See the Abramowitz-Stegun reference. - Wolfdieter Lang, Dec 10 2015
a(n) = (n+1)^2 - A000124(n). - Anton Zakharov, Jun 29 2016
Dirichlet g.f.: (zeta(s-2) + 3*zeta(s-1))/2. - Ilya Gutkovskiy, Jun 30 2016
a(n) = 2*A000290(n+3) - 3*A000217(n+3). - J. M. Bergot, Apr 04 2018
a(n) = Stirling2(n+2, n+1) - 1. - Peter Luschny, Jan 05 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/3 - 5/9. - Amiram Eldar, Jan 10 2021
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = 3.
Product_{n>=1} (1 - 1/a(n)) = 3*cos(sqrt(17)*Pi/2)/(4*Pi). (End)
Product_{n>=0} a(4*n+1)*a(4*n+4)/(a(4*n+2)*a(4*n+3)) = Pi/6. - Michael Jodl, Apr 05 2025

A000346 a(n) = 2^(2*n+1) - binomial(2*n+1, n+1).

Original entry on oeis.org

1, 5, 22, 93, 386, 1586, 6476, 26333, 106762, 431910, 1744436, 7036530, 28354132, 114159428, 459312152, 1846943453, 7423131482, 29822170718, 119766321572, 480832549478, 1929894318332, 7744043540348, 31067656725032, 124613686513778, 499744650202436
Offset: 0

Views

Author

Keywords

Comments

Also a(n) = 2nd elementary symmetric function of binomial(n,0), binomial(n,1), ..., binomial(n,n).
Also a(n) = one half the sum of the heights, over all Dyck (n+2)-paths, of the vertices that are at even height and terminate an upstep. For example with n=1, these vertices are indicated by asterisks in the 5 Dyck 3-paths: UU*UDDD, UU*DU*DD, UDUU*DD, UDUDUD, UU*DDUD, yielding a(1)=(2+4+2+0+2)/2=5. - David Callan, Jul 14 2006
Hankel transform is (-1)^n*(2n+1); the Hankel transform of sum(k=0..n, C(2*n,k) ) - C(2n,n) is (-1)^n*n. - Paul Barry, Jan 21 2007
Row sums of the Riordan matrix (1/(1-4x),(1-sqrt(1-4x))/2) (A187926). - Emanuele Munarini, Mar 16 2011
From Gus Wiseman, Jul 19 2021: (Start)
For n > 0, a(n-1) is also the number of integer compositions of 2n with nonzero alternating sum, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. These compositions are ranked by A053754 /\ A345921. For example, the a(3-1) = 22 compositions of 6 are:
(6) (1,5) (1,1,4) (1,1,1,3) (1,1,1,1,2)
(2,4) (1,2,3) (1,1,3,1) (1,1,2,1,1)
(4,2) (1,4,1) (1,2,1,2) (2,1,1,1,1)
(5,1) (2,1,3) (1,3,1,1)
(2,2,2) (2,1,2,1)
(3,1,2) (3,1,1,1)
(3,2,1)
(4,1,1)
(End)

Examples

			G.f. = 1 + 5*x + 22*x^2 + 93*x^3 + 386*x^4 + 1586*x^5 + 6476*x^6 + ...
		

References

  • T. Myers and L. Shapiro, Some applications of the sequence 1, 5, 22, 93, 386, ... to Dyck paths and ordered trees, Congressus Numerant., 204 (2010), 93-104.
  • D. Phulara and L. W. Shapiro, Descendants in ordered trees with a marked vertex, Congressus Numerantium, 205 (2011), 121-128.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000108, A014137, A014318. A column of A058893. Subdiagonal of A053979.
Bisection of A058622 and (possibly) A007008.
Even bisection of A294175 (without the first two terms).
The following relate to compositions of 2n with alternating sum k.
- The k > 0 case is counted by A000302.
- The k <= 0 case is counted by A000302.
- The k != 0 case is counted by A000346 (this sequence).
- The k = 0 case is counted by A001700 or A088218.
- The k < 0 case is counted by A008549.
- The k >= 0 case is counted by A114121.
A011782 counts compositions.
A086543 counts partitions with nonzero alternating sum (bisection: A182616).
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A345197 counts compositions by length and alternating sum.

Programs

  • Magma
    [2^(2*n+1) - Binomial(2*n+1,n+1): n in [0..30]]; // Vincenzo Librandi, Jun 07 2011
  • Maple
    seq(2^(2*n+1)-binomial(2*n,n)*(2*n+1)/(n+1), n=0..12); # Emanuele Munarini, Mar 16 2011
  • Mathematica
    Table[2^(2n+1)-Binomial[2n,n](2n+1)/(n+1),{n,0,20}] (* Emanuele Munarini, Mar 16 2011 *)
    a[ n_] := If[ n<-4, 0, (4^(n + 1) - Binomial[2 n + 2, n + 1]) / 2]; (* Michael Somos, Jan 25 2014 *)
  • Maxima
    makelist(2^(2*n+1)-binomial(2*n,n)*(2*n+1)/(n+1),n,0,12); /* Emanuele Munarini, Mar 16 2011 */
    
  • PARI
    {a(n) = if( n<-4, 0, n++; (2^(2*n) - binomial(2*n, n)) / 2)}; /* Michael Somos, Jan 25 2014 */
    

Formula

G.f.: c(x)/(1-4x), c(x) = g.f. of Catalan numbers.
Convolution of Catalan numbers and powers of 4.
Also one half of convolution of central binomial coeffs. A000984(n), n=0, 1, 2, ... with shifted central binomial coeffs. A000984(n), n=1, 2, 3, ...
a(n) = A045621(2n+1) = (1/2)*A068551(n+1).
a(n) = Sum_{k=0..n} A000984(k)*A001700(n-k). - Philippe Deléham, Jan 22 2004
a(n) = Sum_{k=0..n+1} binomial(n+k, k-1)2^(n-k+1). - Paul Barry, Nov 13 2004
a(n) = Sum_{i=0..n} binomial(2n+2, i). See A008949. - Ed Catmur (ed(AT)catmur.co.uk), Dec 09 2006
a(n) = Sum_{k=0..n+1, C(2n+2,k)} - C(2n+2,n+1). - Paul Barry, Jan 21 2007
Logarithm g.f. log(1/(2-C(x)))=sum(n>0, a(n)/n*x^n), C(x)=(1-sqrt(1-4*x))/2x (A000108). - Vladimir Kruchinin, Aug 10 2010
D-finite with recurrence: (n+3) a(n+2) - 2(4n+9) a(n+1) + 8(2n+3) a(n) = 0. - Emanuele Munarini, Mar 16 2011
E.g.f.:exp(2*x)*(2*exp(2*x) - BesselI(0,2*x) - BesselI(1,2*x)).
This is the first derivative of exp(2*x)*(exp(2*x) - BesselI(0,2*x))/2. See the e.g.f. of A032443 (which has a plus sign) and the remarks given there. - Wolfdieter Lang, Jan 16 2012
a(n) = Sum_{0<=iMircea Merca, Apr 05 2012
A000346 = A004171 - A001700. See A032443 for the sum. - M. F. Hasler, Jan 02 2014
0 = a(n) * (256*a(n+1) - 224*a(n+2) + 40*a(n+3)) + a(n+1) * (-32*a(n+1) + 56*a(n+2) - 14*a(n+3)) + a(n+2) * (-2*a(n+2) + a(n+3)) if n>-5. - Michael Somos, Jan 25 2014
REVERT transform is [1,-5,28,-168,1056,...] = alternating signed version of A069731. - Michael Somos, Jan 25 2014
Convolution square is A038806. - Michael Somos, Jan 25 2014
BINOMIAL transform of A055217(n-1) is a(n-1). - Michael Somos, Jan 25 2014
(n+1) * a(n) = A033504(n). - Michael Somos, Jan 25 2014
Recurrence: (n+1)*a(n) = 512*(2*n-7)*a(n-5) + 256*(13-5*n)*a(n-4) + 64*(10*n-17)*a(n-3) + 32*(4-5*n)*a(n-2) + 2*(10*n+1)*a(n-1), n>=5. - Fung Lam, Mar 21 2014
Asymptotic approximation: a(n) ~ 2^(2n+1)*(1-1/sqrt(n*Pi)). - Fung Lam, Mar 21 2014
a(n) = Sum_{m = n+2..2*(n+1)} binomial(2*(n+1), m), n >= 0. - Wolfdieter Lang, May 22 2015
a(n) = A000302(n) + A008549(n). - Gus Wiseman, Jul 19 2021
a(n) = Sum_{j=1..n+1} Sum_{k=1..j} 2^(j-k)*binomial(n+k-1, n). - Fabio Visonà, May 04 2022
a(n) = (1/2)*(-1)^n*binomial(-(n+1), n+2)*hypergeom([1, 2*n + 3], [n + 3], 1/2). - Peter Luschny, Nov 29 2023

Extensions

Corrected by Christian G. Bower

A023532 a(n) = 0 if n is of the form m*(m+3)/2, otherwise 1.

Original entry on oeis.org

0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Keywords

Comments

From Stark: "alpha = 0.101101110111101111101111110 ... is irrational. For if alpha were rational, its decimal expansion would be periodic and have a period of length r starting with the k-th digit of the expansion.
"But by the very nature of alpha, there will be blocks of r digits, all 1, in this expansion after the k-th digit and the periodicity would then guarantee that everything after such a block of r digits would also be all ones.
"This contradicts the fact that there will always be zeros occurring after any given point in the expansion of alpha. Hence alpha is irrational."
a(A000096(n)) = 0; a(A007401(n)) = 1. - Reinhard Zumkeller, Dec 04 2012
Sequence B is called a reverse reluctant sequence of sequence A, if B is triangle array read by rows: row number k lists first k elements of the sequence A in reverse order. A023532 is reverse reluctant sequence of sequence A211666. - Boris Putievskiy, Jan 11 2013
An example of a sequence with infinite critical exponent [Vaslet]. - N. J. A. Sloane, May 05 2013

Examples

			From _Boris Putievskiy_, Jan 11 2013: (Start)
As a triangular array written by rows, the sequence begins:
  0;
  1, 0;
  1, 1, 0;
  1, 1, 1, 0;
  1, 1, 1, 1, 0;
  1, 1, 1, 1, 1, 0;
  1, 1, 1, 1, 1, 1, 0;
  ...
(End)
		

References

  • Harold M. Stark, An Introduction to Number Theory, The MIT Press, Cambridge, Mass, eighth printing 1994, page 170.

Crossrefs

Essentially the same sequence as A114607 and A123110. - N. J. A. Sloane, Feb 07 2020

Programs

  • Haskell
    a023532 = (1 -) . a010052 . (+ 9) . (* 8)
    a023532_list = concat $ iterate (\rs -> 1 : rs) [0]
    -- Reinhard Zumkeller, Dec 04 2012
    
  • Maple
    A023532 := proc(n)
        option remember ;
        local m,t ;
        for m from 0 do
            t := m*(m+3)/2 ;
            if t > n then
                return 1 ;
            elif t = n then
                return 0 ;
            end if;
        end do:
    end proc:
    seq(A023532(n),n=0..40) ; # R. J. Mathar, May 15 2025
  • Mathematica
    a = {}; Do[a = Append[a, Join[ {0}, Table[1, {n} ] ] ], {n, 1, 13} ]; a = Flatten[a]
    Table[PadLeft[{0},n,1],{n,0,20}]//Flatten (* Harvey P. Dale, Jul 10 2019 *)
  • PARI
    for(n=1,9,print1("0, ");for(i=1,n,print1("1, "))) \\ Charles R Greathouse IV, Jun 16 2011
    
  • PARI
    a(n)=!issquare(8*n+9) \\ Charles R Greathouse IV, Jun 16 2011
    
  • Python
    from sympy.ntheory.primetest import is_square
    def A023532(n): return bool(is_square((n<<3)+9))^1 # Chai Wah Wu, Feb 10 2023

Formula

a(n) = 0 if and only if 8n+9 is a square. - Charles R Greathouse IV, Jun 16 2011
Blocks of lengths 1, 2, 3, 4, ... of ones separated by a single zero.
a(n) = 1 - floor((sqrt(9+8n)-1)/2) + floor((sqrt(1+8n)-1)/2). - Paul Barry, May 25 2004
a(n) = A211666(m), where m = (t^2 + 3*t + 4)/2n - n, t = floor((-1 + sqrt(8*n-7))/2). - Boris Putievskiy, Jan 11 2013
a(n) = [A002262(n) < A003056(n)]. - Yuchun Ji, May 18 2020
a(n) = 1-A023531(n). - R. J. Mathar, May 15 2025

Extensions

Additional comments from Robert G. Wilson v, Nov 06 2000

A034781 Triangle of number of rooted trees with n >= 2 nodes and height h >= 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 6, 8, 4, 1, 1, 10, 18, 13, 5, 1, 1, 14, 38, 36, 19, 6, 1, 1, 21, 76, 93, 61, 26, 7, 1, 1, 29, 147, 225, 180, 94, 34, 8, 1, 1, 41, 277, 528, 498, 308, 136, 43, 9, 1, 1, 55, 509, 1198, 1323, 941, 487, 188, 53, 10, 1
Offset: 2

Views

Author

Keywords

Examples

			Triangle begins:
  1;
  1  1;
  1  2  1;
  1  4  3  1;
  1  6  8  4  1;
  1 10 18 13  5  1;
  1 14 38 36 19  6 1;
thus there are 10 trees with 7 nodes and height 2.
		

Crossrefs

T(2n,n) = A245102(n), T(2n+1,n) = A245103(n).
Row sums give A000081.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1 or k<1, 0,
         add(binomial(b((i-1)$2, k-1)+j-1, j)*b(n-i*j, i-1, k), j=0..n/i)))
        end:
    T:= (n, k)-> b((n-1)$2, k) -b((n-1)$2, k-1):
    seq(seq(T(n, k), k=1..n-1), n=2..16);  # Alois P. Heinz, Jul 31 2013
  • Mathematica
    Drop[Map[Select[#, # > 0 &] &,
       Transpose[
        Prepend[Table[
          f[n_] :=
           Nest[CoefficientList[
              Series[Product[1/(1 - x^i)^#[[i]], {i, 1, Length[#]}], {x,
                0, 10}], x] &, {1}, n]; f[m] - f[m - 1], {m, 2, 10}],
    Prepend[Table[1, {10}], 0]]]], 1] // Grid (* Geoffrey Critzer, Aug 01 2013 *)
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i<1 || k<1, 0, Sum[Binomial[b[i-1, i-1, k-1]+j-1, j]*b[n-i*j, i-1, k], {j, 0, n/i}]]]; T[n_, k_] := b[n-1, n-1, k]-b[n-1, n-1, k-1]; Table[T[n, k], {n, 2, 16}, {k, 1, n-1}] // Flatten (* Jean-François Alcover, Feb 11 2014, after Alois P. Heinz *)
  • Python
    def A034781(n, k): return A375467(n, k) - A375467(n, k - 1)
    for n in range(2, 10): print([A034781(n, k) for k in range(2, n + 1)])
    # Peter Luschny, Aug 30 2024

Formula

Reference gives recurrence.
T(n, k) = A375467(n, k) - A375467(n, k - 1). - Peter Luschny, Aug 30 2024

Extensions

More terms from Victoria A Sapko (vsapko(AT)canes.gsw.edu), Sep 19 2003

A014132 Complement of triangular numbers (A000217); also array T(n,k) = ((n+k)^2 + n-k)/2, n, k > 0, read by antidiagonals.

Original entry on oeis.org

2, 4, 5, 7, 8, 9, 11, 12, 13, 14, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 46, 47, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 79
Offset: 1

Views

Author

Keywords

Comments

Numbers that are not triangular (nontriangular numbers).
Also definable as follows: a(1)=2; for n>1, a(n) is smallest integer greater than a(n-1) such that the condition "n and a(a(n)) have opposite parities" can always be satisfied. - Benoit Cloitre and Matthew Vandermast, Mar 10 2003
Record values in A256188 that are greater than 1. - Reinhard Zumkeller, Mar 26 2015
From Daniel Forgues, Apr 10 2015: (Start)
With n >= 1, k >= 1:
t(n+k) - k, 1 <= k <= n+k-1, n >= 1;
t(n+k-1) + n, 1 <= n <= n+k-1, k >= 1;
where t(n+k) = t(n+k-1) + (n+k) is the (n+k)-th triangular number, while the number of compositions of n+k into 2 parts is C(n+k-1, 2-1) = n+k-1, the number of nontriangular numbers between t(n+k-1) and t(n+k), just right!
Related to Hilbert's Infinite Hotel:
0) All rooms, numbered through the positive integers, are full;
1) An infinite number of trains, each containing an infinite number of passengers, arrives: i.e., a 2-D lattice of pairs of positive integers;
2) Move occupant of room m, m >= 1, to room t(m) = m*(m+1)/2, where t(m) is the m-th triangular number;
3) Assign n-th passenger from k-th train to room t(n+k-1) + n, 1 <= n <= n+k-1, k >= 1;
4) Everybody has his or her own room, no room is empty, for m >= 1.
If situation 1 happens again, repeat steps 2 and 3, you're back to 4.
(End)
1711 + 2*a(n)*(58 + a(n)) is prime for n<=21. The terms that do not have this property start 29,32,34,43,47,58,59,60,62,63,65,68,70,73,... - Benedict W. J. Irwin, Nov 22 2016
Also numbers k with the property that in the symmetric representation of sigma(k) both Dyck paths have a central peak or both Dyck paths have a central valley. (Cf. A237593.) - Omar E. Pol, Aug 28 2018

Examples

			From _Boris Putievskiy_, Jan 14 2013: (Start)
Start of the sequence as a table (read by antidiagonals, right to left), where the k-th row corresponds to the k-th column of the triangle (shown thereafter):
   2,  4,  7, 11, 16, 22, 29, ...
   5,  8, 12, 17, 23, 30, 38, ...
   9, 13, 18, 24, 31, 39, 48, ...
  14, 19, 25, 32, 40, 49, 59, ...
  20, 26, 33, 41, 50, 60, 71, ...
  27, 34, 42, 51, 61, 72, 84, ...
  35, 43, 52, 62, 73, 85, 98, ...
  (...)
Start of the sequence as a triangle (read by rows), where the i elements of the i-th row are t(i) + 1 up to t(i+1) - 1, i >= 1:
   2;
   4,  5;
   7,  8,  9;
  11, 12, 13, 14;
  16, 17, 18, 19, 20;
  22, 23, 24, 25, 26, 27;
  29, 30, 31, 32, 33, 34, 35;
  (...)
Row number i contains i numbers, where t(i) = i*(i+1)/2:
  t(i) + 1, t(i) + 2, ..., t(i) + i = t(i+1) - 1
(End) [Edited by _Daniel Forgues_, Apr 11 2015]
		

Crossrefs

Cf. A000124 (left edge: quasi-triangular numbers), A000096 (right edge: almost-triangular numbers), A006002 (row sums), A001105 (central terms).
Cf. A242401 (subsequence).
Cf. A145397 (the non-tetrahedral numbers).

Programs

  • Haskell
    a014132 n = n + round (sqrt $ 2 * fromInteger n)
    a014132_list = filter ((== 0) . a010054) [0..]
    -- Reinhard Zumkeller, Dec 12 2012
    
  • Magma
    IsTriangular:=func< n | exists{ k: k in [1..Isqrt(2*n)] | n eq (k*(k+1) div 2)} >; [ n: n in [1..90] | not IsTriangular(n) ]; // Klaus Brockhaus, Jan 04 2011
    
  • Mathematica
    f[n_] := n + Round[Sqrt[2n]]; Array[f, 71] (* or *)
    Complement[ Range[83], Array[ #(# + 1)/2 &, 13]] (* Robert G. Wilson v, Oct 21 2005 *)
    DeleteCases[Range[80],?(OddQ[Sqrt[8#+1]]&)] (* _Harvey P. Dale, Jul 24 2021 *)
  • PARI
    a(n)=if(n<1,0,n+(sqrtint(8*n-7)+1)\2)
    
  • PARI
    isok(n) = !ispolygonal(n,3); \\ Michel Marcus, Mar 01 2016
    
  • Python
    from math import isqrt
    def A014132(n): return n+(isqrt((n<<3)-7)+1>>1) # Chai Wah Wu, Jun 17 2024

Formula

a(n) = n + round(sqrt(2*n)).
a(a(n)) = n + 2*floor(1/2 + sqrt(2n)) + 1.
a(n) = a(n-1) + A035214(n), a(1)=2.
a(n) = A080036(n) - 1.
a(n) = n + A002024(n). - Vincenzo Librandi, Jul 08 2010
A010054(a(n)) = 0. - Reinhard Zumkeller, Dec 10 2012
From Boris Putievskiy, Jan 14 2013: (Start)
a(n) = A007401(n)+1.
a(n) = A003057(n)^2 - A114327(n).
a(n) = ((t+2)^2 + i - j)/2, where
i = n-t*(t+1)/2,
j = (t*t+3*t+4)/2-n,
t = floor((-1+sqrt(8*n-7))/2). (End)
A248952(a(n)) < 0. - Reinhard Zumkeller, Oct 20 2014
a(n) = A256188(A004202(n)). - Reinhard Zumkeller, Mar 26 2015
From Robert Israel, Apr 20 2015 (Start):
a(n) = A118011(n) - n.
G.f.: x/(1-x)^2 + x/(1-x) * Sum(j>=0, x^(j*(j+1)/2)) = x/(1-x)^2 + x^(7/8)/(2-2*x) * Theta2(0,sqrt(x)), where Theta2 is a Jacobi theta function. (End)
G.f. as array: x*y*(2 - 2*y + x^2*y + y^2 - x*(1 + y))/((1 - x)^3*(1 - y)^3). - Stefano Spezia, Apr 22 2024

Extensions

Following Alford Arnold's comment: keyword tabl and correspondent crossrefs added by Reinhard Zumkeller, Dec 12 2012
I restored the original definition. - N. J. A. Sloane, Jan 27 2019

A002802 a(n) = (2*n+3)!/(6*n!*(n+1)!).

Original entry on oeis.org

1, 10, 70, 420, 2310, 12012, 60060, 291720, 1385670, 6466460, 29745716, 135207800, 608435100, 2714556600, 12021607800, 52895074320, 231415950150, 1007340018300, 4365140079300, 18839025605400, 81007810103220, 347176329013800, 1483389769422600
Offset: 0

Views

Author

Keywords

Comments

For n >= 1 a(n) is also the number of rooted bicolored unicellular maps of genus 1 on n+2 edges. - Ahmed Fares (ahmedfares(AT)my-deja.com), Aug 20 2001
a(n) is half the number of (n+2) X 2 Young tableaux with a three horizontal walls between the first and second column. If there is a wall between two cells, the entries may be decreasing; see [Banderier, Wallner 2021], A000984 for one horizontal wall, and A002457 for two. - Michael Wallner, Jan 31 2022
From Robert Coquereaux, Feb 12 2024: (Start)
Call B(p,g) the number of genus g partitions of a set with p elements (genus-dependent Bell number). Up to an appropriate shift the given sequence counts the genus 1 partitions of a set: we have a(n) = B(n+4,1), with a(0)= B(4,1)=1.
When shifted with an offset 4 (i.e., defining b(p)=a(p-4), which starts with 0,0,0,1,10,70, etc., and b(4)=1), the given sequence reads b(p) = (1/( 2^4 3 )) * (1/( (2 p - 1) (2 p - 3))) * (1/(p - 4)!) * (2p)!/p!. In this form it appears as a generalization of Catalan numbers (that indeed count the genus 0 partitions).
Call C[p, [alpha], g] the number of partitions of a set with p elements, of cyclic type [alpha], and of genus g (genus g Faa di Bruno coefficients of type [alpha]). Up to an appropriate shift the given sequence also counts the genus 1 partitions of p=2k into k parts of length 2, which is then called C[2k, [2^k], 1], and we have a(n) = C[2k, [2^k], 1] for k=n+2.
The two previous interpretations of this sequence, leading to a(n) = B(n+4, 1) and to a(n) = C[2(n+2), [2^(n+2)], 1] are not related in any obvious way. (End)

Examples

			G.f. = 1 + 10*x + 70*x^2 + 420*x^3 + 2310*x^4 + 12012*x^5 + 60060*x^6 + ...
		

References

  • C. Jordan, Calculus of Finite Differences. Röttig and Romwalter, Budapest, 1939; Chelsea, NY, 1965, p. 449.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A035309, A000108 (for genus 0 maps), A046521 (third column).
Column g=1 of A370235.

Programs

  • GAP
    F:=Factorial;; List([0..25], n-> F(2*n+3)/(6*F(n)*F(n+1)) ); # G. C. Greubel, Jul 20 2019
  • Magma
    F:=Factorial; [F(2*n+3)/(6*F(n)*F(n+1)): n in [0..25]]; // G. C. Greubel, Jul 20 2019
    
  • Maple
    seq(simplify(4^n*hypergeom([-n,-3/2], [1], 1)),n=0..25); # Peter Luschny, Apr 26 2016
  • Mathematica
    Table[(2*n+3)!/(6*n!*(n+1)!), {n, 0, 25}] (* Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
  • PARI
    {a(n) = if( n<0, 0, (2*n + 3)! / (6 * n! * (n+1)!))}; /* Michael Somos, Sep 16 2013 */
    
  • PARI
    {a(n) = 2^(n+3) * polcoeff( pollegendre(n+4), n) / 3}; /* Michael Somos, Sep 16 2013 */
    
  • Sage
    f=factorial; [f(2*n+3)/(6*f(n)*f(n+1)) for n in (0..25)] # G. C. Greubel, Jul 20 2019
    

Formula

G.f.: (1 - 4*x)^(-5/2) = 1F0(5/2;;4x).
Asymptotic expression for a(n) is a(n) ~ (n+2)^(3/2) * 4^(n+2) / (sqrt(Pi) * 48).
a(n) = Sum_{a+b+c+d+e=n} f(a)*f(b)*f(c)*f(d)*f(e) with f(n) = binomial(2n, n) = A000984(n). - Philippe Deléham, Jan 22 2004
a(n-1) = (1/4)*Sum_{k=1..n} k*(k+1)*binomial(2*k, k). - Benoit Cloitre, Mar 20 2004
a(n) = A051133(n+1)/3 = A000911(n)/6. - Zerinvary Lajos, Jun 02 2007
From Rui Duarte, Oct 08 2011: (Start)
Also convolution of A000984 with A002697, also convolution of A000302 with A002457.
a(n) = ((2n+3)(2n+1)/(3*1)) * binomial(2n, n).
a(n) = binomial(2n+4, 4) * binomial(2n, n) / binomial(n+2, 2).
a(n) = binomial(n+2, 2) * binomial(2n+4, n+2) / binomial(4, 2).
a(n) = binomial(2n+4, n+2) * (n+2)*(n+1) / 12. (End)
D-finite with recurrence: n*a(n) - 2*(2*n+3)*a(n-1) = 0. - R. J. Mathar, Jan 31 2014
a(n) = 4^n*hypergeom([-n,-3/2], [1], 1). - Peter Luschny, Apr 26 2016
Boas-Buck recurrence: a(n) = (10/n)*Sum_{k=0..n-1} 4^(n-k-1)*a(k), n >= 1, a(0) = 1. Proof from a(n) = A046521(n+2, 2). See a comment there. - Wolfdieter Lang, Aug 10 2017
a(n) = (-4)^n*binomial(-5/2, n). - Peter Luschny, Oct 23 2018
Sum_{n>=0} 1/a(n) = 12 - 2*sqrt(3)*Pi. - Amiram Eldar, Oct 13 2020
E.g.f.: (1/12) exp(2 x) x^2 BesselI[2, 2 x]. - Robert Coquereaux, Feb 12 2024

A000065 -1 + number of partitions of n.

Original entry on oeis.org

0, 0, 1, 2, 4, 6, 10, 14, 21, 29, 41, 55, 76, 100, 134, 175, 230, 296, 384, 489, 626, 791, 1001, 1254, 1574, 1957, 2435, 3009, 3717, 4564, 5603, 6841, 8348, 10142, 12309, 14882, 17976, 21636, 26014, 31184, 37337, 44582, 53173, 63260, 75174, 89133, 105557, 124753
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the number of noncongruent n-dimensional integer-sided simplices with diameter n. - Sascha Kurz, Jul 26 2004
Also, the number of partitions of n into parts each less than n.
Also, the number of distinct types of equation which can be derived from the equation [n,0,0] not including itself. (Ince)
Also, the number of rooted trees on n+1 nodes with height exactly 2.
Also, the number of partitions (of any positive integer) whose sum + length is <= n. Example: a(5) = 6 counts 4, 3, 21, 2, 11, 1. Proof: Given a partition of n other than the all 1s partition, subtract 1 from each part and then drop the zeros. This is a bijection to the partitions with sum + length <= n. - David Callan, Nov 29 2007
Number of graphs with n vertices of treewidth n-2. Reason: The complement of a graph with n vertices and treewidth >= n-2 cannot have P3 or K3 as a subgraph (Chlebı́ková 2002, Theorem 10), so every component of it is a star. - Martín Muñoz, Dec 31 2023

Examples

			G.f. = x^2 + 2*x^3 + 4*x^4 + 6*x^5 + 10*x^6 + 14*x^7 + 21*x^8 + 29*x^9 + ...
		

References

  • E. L. Ince, Ordinary Differential Equations, Dover Publications, New York, 1944, p. 498; MR0010757.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A000041 - 1. A column of A058716. A diagonal of A263294.
Column h=2 of A034781.

Programs

  • Magma
    [NumberOfPartitions(n)-1: n in [0..50]]; // Vincenzo Librandi, Aug 25 2013
  • Maple
    with (combstruct):ZL:=proc(m) local i; [T0,{seq(T.i=Prod(Z,Set(T.(i+1))),i=0..m-1), T.m=Z}, unlabeled] end:A:=n -> count(ZL(2),size=n)-count(ZL(1),size=n): seq(A(n),n=1..46); # Zerinvary Lajos, Dec 05 2007
    ZL :=[S, {S = Set(Cycle(Z),1 < card)}, unlabelled]: seq(combstruct[count](ZL, size=n), n=0..45); # Zerinvary Lajos, Mar 25 2008
  • Mathematica
    nn=40;CoefficientList[Series[Product[1/(1-x^i),{i,1,nn}]-1/(1-x),{x,0,nn}],x]  (* Geoffrey Critzer, Oct 28 2012 *)
    PartitionsP[Range[0,50]]-1 (* Harvey P. Dale, Aug 24 2013 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / eta(x + x*O(x^n)), n) - 1)};
    
  • PARI
    {a(n) = if( n<0, 0, numbpart(n) - 1)};
    

Formula

a(n) = A026820(n,n-1) for n>1. - Reinhard Zumkeller, Jan 21 2010
G.f.: x*G(0)/(x-1) where G(k) = 1 - 1/(1-x^(k+1))/(1-x/(x-1/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 23 2013
G.f.: Sum_{k>=2} x^k / Product_{j=1..k} (1 - x^j). - Ilya Gutkovskiy, Sep 07 2021

A002840 Number of polyhedral graphs with n edges.

Original entry on oeis.org

1, 0, 1, 2, 2, 4, 12, 22, 58, 158, 448, 1342, 4199, 13384, 43708, 144810, 485704, 1645576, 5623571, 19358410, 67078828, 233800162, 819267086, 2884908430, 10204782956, 36249143676, 129267865144, 462669746182, 1661652306539, 5986979643542
Offset: 6

Views

Author

Keywords

References

  • M. B. Dillencourt, Polyhedra of small orders and their Hamiltonian properties. Tech. Rep. 92-91, Info. and Comp. Sci. Dept., Univ. Calif. Irvine, 1992.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • T. R. S. Walsh, personal communication.

Crossrefs

Programs

  • PARI
    \\ It is assumed that the 3cp.gp file (from the linked zip archive) has been read before, i.e., \r [path]3cp.gp
    for(k=6,#ThreeConnectedData,print1(#ThreeConnectedData[k],", "));
    \\ printing of the edge lists of the graphs for n <= 11
    print(ThreeConnectedData[6..11]) \\ Hugo Pfoertner, Feb 14 2021

Extensions

a(30)-a(35) from the Numericana link added by Andrey Zabolotskiy, Jun 13 2020

A006384 Number of sensed planar maps with n edges.

Original entry on oeis.org

1, 2, 4, 14, 57, 312, 2071, 15030, 117735, 967850, 8268816, 72833730, 658049140, 6074058060, 57106433817, 545532037612, 5284835906037, 51833908183164, 514019531037910, 5147924676612282, 52017438279806634, 529867070532745464
Offset: 0

Views

Author

Keywords

Comments

The planar maps considered are connected and may contain loops and parallel edges. - Andrew Howroyd, Jan 13 2025

References

  • V. A. Liskovets, A census of nonisomorphic planar maps, in Algebraic Methods in Graph Theory, Vol. II, ed. L. Lovasz and V. T. Sos, North-Holland, 1981.
  • V. A. Liskovets, Enumeration of nonisomorphic planar maps, Selecta Math. Sovietica, 4 (No. 4, 1985), 303-323.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • T. R. S. Walsh, personal communication.

Crossrefs

Antidiagonal sums of A379430.
Cf. A000168 (rooted), A006385 (unsensed), A006443 (achiral), A006402 (2-connected).

Programs

  • Maple
    with(numtheory): a:= n-> `if` (n=0, 1, floor (2*3^n /(n+1)/(n+2) *binomial(2*n, n) +add (phi(n/t) *3^t *binomial(2*t, t), t=divisors(n) minus {n}))/2/n +`if` (irem(n,2)=1, 2*3^((n-1)/2) /(n+1) *binomial(n-1, (n-1)/2), 2*(n-1) *3^((n-2)/2) /n/(n+2) *binomial(n-2, (n-2)/2))): seq (a(n), n=0..30); # Alois P. Heinz, Apr 24 2009
  • Mathematica
    a[0] = 1; a[n_] := (1/(2n))*(2*(3^n/((n+1)*(n+2)))*Binomial[2n, n] + Sum[ EulerPhi[n/k]*3^k*Binomial[ 2k, k], {k, Most[ Divisors[n]]}]) + q[n]; q[n_?OddQ] := 2*(3^((n-1)/2)/(n+1))*Binomial[ n-1, (n-1)/2]; q[n_?EvenQ] := 2*(n-1)*(3^((n-2)/2)/(n*(n+2)))*Binomial[ n-2, (n-2)/2]; Table[ a[n], {n, 0, 21}] (* Jean-François Alcover, after Valery A. Liskovets *)

Formula

For n>0, a(n) = (1/2n)[A'(n)+sum_{kA000010, q(n)=(n+3) A'(n-1/2)/4 if n is odd and q(n) = (n-1)A'(n-2/2)/4 if n is even, where A'(n)=A000168(n), the number of rooted maps. - Valery A. Liskovets, May 27 2006
Equivalently, a(n) = (1/2n)[2*3^n/((n+1)(n+2))*binomial(2n,n) +sum_{kValery A. Liskovets, May 27 2006
a(n) ~ 12^n / (sqrt(Pi) * n^(7/2)). - Vaclav Kotesovec, Sep 12 2014

Extensions

More terms from Alois P. Heinz, Apr 24 2009

A000235 Number of n-node rooted trees of height 3.

Original entry on oeis.org

0, 0, 0, 1, 3, 8, 18, 38, 76, 147, 277, 509, 924, 1648, 2912, 5088, 8823, 15170, 25935, 44042, 74427, 125112, 209411, 348960, 579326, 958077, 1579098, 2593903, 4247768, 6935070, 11290627, 18330973, 29684082, 47946852, 77258764, 124198083
Offset: 1

Views

Author

Keywords

Comments

(1, 1, 2, 3, 5, 8, ...) convolved with (0, 0, 1, 2, 4, 7, ...) = (0, 0, 1, 3, 8, ...). - Gary W. Adamson, Aug 14 2010

References

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

Crossrefs

Column h=3 of A034781.

Programs

  • Maple
    # For Maple program see link.
    with(combstruct):
    ZL:= proc(m) local i; [T0, {seq(T||i=Prod(Z, Set(T||(i+1))), i=0..m-1), T||m=Z}, unlabeled] end: A000235:= n-> count(ZL(3), size=n)-count(ZL(2), size=n): seq(A000235(n), n=1..36); # Zerinvary Lajos, Sep 23 2007
  • Mathematica
    m = 36; Rest @ CoefficientList[ Series[x*Product[(1-x^k)^(-PartitionsP[k-1]), {k, 1, m}], {x, 0, m}], x] - PartitionsP[Range[0, m-1]] (* Jean-François Alcover, Jul 05 2011, after Christian G. Bower *)

Formula

a(n) = A001383(n) - A000041(n-1). - Christian G. Bower
Showing 1-10 of 47 results. Next