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 13 results. Next

A001764 a(n) = binomial(3*n,n)/(2*n+1) (enumerates ternary trees and also noncrossing trees).

Original entry on oeis.org

1, 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715, 8414640, 50067108, 300830572, 1822766520, 11124755664, 68328754959, 422030545335, 2619631042665, 16332922290300, 102240109897695, 642312451217745, 4048514844039120, 25594403741131680, 162250238001816900
Offset: 0

Views

Author

Keywords

Comments

Smallest number of straight line crossing-free spanning trees on n points in the plane.
Number of dissections of some convex polygon by nonintersecting diagonals into polygons with an odd number of sides and having a total number of 2n+1 edges (sides and diagonals). - Emeric Deutsch, Mar 06 2002
Number of lattice paths of n East steps and 2n North steps from (0,0) to (n,2n) and lying weakly below the line y=2x. - David Callan, Mar 14 2004
With interpolated zeros, this has g.f. 2*sqrt(3)*sin(arcsin(3*sqrt(3)*x/2)/3)/(3*x) and a(n) = C(n+floor(n/2),floor(n/2))*C(floor(n/2),n-floor(n/2))/(n+1). This is the first column of the inverse of the Riordan array (1-x^2,x(1-x^2)) (essentially reversion of y-y^3). - Paul Barry, Feb 02 2005
Number of 12312-avoiding matchings on [2n].
Number of complete ternary trees with n internal nodes, or 3n edges.
Number of rooted plane trees with 2n edges, where every vertex has even outdegree ("even trees").
a(n) is the number of noncrossing partitions of [2n] with all blocks of even size. E.g.: a(2)=3 counts 12-34, 14-23, 1234. - David Callan, Mar 30 2007
Pfaff-Fuss-Catalan sequence C^{m}_n for m=3, see the Graham et al. reference, p. 347. eq. 7.66.
Also 3-Raney sequence, see the Graham et al. reference, p. 346-7.
The number of lattice paths from (0,0) to (2n,0) using an Up-step=(1,1) and a Down-step=(0,-2) and staying above the x-axis. E.g., a(2) = 3; UUUUDD, UUUDUD, UUDUUD. - Charles Moore (chamoore(AT)howard.edu), Jan 09 2008
a(n) is (conjecturally) the number of permutations of [n+1] that avoid the patterns 4-2-3-1 and 4-2-5-1-3 and end with an ascent. For example, a(4)=55 counts all 60 permutations of [5] that end with an ascent except 42315, 52314, 52413, 53412, all of which contain a 4-2-3-1 pattern and 42513. - David Callan, Jul 22 2008
Central terms of pendular triangle A167763. - Philippe Deléham, Nov 12 2009
With B(x,t)=x+t*x^3, the comp. inverse in x about 0 is A(x,t) = Sum_{j>=0} a(j) (-t)^j x^(2j+1). Let U(x,t)=(x-A(x,t))/t. Then DU(x,t)/Dt=dU/dt+U*dU/dx=0 and U(x,0)=x^3, i.e., U is a solution of the inviscid Burgers's, or Hopf, equation. Also U(x,t)=U(x-t*U(x,t),0) and dB(x,t)/dt = U(B(x,t),t) = x^3 = U(x,0). The characteristics for the Hopf equation are x(t) = x(0) + t*U(x(t),t) = x(0) + t*U(x(0),0) = x(0) + t*x(0)^3 = B(x(0),t). These results apply to all the Fuss-Catalan sequences with 3 replaced by n>0 and 2 by n-1 (e.g., A000108 with n=2 and A002293 with n=4), see also A086810, which can be generalized to A133437, for associahedra. - Tom Copeland, Feb 15 2014
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Kreweras lattice (noncrossing partitions ordered by refinement) of size n, see the Bernardi & Bonichon (2009) and Kreweras (1972) references. - Noam Zeilberger, Jun 01 2016
Number of sum-indecomposable (4231,42513)-avoiding permutations. Conjecturally, number of sum-indecomposable (2431,45231)-avoiding permutations. - Alexander Burstein, Oct 19 2017
a(n) is the number of topologically distinct endstates for the game Planted Brussels Sprouts on n vertices, see Ji and Propp link. - Caleb Ji, May 14 2018
Number of complete quadrillages of 2n+2-gons. See Baryshnikov p. 12. See also Nov 10 2014 comments in A134264. - Tom Copeland, Jun 04 2018
a(n) is the number of 2-regular words on the alphabet [n] that avoid the patterns 231 and 221. Equivalently, this is the number of 2-regular tortoise-sortable words on the alphabet [n] (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
a(n) is the number of Motzkin paths of length 3n with n steps of each type, with the condition that (1, 0) and (1, 1) steps alternate (starting with (1, 0)). - Helmut Prodinger, Apr 08 2019
a(n) is the number of uniquely sorted permutations of length 2n+1 that avoid the patterns 312 and 1342. - Colin Defant, Jun 08 2019
The compositional inverse o.g.f. pair in Copeland's comment above are related to a pair of quantum fields in Balduf's thesis by Theorem 4.2 on p. 92. - Tom Copeland, Dec 13 2019
The sequences of Fuss-Catalan numbers, of which this is the first after the Catalan numbers A000108 (the next is A002293), appear in articles on random matrices and quantum physics. See Banica et al., Collins et al., and Mlotkowski et al. Interpretations of these sequences in terms of the cardinality of specific sets of noncrossing partitions are provided by A134264. - Tom Copeland, Dec 21 2019
Call C(p, [alpha], g) the number of partitions of a cyclically ordered set with p elements, of cyclic type [alpha], and of genus g (the genus g Faa di Bruno coefficients of type [alpha]). This sequence counts the genus 0 partitions (non-crossing, or planar, partitions) of p = 3n into n parts of length 3: a(n) = C(3n, [3^n], 0). For genus 1 see A371250, for genus 2 see A371251. - Robert Coquereaux, Mar 16 2024
a(n) is the total number of down steps before the first up step in all 2_1-Dyck paths of length 3*n for n > 0. A 2_1-Dyck path is a lattice path with steps (1,2), (1,-1) that starts and ends at y = 0 and does not go below the line y = -1. - Sarah Selkirk, May 10 2020
a(n) is the number of pairs (A<=B) of noncrossing partitions of [n]. - Francesca Aicardi, May 28 2022
a(n) is the number of parking functions of size n avoiding the patterns 231 and 321. - Lara Pudwell, Apr 10 2023
Number of rooted polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {4,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 3 of the family {C(k, n)}A130564.%20-%20_Wolfdieter%20Lang">{n>=0} given in a comment in A130564. - _Wolfdieter Lang, Feb 05 2024
The number of Apollonian networks (planar 3-trees) with n+3 vertices with a given base triangle. - Allan Bickle, Feb 20 2024
Number of rooted polyominoes composed of n tetrahedral cells of the hyperbolic regular tiling with Schläfli symbol {3,3,oo}. A rooted polyomino has one external face identified, and chiral pairs are counted as two. a(n) = T(n) in the second Beineke and Pippert link. - Robert A. Russell, Mar 20 2024

Examples

			a(2) = 3 because the only dissections with 5 edges are given by a square dissected by any of the two diagonals and the pentagon with no dissecting diagonal.
G.f. = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 + 7752*x^7 + 43263*x^8 + ...
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • I. M. H. Etherington, On non-associative combinations, Proc. Royal Soc. Edinburgh, 59 (Part 2, 1938-39), 153-162.
  • I. M. H. Etherington, Some problems of non-associative combinations (I), Edinburgh Math. Notes, 32 (1940), pp. i-vi. Part II is by A. Erdelyi and I. M. H. Etherington, and is on pages vii-xiv of the same issue.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347. See also the Pólya-Szegő reference.
  • W. Kuich, Languages and the enumeration of planted plane trees. Nederl. Akad. Wetensch. Proc. Ser. A 73 = Indag. Math. 32, (1970), 268-280.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, p. 98.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, New York, Heidelberg, Berlin, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • 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. A001762, A001763, A002294 - A002296, A006013, A025174, A063548, A064017, A072247, A072248, A134264, A143603, A258708, A256311, A188687 (binomial transform), A346628 (inverse binomial transform).
A column of triangle A102537.
Bisection of A047749 and A047761.
Row sums of triangles A108410 and A108767.
Second column of triangle A062993.
Mod 3 = A113047.
2D Polyominoes: A005034 (oriented), A005036 (unoriented), A369315 (chiral), A047749 (achiral), A000108 {3,oo}, A002293 {5,oo}.
3D Polyominoes: A007173 (oriented), A027610 (unoriented), A371350 (chiral), A371351 (achiral).
Cf. A130564 (for C(k, n) cases).

Programs

  • GAP
    List([0..25],n->Binomial(3*n,n)/(2*n+1)); # Muniru A Asiru, Oct 31 2018
    
  • Haskell
    a001764 n = a001764_list !! n
    a001764_list = 1 : [a258708 (2 * n) n | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [Binomial(3*n,n)/(2*n+1): n in [0..30]]; // Vincenzo Librandi, Sep 04 2014
    
  • Maple
    A001764 := n->binomial(3*n,n)/(2*n+1): seq(A001764(n), n=0..25);
    with(combstruct): BB:=[T,{T=Prod(Z,F),F=Sequence(B),B=Prod(F,Z,F)}, unlabeled]:seq(count(BB,size=i),i=0..22); # Zerinvary Lajos, Apr 22 2007
    with(combstruct):BB:=[S, {B = Prod(S,S,Z), S = Sequence(B)}, labelled]: seq(count(BB, size=n)/n!, n=0..21); # Zerinvary Lajos, Apr 25 2008
    n:=30:G:=series(RootOf(g = 1+x*g^3, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
    alias(PS=ListTools:-PartialSums): A001764List := proc(m) local A, P, n;
    A := [1,1]; P := [1]; for n from 1 to m - 2 do P := PS(PS([op(P), P[-1]]));
    A := [op(A), P[-1]] od; A end: A001764List(25); # Peter Luschny, Mar 26 2022
  • Mathematica
    InverseSeries[Series[y-y^3, {y, 0, 24}], x] (* then a(n)=y(2n+1)=ways to place non-crossing diagonals in convex (2n+4)-gon so as to create only quadrilateral tiles *) (* Len Smiley, Apr 08 2000 *)
    Table[Binomial[3n,n]/(2n+1),{n,0,25}] (* Harvey P. Dale, Jul 24 2011 *)
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n! / (2*n + 1)!)};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( serreverse( x - x^3 + O(x^(2*n + 2))), 2*n + 1))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( m=1, n, A = 1 + x * A^3); polcoeff(A, n))};
    
  • PARI
    b=vector(22);b[1]=1;for(n=2,22,for(i=1,n-1,for(j=1,n-1,for(k=1,n-1,if((i-1)+(j-1)+(k-1)-(n-2),NULL,b[n]=b[n]+b[i]*b[j]*b[k])))));a(n)=b[n+1]; print1(a(0));for(n=1,21,print1(", ",a(n))) \\ Gerald McGarvey, Oct 08 2008
    
  • PARI
    Vec(1 + serreverse(x / (1+x)^3 + O(x^30))) \\ Gheorghe Coserea, Aug 05 2015
    
  • Python
    from math import comb
    def A001764(n): return comb(3*n,n)//(2*n+1) # Chai Wah Wu, Nov 10 2022
  • Sage
    def A001764_list(n) :
        D = [0]*(n+1); D[1] = 1
        R = []; b = false; h = 1
        for i in range(2*n) :
            for k in (1..h) : D[k] += D[k-1]
            if not b : R.append(D[h])
            else : h += 1
            b = not b
        return R
    A001764_list(22) # Peter Luschny, May 03 2012
    

