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

A001147 Double factorial of odd numbers: a(n) = (2*n-1)!! = 1*3*5*...*(2*n-1).

Original entry on oeis.org

1, 1, 3, 15, 105, 945, 10395, 135135, 2027025, 34459425, 654729075, 13749310575, 316234143225, 7905853580625, 213458046676875, 6190283353629375, 191898783962510625, 6332659870762850625, 221643095476699771875, 8200794532637891559375, 319830986772877770815625
Offset: 0

Views

Author

Keywords

Comments

The solution to Schröder's third problem.
Number of fixed-point-free involutions in symmetric group S_{2n} (cf. A000085).
a(n-2) is the number of full Steiner topologies on n points with n-2 Steiner points. [corrected by Lyle Ramshaw, Jul 20 2022]
a(n) is also the number of perfect matchings in the complete graph K(2n). - Ola Veshta (olaveshta(AT)my-deja.com), Mar 25 2001
Number of ways to choose n disjoint pairs of items from 2*n items. - Ron Zeno (rzeno(AT)hotmail.com), Feb 06 2002
Number of ways to choose n-1 disjoint pairs of items from 2*n-1 items (one item remains unpaired). - Bartosz Zoltak, Oct 16 2012
For n >= 1 a(n) is the number of permutations in the symmetric group S_(2n) whose cycle decomposition is a product of n disjoint transpositions. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 21 2001
a(n) is the number of distinct products of n+1 variables with commutative, nonassociative multiplication. - Andrew Walters (awalters3(AT)yahoo.com), Jan 17 2004. For example, a(3)=15 because the product of the four variables w, x, y and z can be constructed in exactly 15 ways, assuming commutativity but not associativity: 1. w(x(yz)) 2. w(y(xz)) 3. w(z(xy)) 4. x(w(yz)) 5. x(y(wz)) 6. x(z(wy)) 7. y(w(xz)) 8. y(x(wz)) 9. y(z(wx)) 10. z(w(xy)) 11. z(x(wy)) 12. z(y(wx)) 13. (wx)(yz) 14. (wy)(xz) 15. (wz)(xy).
a(n) = E(X^(2n)), where X is a standard normal random variable (i.e., X is normal with mean = 0, variance = 1). So for instance a(3) = E(X^6) = 15, etc. See Abramowitz and Stegun or Hoel, Port and Stone. - Jerome Coleman, Apr 06 2004
Second Eulerian transform of 1,1,1,1,1,1,... The second Eulerian transform transforms a sequence s to a sequence t by the formula t(n) = Sum_{k=0..n} E(n,k)s(k), where E(n,k) is a second-order Eulerian number (A008517). - Ross La Haye, Feb 13 2005
Integral representation as n-th moment of a positive function on the positive axis: a(n) = Integral_{x=0..oo} x^n*exp(-x/2)/sqrt(2*Pi*x) dx, n >= 0. - Karol A. Penson, Oct 10 2005
a(n) is the number of binary total partitions of n+1 (each non-singleton block must be partitioned into exactly two blocks) or, equivalently, the number of unordered full binary trees with n+1 labeled leaves (Stanley, ex 5.2.6). - Mitch Harris, Aug 01 2006
a(n) is the Pfaffian of the skew-symmetric 2n X 2n matrix whose (i,j) entry is i for iDavid Callan, Sep 25 2006
a(n) is the number of increasing ordered rooted trees on n+1 vertices where "increasing" means the vertices are labeled 0,1,2,...,n so that each path from the root has increasing labels. Increasing unordered rooted trees are counted by the factorial numbers A000142. - David Callan, Oct 26 2006
Number of perfect multi Skolem-type sequences of order n. - Emeric Deutsch, Nov 24 2006
a(n) = total weight of all Dyck n-paths (A000108) when each path is weighted with the product of the heights of the terminal points of its upsteps. For example with n=3, the 5 Dyck 3-paths UUUDDD, UUDUDD, UUDDUD, UDUUDD, UDUDUD have weights 1*2*3=6, 1*2*2=4, 1*2*1=2, 1*1*2=2, 1*1*1=1 respectively and 6+4+2+2+1=15. Counting weights by height of last upstep yields A102625. - David Callan, Dec 29 2006
a(n) is the number of increasing ternary trees on n vertices. Increasing binary trees are counted by ordinary factorials (A000142) and increasing quaternary trees by triple factorials (A007559). - David Callan, Mar 30 2007
From Tom Copeland, Nov 13 2007, clarified in first and extended in second paragraph, Jun 12 2021: (Start)
a(n) has the e.g.f. (1-2x)^(-1/2) = 1 + x + 3*x^2/2! + ..., whose reciprocal is (1-2x)^(1/2) = 1 - x - x^2/2! - 3*x^3/3! - ... = b(0) - b(1)*x - b(2)*x^2/2! - ... with b(0) = 1 and b(n+1) = -a(n) otherwise. By the formalism of A133314, Sum_{k=0..n} binomial(n,k)*b(k)*a(n-k) = 0^n where 0^0 := 1. In this sense, the sequence a(n) is essentially self-inverse. See A132382 for an extension of this result. See A094638 for interpretations.
This sequence aerated has the e.g.f. e^(t^2/2) = 1 + t^2/2! + 3*t^4/4! + ... = c(0) + c(1)*t + c(2)*t^2/2! + ... and the reciprocal e^(-t^2/2); therefore, Sum_{k=0..n} cos(Pi k/2)*binomial(n,k)*c(k)*c(n-k) = 0^n; i.e., the aerated sequence is essentially self-inverse. Consequently, Sum_{k=0..n} (-1)^k*binomial(2n,2k)*a(k)*a(n-k) = 0^n. (End)
From Ross Drewe, Mar 16 2008: (Start)
This is also the number of ways of arranging the elements of n distinct pairs, assuming the order of elements is significant but the pairs are not distinguishable, i.e., arrangements which are the same after permutations of the labels are equivalent.
If this sequence and A000680 are denoted by a(n) and b(n) respectively, then a(n) = b(n)/n! where n! = the number of ways of permuting the pair labels.
For example, there are 90 ways of arranging the elements of 3 pairs [1 1], [2 2], [3 3] when the pairs are distinguishable: A = { [112233], [112323], ..., [332211] }.
By applying the 6 relabeling permutations to A, we can partition A into 90/6 = 15 subsets: B = { {[112233], [113322], [221133], [223311], [331122], [332211]}, {[112323], [113232], [221313], [223131], [331212], [332121]}, ....}
Each subset or equivalence class in B represents a unique pattern of pair relationships. For example, subset B1 above represents {3 disjoint pairs} and subset B2 represents {1 disjoint pair + 2 interleaved pairs}, with the order being significant (contrast A132101). (End)
A139541(n) = a(n) * a(2*n). - Reinhard Zumkeller, Apr 25 2008
a(n+1) = Sum_{j=0..n} A074060(n,j) * 2^j. - Tom Copeland, Sep 01 2008
From Emeric Deutsch, Jun 05 2009: (Start)
a(n) is the number of adjacent transpositions in all fixed-point-free involutions of {1,2,...,2n}. Example: a(2)=3 because in 2143=(12)(34), 3412=(13)(24), and 4321=(14)(23) we have 2 + 0 + 1 adjacent transpositions.
a(n) = Sum_{k>=0} k*A079267(n,k).
(End)
Hankel transform is A137592. - Paul Barry, Sep 18 2009
(1, 3, 15, 105, ...) = INVERT transform of A000698 starting (1, 2, 10, 74, ...). - Gary W. Adamson, Oct 21 2009
a(n) = (-1)^(n+1)*H(2*n,0), where H(n,x) is the probabilists' Hermite polynomial. The generating function for the probabilists' Hermite polynomials is as follows: exp(x*t-t^2/2) = Sum_{i>=0} H(i,x)*t^i/i!. - Leonid Bedratyuk, Oct 31 2009
The Hankel transform of a(n+1) is A168467. - Paul Barry, Dec 04 2009
Partial products of odd numbers. - Juri-Stepan Gerasimov, Oct 17 2010
See A094638 for connections to differential operators. - Tom Copeland, Sep 20 2011
a(n) is the number of subsets of {1,...,n^2} that contain exactly k elements from {1,...,k^2} for k=1,...,n. For example, a(3)=15 since there are 15 subsets of {1,2,...,9} that satisfy the conditions, namely, {1,2,5}, {1,2,6}, {1,2,7}, {1,2,8}, {1,2,9}, {1,3,5}, {1,3,6}, {1,3,7}, {1,3,8}, {1,3,9}, {1,4,5}, {1,4,6}, {1,4,7}, {1,4,8}, and {1,4,9}. - Dennis P. Walsh, Dec 02 2011
a(n) is the leading coefficient of the Bessel polynomial y_n(x) (cf. A001498). - Leonid Bedratyuk, Jun 01 2012
For n>0: a(n) is also the determinant of the symmetric n X n matrix M defined by M(i,j) = min(i,j)^2 for 1 <= i,j <= n. - Enrique Pérez Herrero, Jan 14 2013
a(n) is also the numerator of the mean value from 0 to Pi/2 of sin(x)^(2n). - Jean-François Alcover, Jun 13 2013
a(n) is the size of the Brauer monoid on 2n points (see A227545). - James Mitchell, Jul 28 2013
For n>1: a(n) is the numerator of M(n)/M(1) where the numbers M(i) have the property that M(n+1)/M(n) ~ n-1/2 (for example, large Kendell-Mann numbers, see A000140 or A181609, as n --> infinity). - Mikhail Gaichenkov, Jan 14 2014
a(n) = the number of upper-triangular matrix representations required for the symbolic representation of a first order central moment of the multivariate normal distribution of dimension 2(n-1), i.e., E[X_1*X_2...*X_(2n-2)|mu=0, Sigma]. See vignette for symmoments R package on CRAN and Phillips reference below. - Kem Phillips, Aug 10 2014
For n>1: a(n) is the number of Feynman diagrams of order 2n (number of internal vertices) for the vacuum polarization with one charged loop only, in quantum electrodynamics. - Robert Coquereaux, Sep 15 2014
Aerated with intervening zeros (1,0,1,0,3,...) = a(n) (cf. A123023), the e.g.f. is e^(t^2/2), so this is the base for the Appell sequence A099174 with e.g.f. e^(t^2/2) e^(x*t) = exp(P(.,x),t) = unsigned A066325(x,t), the probabilist's (or normalized) Hermite polynomials. P(n,x) = (a. + x)^n with (a.)^n = a_n and comprise the umbral compositional inverses for A066325(x,t) = exp(UP(.,x),t), i.e., UP(n,P(.,t)) = x^n = P(n,UP(.,t)), where UP(n,t) are the polynomials of A066325 and, e.g., (P(.,t))^n = P(n,t). - Tom Copeland, Nov 15 2014
a(n) = the number of relaxed compacted binary trees of right height at most one of size n. A relaxed compacted binary tree of size n is a directed acyclic graph consisting of a binary tree with n internal nodes, one leaf, and n pointers. It is constructed from a binary tree of size n, where the first leaf in a post-order traversal is kept and all other leaves are replaced by pointers. These links may point to any node that has already been visited by the post-order traversal. The right height is the maximal number of right-edges (or right children) on all paths from the root to any leaf after deleting all pointers. The number of unbounded relaxed compacted binary trees of size n is A082161(n). See the Genitrini et al. link. - Michael Wallner, Jun 20 2017
Also the number of distinct adjacency matrices in the n-ladder rung graph. - Eric W. Weisstein, Jul 22 2017
From Christopher J. Smyth, Jan 26 2018: (Start)
a(n) = the number of essentially different ways of writing a probability distribution taking n+1 values as a sum of products of binary probability distributions. See comment of Mitch Harris above. This is because each such way corresponds to a full binary tree with n+1 leaves, with the leaves labeled by the values. (This comment is due to Niko Brummer.)
Also the number of binary trees with root labeled by an (n+1)-set S, its n+1 leaves by the singleton subsets of S, and other nodes labeled by subsets T of S so that the two daughter nodes of the node labeled by T are labeled by the two parts of a 2-partition of T. This also follows from Mitch Harris' comment above, since the leaf labels determine the labels of the other vertices of the tree.
(End)
a(n) is the n-th moment of the chi-squared distribution with one degree of freedom (equivalent to Coleman's Apr 06 2004 comment). - Bryan R. Gillespie, Mar 07 2021
Let b(n) = 0 for n odd and b(2k) = a(k); i.e., let the sequence b(n) be an aerated version of this entry. After expanding the differential operator (x + D)^n and normal ordering the resulting terms, the integer coefficient of the term x^k D^m is n! b(n-k-m) / [(n-k-m)! k! m!] with 0 <= k,m <= n and (k+m) <= n. E.g., (x+D)^2 = x^2 + 2xD + D^2 + 1 with D = d/dx. The result generalizes to the raising (R) and lowering (L) operators of any Sheffer polynomial sequence by replacing x by R and D by L and follows from the disentangling relation e^{t(L+R)} = e^{t^2/2} e^{tR} e^{tL}. Consequently, these are also the coefficients of the reordered 2^n permutations of the binary symbols L and R under the condition LR = RL + 1. E.g., (L+R)^2 = LL + LR + RL + RR = LL + 2RL + RR + 1. (Cf. A344678.) - Tom Copeland, May 25 2021
From Tom Copeland, Jun 14 2021: (Start)
Lando and Zvonkin present several scenarios in which the double factorials occur in their role of enumerating perfect matchings (pairings) and as the nonzero moments of the Gaussian e^(x^2/2).
Speyer and Sturmfels (p. 6) state that the number of facets of the abstract simplicial complex known as the tropical Grassmannian G'''(2,n), the space of phylogenetic T_n trees (see A134991), or Whitehouse complex is a shifted double factorial.
These are also the unsigned coefficients of the x[2]^m terms in the partition polynomials of A134685 for compositional inversion of e.g.f.s, a refinement of A134991.
a(n)*2^n = A001813(n) and A001813(n)/(n+1)! = A000108(n), the Catalan numbers, the unsigned coefficients of the x[2]^m terms in the partition polynomials A133437 for compositional inversion of o.g.f.s, a refinement of A033282, A126216, and A086810. Then the double factorials inherit a multitude of analytic and combinatoric interpretations from those of the Catalan numbers, associahedra, and the noncrossing partitions of A134264 with the Catalan numbers as unsigned-row sums. (End)
Connections among the Catalan numbers A000108, the odd double factorials, values of the Riemann zeta function and its derivative for integer arguments, and series expansions of the reduced action for the simple harmonic oscillator and the arc length of the spiral of Archimedes are given in the MathOverflow post on the Riemann zeta function. - Tom Copeland, Oct 02 2021
b(n) = a(n) / (n! 2^n) = Sum_{k = 0..n} (-1)^n binomial(n,k) (-1)^k a(k) / (k! 2^k) = (1-b.)^n, umbrally; i.e., the normalized double factorial a(n) is self-inverse under the binomial transform. This can be proved by applying the Euler binomial transformation for o.g.f.s Sum_{n >= 0} (1-b.)^n x^n = (1/(1-x)) Sum_{n >= 0} b_n (x / (x-1))^n to the o.g.f. (1-x)^{-1/2} = Sum_{n >= 0} b_n x^n. Other proofs are suggested by the discussion in Watson on pages 104-5 of transformations of the Bessel functions of the first kind with b(n) = (-1)^n binomial(-1/2,n) = binomial(n-1/2,n) = (2n)! / (n! 2^n)^2. - Tom Copeland, Dec 10 2022

Examples

			a(3) = 1*3*5 = 15.
From _Joerg Arndt_, Sep 10 2013: (Start)
There are a(3)=15 involutions of 6 elements without fixed points:
  #:    permutation           transpositions
  01:  [ 1 0 3 2 5 4 ]      (0, 1) (2, 3) (4, 5)
  02:  [ 1 0 4 5 2 3 ]      (0, 1) (2, 4) (3, 5)
  03:  [ 1 0 5 4 3 2 ]      (0, 1) (2, 5) (3, 4)
  04:  [ 2 3 0 1 5 4 ]      (0, 2) (1, 3) (4, 5)
  05:  [ 2 4 0 5 1 3 ]      (0, 2) (1, 4) (3, 5)
  06:  [ 2 5 0 4 3 1 ]      (0, 2) (1, 5) (3, 4)
  07:  [ 3 2 1 0 5 4 ]      (0, 3) (1, 2) (4, 5)
  08:  [ 3 4 5 0 1 2 ]      (0, 3) (1, 4) (2, 5)
  09:  [ 3 5 4 0 2 1 ]      (0, 3) (1, 5) (2, 4)
  10:  [ 4 2 1 5 0 3 ]      (0, 4) (1, 2) (3, 5)
  11:  [ 4 3 5 1 0 2 ]      (0, 4) (1, 3) (2, 5)
  12:  [ 4 5 3 2 0 1 ]      (0, 4) (1, 5) (2, 3)
  13:  [ 5 2 1 4 3 0 ]      (0, 5) (1, 2) (3, 4)
  14:  [ 5 3 4 1 2 0 ]      (0, 5) (1, 3) (2, 4)
  15:  [ 5 4 3 2 1 0 ]      (0, 5) (1, 4) (2, 3)
(End)
G.f. = 1 + x + 3*x^2 + 15*x^3 + 105*x^4 + 945*x^5 + 10395*x^6 + 135135*x^7 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, (26.2.28).
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 317.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 228, #19.
  • Hoel, Port and Stone, Introduction to Probability Theory, Section 7.3.
  • F. K. Hwang, D. S. Richards and P. Winter, The Steiner Tree Problem, North-Holland, 1992, see p. 14.
  • C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill, 1980, pages 466-467.
  • 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).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.6 and also p. 178.
  • R. Vein and P. Dale, Determinants and Their Applications in Mathematical Physics, Springer-Verlag, New York, 1999, p. 73.
  • G. Watson, The Theory of Bessel Functions, Cambridge Univ. Press, 1922.

Crossrefs

Cf. A086677; A055142 (for this sequence, |a(n+1)| + 1 is the number of distinct products which can be formed using commutative, nonassociative multiplication and a nonempty subset of n given variables).
Constant terms of polynomials in A098503. First row of array A099020.
Subsequence of A248652.
Cf. A082161 (relaxed compacted binary trees of unbounded right height).
Cf. A053871 (binomial transform).

Programs

  • GAP
    A001147 := function(n) local i, s, t; t := 1; i := 0; Print(t, ", "); for i in [1 .. n] do t := t*(2*i-1); Print(t, ", "); od; end; A001147(100); # Stefano Spezia, Nov 13 2018
    
  • Haskell
    a001147 n = product [1, 3 .. 2 * n - 1]
    a001147_list = 1 : zipWith (*) [1, 3 ..] a001147_list
    -- Reinhard Zumkeller, Feb 15 2015, Dec 03 2011
    
  • Magma
    A001147:=func< n | n eq 0 select 1 else &*[ k: k in [1..2*n-1 by 2] ] >; [ A001147(n): n in [0..20] ]; // Klaus Brockhaus, Jun 22 2011
    
  • Magma
    I:=[1,3]; [1] cat [n le 2 select I[n]  else (3*n-2)*Self(n-1)-(n-1)*(2*n-3)*Self(n-2): n in [1..25] ]; // Vincenzo Librandi, Feb 19 2015
    
  • Maple
    f := n->(2*n)!/(n!*2^n);
    A001147 := proc(n) doublefactorial(2*n-1); end: # R. J. Mathar, Jul 04 2009
    A001147 := n -> 2^n*pochhammer(1/2, n); # Peter Luschny, Aug 09 2009
    G(x):=(1-2*x)^(-1/2): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..19); # Zerinvary Lajos, Apr 03 2009; aligned with offset by Johannes W. Meijer, Aug 11 2009
    series(hypergeom([1,1/2],[],2*x),x=0,20); # Mark van Hoeij, Apr 07 2013
  • Mathematica
    Table[(2 n - 1)!!, {n, 0, 19}] (* Robert G. Wilson v, Oct 12 2005 *)
    a[ n_] := 2^n Gamma[n + 1/2] / Gamma[1/2]; (* Michael Somos, Sep 18 2014 *)
    Join[{1}, Range[1, 41, 2]!!] (* Harvey P. Dale, Jan 28 2017 *)
    a[ n_] := If[ n < 0, (-1)^n / a[-n], SeriesCoefficient[ Product[1 - (1 - x)^(2 k - 1), {k, n}], {x, 0, n}]]; (* Michael Somos, Jun 27 2017 *)
    (2 Range[0, 20] - 1)!! (* Eric W. Weisstein, Jul 22 2017 *)
  • Maxima
    a(n):=if n=0 then 1 else sum(sum(binomial(n-1,i)*binomial(n-i-1,j)*a(i)*a(j)*a(n-i-j-1),j,0,n-i-1),i,0,n-1); /* Vladimir Kruchinin, May 06 2020 */
  • PARI
    {a(n) = if( n<0, (-1)^n / a(-n), (2*n)! / n! / 2^n)}; /* Michael Somos, Sep 18 2014 */
    
  • PARI
    x='x+O('x^33); Vec(serlaplace((1-2*x)^(-1/2))) \\ Joerg Arndt, Apr 24 2011
    
  • Python
    from sympy import factorial2
    def a(n): return factorial2(2 * n - 1)
    print([a(n) for n in range(101)])  # Indranil Ghosh, Jul 22 2017
    
  • Sage
    [rising_factorial(n+1,n)/2^n for n in (0..15)] # Peter Luschny, Jun 26 2012
    

Formula

E.g.f.: 1 / sqrt(1 - 2*x).
D-finite with recurrence: a(n) = a(n-1)*(2*n-1) = (2*n)!/(n!*2^n) = A010050(n)/A000165(n).
a(n) ~ sqrt(2) * 2^n * (n/e)^n.
Rational part of numerator of Gamma(n+1/2): a(n) * sqrt(Pi) / 2^n = Gamma(n+1/2). - Yuriy Brun, Ewa Dominowska (brun(AT)mit.edu), May 12 2001
With interpolated zeros, the sequence has e.g.f. exp(x^2/2). - Paul Barry, Jun 27 2003
The Ramanujan polynomial psi(n+1, n) has value a(n). - Ralf Stephan, Apr 16 2004
a(n) = Sum_{k=0..n} (-2)^(n-k)*A048994(n, k). - Philippe Deléham, Oct 29 2005
Log(1 + x + 3*x^2 + 15*x^3 + 105*x^4 + 945*x^5 + 10395*x^6 + ...) = x + 5/2*x^2 + 37/3*x^3 + 353/4*x^4 + 4081/5*x^5 + 55205/6*x^6 + ..., where [1, 5, 37, 353, 4081, 55205, ...] = A004208. - Philippe Deléham, Jun 20 2006
1/3 + 2/15 + 3/105 + ... = 1/2. [Jolley eq. 216]
Sum_{j=1..n} j/a(j+1) = (1 - 1/a(n+1))/2. [Jolley eq. 216]
1/1 + 1/3 + 2/15 + 6/105 + 24/945 + ... = Pi/2. - Gary W. Adamson, Dec 21 2006
a(n) = (1/sqrt(2*Pi))*Integral_{x>=0} x^n*exp(-x/2)/sqrt(x). - Paul Barry, Jan 28 2008
a(n) = A006882(2n-1). - R. J. Mathar, Jul 04 2009
G.f.: 1/(1-x-2x^2/(1-5x-12x^2/(1-9x-30x^2/(1-13x-56x^2/(1- ... (continued fraction). - Paul Barry, Sep 18 2009
a(n) = (-1)^n*subs({log(e)=1,x=0},coeff(simplify(series(e^(x*t-t^2/2),t,2*n+1)),t^(2*n))*(2*n)!). - Leonid Bedratyuk, Oct 31 2009
a(n) = 2^n*gamma(n+1/2)/gamma(1/2). - Jaume Oliver Lafont, Nov 09 2009
G.f.: 1/(1-x/(1-2x/(1-3x/(1-4x/(1-5x/(1- ...(continued fraction). - Aoife Hennessy (aoife.hennessy(AT)gmail.com), Dec 02 2009
The g.f. of a(n+1) is 1/(1-3x/(1-2x/(1-5x/(1-4x/(1-7x/(1-6x/(1-.... (continued fraction). - Paul Barry, Dec 04 2009
a(n) = Sum_{i=1..n} binomial(n,i)*a(i-1)*a(n-i). - Vladimir Shevelev, Sep 30 2010
E.g.f.: A(x) = 1 - sqrt(1-2*x) satisfies the differential equation A'(x) - A'(x)*A(x) - 1 = 0. - Vladimir Kruchinin, Jan 17 2011
a(n) = A123023(2*n). - Michael Somos, Jul 24 2011
a(n) = (1/2)*Sum_{i=1..n} binomial(n+1,i)*a(i-1)*a(n-i). See link above. - Dennis P. Walsh, Dec 02 2011
a(n) = Sum_{k=0..n} (-1)^k*binomial(2*n,n+k)*Stirling_1(n+k,k) [Kauers and Ko].
a(n) = A035342(n, 1), n >= 1 (first column of triangle).
a(n) = A001497(n, 0) = A001498(n, n), first column, resp. main diagonal, of Bessel triangle.
From Gary W. Adamson, Jul 19 2011: (Start)
a(n) = upper left term of M^n and sum of top row terms of M^(n-1), where M = a variant of the (1,2) Pascal triangle (Cf. A029635) as the following production matrix:
1, 2, 0, 0, 0, ...
1, 3, 2, 0, 0, ...
1, 4, 5, 2, 0, ...
1, 5, 9, 7, 2, ...
...
For example, a(3) = 15 is the left term in top row of M^3: (15, 46, 36, 8) and a(4) = 105 = (15 + 46 + 36 + 8).
(End)
G.f.: A(x) = 1 + x/(W(0) - x); W(k) = 1 + x + x*2*k - x*(2*k + 3)/W(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 17 2011
a(n) = Sum_{i=1..n} binomial(n,i-1)*a(i-1)*a(n-i). - Dennis P. Walsh, Dec 02 2011
a(n) = A009445(n) / A014481(n). - Reinhard Zumkeller, Dec 03 2011
a(n) = (-1)^n*Sum_{k=0..n} 2^(n-k)*s(n+1,k+1), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n) = (2*n)4! = Gauss_factorial(2*n,4) = Product{j=1..2*n, gcd(j,4)=1} j. - Peter Luschny, Oct 01 2012
G.f.: (1 - 1/Q(0))/x where Q(k) = 1 - x*(2*k - 1)/(1 - x*(2*k + 2)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 19 2013
G.f.: 1 + x/Q(0), where Q(k) = 1 + (2*k - 1)*x - 2*x*(k + 1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 01 2013
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - 2*x*(2*k + 1)/(2*x*(2*k + 1) - 1 + 2*x*(2*k + 2)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 31 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x/(x + 1/(2*k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013
G.f.: G(0), where G(k) = 1 + 2*x*(4*k + 1)/(4*k + 2 - 2*x*(2*k + 1)*(4*k + 3)/(x*(4*k + 3) + 2*(k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 22 2013
a(n) = (2*n - 3)*a(n-2) + (2*n - 2)*a(n-1), n > 1. - Ivan N. Ianakiev, Jul 08 2013
G.f.: G(0), where G(k) = 1 - x*(k+1)/(x*(k+1) - 1/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Aug 04 2013
a(n) = 2*a(n-1) + (2n-3)^2*a(n-2), a(0) = a(1) = 1. - Philippe Deléham, Oct 27 2013
G.f. of reciprocals: Sum_{n>=0} x^n/a(n) = 1F1(1; 1/2; x/2), confluent hypergeometric Function. - R. J. Mathar, Jul 25 2014
0 = a(n)*(+2*a(n+1) - a(n+2)) + a(n+1)*(+a(n+1)) for all n in Z. - Michael Somos, Sep 18 2014
a(n) = (-1)^n / a(-n) = 2*a(n-1) + a(n-1)^2 / a(n-2) for all n in Z. - Michael Somos, Sep 18 2014
From Peter Bala, Feb 18 2015: (Start)
Recurrence equation: a(n) = (3*n - 2)*a(n-1) - (n - 1)*(2*n - 3)*a(n-2) with a(1) = 1 and a(2) = 3.
The sequence b(n) = A087547(n), beginning [1, 4, 52, 608, 12624, ... ], satisfies the same second-order recurrence equation. This leads to the generalized continued fraction expansion lim_{n -> infinity} b(n)/a(n) = Pi/2 = 1 + 1/(3 - 6/(7 - 15/(10 - ... - n*(2*n - 1)/((3*n + 1) - ... )))). (End)
E.g.f of the sequence whose n-th element (n = 1,2,...) equals a(n-1) is 1-sqrt(1-2*x). - Stanislav Sykora, Jan 06 2017
Sum_{n >= 1} a(n)/(2*n-1)! = exp(1/2). - Daniel Suteu, Feb 06 2017
a(n) = A028338(n, 0), n >= 0. - Wolfdieter Lang, May 27 2017
a(n) = (Product_{k=0..n-2} binomial(2*(n-k),2))/n!. - Stefano Spezia, Nov 13 2018
a(n) = Sum_{i=0..n-1} Sum_{j=0..n-i-1} C(n-1,i)*C(n-i-1,j)*a(i)*a(j)*a(n-i-j-1), a(0)=1, - Vladimir Kruchinin, May 06 2020
From Amiram Eldar, Jun 29 2020: (Start)
Sum_{n>=1} 1/a(n) = sqrt(e*Pi/2)*erf(1/sqrt(2)), where erf is the error function.
Sum_{n>=1} (-1)^(n+1)/a(n) = sqrt(Pi/(2*e))*erfi(1/sqrt(2)), where erfi is the imaginary error function. (End)
G.f. of reciprocals: R(x) = Sum_{n>=0} x^n/a(n) satisfies (1 + x)*R(x) = 1 + 2*x*R'(x). - Werner Schulte, Nov 04 2024

Extensions

Removed erroneous comments: neither the number of n X n binary matrices A such that A^2 = 0 nor the number of simple directed graphs on n vertices with no directed path of length two are counted by this sequence (for n = 3, both are 13). - Dan Drake, Jun 02 2009

A000680 a(n) = (2n)!/2^n.

Original entry on oeis.org

1, 1, 6, 90, 2520, 113400, 7484400, 681080400, 81729648000, 12504636144000, 2375880867360000, 548828480360160000, 151476660579404160000, 49229914688306352000000, 18608907752179801056000000, 8094874872198213459360000000, 4015057936610313875842560000000
Offset: 0

Views

Author

Keywords

Comments

Denominators in the expansion of cos(sqrt(2)*x) = 1 - (sqrt(2)*x)^2/2! + (sqrt(2)*x)^4/4! - (sqrt(2)*x)^6/6! + ... = 1 - x^2 + x^4/6 - x^6/90 + ... By Stirling's formula in A000142: a(n) ~ 2^(n+1) * (n/e)^(2n) * sqrt(Pi*n) - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 20 2001
a(n) is also the constant term in the product: Product_{1<=i, j<=n, i!=j} (1 - x_i/x_j)^2. - Sharon Sela (sharonsela(AT)hotmail.com), Feb 12 2002
a(n) is also the number of lattice paths in the n-dimensional lattice [0..2]^n. - T. D. Noe, Jun 06 2002
Representation as the n-th moment of a positive function on the positive half-axis: a(n) = Integral_{x>=0} (x^n*exp(-sqrt(2*x))/sqrt(2*x)), n=0,1,... - Karol A. Penson, Mar 10 2003
Number of permutations of [2n] with no increasing runs of odd length. Example: a(2) = 6 because we have 1234, 13/24, 14/23, 23/14, 24/13 and 34/12 (runs separated by slashes). - Emeric Deutsch, Aug 29 2004
This is also the number of ways of arranging the elements of n distinct pairs, assuming the order of elements is significant and the pairs are distinguishable. When the pairs are not distinguishable, see A001147 and A132101. For example, there are 6 ways of arranging 2 pairs [1,1], [2,2]: {[1122], [1212], [1221], [2211], [2121], [2112]}. - Ross Drewe, Mar 16 2008
n married couples are seated in a row so that every wife is to the left of her husband. The recurrence a(n+1) = a(n)*((2*n + 1) + binomial(2*n+1, 2)) conditions on whether the (n+1)st couple is seated together or separated by at least one other person. - Geoffrey Critzer, Jun 10 2009
a(n) is the number of functions f:[2n]->[n] such that the preimage of {y} has cardinality 2 for every y in [n]. Note that [k] denotes the set {1,2,...,k} and [0] denotes the empty set. - Dennis P. Walsh, Nov 17 2009
a(n) is also the number of n X 2n (0,1)-matrices with row sum 2 and column sum 1. - Shanzhen Gao, Feb 12 2010
Number of ways that 2n people of different heights can be arranged (for a photograph) in two rows of equal length so that every person in the front row is shorter than the person immediately behind them in the back row.
a(n) is the number of functions f:[n]->[n^2] such that, if floor((f(x))^.5) = floor((f(y))^.5), then x = y. For example, with n = 4, the range of f consists of one element from each of the four sets {1,2,3}, {4,5,6,7,8}, {9,10,11,12,13,14,15}, and {16}. Hence there are 1*3*5*7 = 105 ways to choose the range for f, and there are 4! ways to injectively map {1,2,3,4} to the four elements of the range. Thus there are 105*24 = 2520 such functions. Note also that a(n) = n!*(product of the first n odd numbers). - Dennis P. Walsh, Nov 28 2012
a(n) is also the 2*n th difference of n-powers of A000217 (triangular numbers). For example a(2) is the 4th difference of the squares of triangular numbers. - Enric Reverter i Bigas, Jun 24 2013
a(n) is the multinomial coefficient (2*n) over (2, 2, 2, ..., 2) where there are n 2's in the last parenthesis. It is therefore also the number of words of length 2n obtained with n letters, each letter appearing twice. - Robert FERREOL, Jan 14 2018
Number of ways to put socks and shoes on an n-legged animal, if a sock must be put on before a shoe. - Daniel Bishop, Jan 29 2018

Examples

			For n = 2, a(2) = 6 since there are 6 functions f:[4]->[2] with size 2 preimages for both {1} and {2}. In this case, there are binomial(4, 2) = 6 ways to choose the 2 elements of [4] f maps to {1} and the 2 elements of [4] that f maps to {2}. - _Dennis P. Walsh_, Nov 17 2009
		

References

  • G. E. Andrews, R. Askey and R. Roy, Special Functions, Cambridge University Press, 1998.
  • H. T. Davis, Tables of the Mathematical Functions. Vols. 1 and 2, 2nd ed., 1963, Vol. 3 (with V. J. Fisher), 1962; Principia Press of Trinity Univ., San Antonio, TX, Vol. 2, p. 283.
  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 112.
  • Shanzhen Gao and Kenneth Matheis, Closed formulas and integer sequences arising from the enumeration of (0,1)-matrices with row sum two and some constant column sums. In Proceedings of the Forty-First Southeastern International Conference on Combinatorics, Graph Theory and Computing. Congr. Numer. 202 (2010), 45-53.
  • 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).
  • C. B. Tompkins, Methods of successive restrictions in computational problems involving discrete variables. 1963, Proc. Sympos. Appl. Math., Vol. XV pp. 95-106; Amer. Math. Soc., Providence, R.I.

Crossrefs

A diagonal of the triangle in A241171.
Main diagonal of A267479, row sums of A267480.
Row n=2 of A089759.
Column n=2 of A187783.
Even bisection of column k=0 of A097591.

Programs

  • Maple
    A000680 := n->(2*n)!/(2^n);
    a[0]:=1:a[1]:=1:for n from 2 to 50 do a[n]:=a[n-1]*(2*n-1)*n od: seq(a[n], n=0..16); # Zerinvary Lajos, Mar 08 2008
    seq(product(binomial(2*n-2*k,2),k=0..n-1),n=0..16); # Dennis P. Walsh, Nov 17 2009
  • Mathematica
    Table[Product[Binomial[2 i, 2], {i, 1, n}], {n, 0, 16}]
    polygorial[k_, n_] := FullSimplify[ n!/2^n (k -2)^n*Pochhammer[2/(k -2), n]]; Array[ polygorial[6, #] &, 17, 0] (* Robert G. Wilson v, Dec 26 2016 *)
    Table[(2n)!/2^n,{n,0,20}] (* Harvey P. Dale, Sep 21 2020 *)
  • PARI
    a(n) = (2*n)! / 2^n

Formula

E.g.f.: 1/(1 - x^2/2) (with interpolating zeros). - Paul Barry, May 26 2003
a(n) = polygorial(n, 6) = (A000142(n)/A000079(n))*A001813(n) = (n!/2^n)*Product_{i=0..n-1} (4*i + 2) = (n!/2^n)*4^n*Pochhammer(1/2, n) = gamma(2*n+1)/2^n. - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
a(n) = A087127(n,2*n) = Sum_{i=0..2*n} (-1)^(2*n-i)*binomial(2*n, i)*binomial(i+2, 2)^n. Let T(n,k,j) = ((n - k + j)*(2*n - 2*k + 1))^n*binomial(2*n, 2*k-j+1) then a(n) = Sum{k=0..n} (T(n,k,1) - T(n,k,0)). For example a(12) = A087127(12,24) = Sum_{k=0..12} (T(12,k,1) - T(12,k,0)) = 24!/2^12. - André F. Labossière, Mar 29 2004 [Corrected by Jianing Song, Jan 08 2019]
For even n, a(n) = binomial(2n, n)*(a(n/2))^2. For odd n, a(n) = binomial(2n, n+1)*a((n+1)/2)*a((n-1)/2). For positive n, a(n) = binomial(2n, 2)*a(n-1) with a(0) = 1. - Dennis P. Walsh, Nov 17 2009
a(n) = Product_{i=1..n} binomial(2i, 2).
a(n) = a(n-1)*binomial(2n, 2).
From Peter Bala, Feb 21 2011: (Start)
a(n) = Product_{k = 0..n-1} (T(n) - T(k)), where T(n) = n*(n + 1)/2 is the n-th triangular number.
Compare with n! = Product_{k = 0..n-1} (n - k).
Thus we may view a(n) as a generalized factorial function associated with the triangular numbers A000217. Cf. A010050. The corresponding generalized binomial coefficients a(n)/(a(k)*a(n-k)) are triangle A086645. Also cf. A186432.
a(n) = n*(n + n-1)*(n + n-1 + n-2)*...*(n + n-1 + n-2 + ... + 1).
For example, a(5) = 5*(5+4)*(5+4+3)*(5+4+3+2)*(5+4+3+2+1) = 113400. (End).
G.f.: 1/U(0) where U(k)= x*(2*k - 1)*k + 1 - x*(2*k + 1)*(k + 1)/U(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Oct 28 2012
a(n) = n!*(product of the first n odd integers). - Dennis P. Walsh, Nov 28 2012
a(0) = 1, a(n) = a(n-1)*T(2*n-1), where T(n) is the n-th triangular number. For example: a(4) = a(3)*T(7) = 90*28 = 2520. - Enric Reverter i Bigas, Jun 24 2013
E.g.f.: 1/(1 - x/(1 - 2*x/(1 - 3*x/(1 - 4*x/(1 - 5*x/(1 - ...)))))), a continued fraction. - Ilya Gutkovskiy, May 10 2017
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = cosh(sqrt(2)).
Sum_{n>=0} (-1)^n/a(n) = cos(sqrt(2)). (End)
D-finite with recurrence a(n) -n*(2*n-1)*a(n-1)=0. - R. J. Mathar, Jan 28 2022
a(n) = n *A007019(n-1), n>0. - R. J. Mathar, Jan 28 2022

A047974 a(n) = a(n-1) + 2*(n-1)*a(n-2).

Original entry on oeis.org

1, 1, 3, 7, 25, 81, 331, 1303, 5937, 26785, 133651, 669351, 3609673, 19674097, 113525595, 664400311, 4070168161, 25330978113, 163716695587, 1075631907655, 7296866339961, 50322142646161, 356790528924523, 2570964805355607, 18983329135883665, 142389639792952801, 1091556096587136051
Offset: 0

Views

Author

Keywords

Comments

Related to partially ordered sets. - Detlef Pauly (dettodet(AT)yahoo.de), Sep 25 2003
The number of partial permutation matrices P in GL_n with P^2=0. Alternatively, the number of orbits of the Borel group of upper triangular matrices acting by conjugation on the set of matrices M in GL_n with M^2=0. - Brian Rothbach (rothbach(AT)math.berkeley.edu), Apr 16 2004
Number of ways to use the elements of {1..n} once each to form a collection of sequences, each having length 1 or 2. - Bob Proctor, Apr 18 2005
Hankel transform is A108400. - Paul Barry, Feb 11 2008
This is also the number of subsets of equivalent ways to arrange the elements of n pairs, when equivalence is defined under the joint operation of (optional) reversal of elements combined with permutation of the labels and the subset maps to itself. - Ross Drewe, Mar 16 2008
Equals inverse binomial transform of A000898. - Gary W. Adamson, Oct 06 2008
a(n) is also the moment of order n for the measure of density exp(-(x-1)^2/4)/(2*sqrt(Pi)) over the interval -oo..oo. - Groux Roland, Mar 26 2011
The n-th term gives the number of fixed-point-free involutions in S_n^B, the group of permutations on the set {-n,...,-1,1,2,...,n}. - Matt Watson, Jul 26 2012
From Peter Bala, Dec 03 2017: (Start)
a(n+k) == a(n) (mod k) for all n and k. Hence for each k, the sequence a(n) taken modulo k is a periodic sequence and the exact period divides k. Cf. A115329.
More generally, the same divisibility property holds for any sequence with an e.g.f. of the form F(x)*exp(x*G(x)), where F(x) and G(x) are power series with integer coefficients and G(0) = 1. See the Bala link for a proof. (End)

Crossrefs

Row sums of A067147.
Column k=2 of A359762.
Sequences with e.g.f = exp(x + q*x^2): A158968 (q=-9), A158954 (q=-4), A362177 (q=-3), A362176 (q=-2), A293604 (q=-1), A000012 (q=0), this sequence (q=1), A115329 (q=2), A293720 (q=4).

Programs

  • MATLAB
    N = 18; A = zeros(N,1); for n = 1:N; a = factorial(n); s = 0; k = 0; while k <= floor(n/2); b = factorial(n - 2*k); c = factorial(k); s = s + a/(b*c); k = k+1; end; A(n) = s; end; disp(A); % Ross Drewe, Mar 16 2008
    
  • Magma
    [n le 2 select 1 else Self(n-1) + 2*(n-2)*Self(n-2): n in [1..40]]; // G. C. Greubel, Jul 12 2024
    
  • Maple
    seq( add(n!/((n-2*k)!*k!), k=0..floor(n/2)), n=0..30 ); # Detlef Pauly (dettodet(AT)yahoo.de), Nov 15 2001
    with(combstruct):seq(count(([S,{S=Set(Union(Z,Prod(Z,Z)))},labeled],size=n)),n=0..30); # Detlef Pauly (dettodet(AT)yahoo.de), Sep 25 2003
    A047974 := n -> I^(-n)*orthopoly[H](n, I/2):
    seq(A047974(n), n=0..26); # Peter Luschny, Nov 29 2017
  • Mathematica
    Range[0, 23]!*CoefficientList[ Series[ Exp[x*(1-x^2)/(1 - x)], {x, 0,23 }], x] - (* Zerinvary Lajos, Mar 23 2007 *)
    Table[I^(-n)*HermiteH[n, I/2], {n, 0, 23}] - (* Alyssa Byrnes and C. Vignat, Jan 31 2013 *)
  • PARI
    my(x='x+O('x^66)); Vec(serlaplace(exp(x^2+x))) \\ Joerg Arndt, May 04 2013
    
  • SageMath
    [(-i)^n*hermite(n,i/2) for n in range(41)] # G. C. Greubel, Jul 12 2024

Formula

E.g.f.: exp(x^2+x). - Len Smiley, Dec 11 2001
Binomial transform of A001813 (with interpolated zeros). - Paul Barry, May 09 2003
a(n) = Sum_{k=0..n} C(k,n-k)*n!/k!. - Paul Barry, Mar 29 2007
a(n) = Sum_{k=0..floor(n/2)} C(n,2k)*(2k)!/k!; - Paul Barry, Feb 11 2008
G.f.: 1/(1-x-2*x^2/(1-x-4*x^2/(1-x-6*x^2/(1-x-8*x^2/(1-... (continued fraction). -Paul Barry, Apr 10 2009
E.g.f.: Q(0); Q(k) = 1+(x^2+x)/(2*k+1-(x^2+x)*(2*k+1)/((x^2+x)+(2*k+2)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Nov 24 2011
a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator sqrt(1+4*x)*d/dx. Cf. A000085 and A115329. - Peter Bala, Dec 07 2011
a(n) ~ 2^(n/2 - 1/2)*exp(sqrt(n/2) - n/2 - 1/8)*n^(n/2). - Vaclav Kotesovec, Oct 08 2012
E.g.f.: 1 + x*(E(0)-1)/(x+1) where E(k) = 1 + (1+x)/(k+1)/(1-x/(x+1/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
a(n) = i^(-n)*H_{n}(i/2) with i the imaginary unit and H_{n} the Hermite polynomial of degree n. - Alyssa Byrnes and C. Vignat, Jan 31 2013
E.g.f.: -Q(0)/x where Q(k) = 1 - (1+x)/(1 - x/(x - (k+1)/Q(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Mar 06 2013
G.f.: 1/Q(0), where Q(k) = 1 + x*2*k - x/(1 - x*(2*k+2)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, Apr 17 2013
E.g.f.: E(0)-1-x-x^2, where E(k) = 2 + 2*x*(1+x) - 8*k^2 + x^2*(1+x)^2*(2*k+3)*(2*k-1)/E(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 21 2013
E.g.f.: Product_{k>=1} 1/(1 + (-x)^k)^(mu(k)/k). - Ilya Gutkovskiy, May 26 2019
a(n) = Sum_{k=0..floor(n/2)} 2^k*B(n, k), where B are the Bessel numbers A100861. - Peter Luschny, Jun 04 2021

A250001 Number of arrangements of n circles in the affine plane.

Original entry on oeis.org

1, 1, 3, 14, 173, 16951
Offset: 0

Views

Author

Jon Wild, May 16 2014

Keywords

Comments

Two circles are either disjoint or meet in two points. Tangential contacts are not allowed. A point belongs to exactly one or two circles. Three circles may not meet at a point. The circles may have different radii.
This is in the affine plane, rather than the projective plane.
Two arrangements are considered the same if one can be continuously changed to the other while keeping all circles circular (although the radii may be continuously changed), without changing the multiplicity of intersection points, and without a circle passing through an intersection point. Turning the whole configuration over is allowed.
Several variations are possible:
- straight lines instead of circles (see A241600).
- straight lines in general position (see A090338).
- curved lines in general position (see A090339).
- allow circles to meet tangentially but without multiple intersection points (begins 1, 5, ...); more terms are needed.
- again use circles, but allow multiple intersection points (also begins 1, 5, ...); more terms are needed.
- use ellipses rather than circles.
- a question from Walter D. Wallis: what if the circles must all have the same radius?
a(1)-a(5) computed by Jon Wild.
a(n) >= A000081(n+1) - Benoit Jubin, Dec 21 2014. More precisely, there are A000081(n+1) ways to arrange n circles if no two of them meet. - N. J. A. Sloane, May 16 2017
From Daniel Forgues, Aug 08-09 2015: (Start)
A representation for the diagrams in a250001.jpg (in the same order):
a(1) = 1: {{2}};
a(2) = 3: {{2, 3}, {2, 4}, {4, 6}};
a(3) = 14: {{2, 4, 8}, {2, 3, 6}, {2, 3, 4}, {2, 3, 5}, {4, 6, 5},
{4, 6, 15}, {2, 6, 9}, {4, 6, 12}, {2, 8, 12}, {30, 42, 70},
{?, ?, ?}, {?, ?, ?}, {15, 21, 35}, {?, ?, ?}}.
In lexicographic order:
a(3) = 14: {{2, 3, 4}, {2, 3, 5}, {2, 3, 6}, {2, 4, 8}, {2, 6, 9},
{2, 8, 12}, {4, 6, 5}, {4, 6, 12}, {4, 6, 15}, {15, 21, 35},
{30, 42, 70}, {?, ?, ?}, {?, ?, ?}, {?, ?, ?}}.
The smallest integers greater than 1 are used for the representation:
(p_1)^(a_1)*...*(p_m)^(a_m), where
0 <= a_i <= n, for 1 <= i <= m;
(a_1)+...+(a_m) > 0.
Could the Venn diagram interpretation (of the k-wise, 1 <= k <= n, common divisors of k numbers from each subset) reveal a pattern?
Does this representation work for more complex diagrams? (End)
Once you get to n=5, geometric concerns mean that not all topologically-conceivable arrangements are actually circle-drawable. My program enumerated 16977 conceivable arrangements of 5 pseudo-circles, and Christopher Jones and I together have figured out how to show that 26 of these arrangements are not actually circle-drawable. So it seems that a(5) = 16951. This entry will be updated soon, and there will be a new sequence for the number of topologically-conceivable arrangements. - Jon Wild, Aug 25 2016 [The counts in this comment were amended by Jon Wild on Aug 30 2016. I apologize for taking so long to make the corrections here. - N. J. A. Sloane, Jun 11 2017]
a(n) <= 7*13^(binomial(n,3) + binomial(n,2) + 3n - 1) is a (loose) upper bound, see Reddit link. I believe XkF21WNJ's reply shaves off a factor of 13^3 from this bound for all n > 1. - Linus Hamilton, Apr 14 2019
A good upper bound for a(6) is given in sequence A288559, which counts the arrangements of pseudo-circles, i.e. the topologically conceivable arrangements mentioned above, which are not all necessarily realizable with true circles. The number of arrangements of 6 pseudo-circles was found by Andrii Shportko and Jon Wild to be 17,552,169. - Jon Wild, Jun 03 2025
In A288559, a(5) included 26 non-circularizable pseudocircle arrangements, which generated in turn 132,546 6-pseudocircle descendants. These descendants must be excluded from A250001, which means that a tighter upper bound for A250001(6) is 17,419,623. - Andrii Shportko, Jun 06 2025

Examples

			a(2) = 3, because two circles can either be next to each other, overlap with two intersection points, or one may be located within the other (of larger radius). (As per the first comment, the limiting case where they touch in one point is [somewhat arbitrarily] excluded. This would add two more independent configurations, where one touched the other "from inside" or "from outside".) - _M. F. Hasler_, May 03 2025
		

References

  • Jon Wild, Posting to Sequence Fans Mailing List, May 15 2014.

Crossrefs

Row sums of A261070.
Apart from first term, row sums of triangles A249752, A252158, A285996, A274776, A274777.
See A275923 and A275924 for the connected arrangements. See also A288554-A288568.
Cf. A132101 (one-dimensional analog).

Extensions

a(4) is 173, not 168. Corrected by Jon Wild, Aug 08 2015
A duplicate pair of configurations in an older file was spotted by Manfred Scheucher, Aug 13 2016. The pdf and svg files here are now correct.

A132100 Number of distinct Tsuro tiles which are square and have n points per side.

Original entry on oeis.org

1, 2, 35, 2688, 508277, 163715822, 79059439095, 53364540054860, 47974697008198313, 55410773910104281242, 79957746695043660483467, 140965507420235075126987480, 298142048193613276717321211805, 745056978435827991570581878537478
Offset: 0

Views

Author

Keith F. Lynch, Oct 31 2007

Keywords

Comments

Turning over is not allowed, but rotation of the tile is allowed.
In the original Tsuro game the tiles are square and have two points on each side, one third and two thirds of the way along the side and arcs connecting these eight points in various ways.
The shapes of the arcs aren't significant, only which two points they connect is.
Each point is connected to exactly one other point.
There are 35 tiles, agreeing with the entry a(4) = 35 here.
If we vary the shape of the tile and the number of points per side (pps), we get the following table.
....pps:..0....1......2......3......4......5......6......7......8......9.....10
-------------------------------------------------------------------------------
circle....1....0......1......0......2......0......5......0.....18......0....105 (A007769)
monogon...1....0......1......0......3......0.....15......0....105......0....945 (A001147)
digon.....1....1......3.....11.....65....513...5363..68219 .................... (A132101)
triangle..1....0......7......0...3483......0.............0.............0
square....1....2.....35...2688.508277 ......................................... (this entry)
pentagon..1....0....193......0.............0.............0.............0
hexagon...1....5...1799
heptagon..1....0..19311......0.............0.............0.............0
octagon...1...18.254143
9-gon.....1....0.............0.............0.............0.............0
10-gon....1..105
The pps = 2 column is A132102. Blank entries all represent numbers greater than one million.
A monogon is distinct from a circle in that a monogon has not just one side, but also one vertex. Monogons and digons can't exist with straight sides, of course, at least not on a flat plane, but there's no rule that says these tiles have to have straight sides.
If we allow reflections the numbers are smaller (this would be appropriate for a game where the tiles were transparent and could be flipped over):
....pps:..0....1......2......3......4......5......6......7......8......9.....10
-------------------------------------------------------------------------------
circle....1....0......1......0......2......0......5......0.....17......0.....79 (A054499)
monogon...1....0......1......0......3......0.....11......0.....65......0....513 (A132101)
digon.....1....1......3......8.....45....283...2847..34518.511209 ............. (A132103)
triangle..1....0......7......0...1907......0.............0.............0
square....1....2.....30...1447.257107 ......................................... (A132104)
pentagon..1....0....137......0.............0.............0.............0
hexagon...1....5...1065
heptagon..1....0..10307......0.............0.............0.............0
octagon...1...17.130040
9-gon.....1....0.............0.............0.............0.............0
10-gon....1...79
The pps = 2 column is A132105.

Crossrefs

Programs

  • Maple
    # A(n,m) gives the number of n-sided tiles with m points per side (cf. comments)
    # B(n,m) enumerates these tiles, also allowing reflections
    with(numtheory): a:=(p,r)->piecewise(p mod 2 = 1,p^(r/2)*doublefactorial(r-1), sum(p^j*binomial(r, 2*j)*doublefactorial(2*j - 1), j = 0 .. floor(r/2)));
    A := (n,m)->piecewise(n*m mod 2=1,0,add(phi(p)*a(p,m*n/p),p in divisors(n))/n);
    B := (n,m)->A(n,m)/2+piecewise(n*m mod 2=0,piecewise(m mod 2=0,a(2,m*n/2)*2, a(2,m*n/2)+a(2,m*n/2-1))/4,0);
    A132100 := m -> A(4,m);[seq(A132100(m),m=1..15)]; # Laurent Tournier, Jul 09 2014

Formula

From Laurent Tournier, Jul 09 2014: (Start)
a(m) = ((4m-1)!! + sum_{j=0..m} 2^j binomial(2m,2j) (2j-1)!! + 2 sum_{0<=2j<=m} 4^j binomial(m, 2j) (2j-1)!!)/4
More generally, if A(n,m) is the number of n-sided tiles with m points per side (with nm even),
A(n,m) = 1/n sum_{n=pq} phi(p)*alpha(p,mq), phi = Euler's totient function,
alpha(p,r) = sum_{0 <= 2j <= r} p^j binomial(r,2j) (2j-1)!! if p even,
= p^(r/2) (r-1)!! if p odd.
If B(n,m) is the number of n-sided tiles with m points per side (with nm even), allowing reflections,
B(n,m) = (A(n,m) + alpha(2,mn/2))/2 if m even,
= (A(n,m) + alpha(2,mn/2)/2 + alpha(2,mn/2-1)/2)/2 if m odd.
(End)

Extensions

More terms from Laurent Tournier, Jul 09 2014

A210277 a(n) = (3*n)!/3^n.

Original entry on oeis.org

1, 2, 80, 13440, 5913600, 5381376000, 8782405632000, 23361198981120000, 94566133475573760000, 553211880832106496000000, 4492080472356704747520000000, 49017582114356362204938240000000, 699971072593008852286518067200000000
Offset: 0

Views

Author

Mohammad K. Azarian, Mar 20 2012

Keywords

Crossrefs

Programs

Formula

E.g.f.: 1/(1-x^3/3).
a(n) = Product_{i=1..n} (2*binomial(3i,3)). - James Mahoney, Apr 04 2012
From Amiram Eldar, Jan 18 2021: (Start)
Sum_{n>=0} 1/a(n) = exp(3^(1/3))/3 + (2/3)*exp(-3^(1/3)/2)*cos(3^(5/6)/2).
Sum_{n>=0} (-1)^n/a(n) = exp(-3^(1/3))/3 + (2/3)*exp(3^(1/3)/2)*cos(3^(5/6)/2). (End)

A210278 (5n)!/5^n.

Original entry on oeis.org

1, 24, 145152, 10461394944, 3892643213082624, 4963587213865915514880, 16976183027980227752723742720, 132264293969742655099733137120296960, 2088743125114618199924764850166056689336320, 61246577083125859615725138685776750112964471685120
Offset: 0

Views

Author

Mohammad K. Azarian, Mar 20 2012

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(5*n)/5^n: n in [0..10]]; // Vincenzo Librandi, Feb 15 2013
  • Mathematica
    Table[(5 n)!/5^n, {n, 0, 10}] (* Vincenzo Librandi, Feb 15 2013 *)
    With[{nn=100},Take[CoefficientList[Series[1/(1-x^5/5),{x,0,nn}],x] Range[0,nn]!,{1,-1,5}]] (* Harvey P. Dale, May 27 2025 *)

Formula

E.g.f.: 1/(1-x^5/5).

A132102 Number of distinct Tsuro tiles which are n-gonal in shape and have 2 points per side.

Original entry on oeis.org

1, 1, 3, 7, 35, 193, 1799, 19311, 254143, 3828921, 65486307, 1249937335, 26353147811, 608142583137, 15247011443103, 412685556939751, 11993674252049647, 372509404162520641, 12313505313357313047, 431620764875678503143, 15991549339008732109899
Offset: 0

Views

Author

Keith F. Lynch, Oct 31 2007

Keywords

Comments

Turning over is not allowed.
See A132100 for definition and comments.
Even and odd terms can be computed with the help of Burnside Lemma and recursive sequences. - Lionel RAVEL, Sep 18 2013

Crossrefs

Programs

  • Maple
    with(numtheory): a:=(p,q)->piecewise(p mod 2 = 1, p^q*doublefactorial(2*q - 1), sum(p^j*binomial(2*q, 2*j)*doublefactorial(2*j - 1), j = 0 .. q));
    A132102 := n->add(phi(p)*a(p,n/p),p in divisors(n))/n;
    [seq(A132102(n),n=1..20)]; # Laurent Tournier, Jul 09 2014
  • PARI
    a(n)={if(n<1, n==0, sumdiv(n, d, my(m=n/d); eulerphi(d)*sum(j=0, m, (d%2==0 || m-j==0) * binomial(2*m, 2*j) * d^j * (2*j)! / (j!*2^j) ))/n)} \\ Andrew Howroyd, Jan 26 2020

Formula

a(n) = (1/n)*Sum_{d|n} phi(d)*alpha(d, n/d), where phi = Euler's totient function,
alpha(p,q) = Sum_{j=0..q} p^j * binomial(2q, 2j) * (2j-1)!! if p even,
= p^q * (2q-1)!! if p odd. (cf. also A132100) - Laurent Tournier, Jul 09 2014

Extensions

More terms from Lionel RAVEL, Sep 18 2013
a(9) and a(10) corrected, and addition of more terms using formula given above by Laurent Tournier, Jul 09 2014

A210279 (6n)!/6^n.

Original entry on oeis.org

1, 120, 13305600, 29640619008000, 478741050720092160000, 34111736086958726676480000000, 7973107998754741458076119859200000000, 5019026197962676820927435579005599744000000000
Offset: 0

Views

Author

Mohammad K. Azarian, Apr 12 2012

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(6*n)/6^n: n in [0..10]]; // Vincenzo Librandi, Feb 15 2013
  • Mathematica
    Table[(6 n)!/6^n, {n, 0, 11}] (* Vincenzo Librandi, Feb 15 2013 *)
    With[{nn=50},Take[CoefficientList[Series[1/(1-x^6/6),{x,0,nn}],x] Range[0,nn-2]!,{1,-1,6}]] (* Harvey P. Dale, Sep 25 2023 *)

Formula

E.g.f.: 1/(1-x^6/6).

A210280 (7n)!/7^n.

Original entry on oeis.org

1, 720, 1779148800, 148953184174080000, 126983900296423931904000000, 614812159599342234168301977600000000, 11942354952042770431904585727413846016000000000
Offset: 0

Views

Author

Mohammad K. Azarian, Apr 12 2012

Keywords

Crossrefs

Programs

Formula

E.g.f.: 1/(1-x^7/7).
Showing 1-10 of 12 results. Next