Formula

From Karol A. Penson, Nov 08 2001: (Start)
G.f.: (2/sqrt(3*x))*sin((1/3)*arcsin(sqrt(27*x/4))).
E.g.f.: hypergeom([1/3, 2/3], [1, 3/2], 27/4*x).
Integral representation as n-th moment of a positive function on [0, 27/4]: a(n) = Integral_{x=0..27/4} (x^n*((1/12) * 3^(1/2) * 2^(1/3) * (2^(1/3)*(27 + 3 * sqrt(81 - 12*x))^(2/3) - 6 * x^(1/3))/(Pi * x^(2/3)*(27 + 3 * sqrt(81 - 12*x))^(1/3)))), n >= 0. This representation is unique. (End)
G.f. A(x) satisfies A(x) = 1+x*A(x)^3 = 1/(1-x*A(x)^2) [Cyvin (1998)]. - Ralf Stephan, Jun 30 2003
a(n) = n-th coefficient in expansion of power series P(n), where P(0) = 1, P(k+1) = 1/(1 - x*P(k)^2).
G.f. Rev(x/c(x))/x, where c(x) is the g.f. of A000108 (Rev=reversion of). - Paul Barry, Mar 26 2010
From Gary W. Adamson, Jul 07 2011: (Start)
Let M = the production matrix:
1, 1
2, 2, 1
3, 3, 2, 1
4, 4, 3, 2, 1
5, 5, 4, 3, 2, 1
...
a(n) = upper left term in M^n. Top row terms of M^n = (n+1)-th row of triangle A143603, with top row sums generating A006013: (1, 2, 7, 30, 143, 728, ...). (End)
Recurrence: a(0)=1; a(n) = Sum_{i=0..n-1, j=0..n-1-i} a(i)a(j)a(n-1-i-j) for n >= 1 (counts ternary trees by subtrees of the root). - David Callan, Nov 21 2011
G.f.: 1 + 6*x/(Q(0) - 6*x); Q(k) = 3*x*(3*k + 1)*(3*k + 2) + 2*(2*(k^2) + 5*k +3) - 6*x*(2*(k^2) + 5*k + 3)*(3*k + 4)*(3*k + 5)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 27 2011
D-finite with recurrence: 2*n*(2n+1)*a(n) - 3*(3n-1)*(3n-2)*a(n-1) = 0. - R. J. Mathar, Dec 14 2011
REVERT transform of A115140. BINOMIAL transform is A188687. SUMADJ transform of A188678. HANKEL transform is A051255. INVERT transform of A023053. INVERT transform is A098746. - Michael Somos, Apr 07 2012
(n + 1) * a(n) = A174687(n).
G.f.: F([2/3,4/3], [3/2], 27/4*x) / F([2/3,1/3], [1/2], (27/4)*x) where F() is the hypergeometric function. - Joerg Arndt, Sep 01 2012
a(n) = binomial(3*n+1, n)/(3*n+1) = A062993(n+1,1). - Robert FERREOL, Apr 03 2015
a(n) = A258708(2*n,n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
0 = a(n)*(-3188646*a(n+2) + 20312856*a(n+3) - 11379609*a(n+4) + 1437501*a(n+5)) + a(n+1)*(177147*a(n+2) - 2247831*a(n+3) + 1638648*a(n+4) - 238604*a(n+5)) + a(n+2)*(243*a(n+2) + 31497*a(n+3) - 43732*a(n+4) + 8288*a(n+5)) for all integer n. - Michael Somos, Jun 03 2016
a(n) ~ 3^(3*n + 1/2)/(sqrt(Pi)*4^(n+1)*n^(3/2)). - Ilya Gutkovskiy, Nov 21 2016
Given g.f. A(x), then A(1/8) = -1 + sqrt(5), A(2/27) = (-1 + sqrt(3))*3/2, A(4/27) = 3/2, A(3/64) = -2 + 2*sqrt(7/3), A(5/64) = (-1 + sqrt(5))*2/sqrt(5), etc. A(n^2/(n+1)^3) = (n+1)/n if n > 1. - Michael Somos, Jul 17 2018
From Peter Bala, Sep 14 2021: (Start)
A(x) = exp( Sum_{n >= 1} (1/3)*binomial(3*n,n)*x^n/n ).
The sequence defined by b(n) := [x^n] A(x)^n = A224274(n) for n >= 1 and satisfies the congruence b(p) == b(1) (mod p^3) for prime p >= 3. Cf. A060941. (End)
G.f.: 1/sqrt(B(x)+(1-6*x)/(9*B(x))+1/3), with B(x):=((27*x^2-18*x+2)/54-(x*sqrt((-(4-27*x))*x))/(2*3^(3/2)))^(1/3). - Vladimir Kruchinin, Sep 28 2021
x*A'(x)/A(x) = (A(x) - 1)/(- 2*A(x) + 3) = x + 5*x^2 + 28*x^3 + 165*x^4 + ... is the o.g.f. of A025174. Cf. A002293 - A002296. - Peter Bala, Feb 04 2022
a(n) = hypergeom([1 - n, -2*n], [2], 1). Row sums of A108767. - Peter Bala, Aug 30 2023
G.f.: z*exp(3*z*hypergeom([1, 1, 4/3, 5/3], [3/2, 2, 2], (27*z)/4)) + 1.
- Karol A. Penson, Dec 19 2023
G.f.: hypergeometric([1/3, 2/3], [3/2], (3^3/2^2)*x). See the e.g.f. above. - Wolfdieter Lang, Feb 04 2024
a(n) = (3*n)! / (n!*(2*n+1)!). - Allan Bickle, Feb 20 2024
Sum_{n >= 0} a(n)*x^n/(1 + x)^(3*n+1) = 1. See A316371 and A346627. - Peter Bala, Jun 02 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^5). - Seiichi Manyama, Jun 16 2025

A047749 If n = 2*m then a(n) = binomial(3*m, m)/(2*m+1), if n=2*m+1 then a(n) = binomial(3*m+1, m+1)/(2*m+1).

Original entry on oeis.org

1, 1, 1, 2, 3, 7, 12, 30, 55, 143, 273, 728, 1428, 3876, 7752, 21318, 43263, 120175, 246675, 690690, 1430715, 4032015, 8414640, 23841480, 50067108, 142498692, 300830572, 859515920, 1822766520, 5225264024, 11124755664, 31983672534, 68328754959, 196947587823, 422030545335
Offset: 0

Views

Author

Keywords

Comments

Hankel transform appears to be a signed aerated version of A059492. - Paul Barry, Apr 16 2008
Row sums of inverse Riordan array (1, x*(1-x^2))^(-1). - Paul Barry, Apr 16 2008
a(n) is the number of permutations of length n avoiding 213 in the classical sense which are breadth-first search reading words of increasing unary-binary trees. For more details, see the entry for permutations avoiding 231 at A245898. - Manda Riehl, Aug 05 2014
From David Callan, Aug 22 2014: (Start)
a(n) is the number of ordered trees (A000108) with n vertices in which every non-root non-leaf vertex has exactly one leaf child (no restriction on its non-leaf children). For example, a(4) counts the 3 trees
| |
\/ \|/ \/
(End)
From Emeric Deutsch, Oct 28 2014: (Start)
a(n) is the number of symmetric ternary trees having n internal nodes.
a(n) is the number of symmetric non-crossing rooted trees having n edges.
a(n) is the number of symmetric even trees having 2n edges.
a(n) is the number of symmetric diagonally convex directed polyominoes having n diagonals.
(End)
For the above 4 items see the Deutsch-Feretic-Noy reference.
a(n) is also the number of self-dual labeled non-crossing trees with n edges. See my paper in the links section. - Nikos Apostolakis, Jun 11 2019
Number of achiral polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A stereographic projection of this tiling on the Poincaré disk can be obtained via the Christensson link. An achiral polyomino is identical to its reflection. - Robert A. Russell, Jan 20 2024

Examples

			G.f. = 1 + x + x^2 + 2*x^3 + 3*x^4 + 7*x^5 + 12*x^6 + 30*x^7 + 55*x^8 + ...
		

Crossrefs

Column k=3 of A369929 and k=4 of A370062.
Cf. A006013 is the odd-indexed terms of this sequence.
Polyominoes: A005034 (oriented), A005036 (unoriented), A369315 (chiral), A385149 (asymmetric), A001764 (rooted), A208355(n-1) {3,oo}, A369472 {5,oo}.

Programs

  • Magma
    G:=Gamma; [Round((1+(-1)^n)*G(3*n/2+1)/(G(n/2+1)*Factorial(n+1)) + (1-(-1)^n)*G((3*n+1)/2)/(G((n+3)/2)*Factorial(n)))/2: n in [0..35]]; // G. C. Greubel, Jul 07 2019
    
  • Maple
    A047749 := proc(m) if m mod 2 = 1 then x := (m-1)/2; RETURN((3*x+1)!/((x+1)!*(2*x+1)!)); fi; x := m/2; RETURN((3*x)!/(x!*(2*x+1)!)); end;
    A047749 := proc(m) local x; if m mod 2 = 1 then x := (m-1)/2; RETURN((3*x+1)!/((x+1)!*(2*x+1)!)); fi; RETURN(A001764(m/2)); end;
  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], SeriesCoefficient[ InverseSeries[ Series[ (x + 2 x^2) / (1 + x)^3, {x, 0, n}]], {x, 0, n}]]; (* Michael Somos, Oct 29 2014 *)
    Table[If[OddQ[n],2Binomial[(3n-1)/2,(n-1)/2],Binomial[3n/2,n/2]]/(n+1),{n,0,40}] (* Robert A. Russell, Jan 19 2024 *)
  • PARI
    {a(n)=local(A=1+x);for(i=1,n,A=1+x*A^2*subst(A,x,-x+x*O(x^n)));polcoeff(A,n)} \\ Paul D. Hanna, Sep 20 2009
    
  • PARI
    x='x+O('x^66);
    C(x)=serreverse(x-x^3); /* =x+x^3+3*x^5+12*x^7+55*x^9 +..., cf. A001764 */
    s=1/(1-C(x)); /* g.f. */
    Vec(s) /* Joerg Arndt, Apr 16 2011 */
    
  • PARI
    apr(n, p, r) = r*binomial(n*p+r, n)/(n*p+r);
    a(n) = apr(n\2, 3, n%2+1); \\ Seiichi Manyama, Jul 20 2025
    
  • Python
    from math import comb
    def A047749(n): return comb(n+(a:=n>>1),a+(b:=n&1))//(n+1-b) # Chai Wah Wu, Jul 30 2022
  • Sage
    def A047749_list(n) :
        D = [0]*n; D[1] = 1
        R = []; b = False; h = 1
        for i in range(n) :
            for k in (1..h) :
                D[k] = D[k] + D[k-1]
            R.append(D[h])
            if b : h += 1
            b = not b
        return R
    A047749_list(35) # Peter Luschny, May 03 2012
    
  • Sage
    [1]+[((1+(-1)^n)*binomial(3*n/2,n/2)/(n+1) + (1-(-1)^n)* binomial((3*n-1)/2, (n+1)/2)/n)/2 for n in (1..35)] # G. C. Greubel, Jul 07 2019
    

Formula

G.f. is 1+Z, where Z satisfies x*Z^3 + (3*x-2)*Z^2 + (3*x-1)*Z + x = 0. Equivalently, the g.f. Y satisfies x*Y^3 - 2*Y^2 + 3*Y - 1 = 0. - Vladeta Jovovic, Dec 06 2002
Reversion of g.f. (x-2*x^2)/(1-x)^3 (ignoring signs). - Ralf Stephan, Mar 22 2004
G.f.: (4/(3*x))*(sin((1/3)*asin(sqrt(27*x^2/4))))^2 +(2/sqrt(3*x^2))*sin((1/3)*asin(sqrt(27*x^2/4))). - Paul Barry, Nov 08 2006
G.f.: 1/(1-2*sin(asin(3*sqrt(3)*x/2)/3)/sqrt(3)). - Paul Barry, Apr 16 2008
From Paul D. Hanna, Sep 20 2009: (Start)
G.f. satisfies: A(x) = 1 + x*A(x)^2*A(-x);
also, A(x)*A(-x) = B(x^2) where B(x) = 1 + x*B(x)^3 = g.f. of A001764. (End)
G.f.: 1/(1-C(x)) where C(x) = Reverse(x-x^3) = x + x^3 + 3*x^5 + 12*x^7 + 55*x^9 + ... (cf. A001764). - Joerg Arndt, Apr 16 2011
G.f.: G(z^2)+z*G(z^2)^2, where G(z) = 1 + z*G(z)^3, the generating function for A001764. - Robert A. Russell, Jan 26 2024
From Gary W. Adamson, Jul 14 2011: (Start)
a(n) is the upper left term in M^n, M = the infinite square production matrix:
1, 1, 0, 0, 0, 0, ...
0, 0, 1, 0, 0, 0, ...
1, 1, 0, 1, 0, 0, ...
0, 0, 1, 0, 1, 0, ...
1, 1, 0, 1, 0, 1, ...
... (End)
Conjecture D-finite with recurrence: 8*n*(n+1)*a(n) + 36*n*(n-2)*a(n-1) - 6*(9*n^2-18*n+14)*a(n-2) - 27*(3*n-7)*(3*n-8)*a(n-3) = 0. - R. J. Mathar, Dec 19 2011
0 = a(n)*(+7308954*a(n+2) - 16659999*a(n+3) - 4854519*a(n+4) + 6201838*a(n+5)) + a(n+1)*(-406053*a(n+2) - 1627560*a(n+3) + 1683538*a(n+4) - 245747*a(n+5)) + a(n+2)*(+45117*a(n+2) + 235870*a(n+3) + 173953*a(n+4) - 484295*a(n+5)) + a(n+3)*(-41820*a(n+3) - 50184*a(n+4) + 22304*a(n+5)) for all n in Z if a(-1) = -2/3. - Michael Somos, Oct 29 2014
a(0) = 1; a(n) = Sum_{i=0..n-1} Sum_{j=0..n-i-1} (-1)^i * a(i) * a(j) * a(n-i-j-1). - Ilya Gutkovskiy, Jul 28 2021
a(n) = binomial(A032766(n),floor((n+1)/2))/(2*floor(n/2)+1). - Miko Labalan, Nov 28 2023
a(n) = 2*A005036(n) - A005034(n) = A005034(n) - 2*A369315(n) = A005036(n) - A369315(n). - Robert A. Russell, Jan 20 2024
From Robert A. Russell, Mar 20 2024: (Start)
a(n) = U(n) in the Beineke and Pippert link.
G.f.: E(1)(t*E(3)(t^2)) (second entry in Table 1), where E(d)(t) is defined in formula 3 of Hering link. (End)
From Robert A. Russell, Jul 15 2024: (Start)
a(2m) = A001764(m) ~ (3^3/2^2)^m*sqrt(3/(2*Pi*(2*m)^3)).
a(n+2)/a(n) ~ 27/4; a(2m+1)/a(2m) ~ 3; a(2m)/a(2m-1) ~ 9/4. (End)
a(n) ~ 3^((6n+3)/4)/(sqrt(Pi)*2^((2n-1)/2)*(2n+1)^(3/2)). - Miko Labalan, Dec 05 2024
a(0) = 1; a(n) = Sum_{k=0..floor((n-1)/2)} a(2*k) * a(n-1-2*k). - Seiichi Manyama, Jul 07 2025

A001683 Number of one-sided triangulations of the disk; or flexagons of order n; or unlabeled plane trivalent trees (n-2 internal vertices, all of degree 3 and hence n leaves).

Original entry on oeis.org

1, 1, 1, 1, 4, 6, 19, 49, 150, 442, 1424, 4522, 14924, 49536, 167367, 570285, 1965058, 6823410, 23884366, 84155478, 298377508, 1063750740, 3811803164, 13722384546, 49611801980, 180072089896, 655977266884, 2397708652276, 8791599732140, 32330394085528
Offset: 2

Views

Author

Keywords

Comments

a(n) is the number of triangulations of an n-gon (equivalently, the number of vertices of the (n - 3)-dimensional associahedron) modulo the cyclic action [Bowman and Regev]. - N. J. A. Sloane, Dec 29 2012
a(n) is also the number of non-isomorphic cluster-tilted algebras of type A_(n-3), for n greater than or equal to 5. Equivalently it is the number of non-isomorphic quivers in the mutation class of any quiver with underlying graph A_(n-3) for n greater than or equal to 5. - Hermund A. Torkildsen (hermunda(AT)math.ntnu.no), Aug 06 2008
Number of oriented polyominoes composed of n-2 triangular cells of the hyperbolic regular tiling with Schläfli symbol {3,oo}. A stereographic projection of this tiling on the Poincaré disk can be obtained via the Christensson link. For oriented polyominoes, chiral pairs are counted as two. - Robert A. Russell, Jan 20 2024

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 k=3 of A295224.
A row or column of the array in A262586.
Polyominoes: A000207 (unoriented), A369314 (chiral), A208355(n-1) (achiral), A005034 {4,oo}, A007173 {3,3,oo}.

Programs

  • Maple
    C := n->binomial(2*n,n)/(n+1); c := x->if whattype(x) = integer then C(x) else 0; fi; A001683 := n->C(n-2)/n + c(n/2-1)/2+(2/3)*c(n/3-1);
  • Mathematica
    p=3; Table[Binomial[(p-1)n, n]/(((p-2)n+1)((p-2)n+2)) +If[OddQ[n], 0, Binomial[(p-1)n/2, n/2]/((p-2)n+2)]+Plus @@ Map[EulerPhi[ # ]Binomial[((p-1)n+1)/#, (n-1)/# ]/((p-1)n+1)&, Complement[Divisors[GCD[p, n-1]], {1}]], {n, 0, 20}] (* Robert A. Russell, Dec 11 2004 *)
    Rest[Rest[CoefficientList[Series[(6 + (1 - 4 x)^(3/2) + 6 x - 3(1 - 4 x^2)^(1/2) - 4 (1 - 4 x^3)^(1/2))/12, {x, 0, 33}], x]]] (* Vincenzo Librandi, Nov 25 2015 *)
  • PARI
    Cat(n)=if(n==floor(n),return(binomial(2*n,n)/(n+1)));0
    for(n=2,100,print1(Cat(n-2)/n+Cat(n/2-1)/2+(2/3)*Cat(n/3-1),", ")) \\ Derek Orr, Feb 26 2017

Formula

a(n) = C(n-2)/n + C(n/2-1)/2 + (2/3)*C(n/3-1), where C(n) = Catalan(n) (A000108) and terms are omitted if their subscripts are not integers.
G.f.: (6 + (1 - 4*x)^(3/2) + 6*x - 3*(1 - 4*x^2)^(1/2) - 4*(1 - 4*x^3)^(1/2))/12. - David Callan, Aug 01 2004
a(n) ~ 2^(2*n-4) / (sqrt(Pi) * n^(5/2)). - Vaclav Kotesovec, Mar 13 2016
a(n+2) = A000207(n) + A369314(n) = 2*A000207(n) - A208355(n-1) = 2*A369314(n) + A208355(n-1). - Robert A. Russell, Jan 19 2024
G.f.: z^2 * (4*G(z) - G(z)^2 + 3*G(z^2) + 4*z*G(z^3)) / 6, where G(z) = 1 + z*G(z)^2 is the g.f. for A000108. - Robert A. Russell, Apr 06 2024

A005036 Number of nonequivalent dissections of a polygon into n quadrilaterals by nonintersecting diagonals up to rotation and reflection.

Original entry on oeis.org

1, 1, 2, 5, 16, 60, 261, 1243, 6257, 32721, 175760, 963900, 5374400, 30385256, 173837631, 1004867079, 5861610475, 34469014515, 204161960310, 1217145238485, 7299007647552, 44005602441840
Offset: 1

Views

Author

Keywords

Comments

Closed formula is given in my paper linked below. - Nikos Apostolakis, Aug 01 2018
Number of unoriented polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A stereographic projection of this tiling on the Poincaré disk can be obtained via the Christensson link. For unoriented polyominoes, chiral pairs are counted as one. - Robert A. Russell, Jan 20 2024

References

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

Crossrefs

Column k=4 of A295260.
Polyominoes: A005034 (oriented), A369315 (chiral), A047749 (achiral), A385149 (asymmetric), A001764 (rooted), A000207 {3,oo}, A005040 {5,oo}, A005038 {5,oo} (oriented).

Programs

  • Mathematica
    p=4; Table[(Binomial[(p-1)n, n]/(((p-2)n+1)((p-2)n+2)) + If[OddQ[n], If[OddQ[p], Binomial[(p-1)n/2, (n-1)/2]/n, (p+1)Binomial[((p-1)n-1)/2, (n-1)/2]/((p-2)n+2)], 3Binomial[(p-1)n/2, n/2]/((p-2)n+2)]+Plus @@ Map[EulerPhi[ # ]Binomial[((p-1)n+1)/#, (n-1)/# ]/((p-1)n+1)&, Complement[Divisors[GCD[p, n-1]], {1, 2}]])/2, {n, 1, 20}] (* Robert A. Russell, Dec 11 2004 *)
    Table[(3Binomial[3n,n]/(2n+1)-Binomial[3n+1,n]/(n+1)-If[OddQ[n],-10Binomial[(3n-1)/2,(n-1)/2]-If[1==Mod[n,4],4Binomial[(3n-3)/4,(n-1)/4],0],-6Binomial[3n/2,n/2]]/(n+1))/8,{n,0,30}] (* Robert A. Russell, Jun 19 2025 *)

Formula

a(n) ~ 3^(3*n + 1/2) / (sqrt(Pi) * n^(5/2) * 2^(2*n + 4)). - Vaclav Kotesovec, Mar 13 2016
a(n) = A005034(n) - A369315(n) = (A005034(n) + A047749(n)) / 2 = A369315(n) + A047749(n). - Robert A. Russell, Jan 19 2024
G.f.: (3*G(z) - G(z)^2 + 6*G(z^2) + 5z*G(z^2)^2 + 2z*G(z^4)) / 8, where G(z)=1+z*G(z)^3 is the g.f. for A001764. - Robert A. Russell, Jun 19 2025

Extensions

More terms from Sascha Kurz, Oct 13 2001
Name edited by Andrew Howroyd, Nov 20 2017

A216371 Odd primes with one coach: primes p such that A135303((p-1)/2) = 1.

Original entry on oeis.org

3, 5, 7, 11, 13, 19, 23, 29, 37, 47, 53, 59, 61, 67, 71, 79, 83, 101, 103, 107, 131, 139, 149, 163, 167, 173, 179, 181, 191, 197, 199, 211, 227, 239, 263, 269, 271, 293, 311, 317, 347, 349, 359, 367, 373, 379, 383, 389, 419, 421, 443, 461, 463, 467, 479, 487
Offset: 1

Views

Author

Gary W. Adamson, Sep 05 2012

Keywords

Comments

Given that prime p has only one coach, the corresponding value of k in A003558 must be (p-1)/2, and vice versa. Using the Coach theorem of Jean Pedersen et al., phi(b) = 2 * c * k, with b odd. Let b = p, prime. Then phi(p) = (p-1), and k must be (p-1)/2 iff c = 1. Or, phi(p) = (p-1) = 2 * 1 * (p-1)/2.
Conjecture relating to odd integers: iff an integer is in the set A216371 and is either of the form 4q - 1 or 4q + 1, (q>0); then the top row of its coach (cf. A003558) is composed of a permutation of the first q odd integers. Examples: 11 is of the form 4q - 1, q = 3; with the top row of its coach [1, 5, 3]. 13 is of the form 4q + 1, q = 3; so has a coach of [1, 3, 5]. 37 is of the form 4q + 1, q = 9; so has a coach with the top row composed of a permutation of the first 9 odd integers: [1, 9, 7, 15, 11, 13, 3, 17, 5]. - Gary W. Adamson, Sep 08 2012
Odd primes p such that 2^m is not congruent to 1 or -1 (mod p) for 0 < m < (p-1)/2. - Charles R Greathouse IV, Sep 15 2012
These are also the odd primes a(n) for which there is only one periodic Schick sequence (see the reference, and also the Brändli and Beyne link, eq. (2) for the recurrence but using various inputs. See also a comment in A332439). This sequence has primitive period length (named pes in Schick's book) A003558((a(n)-1)/2) = A005034(a(n)) = A000010(a(n))/2 = (a(n) - 1)/2, for n >= 1. - Wolfdieter Lang, Apr 09 2020
From Jianing Song, Dec 24 2022: (Start)
Primes p such that the multiplicative order of 4 modulo p is (p-1)/2. Proof of equivalence: let ord(a,k) be the multiplicative of a modulo k.
If 2^m is not 1 or -1 (mod p) for 0 < m < (p-1)/2, then ord(2,p) is either p-1 or (p-1)/2. If ord(2,p) = p-1, then ord(4,p) = (p-1)/2. If ord(2,p) = (p-1)/2, then p == 3 (mod 4), otherwise 2^((p-1)/4) == -1 (mod p), so ord(4,p) = (p-1)/2.
Conversely, if ord(4,p) = (p-1)/2, then ord(2,p) = p-1, or ord(2,p) = (p-1)/2 and p == 3 (mod 4) (otherwise ord(4,p) = (p-1)/4). In the first case, (p-1)/2 is the smallest m > 0 such that 2^m == +-1 (mod p); in the second case, since (p-1)/2 is odd, 2^m == -1 (mod p) has no solution. In either case, so 2^m is not 1 or -1 (mod p) for 0 < m < (p-1)/2.
{(a(n)-1)/2} is the sequence of indices of fixed points of A053447.
A prime p is a term if and only if one of the two following conditions holds: (a) 2 is a primitive root modulo p; (b) p == 3 (mod 4), and the multiplicative order of 2 modulo p is (p-1)/2 (in this case, we have p == 7 (mod 8) since 2 is a quadratic residue modulo p). (End)
From Jianing Song, Aug 11 2023: (Start)
Primes p such that 2 or -2 (or both) is a primitive root modulo p. Proof of equivalence: if ord(2,p) = p-1, then clearly ord(4,p) = (p-1)/2. If ord(-2,p) = p-1, then we also have ord(4,p) = (p-1)/2. Conversely, suppose that ord(4,p) = (p-1)/2, then ord(2,p) = p-1 or (p-1)/2, and ord(-2,p) = p-1 or (p-1)/2. If ord(2,p) = ord(-2,p) = (p-1)/2, then we have that (p-1)/2 is odd and (-1)^((p-1)/2) == 1 (mod p), a contradiction.
A prime p is a term if and only if one of the two following conditions holds: (a) -2 is a primitive root modulo p; (b) p == 3 (mod 4), and the multiplicative order of -2 modulo p is (p-1)/2 (in this case, we have p == 3 (mod 8) since -2 is a quadratic residue modulo p). (End)
No terms are congruent to 1 modulo 8, since otherwise we would have 4^((p-1)/4) = (+-2)^((p-1)/2) == 1 (mod p). - Jianing Song, May 14 2024
The n-th prime A000040(n) is a term iff A376010(n) = 2. - Max Alekseyev, Sep 05 2024

Examples

			Prime 23 has a k value of 11 = (23 - 1)/2 (Cf. A003558(11)). It follows that 23 has only one coach (A135303(11) = 1). 23 is thus in the set. On the other hand 31 is not in the set since A135303(15) shows 3 coaches, with A003558(15) = 5.
13 is in the set since A135303(6) = 1; but 17 isn't since A135303(8) = 2.
		

References

  • P. Hilton and J. Pedersen, A Mathematical Tapestry, Demonstrating the Beautiful Unity of Mathematics, 2010, Cambridge University Press, pages 260-264.
  • Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6). Tables 3.1 to 3.10, for odd p = 3..113 (with gaps), pp. 158-166.

Crossrefs

Union of A001122 and A105874.
A105876 is the subsequence of terms congruent to 3 modulo 4.
Complement of A268923 in the set of odd primes.
Cf. A082654 (order of 4 mod n-th prime), A000010, A000040, A003558, A005034, A053447, A054639, A135303, A364867, A376010.

Programs

  • Maple
    isA216371 := proc(n)
        if isprime(n) then
            if A135303((n-1)/2) = 1 then
                true;
            else
                false;
            end if;
        else
            false;
        end if;
    end proc:
    A216371 := proc(n)
        local p;
        if n = 1 then
            3;
        else
            p := nextprime(procname(n-1)) ;
            while true do
                if isA216371(p) then
                    return p;
                end if;
                p := nextprime(p) ;
            end do:
        end if;
    end proc:
    seq(A216371(n),n=1..40) ; # R. J. Mathar, Dec 01 2014
  • Mathematica
    Suborder[a_, n_] := If[n > 1 && GCD[a, n] == 1, Min[MultiplicativeOrder[a, n, {-1, 1}]], 0]; nn = 150; Select[Prime[Range[2, nn]], EulerPhi[#]/(2*Suborder[2, #]) == 1 &] (* T. D. Noe, Sep 18 2012 *)
    f[p_] := Sum[Cos[2^n Pi/((2 p + 1))], {n, p}]; 1 + 2 * Select[Range[500], Reduce[f[#] == -1/2, Rationals] &]; (* Gerry Martens, May 01 2016 *)
  • PARI
    is(p)=for(m=1,p\2-1, if(abs(centerlift(Mod(2,p)^m))==1, return(0))); p>2 && isprime(p) \\ Charles R Greathouse IV, Sep 18 2012
    
  • PARI
    is(p) = isprime(p) && (p>2) && znorder(Mod(4,p)) == (p-1)/2 \\ Jianing Song, Dec 24 2022

Formula

a(n) = 2*A054639(n) + 1. - L. Edson Jeffery, Dec 18 2012

A295224 Array read by antidiagonals: T(n,k) = number of nonequivalent dissections of a polygon into n k-gons by nonintersecting diagonals up to rotation (k >= 3).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 4, 1, 1, 2, 7, 6, 1, 1, 3, 12, 25, 19, 1, 1, 3, 19, 57, 108, 49, 1, 1, 4, 26, 118, 366, 492, 150, 1, 1, 4, 35, 203, 931, 2340, 2431, 442, 1, 1, 5, 46, 332, 1989, 7756, 16252, 12371, 1424, 1, 1, 5, 57, 494, 3766, 20254, 68685, 115940, 65169, 4522
Offset: 1

Views

Author

Andrew Howroyd, Nov 17 2017

Keywords

Comments

The polygon prior to dissection will have n*(k-2)+2 sides.
In the Harary, Palmer and Read reference these are the sequences called H.
T(n,k) is the number of oriented polyominoes containing n k-gonal tiles of the hyperbolic regular tiling with Schläfli symbol {k,oo}. Stereographic projections of several of these tilings on the Poincaré disk can be obtained via the Christensson link. For oriented polyominoes, chiral pairs are counted as two. T(n,2) could represent polyominoes of the Euclidean regular tiling with Schläfli symbol {2,oo}; T(n,2) = 1. - Robert A. Russell, Jan 21 2024

Examples

			Array begins:
  =====================================================
  n\k|    3     4      5       6        7         8
  ---|-------------------------------------------------
   1 |    1     1      1       1        1         1 ...
   2 |    1     1      1       1        1         1 ...
   3 |    1     2      2       3        3         4 ...
   4 |    4     7     12      19       26        35 ...
   5 |    6    25     57     118      203       332 ...
   6 |   19   108    366     931     1989      3766 ...
   7 |   49   492   2340    7756    20254     45448 ...
   8 |  150  2431  16252   68685   219388    580203 ...
   9 |  442 12371 115940  630465  2459730   7684881 ...
  10 | 1424 65169 854981 5966610 28431861 104898024 ...
  ...
		

Crossrefs

Columns k=3..6 are A001683(n+2), A005034, A005038, A221184(n-1).
Polyominoes: A295260 (unoriented), A070914 (rooted).

Programs

  • Mathematica
    u[n_, k_, r_] := r*Binomial[(k - 1)*n + r, n]/((k - 1)*n + r);
    T[n_, k_] := u[n, k, 1] + (If[EvenQ[n], u[n/2, k, 1], 0] - u[n, k, 2])/2 + DivisorSum[GCD[n - 1, k], EulerPhi[#]*u[(n - 1)/#, k, k/#]&]/k;
    Table[T[n - k + 1, k], {n, 1, 13}, {k, n, 3, -1}] // Flatten (* Jean-François Alcover, Nov 21 2017, after Andrew Howroyd *)
  • PARI
    \\ here u is Fuss-Catalan sequence with p = k+1.
    u(n, k, r)={r*binomial((k - 1)*n + r, n)/((k - 1)*n + r)}
    T(n,k) = u(n,k,1) + (if(n%2==0, u(n/2,k,1))-u(n,k,2))/2 + sumdiv(gcd(n-1,k), d, eulerphi(d)*u((n-1)/d,k,k/d))/k;
    for(n=1, 10, for(k=3, 8, print1(T(n, k), ", ")); print);
    
  • Python
    from sympy import binomial, gcd, totient, divisors
    def u(n, k, r): return r*binomial((k - 1)*n + r, n)//((k - 1)*n + r)
    def T(n, k): return u(n, k, 1) + ((u(n//2, k, 1) if n%2==0 else 0) - u(n, k, 2))//2 + sum([totient(d)*u((n - 1)//d, k, k//d) for d in divisors(gcd(n - 1, k))])//k
    for n in range(1, 11): print([T(n, k) for k in range(3, 9)]) # Indranil Ghosh, Dec 13 2017, after PARI code

Formula

T(n,k) ~ A295222(n,k)/n for fixed k.

A005038 Number of nonequivalent dissections of a polygon into n pentagons by nonintersecting diagonals up to rotation.

Original entry on oeis.org

1, 1, 2, 12, 57, 366, 2340, 16252, 115940, 854981, 6444826, 49554420, 387203390, 3068067060, 24604111560, 199398960212, 1631041938108, 13451978877748, 111765327780200, 934774244822704, 7865200653146905
Offset: 1

Views

Author

Keywords

Comments

Also, with a different offset, number of colored quivers in the 3-mutation class of a quiver of Dynkin type A_n. - N. J. A. Sloane, Jan 22 2013
Number of oriented polyominoes composed of n pentagonal cells of the hyperbolic regular tiling with Schläfli symbol {5,oo}. A stereographic projection of this tiling on the Poincaré disk can be obtained via the Christensson link. For oriented polyominoes, chiral pairs are counted as two. - Robert A. Russell, Jan 23 2024

References

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

Crossrefs

Column k=5 of A295224.
Polyominoes: A005040 (unoriented), A369471 (chiral), A369472 (achiral), A001683(n+2) {3,oo}, A005034 {4,oo}, A221184{n-1} {6,oo}.

Programs

  • Mathematica
    p=5; Table[Binomial[(p-1)n, n]/(((p-2)n+1)((p-2)n+2)) +If[OddQ[n], 0, Binomial[(p-1)n/2, n/2]/((p-2)n+2)]+Plus @@ Map[EulerPhi[ # ]Binomial[((p-1)n+1)/#, (n-1)/# ]/((p-1)n+1)&, Complement[Divisors[GCD[p, n-1]], {1}]], {n, 1, 20}] (* Robert A. Russell, Dec 11 2004 *)

Formula

a(n) ~ 2^(8*n + 1/2) / (sqrt(Pi) * n^(5/2) * 3^(3*n + 5/2)). - Vaclav Kotesovec, Mar 13 2016
a(n) = A005040(n) + A369471(n) = 2*A005040(n) - A369472(n) = 2*A369471(n) + A369472(n). - Robert A. Russell, Jan 23 2024

Extensions

a(21) corrected by Sean A. Irvine, Mar 11 2016
Name edited by Andrew Howroyd, Nov 20 2017

A003455 Number of nonequivalent dissections of an n-gon by nonintersecting diagonals up to rotation.

Original entry on oeis.org

1, 2, 3, 11, 29, 122, 479, 2113, 9369, 43392, 203595, 975563, 4736005, 23296394, 115811855, 581324861, 2942579633, 15008044522, 77064865555, 398150807179, 2068470765261, 10800665952376, 56658467018647, 298489772155137, 1578702640556193
Offset: 3

Views

Author

Keywords

Comments

Total number of dissections of an n-gon into polygons without reflection. - Sean A. Irvine, May 15 2015

References

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

Crossrefs

Programs

  • PARI
    \\ See A295495 for DissectionsModCyclic().
    DissectionsModCyclic(apply(v->1, [1..30])) \\ Andrew Howroyd, Nov 22 2017

Extensions

More terms from Sean A. Irvine, May 15 2015
Name clarified by Andrew Howroyd, Nov 22 2017

A221184 Number of colored quivers in the 4-mutation class of a quiver of Dynkin type A_n.

Original entry on oeis.org

1, 1, 3, 19, 118, 931, 7756, 68685, 630465, 5966610, 57805410, 571178751, 5737638778, 58455577800, 602859152496, 6283968796705, 66119469155523, 701526880303315, 7498841128986109, 80696081185766970, 873654669882574580
Offset: 0

Views

Author

N. J. A. Sloane, Jan 22 2013

Keywords

Comments

Also, number of nonequivalent dissections of a polygon into n+1 hexagons by nonintersecting diagonals up to rotation. - Andrew Howroyd, Nov 20 2017
Number of oriented polyominoes composed of n+1 hexagonal cells of the hyperbolic regular tiling with Schläfli symbol {6,oo}. A stereographic projection of this tiling on the Poincaré disk can be obtained via the Christensson link. For oriented polyominoes, chiral pairs are counted as two. - Robert A. Russell, Jan 23 2024

Crossrefs

Column k=6 of A295224.
Polyominoes: A004127 (unoriented), A369473 (chiral), A143546 (achiral), A001683(n+2) {3,oo}, A005034 {4,oo}, A005038 {5,oo}.

Programs

  • Mathematica
    u[n_, k_, r_] := r*Binomial[(k - 1)*n + r, n]/((k - 1)*n + r);
    T[n_, k_] := u[n, k, 1] + (If[EvenQ[n], u[n/2, k, 1], 0] - u[n, k, 2])/2 + DivisorSum[GCD[n - 1, k], EulerPhi[#]*u[(n - 1)/#, k, k/#] &]/k;
    a[n_] := T[n + 1, 6];
    Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jun 14 2018, after Andrew Howroyd *)
    p=6; Table[Binomial[(p-1)n,n]/(((p-2)n+1)((p-2)n+2))+If[OddQ[n],0,Binomial[(p-1)n/2,n/2]/((p-2)n+2)]+DivisorSum[GCD[p,n-1],EulerPhi[#]Binomial[((p-1)n+1)/#,(n-1)/#]/((p-1)n+1)&,#>1&],{n,30}] (* Robert A. Russell, Jan 23 2024 *)

Formula

a(n) ~ 5^(5*n + 11/2) / (sqrt(Pi) * n^(5/2) * 2^(8*n + 27/2)). - Vaclav Kotesovec, Jun 15 2018
a(n-1) = A004127(n) + A369473(n) = 2*A004127(n) - A143546(n) = 2*A369473(n) + A143546(n). - Robert A. Russell, Jan 23 2024

Extensions

a(0)=1 and a(18)-a(20) corrected by Andrew Howroyd, Nov 20 2017

A369315 Number of chiral pairs of polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}.

Original entry on oeis.org

2, 9, 48, 231, 1188, 6114, 32448, 175032, 962472, 5370524, 30377504, 173816313, 1004823816, 5861490300, 34468767840, 204161269620, 1217143807770, 7299003615537, 44005594027200, 266608363362900
Offset: 4

Views

Author

Robert A. Russell, Jan 19 2024

Keywords

Comments

A stereographic projection of the {4,oo} tiling on the Poincaré disk can be obtained via the Christensson link. Each member of a chiral pair is a reflection but not a rotation of the other.

Examples

			 __ __ __    __ __ __     __ __          __ __
|__|__|__|  |__|__|__|   |__|__|__    __|__|__|  a(4) = 2.
      |__|  |__|            |__|__|  |__|__|
		

Crossrefs

Polyominoes: A005034 (oriented), A005036 (unoriented), A047749 (achiral), A385149 (asymmetric), A001764 (rooted), A369314 {3,oo}.

Programs

  • Mathematica
    p=4; Table[(Binomial[(p-1)n,n]/(((p-2)n+1)((p-2)n+2))-If[OddQ[n],If[OddQ[p],Binomial[(p-1)n/2,(n-1)/2]/n,(p+1)Binomial[((p-1)n-1)/2,(n-1)/2]/((p-2)n+2)-Binomial[((p-1)n+1)/2,(n-1)/2]/((p-1)n+1)],Binomial[(p-1)n/2,n/2]/((p-2)n+2)]+DivisorSum[GCD[p,n-1],EulerPhi[#]Binomial[((p-1)n+1)/#,(n-1)/#]/((p-1)n+1)&,#>1&])/2,{n,4,30}]
    Table[(3Binomial[3n,n]/(2n+1)-Binomial[3n+1,n]/(n+1)-If[OddQ[n],6Binomial[(3n-1)/2,(n-1)/2]-If[1==Mod[n,4],4Binomial[(3n-3)/4,(n-1)/4],0],2Binomial[3n/2,n/2]]/(n+1))/8,{n,0,30}] (* Robert A. Russell, Jun 19 2025 *)

Formula

a(n) = A005034(n) - A005036(n) = (A005034(n) - A047749(n)) / 2 = A005036(n) - A047749(n).
G.f.: (3*G(z) - G(z)^2 - 2*G(z^2) - 3z*G(z^2)^2 + 2z*G(z^4)) / 8, where G(z)=1+z*G(z)^3 is the g.f. for A001764. - Robert A. Russell, Jun 19 2025
Showing 1-10 of 13 results. Next